Aca se analiza si la prevalencia de severidad ajustada se asocia con las variables de riesgo
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
# load('C:/Users/usuario/Dropbox/odontologia/maestria licet/octubre
# 2013/datos_licet_03112013.RData')
load("~/Dropbox/odontologia/maestria licet/octubre 2013/datos_licet_03112013.RData")
diseniopost1$variables$Severas.rec <- recode(diseniopost1$variables$Severas,
"0=0;1:4=1")
diseniopost1$variables$Severas.rec <- as.factor(diseniopost1$variables$Severas.rec)
table(diseniopost1$variables$Severas, diseniopost1$variables$Severas.rec)
##
## 0 1
## 0 1089 0
## 1 0 23
## 2 0 16
## 3 0 2
## 4 0 6
PREVALENCIA DE SEVERIDAD AJUSTADA
Severidad6 <- svyby(~Severas.rec, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad6[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-F 0,98 0,02 0,01
## 2-M 0,94 0,06 0,02
round(confint(Severidad6) * 1, 2)
## 2,5 % 97,5 %
## 1-F:Severas.rec0 0,96 0,99
## 2-M:Severas.rec0 0,90 0,97
## 1-F:Severas.rec1 0,01 0,04
## 2-M:Severas.rec1 0,03 0,10
Severidad7 <- svyby(~Severas.rec, ~Nivel.Socieconomico.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad7[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-BAJO 0,96 0,04 0,01
## 2-MEDIO 0,96 0,04 0,01
## 3-ALTO 0,94 0,06 0,02
round(confint(Severidad7) * 1, 2)
## 2,5 % 97,5 %
## 1-BAJO:Severas.rec0 0,94 0,98
## 2-MEDIO:Severas.rec0 0,93 0,99
## 3-ALTO:Severas.rec0 0,90 0,99
## 1-BAJO:Severas.rec1 0,02 0,06
## 2-MEDIO:Severas.rec1 0,01 0,07
## 3-ALTO:Severas.rec1 0,01 0,10
Severidad8 <- svyby(~Severas.rec, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Severidad8[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Incomplete Education 0,87 0,13 0,05
## 2-Basic 0,96 0,04 0,01
## 3-Medium 0,97 0,03 0,01
## 4-High 0,95 0,05 0,02
round(confint(Severidad8) * 1, 2)
## 2,5 % 97,5 %
## 1-Incomplete Education:Severas.rec0 0,77 0,98
## 2-Basic:Severas.rec0 0,94 0,98
## 3-Medium:Severas.rec0 0,94 0,99
## 4-High:Severas.rec0 0,91 0,99
## 1-Incomplete Education:Severas.rec1 0,02 0,23
## 2-Basic:Severas.rec1 0,02 0,06
## 3-Medium:Severas.rec1 0,01 0,06
## 4-High:Severas.rec1 0,01 0,09
Severidad8a <- svyby(~Severas.rec, ~Nive.Educativo.de.la.Madre2.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Severidad8a[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-ELEMENTARY SCHOOL 0,94 0,06 0,02
## 2-HIGH SCHOOL 0,97 0,03 0,01
## 3-COLLEGE-UNIVERSITY 0,95 0,05 0,02
round(confint(Severidad8a) * 1, 2)
## 2,5 % 97,5 %
## 1-ELEMENTARY SCHOOL:Severas.rec0 0,91 0,98
## 2-HIGH SCHOOL:Severas.rec0 0,94 0,99
## 3-COLLEGE-UNIVERSITY:Severas.rec0 0,91 0,98
## 1-ELEMENTARY SCHOOL:Severas.rec1 0,02 0,09
## 2-HIGH SCHOOL:Severas.rec1 0,01 0,06
## 3-COLLEGE-UNIVERSITY:Severas.rec1 0,02 0,09
Severidad9 <- svyby(~Severas.rec, ~Tipo.de.Escuela.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad9[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Public 0,98 0,02 0,01
## 2-Private 0,95 0,05 0,01
round(confint(Severidad9) * 1, 2)
## 2,5 % 97,5 %
## 1-Public:Severas.rec0 0,95 1,00
## 2-Private:Severas.rec0 0,93 0,97
## 1-Public:Severas.rec1 0,00 0,05
## 2-Private:Severas.rec1 0,03 0,07
Severidad10 <- svyby(~Severas.rec, ~FrCepDenti.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad10[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-1 vez al dia o menos 0,96 0,04 0,01
## 2- 2 veces al dia 0,96 0,04 0,02
## 3-3 o mas veces al dia 0,95 0,05 0,01
round(confint(Severidad10) * 1, 2)
## 2,5 % 97,5 %
## 1-1 vez al dia o menos:Severas.rec0 0,93 0,99
## 2- 2 veces al dia:Severas.rec0 0,92 1,00
## 3-3 o mas veces al dia:Severas.rec0 0,92 0,97
## 1-1 vez al dia o menos:Severas.rec1 0,01 0,07
## 2- 2 veces al dia:Severas.rec1 0,00 0,08
## 3-3 o mas veces al dia:Severas.rec1 0,03 0,08
Severidad11 <- svyby(~Severas.rec, ~UsoDentifrico3.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad11[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Yes 0,96 0,04 0,01
## 2-No 0,94 0,06 0,05
round(confint(Severidad11) * 1, 2)
## 2,5 % 97,5 %
## 1-Yes:Severas.rec0 0,94 0,97
## 2-No:Severas.rec0 0,85 1,03
## 1-Yes:Severas.rec1 0,03 0,06
## 2-No:Severas.rec1 -0,03 0,15
Severidad12 <- svyby(~Severas.rec, ~Consitencia_Cepillo.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad12[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Soft 0,92 0,08 0,02
## 2-Medium 0,97 0,03 0,01
## 3-Hard 0,93 0,07 0,04
round(confint(Severidad12) * 1, 2)
## 2,5 % 97,5 %
## 1-Soft:Severas.rec0 0,88 0,96
## 2-Medium:Severas.rec0 0,95 0,99
## 3-Hard:Severas.rec0 0,85 1,00
## 1-Soft:Severas.rec1 0,04 0,12
## 2-Medium:Severas.rec1 0,01 0,05
## 3-Hard:Severas.rec1 0,00 0,15
Severidad13 <- svyby(~Severas.rec, ~IGS.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad13[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Mas de 60 0,99 0,01 0,01
## 2-De 45 a 60 0,96 0,04 0,03
## 3 -Menos o igual a 45 0,95 0,05 0,01
round(confint(Severidad13) * 1, 2)
## 2,5 % 97,5 %
## 1-Mas de 60:Severas.rec0 0,98 1,01
## 2-De 45 a 60:Severas.rec0 0,91 1,00
## 3 -Menos o igual a 45:Severas.rec0 0,93 0,98
## 1-Mas de 60:Severas.rec1 -0,01 0,02
## 2-De 45 a 60:Severas.rec1 0,00 0,09
## 3 -Menos o igual a 45:Severas.rec1 0,02 0,07
Severidad14 <- svyby(~Severas.rec, ~RefrCola.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad14[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Nunca o raramente 0,98 0,02 0,01
## 2-Todos los dias 0,96 0,04 0,01
## 3-Mas de 3 veces al dia 0,93 0,07 0,03
round(confint(Severidad14) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente:Severas.rec0 0,96 1,01
## 2-Todos los dias:Severas.rec0 0,94 0,98
## 3-Mas de 3 veces al dia:Severas.rec0 0,88 0,99
## 1-Nunca o raramente:Severas.rec1 -0,01 0,04
## 2-Todos los dias:Severas.rec1 0,02 0,06
## 3-Mas de 3 veces al dia:Severas.rec1 0,01 0,12
Severidad15 <- svyby(~Severas.rec, ~JugFrutas.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad15[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Nunca o raramente 0,98 0,02 0,02
## 2-Todos los dias 0,96 0,04 0,01
## 3-Mas de 3 veces al dia 0,93 0,07 0,02
round(confint(Severidad15) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente:Severas.rec0 0,93 1,02
## 2-Todos los dias:Severas.rec0 0,94 0,98
## 3-Mas de 3 veces al dia:Severas.rec0 0,88 0,98
## 1-Nunca o raramente:Severas.rec1 -0,02 0,07
## 2-Todos los dias:Severas.rec1 0,02 0,06
## 3-Mas de 3 veces al dia:Severas.rec1 0,02 0,12
Severidad16 <- svyby(~Severas.rec, ~Bebidas_energizantes.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Severidad16[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Nunca o raramente 0,97 0,03 0,01
## 2-Todos los dias 0,96 0,04 0,01
## Mas de 3 veces al dia 0,78 0,22 0,12
round(confint(Severidad16) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente:Severas.rec0 0,95 0,99
## 2-Todos los dias:Severas.rec0 0,93 0,98
## Mas de 3 veces al dia:Severas.rec0 0,55 1,02
## 1-Nunca o raramente:Severas.rec1 0,01 0,05
## 2-Todos los dias:Severas.rec1 0,02 0,07
## Mas de 3 veces al dia:Severas.rec1 -0,02 0,45
Severidad17 <- svyby(~Severas.rec, ~Yogurt.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad17[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Mas de 3 veces al dia 0,90 0,10 0,06
## 2-Todos los dias 0,96 0,04 0,01
## 3-Nunca o raramente 0,96 0,04 0,01
round(confint(Severidad17) * 1, 2)
## 2,5 % 97,5 %
## 1-Mas de 3 veces al dia:Severas.rec0 0,79 1,01
## 2-Todos los dias:Severas.rec0 0,94 0,98
## 3-Nunca o raramente:Severas.rec0 0,93 0,99
## 1-Mas de 3 veces al dia:Severas.rec1 -0,01 0,21
## 2-Todos los dias:Severas.rec1 0,02 0,06
## 3-Nunca o raramente:Severas.rec1 0,01 0,07
Severidad18 <- svyby(~Severas.rec, ~Bruxismo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad18[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 0,96 0,04 0,01
## 2-Yes 0,96 0,04 0,02
round(confint(Severidad18) * 1, 2)
## 2,5 % 97,5 %
## 1-No:Severas.rec0 0,93 0,98
## 2-Yes:Severas.rec0 0,93 0,99
## 1-No:Severas.rec1 0,02 0,07
## 2-Yes:Severas.rec1 0,01 0,07
Severidad19 <- svyby(~Severas.rec, ~bedeportediario.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad19[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Water 0,96 0,04 0,01
## 2-Gatorade 0,90 0,10 0,06
## 3-other 0,99 0,01 0,01
round(confint(Severidad19) * 1, 2)
## 2,5 % 97,5 %
## 1-Water:Severas.rec0 0,93 0,98
## 2-Gatorade:Severas.rec0 0,79 1,02
## 3-other:Severas.rec0 0,97 1,01
## 1-Water:Severas.rec1 0,02 0,07
## 2-Gatorade:Severas.rec1 -0,02 0,21
## 3-other:Severas.rec1 -0,01 0,03
Severidad21 <- svyby(~Severas.rec, ~BuchTragar.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad21[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 0,96 0,04 0,01
## 2-Yes 0,91 0,09 0,02
round(confint(Severidad21) * 1, 2)
## 2,5 % 97,5 %
## 1-No:Severas.rec0 0,94 0,98
## 2-Yes:Severas.rec0 0,87 0,96
## 1-No:Severas.rec1 0,02 0,06
## 2-Yes:Severas.rec1 0,04 0,13
Severidad22 <- svyby(~Severas.rec, ~FormBeber.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad22[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Con sorbete 0,96 0,04 0,03
## 2-Por el pico 0,98 0,02 0,01
## 3-Con vaso 0,95 0,05 0,01
round(confint(Severidad22) * 1, 2)
## 2,5 % 97,5 %
## 1-Con sorbete:Severas.rec0 0,90 1,01
## 2-Por el pico:Severas.rec0 0,96 1,00
## 3-Con vaso:Severas.rec0 0,93 0,97
## 1-Con sorbete:Severas.rec1 -0,01 0,10
## 2-Por el pico:Severas.rec1 0,00 0,04
## 3-Con vaso:Severas.rec1 0,03 0,07
Severidad23 <- svyby(~Severas.rec, ~Natac2vec.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad23[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 0,95 0,05 0,01
## 2-Yes 0,95 0,05 0,02
round(confint(Severidad23) * 1, 2)
## 2,5 % 97,5 %
## 1-No:Severas.rec0 0,93 0,98
## 2-Yes:Severas.rec0 0,92 0,98
## 1-No:Severas.rec1 0,02 0,07
## 2-Yes:Severas.rec1 0,02 0,08
Severidad24 <- svyby(~Severas.rec, ~MedResp.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad24[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 0,95 0,05 0,01
## 2-Yes 0,97 0,03 0,02
round(confint(Severidad24) * 1, 2)
## 2,5 % 97,5 %
## 1-No:Severas.rec0 0,93 0,97
## 2-Yes:Severas.rec0 0,94 1,00
## 1-No:Severas.rec1 0,03 0,07
## 2-Yes:Severas.rec1 0,00 0,06
Severidad25 <- svyby(~Severas.rec, ~AlterGastrica.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad25[, 2:4] * 1, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 0,95 0,05 0,01
## 2-Yes 1,00 0,00 0,00
round(confint(Severidad25) * 1, 2)
## 2,5 % 97,5 %
## 1-No:Severas.rec0 0,93 0,97
## 2-Yes:Severas.rec0 1,00 1,00
## 1-No:Severas.rec1 0,03 0,07
## 2-Yes:Severas.rec1 0,00 0,00