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)
addmargins(table(diseniopost1$variables$Severas, diseniopost1$variables$Severas.rec))
##
## 0 1 Sum
## 0 1089 0 1089
## 1 0 23 23
## 2 0 16 16
## 3 0 2 2
## 4 0 6 6
## Sum 1089 47 1136
PREVALENCIA DE SEVERIDAD AJUSTADA
Severidad6 <- svyby(~Severas.rec, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad6[, 2:4] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-F 97,66 2,34 0,84
## 2-M 93,62 6,38 1,72
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-BAJO 96,04 3,96 1,20
## 2-MEDIO 95,95 4,05 1,36
## 3-ALTO 94,30 5,70 2,25
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Incomplete Education 87,25 12,75 5,34
## 2-Basic 96,18 3,82 1,17
## 3-Medium 96,67 3,33 1,27
## 4-High 95,28 4,72 2,10
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-ELEMENTARY SCHOOL 94,26 5,74 1,74
## 2-HIGH SCHOOL 96,57 3,43 1,14
## 3-COLLEGE-UNIVERSITY 94,67 5,33 1,75
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Public 97,67 2,33 1,37
## 2-Private 94,86 5,14 1,20
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-1 vez al dia o menos 96,16 3,84 1,41
## 2- 2 veces al dia 95,98 4,02 1,82
## 3-3 o mas veces al dia 94,83 5,17 1,32
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Yes 95,59 4,41 0,97
## 2-No 94,29 5,71 4,53
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Soft 92,04 7,96 2,09
## 2-Medium 96,96 3,04 0,84
## 3-Hard 92,59 7,41 3,98
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Mas de 60 99,19 0,81 0,86
## 2-De 45 a 60 95,54 4,46 2,52
## 3 -Menos o igual a 45 95,26 4,74 1,18
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Nunca o raramente 98,38 1,62 1,20
## 2-Todos los dias 96,08 3,92 1,11
## 3-Mas de 3 veces al dia 93,44 6,56 2,61
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Nunca o raramente 97,63 2,37 2,21
## 2-Todos los dias 96,19 3,81 0,90
## 3-Mas de 3 veces al dia 92,74 7,26 2,43
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Nunca o raramente 96,59 3,41 1,00
## 2-Todos los dias 95,54 4,46 1,38
## Mas de 3 veces al dia 78,42 21,58 12,17
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Mas de 3 veces al dia 89,81 10,19 5,62
## 2-Todos los dias 95,94 4,06 1,04
## 3-Nunca o raramente 95,77 4,23 1,42
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 95,51 4,49 1,06
## 2-Yes 95,86 4,14 1,57
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Water 95,77 4,23 1,35
## 2-Gatorade 90,28 9,72 5,95
## 3-other 98,87 1,13 1,19
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 95,83 4,17 1,18
## 2-Yes 91,40 8,60 2,38
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-Con sorbete 95,51 4,49 2,75
## 2-Por el pico 97,88 2,12 1,04
## 3-Con vaso 95,12 4,88 1,21
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 95,44 4,56 1,25
## 2-Yes 95,09 4,91 1,70
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 95,32 4,68 1,03
## 2-Yes 97,00 3,00 1,69
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] * 100, 2)
## Severas.rec0 Severas.rec1 se.Severas.rec0
## 1-No 95,31 4,69 1,01
## 2-Yes 100,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