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
svychisq(~prebewe.rec + Bebidas_energizantes.rec1, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~prebewe.rec + Bebidas_energizantes.rec1, diseniopost1, statistic = "Wald")
## F = 2.688, ndf = 1, ddf = 39, p-value = 0.1091
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
svychisq(~prebewe.rec + RefrCola.rec1, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~prebewe.rec + RefrCola.rec1, diseniopost1, statistic = "Wald")
## F = 1.097, ndf = 1, ddf = 39, p-value = 0.3014
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
svychisq(~prebewe.rec + Yogurt.rec1, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~prebewe.rec + Yogurt.rec1, diseniopost1, statistic = "Wald")
## F = 1.061, ndf = 1, ddf = 39, p-value = 0.3092
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
svyttest(Bewes.ajust ~ Bebidas_energizantes.rec1, diseniopost1, na.rm = TRUE,
deff = TRUE)
##
## Design-based t-test
##
## data: Bewes.ajust ~ Bebidas_energizantes.rec1
## t = -2.099, df = 38, p-value = 0.04256
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -3.314
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
svyttest(Bewes.ajust ~ RefrCola.rec1, diseniopost1, na.rm = TRUE, deff = TRUE)
##
## Design-based t-test
##
## data: Bewes.ajust ~ RefrCola.rec1
## t = 1.117, df = 38, p-value = 0.2709
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 1.151
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
svyttest(Bewes.ajust ~ Yogurt.rec1, diseniopost1, na.rm = TRUE, deff = TRUE)
##
## Design-based t-test
##
## data: Bewes.ajust ~ Yogurt.rec1
## t = 0.566, df = 38, p-value = 0.5747
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.9288
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
svyttest(Severas ~ Bebidas_energizantes.rec1, diseniopost1.sub, na.rm = TRUE,
deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Bebidas_energizantes.rec1
## t = 1.136, df = 18, p-value = 0.2708
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.5322
tabla11 <- svyby(~Severas, ~RefrCola.rec1, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
tabla11
## RefrCola.rec1 Severas se DEff.Severas
## 1-Menos de 3 veces 1-Menos de 3 veces 2.075 0.2255 1.274
## 2-Mas de 3 veces 2-Mas de 3 veces 1.665 0.4184 2.268
svyttest(Severas ~ RefrCola.rec1, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ RefrCola.rec1
## t = -0.9545, df = 18, p-value = 0.3525
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -0.4102
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
svyttest(Severas ~ Yogurt.rec1, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Yogurt.rec1
## t = -0.3463, df = 18, p-value = 0.7331
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -0.2622
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
svyttest(Severas ~ Sexo.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Sexo.rec
## t = 1.512, df = 18, p-value = 0.1478
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.7166
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
svyttest(Severas ~ Nivel.Socieconomico.rec, diseniopost1.sub, na.rm = TRUE,
deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Nivel.Socieconomico.rec
## t = 0.2211, df = 18, p-value = 0.8275
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.07332
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
svyttest(Severas ~ Nivel.Educativo.de.la.Madre1.rec, diseniopost1.sub, na.rm = TRUE,
deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Nivel.Educativo.de.la.Madre1.rec
## t = -0.8409, df = 18, p-value = 0.4114
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -0.2806
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
svyttest(Severas ~ Nive.Educativo.de.la.Madre2.rec, diseniopost1.sub, na.rm = TRUE,
deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Nive.Educativo.de.la.Madre2.rec
## t = 0.2212, df = 18, p-value = 0.8275
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.06634
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
svyttest(Severas ~ Tipo.de.Escuela.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Tipo.de.Escuela.rec
## t = -1.993, df = 18, p-value = 0.06169
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -0.8961
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
svyttest(Severas ~ FrCepDenti.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ FrCepDenti.rec
## t = 1.561, df = 18, p-value = 0.136
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.8271
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
svyttest(Severas ~ UsoDentifrico3.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ UsoDentifrico3.rec
## t = 1.462, df = 18, p-value = 0.1611
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 1.212
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
svyttest(Severas ~ Consitencia_Cepillo.rec, diseniopost1.sub, na.rm = TRUE,
deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Consitencia_Cepillo.rec
## t = 1.201, df = 18, p-value = 0.2453
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.1834
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
svyttest(Severas ~ IGS.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ IGS.rec
## t = 2.026, df = 18, p-value = 0.0578
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.9036
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
svyttest(Severas ~ JugFrutas.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ JugFrutas.rec
## t = 3.051, df = 18, p-value = 0.006872
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 1.088
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
svyttest(Severas ~ Bruxismo.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ Bruxismo.rec
## t = -1.512, df = 18, p-value = 0.1478
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -0.5555
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
svyttest(Severas ~ bedeportediario.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ bedeportediario.rec
## t = -1.966, df = 18, p-value = 0.06488
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -0.7369
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
svyttest(Severas ~ BuchTragar.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ BuchTragar.rec
## t = -2.576, df = 18, p-value = 0.01903
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -0.7927
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
svyttest(Severas ~ FormBeber.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ FormBeber.rec
## t = 0.6727, df = 18, p-value = 0.5097
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## 0.4504
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
svyttest(Severas ~ Natac2vec.rec, , diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Error: argument "design" is missing, with no default
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
svyttest(Severas ~ MedResp.rec, diseniopost1.sub, na.rm = TRUE, deff = TRUE)
## Warning: observations with zero weight not used for calculating dispersion
## Warning: observations with zero weight not used for calculating dispersion
##
## Design-based t-test
##
## data: Severas ~ MedResp.rec
## t = -1.716, df = 18, p-value = 0.1033
## alternative hypothesis: true difference in mean is not equal to 0
## sample estimates:
## difference in mean
## -0.4821