Me falta aun poder hacer la prevalencia para tabla 2 (no me doy cuenta sobre que hacerla)
load("~/Dropbox/odontologia/maestria licet/diciembre2013/datos_licet_23122013.RData")
library(survey)
## Attaching package: 'survey'
## The following object(s) are masked from 'package:graphics':
##
## dotchart
Prevalencia Tabla 1
tabla1 <- svyby(~prebewe.rec, ~Bebidas_energizantes.rec1, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
tabla1
## Bebidas_energizantes.rec1 prebewe.rec se
## 1-Menos de 3 veces 1-Menos de 3 veces 0.5329 0.03407
## Mas de 3 veces al dia Mas de 3 veces al dia 0.2158 0.12168
## DEff.prebewe.rec
## 1-Menos de 3 veces 5.257
## Mas de 3 veces al dia 1.086
tabla2 <- svyby(~prebewe.rec, ~RefrCola.rec1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
tabla2
## RefrCola.rec1 prebewe.rec se DEff.prebewe.rec
## 1-Menos de 3 veces 1-Menos de 3 veces 0.5184 0.03586 4.537
## 2-Mas de 3 veces 2-Mas de 3 veces 0.5609 0.05083 3.036
tabla3 <- svyby(~prebewe.rec, ~Yogurt.rec1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
tabla3
## Yogurt.rec1 prebewe.rec se DEff.prebewe.rec
## 1-Menos de 3 veces 1-Menos de 3 veces 0.5276 0.03762 6.305
## 2-Mas de 3 veces 2-Mas de 3 veces 0.6365 0.08815 1.513
Extension Tabla 1
tabla4 <- svyby(~Bewes.ajust, ~Bebidas_energizantes.rec1, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
tabla4
## Bebidas_energizantes.rec1 Bewes.ajust se
## 1-Menos de 3 veces 1-Menos de 3 veces 5.895 0.7579
## Mas de 3 veces al dia Mas de 3 veces al dia 2.580 1.5209
## DEff.Bewes.ajust
## 1-Menos de 3 veces 13.006
## Mas de 3 veces al dia 1.102
tabla5 <- svyby(~Bewes.ajust, ~RefrCola.rec1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
tabla5
## RefrCola.rec1 Bewes.ajust se DEff.Bewes.ajust
## 1-Menos de 3 veces 1-Menos de 3 veces 5.591 0.5931 6.654
## 2-Mas de 3 veces 2-Mas de 3 veces 6.742 1.4179 9.837
tabla6 <- svyby(~Bewes.ajust, ~Yogurt.rec1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
tabla6
## Yogurt.rec1 Bewes.ajust se DEff.Bewes.ajust
## 1-Menos de 3 veces 1-Menos de 3 veces 5.853 0.7945 13.955
## 2-Mas de 3 veces 2-Mas de 3 veces 6.782 1.2557 1.885
Extension Tabla 2
tabla10 <- svyby(~Severas, ~Bebidas_energizantes.rec1, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
tabla10
## Bebidas_energizantes.rec1 Severas se
## 1-Menos de 3 veces 1-Menos de 3 veces 1.829 0.2983
## Mas de 3 veces al dia Mas de 3 veces al dia 2.362 0.3565
## DEff.Severas
## 1-Menos de 3 veces 2.9499
## Mas de 3 veces al dia 0.6754
tabla11 <- svyby(~Severas, ~RefrCola.rec1, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
## Error: objeto 'RefrCola.rec1' no encontrado
tabla11
## RefrColalight.rec1 Severas se DEff.Severas
## 1-Menos de 3 veces 1-Menos de 3 veces 1.937 0.2646 2.589
## 2-Mas de 3 veces 2-Mas de 3 veces 1.000 0.0000 NaN
tabla12 <- svyby(~Severas, ~Yogurt.rec1, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
tabla12
## Yogurt.rec1 Severas se DEff.Severas
## 1-Menos de 3 veces 1-Menos de 3 veces 1.937 0.3145 3.5790
## 2-Mas de 3 veces 2-Mas de 3 veces 1.675 0.5423 0.5769
SEVERIDAD calculado sobre lo que corresponde
Severidad6 <- svyby(~Severas, ~Sexo.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad6[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-F 1.35 0.17 1.59
## 2-M 2.06 0.37 3.40
round(confint(Severidad6) * 1, 2)
## 2.5 % 97.5 %
## 1-F 1.01 1.68
## 2-M 1.34 2.79
Severidad7 <- svyby(~Severas, ~Nivel.Socieconomico.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad7[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-BAJO 1.64 0.22 1.05
## 2-MEDIO 1.72 0.31 1.29
## 3-ALTO 2.23 0.51 2.46
round(confint(Severidad7) * 1, 2)
## 2.5 % 97.5 %
## 1-BAJO 1.22 2.07
## 2-MEDIO 1.10 2.33
## 3-ALTO 1.22 3.23
Severidad8 <- svyby(~Severas, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1.sub,
svymean, na.rm = TRUE, deff = TRUE)
round(Severidad8[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Incomplete Education 2.04 0.39 1.55
## 2-Basic 1.76 0.30 1.59
## 3-Medium 1.71 0.56 1.24
## 4-High 2.23 0.44 1.27
round(confint(Severidad8) * 1, 2)
## 2.5 % 97.5 %
## 1-Incomplete Education 1.28 2.80
## 2-Basic 1.18 2.34
## 3-Medium 0.62 2.81
## 4-High 1.37 3.09
Severidad8a <- svyby(~Severas, ~Nive.Educativo.de.la.Madre2.rec, diseniopost1.sub,
svymean, na.rm = TRUE, deff = TRUE)
round(Severidad8a[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-ELEMENTARY SCHOOL 1.86 0.28 1.81
## 2-HIGH SCHOOL 1.93 0.35 1.22
## 3-COLLEGE-UNIVERSITY 1.85 0.35 1.16
round(confint(Severidad8a) * 1, 2)
## 2.5 % 97.5 %
## 1-ELEMENTARY SCHOOL 1.31 2.41
## 2-HIGH SCHOOL 1.25 2.61
## 3-COLLEGE-UNIVERSITY 1.17 2.53
Severidad9 <- svyby(~Severas, ~Tipo.de.Escuela.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad9[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Public 2.65 0.36 0.99
## 2-Private 1.76 0.27 2.61
round(confint(Severidad9) * 1, 2)
## 2.5 % 97.5 %
## 1-Public 1.95 3.36
## 2-Private 1.22 2.29
Severidad10 <- svyby(~Severas, ~FrCepDenti.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad10[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-1 vez al dia o menos 1.84 0.42 1.42
## 2- 2 veces al dia 2.67 0.40 2.02
## 3-3 o mas veces al dia 1.46 0.20 1.28
round(confint(Severidad10) * 1, 2)
## 2.5 % 97.5 %
## 1-1 vez al dia o menos 1.01 2.67
## 2- 2 veces al dia 1.88 3.46
## 3-3 o mas veces al dia 1.07 1.85
Severidad11 <- svyby(~Severas, ~UsoDentifrico3.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad11[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Yes 1.85 0.25 2.47
## 2-No 3.06 0.74 0.59
round(confint(Severidad11) * 1, 2)
## 2.5 % 97.5 %
## 1-Yes 1.36 2.34
## 2-No 1.62 4.51
Severidad12 <- svyby(~Severas, ~Consitencia_Cepillo.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad12[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Soft 1.35 0.19 0.99
## 2-Medium 1.54 0.22 1.30
## 3-Hard 2.41 0.60 2.50
round(confint(Severidad12) * 1, 2)
## 2.5 % 97.5 %
## 1-Soft 0.99 1.72
## 2-Medium 1.10 1.97
## 3-Hard 1.24 3.58
Severidad13 <- svyby(~Severas, ~IGS.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad13[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Mas de 60 1.00 0.00 NaN
## 2-De 45 a 60 1.90 0.45 0.76
## 3 -Menos o igual a 45 1.89 0.28 2.76
round(confint(Severidad13) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 60 1.00 1.00
## 2-De 45 a 60 1.03 2.78
## 3 -Menos o igual a 45 1.34 2.44
Severidad14 <- svyby(~Severas, ~RefrCola.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad14[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Nunca o raramente 1.98 0.74 0.77
## 2-Todos los dias 2.08 0.22 1.09
## 3-Mas de 3 veces al dia 1.66 0.42 2.27
round(confint(Severidad14) * 1, 2)
## 2.5 % 97.5 %
## 1-Nunca o raramente 0.54 3.43
## 2-Todos los dias 1.65 2.51
## 3-Mas de 3 veces al dia 0.84 2.48
Severidad15 <- svyby(~Severas, ~JugFrutas.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad15[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Nunca o raramente 1.00 0.00 NaN
## 2-Todos los dias 2.09 0.36 3.36
## 3-Mas de 3 veces al dia 1.57 0.28 1.03
round(confint(Severidad15) * 1, 2)
## 2.5 % 97.5 %
## 1-Nunca o raramente 1.00 1.00
## 2-Todos los dias 1.39 2.79
## 3-Mas de 3 veces al dia 1.02 2.12
Severidad16 <- svyby(~Severas, ~Bebidas_energizantes.rec, diseniopost1.sub,
svymean, na.rm = TRUE, deff = TRUE)
round(Severidad16[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Nunca o raramente 1.63 0.31 1.84
## 2-Todos los dias 2.05 0.38 2.12
## Mas de 3 veces al dia 2.36 0.36 0.68
round(confint(Severidad16) * 1, 2)
## 2.5 % 97.5 %
## 1-Nunca o raramente 1.02 2.24
## 2-Todos los dias 1.30 2.79
## Mas de 3 veces al dia 1.66 3.06
Severidad17 <- svyby(~Severas, ~Yogurt.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad17[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Mas de 3 veces al dia 1.68 0.54 0.58
## 2-Todos los dias 1.94 0.39 3.05
## 3-Nunca o raramente 1.92 0.30 2.13
round(confint(Severidad17) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 0.61 2.74
## 2-Todos los dias 1.18 2.71
## 3-Nunca o raramente 1.32 2.51
Severidad18 <- svyby(~Severas, ~Bruxismo.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad18[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 1.97 0.29 2.33
## 2-Yes 1.42 0.20 1.39
round(confint(Severidad18) * 1, 2)
## 2.5 % 97.5 %
## 1-No 1.41 2.53
## 2-Yes 1.02 1.81
Severidad19 <- svyby(~Severas, ~bedeportediario.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad19[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Water 2.00 0.27 2.13
## 2-Gatorade 1.26 0.24 0.74
## 3-other 2.00 0.00 NaN
round(confint(Severidad19) * 1, 2)
## 2.5 % 97.5 %
## 1-Water 1.47 2.53
## 2-Gatorade 0.80 1.73
## 3-other 2.00 2.00
Severidad21 <- svyby(~Severas, ~BuchTragar.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad21[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 2.02 0.28 2.20
## 2-Yes 1.23 0.17 0.47
round(confint(Severidad21) * 1, 2)
## 2.5 % 97.5 %
## 1-No 1.48 2.57
## 2-Yes 0.90 1.56
Severidad22 <- svyby(~Severas, ~FormBeber.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad22[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Con sorbete 1.45 0.48 0.87
## 2-Por el pico 1.90 0.49 1.20
## 3-Con vaso 1.92 0.28 2.36
round(confint(Severidad22) * 1, 2)
## 2.5 % 97.5 %
## 1-Con sorbete 0.50 2.40
## 2-Por el pico 0.94 2.85
## 3-Con vaso 1.37 2.47
Severidad23 <- svyby(~Severas, ~Natac2vec.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad23[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 1.90 0.34 2.56
## 2-Yes 1.36 0.16 1.33
round(confint(Severidad23) * 1, 2)
## 2.5 % 97.5 %
## 1-No 1.25 2.56
## 2-Yes 1.04 1.68
Severidad24 <- svyby(~Severas, ~MedResp.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad24[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 1.94 0.26 2.23
## 2-Yes 1.46 0.09 0.16
round(confint(Severidad24) * 1, 2)
## 2.5 % 97.5 %
## 1-No 1.42 2.46
## 2-Yes 1.29 1.63
Severidad25 <- svyby(~Severas, ~AlterGastrica.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad25[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 1.89 0.24 2.32
round(confint(Severidad25) * 1, 2)
## 2.5 % 97.5 %
## 1-No 1.42 2.37