Datos de Recuentos
Tipo de datos estadísticos en el que las observaciones pueden tomar solamente valores enteros positivos, por lo general, son aquellos que determinan el número de sucesos o eventos que ocurren en una misma unidad de observación en un intervalo espacial o temporal definido (Lindsey, 2007).
Ejemplo:
Número de personas que visitan el consultorio de un médico por mes.
Número de actos agresivos de los niños durante un período de recreo.
Regresión de Poisson
library(readxl)
library(DT)
<- read_excel("C:/Users/wsand/Dropbox/Centro Oriente/Diario de Sintomas/Regresiones Poisson/Tos.xlsx",
Tos col_types = c("text", "numeric", "text",
"numeric", "numeric", "numeric",
"text", "text", "text", "text", "text",
"text", "text", "text", "text", "text",
"text", "text"))
::datatable(Tos) DT
colnames(Tos)
## [1] "Nombre" "cod_localidad"
## [3] "localidad" "casos"
## [5] "nulos" "Dias de participacion"
## [7] "Exposición" "SEXO"
## [9] "Estrato" "PESO MENOS 2,5 KG"
## [11] "LACTANCIA MATERNA" "DX MALNUTRICIÓN"
## [13] "MADRE FUMÓ EMBARAZO" "MADRE CONTACTO FUMADORES"
## [15] "MADRE FUMÓ DOS AÑOS NIÑO(A)" "NIÑO(A) CONTACTO FUMADORES"
## [17] "DX PADRES ENF.RESPIRATORIA" "TIPO VIVIENDA NIÑO"
ls.str(Tos)
## casos : num [1:621] 0 0 0 0 0 1 0 1 1 3 ...
## cod_localidad : num [1:621] 6 14 18 15 13 7 1 5 9 3 ...
## Dias de participacion : num [1:621] 304 304 304 304 273 304 59 304 304 304 ...
## DX MALNUTRICIÓN : chr [1:621] "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" ...
## DX PADRES ENF.RESPIRATORIA : chr [1:621] "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" ...
## Estrato : chr [1:621] "2" "3" "3" "3" "3" "2" "2" "1" "2" "2" "1" "3" "3" "2" "4" ...
## Exposición : chr [1:621] "Riesgo alto" "Riesgo bajo" "Riesgo bajo" "Riesgo alto" ...
## LACTANCIA MATERNA : chr [1:621] "1" "1" "1" "1" "1" "1" "1" "1" "1" "2" "1" "2" "1" "1" "1" ...
## localidad : chr [1:621] "Tunjuelito" "Mártires" "Rafael_Uribe" "Antonio_Nariño" ...
## MADRE CONTACTO FUMADORES : chr [1:621] "2" "2" "2" "2" "2" "2" "2" "1" "2" "2" "1" "1" "2" "2" "2" ...
## MADRE FUMÓ DOS AÑOS NIÑO(A) : chr [1:621] "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "1" "2" "2" "2" "2" ...
## MADRE FUMÓ EMBARAZO : chr [1:621] "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" ...
## NIÑO(A) CONTACTO FUMADORES : chr [1:621] "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "1" "1" "2" "2" "2" ...
## Nombre : chr [1:621] "AARON MATIAS MONROY ORTIZ" "ABBY MICHELLE DELGADO RODRIGUEZ" ...
## nulos : num [1:621] 0 0 0 0 31 0 245 0 0 0 ...
## PESO MENOS 2,5 KG : chr [1:621] "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" "2" ...
## SEXO : chr [1:621] "Masculino" "Femenino" "Femenino" "Femenino" "Femenino" ...
## TIPO VIVIENDA NIÑO : chr [1:621] "Apartamento (Independiente)" "Apartamento (Independiente)" ...
$Estrato<-as.factor(Tos$Estrato) Tos
<-glm(casos ~ Exposición+ Estrato+ SEXO+ Tos$`TIPO VIVIENDA NIÑO`+Tos$`PESO MENOS 2,5 KG`, data=Tos, family = poisson(link = "log"))
poisson.modelsummary(poisson.model)
##
## Call:
## glm(formula = casos ~ Exposición + Estrato + SEXO + Tos$`TIPO VIVIENDA NIÑO` +
## Tos$`PESO MENOS 2,5 KG`, family = poisson(link = "log"),
## data = Tos)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.4330 -1.1827 -1.0658 0.3793 4.2746
##
## Coefficients:
## Estimate Std. Error
## (Intercept) -0.66532 0.42359
## ExposiciónRiesgo bajo 0.20821 0.09890
## Estrato2 0.16320 0.17252
## Estrato3 0.04735 0.17933
## Estrato4 -0.96403 0.44833
## SEXOMasculino 0.10977 0.09833
## Tos$`TIPO VIVIENDA NIÑO`Apartamento (Propiedad horizontal) -0.12286 0.21556
## Tos$`TIPO VIVIENDA NIÑO`Casa (Independiente) -0.01741 0.10625
## Tos$`TIPO VIVIENDA NIÑO`Habitación -14.89731 519.34306
## Tos$`TIPO VIVIENDA NIÑO`Inquilinato 0.26804 0.45334
## Tos$`PESO MENOS 2,5 KG`2 0.05221 0.38177
## z value Pr(>|z|)
## (Intercept) -1.571 0.1163
## ExposiciónRiesgo bajo 2.105 0.0353 *
## Estrato2 0.946 0.3442
## Estrato3 0.264 0.7918
## Estrato4 -2.150 0.0315 *
## SEXOMasculino 1.116 0.2643
## Tos$`TIPO VIVIENDA NIÑO`Apartamento (Propiedad horizontal) -0.570 0.5687
## Tos$`TIPO VIVIENDA NIÑO`Casa (Independiente) -0.164 0.8698
## Tos$`TIPO VIVIENDA NIÑO`Habitación -0.029 0.9771
## Tos$`TIPO VIVIENDA NIÑO`Inquilinato 0.591 0.5543
## Tos$`PESO MENOS 2,5 KG`2 0.137 0.8912
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for poisson family taken to be 1)
##
## Null deviance: 962.18 on 620 degrees of freedom
## Residual deviance: 936.03 on 610 degrees of freedom
## AIC: 1518.4
##
## Number of Fisher Scoring iterations: 13
anova(poisson.model, test = "Chisq")
## Analysis of Deviance Table
##
## Model: poisson, link: log
##
## Response: casos
##
## Terms added sequentially (first to last)
##
##
## Df Deviance Resid. Df Resid. Dev Pr(>Chi)
## NULL 620 962.18
## Exposición 1 3.7554 619 958.42 0.052638 .
## Estrato 3 12.3305 616 946.09 0.006333 **
## SEXO 1 1.4282 615 944.66 0.232054
## Tos$`TIPO VIVIENDA NIÑO` 4 8.6131 611 936.05 0.071532 .
## Tos$`PESO MENOS 2,5 KG` 1 0.0190 610 936.03 0.890304
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1