# 09 de setiembre de 2013
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
modelomulti_ero.logit1 <- svyglm(erodentina.rec ~ Sexo + Nivel.Educativo.de.la.Madre1 +
Tipo.de.Escuela + IGS.rec + RefrCola.rec + Bebidas_energizantes.rec + Yogurt.rec +
AlterGastrica.rec, design = diseniopost1, family = quasibinomial())
summary(modelomulti_ero.logit1)
##
## Call:
## svyglm(formula = erodentina.rec ~ Sexo + Nivel.Educativo.de.la.Madre1 +
## 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
## (Intercept) -4,408 1,389
## SexoM 1,282 0,495
## Nivel.Educativo.de.la.Madre12-Elementary School -1,302 0,492
## Nivel.Educativo.de.la.Madre13-High School -1,270 0,557
## Nivel.Educativo.de.la.Madre14-College-University -1,367 0,727
## Tipo.de.Escuela2-Privada 0,802 0,628
## IGS.rec2-De 45 a 60 1,368 1,312
## IGS.rec3 -Menos o igual a 45 1,984 1,177
## RefrCola.rec2-Todos los dias 0,317 0,546
## RefrCola.rec3-Mas de 3 veces al dia 0,913 0,794
## Bebidas_energizantes.rec2-Todos los dias 0,331 0,327
## Bebidas_energizantes.recMas de 3 veces al dia 1,805 0,803
## Yogurt.rec2-Todos los dias -1,627 0,771
## Yogurt.rec3-Nunca o raramente -1,476 0,788
## AlterGastrica.rec2-Yes -15,572 0,395
## t value Pr(>|t|)
## (Intercept) -3,17 0,004 **
## SexoM 2,59 0,016 *
## Nivel.Educativo.de.la.Madre12-Elementary School -2,65 0,014 *
## Nivel.Educativo.de.la.Madre13-High School -2,28 0,031 *
## Nivel.Educativo.de.la.Madre14-College-University -1,88 0,072 .
## Tipo.de.Escuela2-Privada 1,28 0,213
## IGS.rec2-De 45 a 60 1,04 0,307
## IGS.rec3 -Menos o igual a 45 1,68 0,105
## RefrCola.rec2-Todos los dias 0,58 0,567
## RefrCola.rec3-Mas de 3 veces al dia 1,15 0,261
## Bebidas_energizantes.rec2-Todos los dias 1,01 0,322
## Bebidas_energizantes.recMas de 3 veces al dia 2,25 0,034 *
## Yogurt.rec2-Todos los dias -2,11 0,045 *
## Yogurt.rec3-Nunca o raramente -1,87 0,073 .
## AlterGastrica.rec2-Yes -39,41 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1
##
## (Dispersion parameter for quasibinomial family taken to be 0,7978)
##
## Number of Fisher Scoring iterations: 17
confint(modelomulti_ero.logit1)
## 2,5 % 97,5 %
## (Intercept) -7,1307 -1,68486
## SexoM 0,3127 2,25120
## Nivel.Educativo.de.la.Madre12-Elementary School -2,2664 -0,33725
## Nivel.Educativo.de.la.Madre13-High School -2,3617 -0,17846
## Nivel.Educativo.de.la.Madre14-College-University -2,7923 0,05863
## Tipo.de.Escuela2-Privada -0,4282 2,03303
## IGS.rec2-De 45 a 60 -1,2030 3,93896
## IGS.rec3 -Menos o igual a 45 -0,3242 4,29124
## RefrCola.rec2-Todos los dias -0,7530 1,38622
## RefrCola.rec3-Mas de 3 veces al dia -0,6439 2,47034
## Bebidas_energizantes.rec2-Todos los dias -0,3108 0,97285
## Bebidas_energizantes.recMas de 3 veces al dia 0,2317 3,37810
## Yogurt.rec2-Todos los dias -3,1391 -0,11530
## Yogurt.rec3-Nunca o raramente -3,0209 0,06959
## AlterGastrica.rec2-Yes -16,3461 -14,79715
reporte <- data.frame(modelomulti_ero.logit1$coefficients, exp(modelomulti_ero.logit1$coefficients))
colnames(reporte) <- c("coef", "exp_coef")
reporte
## coef exp_coef
## (Intercept) -4,4078 1,218e-02
## SexoM 1,2820 3,604e+00
## Nivel.Educativo.de.la.Madre12-Elementary School -1,3018 2,720e-01
## Nivel.Educativo.de.la.Madre13-High School -1,2701 2,808e-01
## Nivel.Educativo.de.la.Madre14-College-University -1,3668 2,549e-01
## Tipo.de.Escuela2-Privada 0,8024 2,231e+00
## IGS.rec2-De 45 a 60 1,3680 3,928e+00
## IGS.rec3 -Menos o igual a 45 1,9835 7,268e+00
## RefrCola.rec2-Todos los dias 0,3166 1,372e+00
## RefrCola.rec3-Mas de 3 veces al dia 0,9132 2,492e+00
## Bebidas_energizantes.rec2-Todos los dias 0,3310 1,392e+00
## Bebidas_energizantes.recMas de 3 veces al dia 1,8049 6,079e+00
## Yogurt.rec2-Todos los dias -1,6272 1,965e-01
## Yogurt.rec3-Nunca o raramente -1,4757 2,286e-01
## AlterGastrica.rec2-Yes -15,5716 1,727e-07