Parcial

setwd("C:/Users/tnadi/OneDrive/Documentos/2024-1/Estadística/data")

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

EXAMEN PARCIAL

#Pregunta 1

reg_ols = lm(buenEstado ~ contribuyentesSunat + peaOcupada, data = data)
summary(reg_ols)

## 
## Call:
## lm(formula = buenEstado ~ contribuyentesSunat + peaOcupada, 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

#Pregunta 2

library(knitr)
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)

h1=formula(peaOcupada~contribuyentesSunat + buenEstado)
    
rp1=glm(h1
        , data = data, 
        offset=log(pobTotal),
        family = poisson(link = "log"))

modelsPois=list('POISSON asegurados (I)'=rp1)
modelsummary(modelsPois, 
             title = "Regresiones Poisson anidadas",
             stars = TRUE,
             output = "kableExtra")

Regresiones Poisson anidadas

	POISSON asegurados (I)
(Intercept)	-1.159***
	(0.001)
contribuyentesSunat	0.000***
	(0.000)
buenEstado	0.008***
	(0.000)
Num.Obs.	25
AIC	55974.5
BIC	55978.2
Log.Lik.	-27984.274
F	68713.131
RMSE	30706.19
+ p < 0.1, * p < 0.05, ** p < 0.01, *** p < 0.001