# 14 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.logit2 <- svyglm(erodentina.rec ~ Sexo + Nivel.Educativo.de.la.Madre1 +
Tipo.de.Escuela + IGS.rec + RefrCola.rec + Bebidas_energizantes.rec + Yogurt.rec +
AlterGastrica.rec + bedeportediario.rec, design = diseniopost1, family = quasibinomial())
summary(modelomulti_ero.logit2)
##
## 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 + bedeportediario.rec, design = diseniopost1,
## family = quasibinomial())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error
## (Intercept) -4,466 1,645
## SexoM 1,081 0,459
## Nivel.Educativo.de.la.Madre12-Elementary School -1,104 0,586
## Nivel.Educativo.de.la.Madre13-High School -1,289 0,695
## Nivel.Educativo.de.la.Madre14-College-University -1,147 0,876
## Tipo.de.Escuela2-Privada 0,802 0,635
## IGS.rec2-De 45 a 60 1,232 1,445
## IGS.rec3 -Menos o igual a 45 1,941 1,234
## RefrCola.rec2-Todos los dias 0,460 0,855
## RefrCola.rec3-Mas de 3 veces al dia 1,070 1,023
## Bebidas_energizantes.rec2-Todos los dias 0,280 0,503
## Bebidas_energizantes.recMas de 3 veces al dia 1,489 0,966
## Yogurt.rec2-Todos los dias -1,615 0,793
## Yogurt.rec3-Nunca o raramente -1,596 0,766
## AlterGastrica.rec2-Yes -15,633 0,449
## bedeportediario.rec2-Gatorade 0,496 0,851
## bedeportediario.rec3-other -1,394 1,218
## t value Pr(>|t|)
## (Intercept) -2,71 0,014 *
## SexoM 2,35 0,030 *
## Nivel.Educativo.de.la.Madre12-Elementary School -1,88 0,076 .
## Nivel.Educativo.de.la.Madre13-High School -1,86 0,080 .
## Nivel.Educativo.de.la.Madre14-College-University -1,31 0,207
## Tipo.de.Escuela2-Privada 1,26 0,223
## IGS.rec2-De 45 a 60 0,85 0,405
## IGS.rec3 -Menos o igual a 45 1,57 0,133
## RefrCola.rec2-Todos los dias 0,54 0,597
## RefrCola.rec3-Mas de 3 veces al dia 1,05 0,310
## Bebidas_energizantes.rec2-Todos los dias 0,56 0,584
## Bebidas_energizantes.recMas de 3 veces al dia 1,54 0,141
## Yogurt.rec2-Todos los dias -2,04 0,057 .
## Yogurt.rec3-Nunca o raramente -2,08 0,052 .
## AlterGastrica.rec2-Yes -34,79 <2e-16 ***
## bedeportediario.rec2-Gatorade 0,58 0,567
## bedeportediario.rec3-other -1,14 0,267
## ---
## Signif. codes: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1
##
## (Dispersion parameter for quasibinomial family taken to be 0,7175)
##
## Number of Fisher Scoring iterations: 17
confint(modelomulti_ero.logit2)
## 2,5 % 97,5 %
## (Intercept) -7,6899 -1,24149
## SexoM 0,1812 1,98015
## Nivel.Educativo.de.la.Madre12-Elementary School -2,2530 0,04577
## Nivel.Educativo.de.la.Madre13-High School -2,6509 0,07292
## Nivel.Educativo.de.la.Madre14-College-University -2,8649 0,57062
## Tipo.de.Escuela2-Privada -0,4433 2,04721
## IGS.rec2-De 45 a 60 -1,6001 4,06429
## IGS.rec3 -Menos o igual a 45 -0,4778 4,36006
## RefrCola.rec2-Todos los dias -1,2167 2,13634
## RefrCola.rec3-Mas de 3 veces al dia -0,9364 3,07564
## Bebidas_energizantes.rec2-Todos los dias -0,7054 1,26589
## Bebidas_energizantes.recMas de 3 veces al dia -0,4048 3,38253
## Yogurt.rec2-Todos los dias -3,1698 -0,06065
## Yogurt.rec3-Nunca o raramente -3,0979 -0,09438
## AlterGastrica.rec2-Yes -16,5139 -14,75230
## bedeportediario.rec2-Gatorade -1,1724 2,16527
## bedeportediario.rec3-other -3,7803 0,99327
reporte <- data.frame(modelomulti_ero.logit2$coefficients, exp(modelomulti_ero.logit2$coefficients))
colnames(reporte) <- c("coef", "exp_coef")
reporte
## coef exp_coef
## (Intercept) -4,4657 1,150e-02
## SexoM 1,0807 2,947e+00
## Nivel.Educativo.de.la.Madre12-Elementary School -1,1036 3,317e-01
## Nivel.Educativo.de.la.Madre13-High School -1,2890 2,755e-01
## Nivel.Educativo.de.la.Madre14-College-University -1,1471 3,175e-01
## Tipo.de.Escuela2-Privada 0,8020 2,230e+00
## IGS.rec2-De 45 a 60 1,2321 3,428e+00
## IGS.rec3 -Menos o igual a 45 1,9411 6,967e+00
## RefrCola.rec2-Todos los dias 0,4598 1,584e+00
## RefrCola.rec3-Mas de 3 veces al dia 1,0696 2,914e+00
## Bebidas_energizantes.rec2-Todos los dias 0,2802 1,323e+00
## Bebidas_energizantes.recMas de 3 veces al dia 1,4889 4,432e+00
## Yogurt.rec2-Todos los dias -1,6152 1,988e-01
## Yogurt.rec3-Nunca o raramente -1,5961 2,027e-01
## AlterGastrica.rec2-Yes -15,6331 1,624e-07
## bedeportediario.rec2-Gatorade 0,4964 1,643e+00
## bedeportediario.rec3-other -1,3935 2,482e-01
modelomulti_ero.logit2
## Stratified 1 - level Cluster Sampling design (with replacement)
## With (44) clusters.
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## 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 + bedeportediario.rec, design = diseniopost1,
## family = quasibinomial())
##
## Coefficients:
## (Intercept)
## -4,466
## SexoM
## 1,081
## Nivel.Educativo.de.la.Madre12-Elementary School
## -1,104
## Nivel.Educativo.de.la.Madre13-High School
## -1,289
## Nivel.Educativo.de.la.Madre14-College-University
## -1,147
## Tipo.de.Escuela2-Privada
## 0,802
## IGS.rec2-De 45 a 60
## 1,232
## IGS.rec3 -Menos o igual a 45
## 1,941
## RefrCola.rec2-Todos los dias
## 0,460
## RefrCola.rec3-Mas de 3 veces al dia
## 1,070
## Bebidas_energizantes.rec2-Todos los dias
## 0,280
## Bebidas_energizantes.recMas de 3 veces al dia
## 1,489
## Yogurt.rec2-Todos los dias
## -1,615
## Yogurt.rec3-Nunca o raramente
## -1,596
## AlterGastrica.rec2-Yes
## -15,633
## bedeportediario.rec2-Gatorade
## 0,496
## bedeportediario.rec3-other
## -1,394
##
## Degrees of Freedom: 876 Total (i.e. Null); 18 Residual
## (277 observations deleted due to missingness)
## Null Deviance: 321
## Residual Deviance: 274 AIC: NA