PARCIAL

library(rio)
data=import("dataPeru.xlsx")

h1=formula(buenEstado~contribuyentesSunat+peaOcupada)
reg1=lm(h1,data=data)
summary(reg1)

## 
## Call:
## lm(formula = h1, data = data)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -10.589  -3.966  -1.347   1.907  21.518 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          1.865e+01  2.694e+00   6.922 5.98e-07 ***
## contribuyentesSunat  1.786e-05  2.060e-05   0.867    0.395    
## peaOcupada          -1.596e-05  2.241e-05  -0.712    0.484    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.925 on 22 degrees of freedom
## Multiple R-squared:  0.1561, Adjusted R-squared:  0.07939 
## F-statistic: 2.035 on 2 and 22 DF,  p-value: 0.1546

library(modelsummary)

## `modelsummary` 2.0.0 now uses `tinytable` as its default table-drawing
##   backend. Learn more at: https://vincentarelbundock.github.io/tinytable/
## 
## Revert to `kableExtra` for one session:
## 
##   options(modelsummary_factory_default = 'kableExtra')
## 
## Change the default backend persistently:
## 
##   config_modelsummary(factory_default = 'gt')
## 
## Silence this message forever:
## 
##   config_modelsummary(startup_message = FALSE)

model1=list('apropiacion (I)'=reg1)
modelsummary(model1, title = "Regresion: modelo 1",
             stars = TRUE,
             output = "kableExtra")

Regresion: modelo 1

	apropiacion (I)
(Intercept)	18.646***
	(2.694)
contribuyentesSunat	0.000
	(0.000)
peaOcupada	0.000
	(0.000)
Num.Obs.	25
R2	0.156
R2 Adj.	0.079
AIC	179.3
BIC	184.1
Log.Lik.	-85.626
F	2.035
RMSE	7.43
+ p < 0.1, * p < 0.05, ** p < 0.01, *** p < 0.001

h2=formula(peaOcupada~contribuyentesSunat+buenEstado)
reg2=lm(h2,data=data)
summary(reg2)

## 
## Call:
## lm(formula = h2, data = data)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -91867 -58573 -11166  46174 155851 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          1.155e+05  3.787e+04   3.049  0.00588 ** 
## contribuyentesSunat  9.206e-01  1.741e-02  52.872  < 2e-16 ***
## buenEstado          -1.412e+03  1.983e+03  -0.712  0.48395    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 74540 on 22 degrees of freedom
## Multiple R-squared:  0.9932, Adjusted R-squared:  0.9926 
## F-statistic:  1603 on 2 and 22 DF,  p-value: < 2.2e-16

library(modelsummary)
model2=list('apropiacion (I)'=reg2)
modelsummary(model2, title = "Regresion: modelo 1",
             stars = TRUE,
             output = "kableExtra")

Regresion: modelo 1

	apropiacion (I)
(Intercept)	115468.390**
	(37869.931)
contribuyentesSunat	0.921***
	(0.017)
buenEstado	-1411.649
	(1982.682)
Num.Obs.	25
R2	0.993
R2 Adj.	0.993
AIC	636.7
BIC	641.6
Log.Lik.	-314.354
F	1603.130
RMSE	69927.80
+ p < 0.1, * p < 0.05, ** p < 0.01, *** p < 0.001