Tablas de variables explicativas en forma univariada
[424] “categor.rec”
[425] “prebewe.rec”
[426] “GolosinasB”
[427] “erodentina.rec”
[428] “riesgoerosion.rec”
[429] “Frecuencia.cepillado.rec2”
[430] “Frecuencia.cepillado.rec1”
[431] “bewes”
[432] “Sexo.rec”
[433] “Nivel.Socieconomico.rec”
[434] “Nivel.Educativo.de.la.Madre1.rec”
[435] “Tipo.de.Escuela.rec”
[436] “FrCepDenti.rec”
[437] “UsoDentifrico3.rec”
[438] “Consitencia_Cepillo.rec”
[439] “IGS.rec”
[440] “RefrColalight.rec”
[441] “JugFrutas.rec”
[442] “Bebidas_energizantes.rec”
[443] “Yogurt.rec”
[444] “Bruxismo.rec”
[445] “BuchTragar.rec”
[446] “FormBeber.rec”
[447] “AlterGastrica.rec”
[448] “MedResp.rec”
[449] “Natac2vec.rec”
[450] “bedeportediario.rec”
[451] “RefrCola.rec”
[452] “Nive.Educativo.de.la.Madre2.rec”
[453] “filtro”
load("~/Dropbox/odontologia/maestria licet/octubre 2013/datos_licet_13102013.RData")
library(survey)
## Attaching package: 'survey'
##
## The following object(s) are masked from 'package:graphics':
##
## dotchart
options(OutDec = ",")
# subset(dclus1, sch.wide=='Yes' & comp.imp=='Yes')
diseniopost1.sub <- subset(diseniopost1, diseniopost1$variables$filtro == 1)
# str(diseniopost1.sub)
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)] * 100, 1)
## [1] 31,1 9,9
tabla6 <- svyby(~erodentina, ~Sexo, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
tabla6[, 2:4] * 1
## erodentina se DEff.erodentina
## F 1,427 0,1758 1,264
## M 2,368 0,4747 3,432
round(confint(tabla6) * 1, 2)
## 2,5 % 97,5 %
## F 1,08 1,77
## M 1,44 3,30
tabla7 <- svyby(~erodentina, ~Nivel.Socieconomico.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla7[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 1-BAJO 1,9 0,3 1,7
## 2-MEDIO 1,9 0,4 0,9
## 3-ALTO 2,5 0,5 1,6
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(~erodentina, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1.sub,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 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(~erodentina, ~Nive.Educativo.de.la.Madre2.rec, diseniopost1.sub,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8a[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 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(~erodentina, ~Tipo.de.Escuela.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla9[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 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(~erodentina, ~FrCepDenti.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla10[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 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(~erodentina, ~IGS.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla11[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 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(tabla11) * 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
tabla12 <- svyby(~erodentina, ~MedResp.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla12[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 1-No 2,2 0,3 2,5
## 2-Yes 1,5 0,1 0,2
round(confint(tabla12) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,52 2,88
## 2-Yes 1,29 1,63
tabla13 <- svyby(~erodentina, ~AlterGastrica.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla13[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 1-No 2,1 0,3 2,6
round(confint(tabla13) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,5 2,76
tabla14 <- svyby(~erodentina, ~RefrCola.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla14[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 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(~erodentina, ~UsoDentifrico3.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla15[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 1-Yes 2,1 0,3 3,0
## 2-No 4,7 1,8 0,6
round(confint(tabla15) * 1, 2)
## 2,5 % 97,5 %
## 1-Yes 1,42 2,69
## 2-No 1,05 8,27
tabla16 <- svyby(~erodentina, ~Consitencia_Cepillo.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla16[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 1 2,7 0,7 2,6
## 2 1,4 0,2 1,0
## 3 1,7 0,3 1,1
round(confint(tabla16) * 1, 2)
## 2,5 % 97,5 %
## 1 1,37 4,07
## 2 0,99 1,72
## 3 1,09 2,25
tabla17 <- svyby(~erodentina, ~Natac2vec.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 1-No 2,1 0,4 2,5
## 2-Yes 1,4 0,2 1,2
round(confint(tabla17) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,28 2,87
## 2-Yes 1,09 1,75
tabla18 <- svyby(~erodentina, ~BuchTragar.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla18[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 1-No 2,2 0,3 2,7
## 2-Yes 1,4 0,3 0,5
round(confint(tabla18) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,50 2,81
## 2-Yes 0,76 2,03
tabla19 <- svyby(~erodentina, ~FormBeber.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla19[, 2:4] * 1, 1)
## erodentina se DEff.erodentina
## 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(tabla19) * 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