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:

Regresión de Poisson

library(readxl)
library(DT)
Tos <- read_excel("C:/Users/wsand/Dropbox/Centro Oriente/Diario de Sintomas/Regresiones Poisson/Tos.xlsx", 
    col_types = c("text", "numeric", "text", 
        "numeric", "numeric", "numeric", 
        "text", "text", "text", "text", "text", 
       "text", "text", "text", "text", "text", 
        "text", "text"))
DT::datatable(Tos)
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)" ...
Tos$Estrato<-as.factor(Tos$Estrato)
poisson.model<-glm(casos ~ Exposición+ Estrato+ SEXO+ Tos$`TIPO VIVIENDA NIÑO`+Tos$`PESO MENOS 2,5 KG`, data=Tos, family = poisson(link = "log"))
summary(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