# 15 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.logit <- svyglm(erodentina.rec ~ Sexo + Nivel.Educativo.de.la.Madre1 +
Nivel.Socieconomico + 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 + Nivel.Educativo.de.la.Madre1 +
## Nivel.Socieconomico + 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,802 1,636
## SexoM 1,279 0,500
## Nivel.Educativo.de.la.Madre12-Elementary School -1,290 0,439
## Nivel.Educativo.de.la.Madre13-High School -1,212 0,430
## Nivel.Educativo.de.la.Madre14-College-University -1,397 0,738
## Nivel.Socieconomico2-MEDIO 0,332 0,728
## Nivel.Socieconomico3-ALTO 0,596 0,335
## Tipo.de.Escuela2-Privada 0,880 0,719
## IGS.rec2-De 45 a 60 1,347 1,309
## IGS.rec3 -Menos o igual a 45 2,031 1,184
## RefrCola.rec2-Todos los dias 0,322 0,550
## RefrCola.rec3-Mas de 3 veces al dia 0,990 0,817
## Bebidas_energizantes.rec2-Todos los dias 0,320 0,322
## Bebidas_energizantes.recMas de 3 veces al dia 1,660 0,943
## Yogurt.rec2-Todos los dias -1,618 0,797
## Yogurt.rec3-Nunca o raramente -1,388 0,778
## AlterGastrica.rec2-Yes -15,583 0,362
## t value Pr(>|t|)
## (Intercept) -2,93 0,0075 **
## SexoM 2,56 0,0176 *
## Nivel.Educativo.de.la.Madre12-Elementary School -2,94 0,0074 **
## Nivel.Educativo.de.la.Madre13-High School -2,82 0,0097 **
## Nivel.Educativo.de.la.Madre14-College-University -1,89 0,0708 .
## Nivel.Socieconomico2-MEDIO 0,46 0,6527
## Nivel.Socieconomico3-ALTO 1,78 0,0888 .
## Tipo.de.Escuela2-Privada 1,22 0,2336
## IGS.rec2-De 45 a 60 1,03 0,3143
## IGS.rec3 -Menos o igual a 45 1,72 0,0997 .
## RefrCola.rec2-Todos los dias 0,58 0,5645
## RefrCola.rec3-Mas de 3 veces al dia 1,21 0,2377
## Bebidas_energizantes.rec2-Todos los dias 0,99 0,3303
## Bebidas_energizantes.recMas de 3 veces al dia 1,76 0,0916 .
## Yogurt.rec2-Todos los dias -2,03 0,0541 .
## Yogurt.rec3-Nunca o raramente -1,78 0,0875 .
## AlterGastrica.rec2-Yes -43,05 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1
##
## (Dispersion parameter for quasibinomial family taken to be 0,8047)
##
## Number of Fisher Scoring iterations: 17
confint(modelomulti_ero.logit)
## 2,5 % 97,5 %
## (Intercept) -8,00882 -1,59431
## SexoM 0,29852 2,25935
## Nivel.Educativo.de.la.Madre12-Elementary School -2,15018 -0,42922
## Nivel.Educativo.de.la.Madre13-High School -2,05439 -0,36945
## Nivel.Educativo.de.la.Madre14-College-University -2,84323 0,04823
## Nivel.Socieconomico2-MEDIO -1,09541 1,75961
## Nivel.Socieconomico3-ALTO -0,06135 1,25300
## Tipo.de.Escuela2-Privada -0,52984 2,28981
## IGS.rec2-De 45 a 60 -1,21933 3,91328
## IGS.rec3 -Menos o igual a 45 -0,28923 4,35148
## RefrCola.rec2-Todos los dias -0,75687 1,40038
## RefrCola.rec3-Mas de 3 veces al dia -0,61050 2,59015
## Bebidas_energizantes.rec2-Todos los dias -0,31053 0,95031
## Bebidas_energizantes.recMas de 3 veces al dia -0,18817 3,50754
## Yogurt.rec2-Todos los dias -3,17959 -0,05589
## Yogurt.rec3-Nunca o raramente -2,91237 0,13614
## AlterGastrica.rec2-Yes -16,29275 -14,87397
reporte <- data.frame(modelomulti_ero.logit$coefficients, exp(modelomulti_ero.logit$coefficients))
colnames(reporte) <- c("coef", "exp_coef")
reporte
## coef exp_coef
## (Intercept) -4,8016 8,217e-03
## SexoM 1,2789 3,593e+00
## Nivel.Educativo.de.la.Madre12-Elementary School -1,2897 2,754e-01
## Nivel.Educativo.de.la.Madre13-High School -1,2119 2,976e-01
## Nivel.Educativo.de.la.Madre14-College-University -1,3975 2,472e-01
## Nivel.Socieconomico2-MEDIO 0,3321 1,394e+00
## Nivel.Socieconomico3-ALTO 0,5958 1,815e+00
## Tipo.de.Escuela2-Privada 0,8800 2,411e+00
## IGS.rec2-De 45 a 60 1,3470 3,846e+00
## IGS.rec3 -Menos o igual a 45 2,0311 7,623e+00
## RefrCola.rec2-Todos los dias 0,3218 1,380e+00
## RefrCola.rec3-Mas de 3 veces al dia 0,9898 2,691e+00
## Bebidas_energizantes.rec2-Todos los dias 0,3199 1,377e+00
## Bebidas_energizantes.recMas de 3 veces al dia 1,6597 5,258e+00
## Yogurt.rec2-Todos los dias -1,6177 1,983e-01
## Yogurt.rec3-Nunca o raramente -1,3881 2,495e-01
## AlterGastrica.rec2-Yes -15,5834 1,707e-07