Nuevos cuadros
library(survey)
## Attaching package: 'survey'
## The following object(s) are masked from 'package:graphics':
##
## dotchart
options(OutDec = ",")
# load('~/Dropbox/odontologia/maestria
# licet/julio_2013/datos_licet_25072013.RData')
# load('C:/Users/usuario/Dropbox/odontologia/maestria
# licet/julio_2013/datos_licet_25072013.RData')
# load('~/Dropbox/odontologia/maestria
# licet/julio_2013/datos_licet_25072013.RData')
load("~/Dropbox/odontologia/maestria licet/setiembre2013/datos_licet_12092013.RData")
library(car)
## Loading required package: MASS
## Loading required package: nnet
# Ramon: estuve hablando con la Tana y ella tb tuvo que reclasificar estas
# variables de acuerdo al comportamiento de caries. Creo que sería
# adecuado reagrupar la variable así: unir incomplete y basic como
# elementary (326); basic y médium como high (577); médium y high como
# college (209).te pido tirar la uni así y si da significativo tiramos la
# multi. De cualquier forma para que la profe se quede contenta te pido se
# tiras la multiple tal como estaba antes pero retirando nivel
# socioeconómico que la incluiste pero no era significativa.Espero que con
# esto se haga la luz. Gracias licet
diseniopost1$variables$Nive.Educativo.de.la.Madre2.rec <- recode(diseniopost1$variables$Nive.Educativo.de.la.Madre2,
"'COLLEGE-UNIVERSITY'='3-COLLEGE-UNIVERSITY';'ELEMENTARY SCHOOL'='1-ELEMENTARY SCHOOL';'HIGH SCHOOL'='2-HIGH SCHOOL';99=NA")
addmargins(table(diseniopost1$variables$Nivel.Educativo.de.la.Madre1.rec, diseniopost1$variables$erodentina.rec))
##
## ausente presente Sum
## 1-Incomplete Education 83 7 90
## 2-Basic 656 25 681
## 3-Medium 189 7 196
## 4-High 119 8 127
## Sum 1047 47 1094
addmargins(table(diseniopost1$variables$Nive.Educativo.de.la.Madre2.rec, diseniopost1$variables$erodentina.rec))
##
## ausente presente Sum
## 1-Elementary School 308 17 325
## 2-High School 555 18 573
## 3-College-University 184 12 196
## Sum 1047 47 1094
modelo_ero_3a.logit <- svyglm(erodentina.rec ~ Nive.Educativo.de.la.Madre2.rec,
design = diseniopost1, family = quasibinomial())
summary(modelo_ero_3a.logit)
##
## Call:
## svyglm(formula = erodentina.rec ~ Nive.Educativo.de.la.Madre2.rec,
## design = diseniopost1, family = quasibinomial())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error
## (Intercept) -2,7780 0,3161
## Nive.Educativo.de.la.Madre2.rec2-High School -0,4878 0,3562
## Nive.Educativo.de.la.Madre2.rec3-College-University -0,0981 0,4784
## t value Pr(>|t|)
## (Intercept) -8,79 1,4e-10 ***
## Nive.Educativo.de.la.Madre2.rec2-High School -1,37 0,18
## Nive.Educativo.de.la.Madre2.rec3-College-University -0,21 0,84
## ---
## Signif. codes: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1
##
## (Dispersion parameter for quasibinomial family taken to be 0,9559)
##
## Number of Fisher Scoring iterations: 5
confint(modelo_ero_3a.logit)
## 2,5 % 97,5 %
## (Intercept) -3,398 -2,1583
## Nive.Educativo.de.la.Madre2.rec2-High School -1,186 0,2103
## Nive.Educativo.de.la.Madre2.rec3-College-University -1,036 0,8395
modelomulti_ero.logit <- svyglm(erodentina.rec ~ Sexo + Tipo.de.Escuela + IGS.rec +
RefrCola.rec + Bebidas_energizantes.rec + Yogurt.rec + AlterGastrica.rec,
design = diseniopost1, family = quasibinomial())
summary(modelomulti_ero.logit)
##
## Call:
## svyglm(formula = erodentina.rec ~ Sexo + Tipo.de.Escuela + IGS.rec +
## RefrCola.rec + Bebidas_energizantes.rec + Yogurt.rec + AlterGastrica.rec,
## design = diseniopost1, family = quasibinomial())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) -5,544 1,287 -4,31
## SexoM 1,199 0,478 2,51
## Tipo.de.Escuela2-Privada 0,892 0,743 1,20
## IGS.rec2-De 45 a 60 1,186 1,291 0,92
## IGS.rec3 -Menos o igual a 45 1,832 1,111 1,65
## RefrCola.rec2-Todos los dias 0,375 0,509 0,74
## RefrCola.rec3-Mas de 3 veces al dia 0,926 0,715 1,29
## Bebidas_energizantes.rec2-Todos los dias 0,401 0,312 1,28
## Bebidas_energizantes.recMas de 3 veces al dia 1,787 0,732 2,44
## Yogurt.rec2-Todos los dias -1,585 0,681 -2,33
## Yogurt.rec3-Nunca o raramente -1,270 0,706 -1,80
## AlterGastrica.rec2-Yes -15,686 0,384 -40,87
## Pr(>|t|)
## (Intercept) 0,00018 ***
## SexoM 0,01814 *
## Tipo.de.Escuela2-Privada 0,24016
## IGS.rec2-De 45 a 60 0,36622
## IGS.rec3 -Menos o igual a 45 0,11051
## RefrCola.rec2-Todos los dias 0,46729
## RefrCola.rec3-Mas de 3 veces al dia 0,20633
## Bebidas_energizantes.rec2-Todos los dias 0,20967
## Bebidas_energizantes.recMas de 3 veces al dia 0,02126 *
## Yogurt.rec2-Todos los dias 0,02731 *
## Yogurt.rec3-Nunca o raramente 0,08271 .
## AlterGastrica.rec2-Yes < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1
##
## (Dispersion parameter for quasibinomial family taken to be 0,8205)
##
## Number of Fisher Scoring iterations: 17
confint(modelomulti_ero.logit)
## 2,5 % 97,5 %
## (Intercept) -8,0673 -3,0214
## SexoM 0,2625 2,1346
## Tipo.de.Escuela2-Privada -0,5646 2,3479
## IGS.rec2-De 45 a 60 -1,3448 3,7167
## IGS.rec3 -Menos o igual a 45 -0,3466 4,0101
## RefrCola.rec2-Todos los dias -0,6226 1,3729
## RefrCola.rec3-Mas de 3 veces al dia -0,4766 2,3278
## Bebidas_energizantes.rec2-Todos los dias -0,2109 1,0120
## Bebidas_energizantes.recMas de 3 veces al dia 0,3517 3,2216
## Yogurt.rec2-Todos los dias -2,9192 -0,2511
## Yogurt.rec3-Nunca o raramente -2,6525 0,1132
## AlterGastrica.rec2-Yes -16,4379 -14,9333
reporte <- data.frame(modelomulti_ero.logit$coefficients, exp(modelomulti_ero.logit$coefficients))
colnames(reporte) <- c("coef", "exp_coef")
reporte
## coef exp_coef
## (Intercept) -5,5444 3,909e-03
## SexoM 1,1986 3,315e+00
## Tipo.de.Escuela2-Privada 0,8917 2,439e+00
## IGS.rec2-De 45 a 60 1,1859 3,274e+00
## IGS.rec3 -Menos o igual a 45 1,8317 6,245e+00
## RefrCola.rec2-Todos los dias 0,3752 1,455e+00
## RefrCola.rec3-Mas de 3 veces al dia 0,9256 2,523e+00
## Bebidas_energizantes.rec2-Todos los dias 0,4006 1,493e+00
## Bebidas_energizantes.recMas de 3 veces al dia 1,7867 5,969e+00
## Yogurt.rec2-Todos los dias -1,5852 2,049e-01
## Yogurt.rec3-Nunca o raramente -1,2697 2,809e-01
## AlterGastrica.rec2-Yes -15,6856 1,541e-07