load("~/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
FRECUENCIAS
table(diseniopost1$variables$RefrCola.rec)
##
## 1-Nunca o raramente 2-Todos los dias 3-Mas de 3 veces al dia
## 141 720 280
table(diseniopost1$variables$JugFrutas.rec)
##
## 1-Nunca o raramente 2-Todos los dias 3-Mas de 3 veces al dia
## 48 889 204
table(diseniopost1$variables$Bebidas_energizantes.rec)
##
## 1-Nunca o raramente 2-Todos los dias Mas de 3 veces al dia
## 665 434 13
table(diseniopost1$variables$Yogurt.rec)
##
## 1-Mas de 3 veces al dia 2-Todos los dias 3-Nunca o raramente
## 44 784 299
table(diseniopost1$variables$Bruxismo.rec)
##
## 1-No 2-Yes
## 954 164
PREBEWE
tabla14a <- svyby(~prebewe.rec, ~RefrCola.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla14a[, 2:4] * 100, 2)
## prebewe.rec se DEff.prebewe.rec
## 1-Nunca o raramente 64,07 5,73 205,5
## 2-Todos los dias 49,45 3,72 407,3
## 3-Mas de 3 veces al dia 56,09 5,08 303,6
round(confint(tabla14a) * 100, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 52,83 75,31
## 2-Todos los dias 42,16 56,74
## 3-Mas de 3 veces al dia 46,13 66,05
tabla15a <- svyby(~prebewe.rec, ~JugFrutas.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla15a[, 2:4] * 100, 2)
## prebewe.rec se DEff.prebewe.rec
## 1-Nunca o raramente 64,80 7,38 116,2
## 2-Todos los dias 51,26 3,51 450,3
## 3-Mas de 3 veces al dia 57,46 4,65 183,7
round(confint(tabla15a) * 100, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 50,33 79,27
## 2-Todos los dias 44,38 58,14
## 3-Mas de 3 veces al dia 48,34 66,58
tabla16a <- svyby(~prebewe.rec, ~Bebidas_energizantes.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla16a[, 2:4] * 100, 2)
## prebewe.rec se DEff.prebewe.rec
## 1-Nunca o raramente 52,00 3,71 375,6
## 2-Todos los dias 55,08 4,63 384,5
## Mas de 3 veces al dia 21,58 12,17 108,6
round(confint(tabla16a) * 100, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 44,74 59,27
## 2-Todos los dias 46,01 64,15
## Mas de 3 veces al dia -2,27 45,43
tabla17a <- svyby(~prebewe.rec, ~Yogurt.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17a[, 2:4] * 100, 2)
## prebewe.rec se DEff.prebewe.rec
## 1-Mas de 3 veces al dia 63,65 8,81 151,3
## 2-Todos los dias 52,09 3,88 487,3
## 3-Nunca o raramente 54,41 4,23 217,6
round(confint(tabla17a) * 100, 2)
## 2,5 % 97,5 %
## 1-Mas de 3 veces al dia 46,37 80,93
## 2-Todos los dias 44,49 59,70
## 3-Nunca o raramente 46,12 62,70
VALORES P
svychisq(~RefrCola.rec + prebewe.rec, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~RefrCola.rec + prebewe.rec, diseniopost1, statistic = "Wald")
## F = 3,558, ndf = 2, ddf = 39, p-value = 0,03806
svychisq(~JugFrutas.rec + prebewe.rec, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~JugFrutas.rec + prebewe.rec, diseniopost1, statistic = "Wald")
## F = 3,329, ndf = 2, ddf = 39, p-value = 0,04625
svychisq(~Bebidas_energizantes.rec + prebewe.rec, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~Bebidas_energizantes.rec + prebewe.rec, diseniopost1, statistic = "Wald")
## F = 1,565, ndf = 2, ddf = 39, p-value = 0,2218
svychisq(~Yogurt.rec + prebewe.rec, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~Yogurt.rec + prebewe.rec, diseniopost1, statistic = "Wald")
## F = 0,857, ndf = 2, ddf = 39, p-value = 0,4323
ERONDENTINA Y VALORES P
tabla16 <- svyby(~erodentina.rec, ~Bebidas_energizantes.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla16[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Nunca o raramente 96,59 3,41
## 2-Todos los dias 94,93 5,07
## Mas de 3 veces al dia 78,42 21,58
## se.erodentina.recausente
## 1-Nunca o raramente 1,00
## 2-Todos los dias 1,43
## Mas de 3 veces al dia 12,17
round(confint(tabla16) * 100, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente:erodentina.recausente 94,63 98,54
## 2-Todos los dias:erodentina.recausente 92,12 97,74
## Mas de 3 veces al dia:erodentina.recausente 54,57 102,27
## 1-Nunca o raramente:erodentina.recpresente 1,46 5,37
## 2-Todos los dias:erodentina.recpresente 2,26 7,88
## Mas de 3 veces al dia:erodentina.recpresente -2,27 45,43
tabla17 <- svyby(~erodentina.rec, ~Yogurt.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Mas de 3 veces al dia 89,81 10,19
## 2-Todos los dias 95,94 4,06
## 3-Nunca o raramente 94,88 5,12
## se.erodentina.recausente
## 1-Mas de 3 veces al dia 5,62
## 2-Todos los dias 1,04
## 3-Nunca o raramente 1,62
round(confint(tabla17) * 100, 2)
## 2,5 % 97,5 %
## 1-Mas de 3 veces al dia:erodentina.recausente 78,80 100,82
## 2-Todos los dias:erodentina.recausente 93,90 97,97
## 3-Nunca o raramente:erodentina.recausente 91,69 98,06
## 1-Mas de 3 veces al dia:erodentina.recpresente -0,82 21,20
## 2-Todos los dias:erodentina.recpresente 2,03 6,10
## 3-Nunca o raramente:erodentina.recpresente 1,94 8,31
svychisq(~Bebidas_energizantes.rec + erodentina.rec, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~Bebidas_energizantes.rec + erodentina.rec, diseniopost1, statistic = "Wald")
## F = 3,658, ndf = 2, ddf = 39, p-value = 0,03498
svychisq(~Yogurt.rec + erodentina.rec, diseniopost1, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~Yogurt.rec + erodentina.rec, diseniopost1, statistic = "Wald")
## F = 0,544, ndf = 2, ddf = 39, p-value = 0,5848
PROMEDIOS DE LESIONES SEVERAS
tabla16c <- svyby(~erodentina2, ~Bebidas_energizantes.rec, diseniopost1.sub,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla16c[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Nunca o raramente 1,8 0,3 1,7
## 2-Todos los dias 2,3 0,5 2,1
## Mas de 3 veces al dia 2,9 0,9 0,7
round(confint(tabla16c) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 1,15 2,36
## 2-Todos los dias 1,35 3,29
## Mas de 3 veces al dia 1,16 4,65
tabla17c <- svyby(~erodentina2, ~Yogurt.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17c[, 2:4] * 1, 1)
## erodentina2 se DEff.erodentina2
## 1-Mas de 3 veces al dia 1,7 0,5 0,6
## 2-Todos los dias 2,2 0,5 2,8
## 3-Nunca o raramente 2,3 0,4 2,3
round(confint(tabla17c) * 1, 2)
## 2,5 % 97,5 %
## 1-Mas de 3 veces al dia 0,61 2,74
## 2-Todos los dias 1,22 3,12
## 3-Nunca o raramente 1,48 3,02
tabla18c <- svyby(~erodentina2, ~Bruxismo.rec, diseniopost1.sub, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla18c[, 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(tabla18c) * 1, 2)
## 2,5 % 97,5 %
## 1-No 1,47 2,84
## 2-Yes 1,37 2,54
tabla19c <- svyby(~erodentina2, ~bedeportediario.rec, diseniopost1.sub, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla19c[, 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(tabla19c) * 1, 2)
## 2,5 % 97,5 %
## 1-Water 1,55 2,83
## 2-Gatorade 0,80 1,73
## 3-other 2,00 2,00