Promedios para niveles de severidad
# load('~/Dropbox/odontologia/maestria licet/octubre
# 2013/datos_licet_13102013.RData')
load("C:/Users/usuario/Dropbox/odontologia/maestria licet/octubre 2013/datos_licet_16102013.RData")
library(survey)
## Attaching package: 'survey'
## The following object(s) are masked from 'package:graphics':
##
## dotchart
options(OutDec = ",")
library(car)
## Loading required package: MASS
## Loading required package: nnet
summary(diseniopost1.sub$variables$erodentina1)
## 0 1 2 3 4 6 7 NA's
## 1086 22 17 3 6 1 1 18
summary(diseniopost1.sub$variables$erodentina2)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 0,00 0,00 0,00 0,09 0,00 7,00 18
tabla1 <- svymean(~Sexo.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
round(tabla1[c(1:2)] * 100, 2)
## Sexo.rec1-F Sexo.rec2-M
## 26,97 73,03
round(confint(tabla1)[c(1:4)] * 100, 2)
## [1] 9,60 55,65 44,35 90,40
tabla2 <- svymean(~Nivel.Educativo.de.la.Madre1.rec, diseniopost1.sub, na.rm = TRUE,
deff = TRUE)
round(tabla2 * 100, 1)
## mean SE
## Nivel.Educativo.de.la.Madre1.rec1-Incomplete Education 18,9000 0,0891
## Nivel.Educativo.de.la.Madre1.rec2-Basic 48,2000 0,0803
## Nivel.Educativo.de.la.Madre1.rec3-Medium 16,7000 0,0731
## Nivel.Educativo.de.la.Madre1.rec4-High 16,2000 0,0769
## DEff
## Nivel.Educativo.de.la.Madre1.rec1-Incomplete Education 2,47
## Nivel.Educativo.de.la.Madre1.rec2-Basic 1,23
## Nivel.Educativo.de.la.Madre1.rec3-Medium 1,84
## Nivel.Educativo.de.la.Madre1.rec4-High 2,09
round(confint(tabla2) * 100, 1)
## 2,5 % 97,5 %
## Nivel.Educativo.de.la.Madre1.rec1-Incomplete Education 1,5 36,4
## Nivel.Educativo.de.la.Madre1.rec2-Basic 32,5 63,9
## Nivel.Educativo.de.la.Madre1.rec3-Medium 2,3 31,0
## Nivel.Educativo.de.la.Madre1.rec4-High 1,1 31,3
tabla3 <- svymean(~Nive.Educativo.de.la.Madre2.rec, diseniopost1.sub, na.rm = TRUE,
deff = TRUE)
round(tabla3[c(1:3)] * 100, 1)
## Nive.Educativo.de.la.Madre2.rec1-ELEMENTARY SCHOOL
## 31,3
## Nive.Educativo.de.la.Madre2.rec2-HIGH SCHOOL
## 40,8
## Nive.Educativo.de.la.Madre2.rec3-COLLEGE-UNIVERSITY
## 27,9
round(confint(tabla3) * 100, 1)
## 2,5 % 97,5 %
## Nive.Educativo.de.la.Madre2.rec1-ELEMENTARY SCHOOL 13,1 49,5
## Nive.Educativo.de.la.Madre2.rec2-HIGH SCHOOL 25,6 56,1
## Nive.Educativo.de.la.Madre2.rec3-COLLEGE-UNIVERSITY 10,9 44,9
tabla4 <- svymean(~Tipo.de.Escuela.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
round(tabla4[c(1:3)] * 100, 1)
## Tipo.de.Escuela.rec1-Public Tipo.de.Escuela.rec2-Private
## 12,9 87,1
## <NA>
## NA
round(confint(tabla4) * 100, 1)
## 2,5 % 97,5 %
## Tipo.de.Escuela.rec1-Public -1,0 26,8
## Tipo.de.Escuela.rec2-Private 73,2 101,0
tabla5 <- svymean(~Nivel.Socieconomico.rec, diseniopost1.sub, na.rm = TRUE,
deff = TRUE)
round(tabla5[c(1:2)] * 100, 1)
## Nivel.Socieconomico.rec1-BAJO Nivel.Socieconomico.rec2-MEDIO
## 45,2 26,8
round(confint(tabla5)[c(1:2)] * 1, 1)
## [1] 0,3 0,1
tabla6 <- svyby(~erodentina2, ~Sexo.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla6[, 2:4] * 1, 2)
## erodentina2 se DEff.erodentina2
## 1-F 1,43 0,18 1,26
## 2-M 2,37 0,47 3,43
round(confint(tabla6) * 1, 2)
## 2,5 % 97,5 %
## 1-F 1,08 1,77
## 2-M 1,44 3,30
tabla7 <- svyby(~erodentina2, ~Nivel.Socieconomico.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla7[, 2:4] * 1, 2)
## erodentina2 se DEff.erodentina2
## 1-BAJO 1,92 0,33 1,69
## 2-MEDIO 1,85 0,38 0,91
## 3-ALTO 2,49 0,54 1,64
round(confint(tabla7) * 1, 2)
## 2,5 % 97,5 %
## 1-BAJO 1,28 2,57
## 2-MEDIO 1,11 2,60
## 3-ALTO 1,42 3,56
tabla8 <- svyby(~erodentina2, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1.sub,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Incomplete Education 2,6 0,7 2,4
## 2-Basic 2,1 0,4 1,3
## 3-Medium 1,7 0,6 1,2
## 4-High 2,2 0,4 1,3
round(confint(tabla8) * 1, 2)
## 2,5 % 97,5 %
## 1-Incomplete Education 1,34 3,94
## 2-Basic 1,35 2,79
## 3-Medium 0,62 2,81
## 4-High 1,37 3,09
tabla8a <- svyby(~erodentina2, ~Nive.Educativo.de.la.Madre2.rec, diseniopost1.sub,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8a[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-ELEMENTARY SCHOOL 2,5 0,6 2,3
## 2-HIGH SCHOOL 2,1 0,4 1,0
## 3-COLLEGE-UNIVERSITY 1,9 0,3 1,2
round(confint(tabla8a) * 1, 2)
## 2,5 % 97,5 %
## 1-ELEMENTARY SCHOOL 1,39 3,62
## 2-HIGH SCHOOL 1,37 2,77
## 3-COLLEGE-UNIVERSITY 1,17 2,53
tabla9 <- svyby(~erodentina2, ~Tipo.de.Escuela.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla9[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Public 2,7 0,4 1,0
## 2-Private 2,0 0,4 2,7
round(confint(tabla9) * 1, 2)
## 2,5 % 97,5 %
## 1-Public 1,95 3,36
## 2-Private 1,32 2,75
tabla10 <- svyby(~erodentina2, ~FrCepDenti.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla10[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-1 vez al dia o menos 1,9 0,4 1,5
## 2- 2 veces al dia 3,3 0,5 1,5
## 3-3 o mas veces al dia 1,6 0,2 1,4
round(confint(tabla10) * 1, 2)
## 2,5 % 97,5 %
## 1-1 vez al dia o menos 1,13 2,60
## 2- 2 veces al dia 2,26 4,28
## 3-3 o mas veces al dia 1,14 2,01
tabla11 <- svyby(~erodentina2, ~UsoDentifrico3.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla11[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Yes 2,1 0,3 3,0
## 2-No 4,7 1,8 0,6
round(confint(tabla11) * 1, 2)
## 2,5 % 97,5 %
## 1-Yes 1,42 2,69
## 2-No 1,05 8,27
tabla12 <- svyby(~erodentina2, ~Consitencia_Cepillo.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla12[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1 2,7 0,7 2,6
## 2 1,4 0,2 1,0
## 3 1,7 0,3 1,1
round(confint(tabla12) * 1, 2)
## 2,5 % 97,5 %
## 1 1,37 4,07
## 2 0,99 1,72
## 3 1,09 2,25
tabla13 <- svyby(~erodentina2, ~IGS.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla13[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Mas de 60 1,0 0,0 NaN
## 2-De 45 a 60 2,8 1,1 0,9
## 3 -Menos o igual a 45 2,0 0,3 3,4
round(confint(tabla13) * 1, 2)
## 2,5 % 97,5 %
## 1-Mas de 60 1,00 1,00
## 2-De 45 a 60 0,65 4,98
## 3 -Menos o igual a 45 1,39 2,66
tabla14 <- svyby(~erodentina2, ~RefrCola.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla14[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Nunca o raramente 2,0 0,5 0,7
## 2-Todos los dias 2,3 0,4 2,0
## 3-Mas de 3 veces al dia 2,0 0,6 2,0
round(confint(tabla14) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 1,09 2,89
## 2-Todos los dias 1,59 2,97
## 3-Mas de 3 veces al dia 0,75 3,19
tabla15 <- svyby(~erodentina2, ~JugFrutas.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla15[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Nunca o raramente 2,1 0,4 0,6
## 2-Todos los dias 2,3 0,5 3,5
## 3-Mas de 3 veces al dia 1,7 0,4 0,8
round(confint(tabla15) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 1,24 2,91
## 2-Todos los dias 1,45 3,22
## 3-Mas de 3 veces al dia 0,95 2,48
tabla16 <- svyby(~erodentina2, ~Bebidas_energizantes_b, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla16[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-No diario 1,9 0,4 3,8
## 2-Diario 3,4 0,8 1,2
round(confint(tabla16) * 1, 2)
## 2,5 % 97,5 %
## 1-No diario 1,17 2,57
## 2-Diario 1,81 4,98
tabla17 <- svyby(~erodentina2, ~Yogurt_b, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Diario 2,1 0,6 1,2
## 2-No diario 2,2 0,4 3,8
round(confint(tabla17) * 1, 2)
## 2,5 % 97,5 %
## 1-Diario 0,85 3,38
## 2-No diario 1,34 2,98
tabla18 <- svyby(~erodentina2, ~Bruxismo.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla18[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-No 2,2 0,3 2,2
## 2-Yes 2,0 0,3 1,3
round(confint(tabla18) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,47 2,84
## 2-Yes 1,37 2,54
tabla19 <- svyby(~erodentina2, ~bedeportediario.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla19[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Water 2,2 0,3 2,1
## 2-Gatorade 1,3 0,2 0,7
## 3-other 2,0 0,0 NaN
round(confint(tabla19) * 1, 2)
## 2,5 % 97,5 %
## 1-Water 1,55 2,83
## 2-Gatorade 0,80 1,73
## 3-other 2,00 2,00
tabla21 <- svyby(~erodentina2, ~BuchTragar.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla21[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-No 2,2 0,3 2,7
## 2-Yes 1,4 0,3 0,5
round(confint(tabla21) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,50 2,81
## 2-Yes 0,76 2,03
tabla22 <- svyby(~erodentina2, ~FormBeber.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla22[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Con sorbete 1,4 0,5 0,9
## 2-Por el pico 2,7 0,8 0,8
## 3-Con vaso 2,1 0,3 2,6
round(confint(tabla22) * 1, 2)
## 2,5 % 97,5 %
## 1-Con sorbete 0,50 2,40
## 2-Por el pico 1,21 4,27
## 3-Con vaso 1,44 2,78
tabla23 <- svyby(~erodentina2, ~Natac2vec.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla23[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-No 2,1 0,4 2,5
## 2-Yes 1,4 0,2 1,2
round(confint(tabla23) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,28 2,87
## 2-Yes 1,09 1,75
tabla24 <- svyby(~erodentina2, ~MedResp.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla24[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-No 2,2 0,3 2,5
## 2-Yes 1,5 0,1 0,2
round(confint(tabla24) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,52 2,88
## 2-Yes 1,29 1,63
tabla25 <- svyby(~erodentina2, ~AlterGastrica.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla25[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-No 2,1 0,3 2,6
round(confint(tabla25) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,5 2,76