Prevalencias para erodentina
# load('~/Dropbox/odontologia/maestria licet/octubre
# 2013/datos_licet_13102013.RData')
load("C:/Users/usuario/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
summary(diseniopost1$variables$erodentina1)
## 0 1 2 3 4 6 7 NA's
## 1086 22 17 3 6 1 1 18
summary(diseniopost1$variables$erodentina.rec)
## ausente presente NA's
## 1086 50 18
tabla1 <- svymean(~Sexo.rec, diseniopost1, na.rm = TRUE, deff = TRUE)
round(tabla1[c(1:2)] * 100, 2)
## Sexo.rec1-F Sexo.rec2-M
## 48,25 51,75
round(confint(tabla1)[c(1:4)] * 100, 2)
## [1] 48,25 51,75 48,25 51,75
tabla2 <- svymean(~Nivel.Educativo.de.la.Madre1.rec, diseniopost1, na.rm = TRUE,
deff = TRUE)
round(tabla2 * 100, 2)
## mean SE
## Nivel.Educativo.de.la.Madre1.rec1-Incomplete Education 6,7400 0,0137
## Nivel.Educativo.de.la.Madre1.rec2-Basic 53,2600 0,0397
## Nivel.Educativo.de.la.Madre1.rec3-Medium 23,3100 0,0387
## Nivel.Educativo.de.la.Madre1.rec4-High 16,7000 0,0356
## DEff
## Nivel.Educativo.de.la.Madre1.rec1-Incomplete Education 3,44
## Nivel.Educativo.de.la.Madre1.rec2-Basic 7,34
## Nivel.Educativo.de.la.Madre1.rec3-Medium 9,68
## Nivel.Educativo.de.la.Madre1.rec4-High 10,58
round(confint(tabla2) * 100, 2)
## 2,5 % 97,5 %
## Nivel.Educativo.de.la.Madre1.rec1-Incomplete Education 4,06 9,41
## Nivel.Educativo.de.la.Madre1.rec2-Basic 45,48 61,05
## Nivel.Educativo.de.la.Madre1.rec3-Medium 15,73 30,88
## Nivel.Educativo.de.la.Madre1.rec4-High 9,71 23,69
tabla3 <- svymean(~Nive.Educativo.de.la.Madre2.rec, diseniopost1, na.rm = TRUE,
deff = TRUE)
round(tabla3[c(1:3)] * 100, 2)
## Nive.Educativo.de.la.Madre2.rec1-ELEMENTARY SCHOOL
## 24,20
## Nive.Educativo.de.la.Madre2.rec2-HIGH SCHOOL
## 50,57
## Nive.Educativo.de.la.Madre2.rec3-COLLEGE-UNIVERSITY
## 25,23
round(confint(tabla3) * 100, 2)
## 2,5 % 97,5 %
## Nive.Educativo.de.la.Madre2.rec1-ELEMENTARY SCHOOL 18,55 29,84
## Nive.Educativo.de.la.Madre2.rec2-HIGH SCHOOL 44,02 57,12
## Nive.Educativo.de.la.Madre2.rec3-COLLEGE-UNIVERSITY 18,64 31,82
tabla4 <- svymean(~Tipo.de.Escuela.rec, diseniopost1, na.rm = TRUE, deff = TRUE)
round(tabla4[c(1:3)] * 100, 2)
## Tipo.de.Escuela.rec1-Public Tipo.de.Escuela.rec2-Private
## 26,66 73,34
## <NA>
## NA
round(confint(tabla4) * 100, 2)
## 2,5 % 97,5 %
## Tipo.de.Escuela.rec1-Public 26,66 26,66
## Tipo.de.Escuela.rec2-Private 73,34 73,34
tabla5 <- svymean(~Nivel.Socieconomico.rec, diseniopost1, na.rm = TRUE, deff = TRUE)
round(tabla5[c(1:2)] * 100, 2)
## Nivel.Socieconomico.rec1-BAJO Nivel.Socieconomico.rec2-MEDIO
## 48,64 31,35
round(confint(tabla5)[c(1:2)] * 100, 2)
## [1] 43,54 25,00
tabla6 <- svyby(~erodentina.rec, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla6[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente se.erodentina.recausente
## 1-F 97,41 2,59 0,87
## 2-M 93,38 6,62 1,78
round(confint(tabla6) * 100, 2)
## 2,5 % 97,5 %
## 1-F:erodentina.recausente 95,71 99,11
## 2-M:erodentina.recausente 89,90 96,87
## 1-F:erodentina.recpresente 0,89 4,29
## 2-M:erodentina.recpresente 3,13 10,10
tabla7 <- svyby(~erodentina.rec, ~Nivel.Socieconomico.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla7[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-BAJO 95,79 4,21
## 2-MEDIO 95,95 4,05
## 3-ALTO 93,70 6,30
## se.erodentina.recausente
## 1-BAJO 1,20
## 2-MEDIO 1,36
## 3-ALTO 2,40
round(confint(tabla7) * 100, 2)
## 2,5 % 97,5 %
## 1-BAJO:erodentina.recausente 93,44 98,13
## 2-MEDIO:erodentina.recausente 93,29 98,60
## 3-ALTO:erodentina.recausente 89,00 98,40
## 1-BAJO:erodentina.recpresente 1,87 6,56
## 2-MEDIO:erodentina.recpresente 1,40 6,71
## 3-ALTO:erodentina.recpresente 1,60 11,00
tabla8 <- svyby(~erodentina.rec, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Incomplete Education 87,25 12,75
## 2-Basic 95,90 4,10
## 3-Medium 96,67 3,33
## 4-High 95,28 4,72
## se.erodentina.recausente
## 1-Incomplete Education 5,34
## 2-Basic 1,20
## 3-Medium 1,27
## 4-High 2,10
round(confint(tabla8) * 100, 2)
## 2,5 % 97,5 %
## 1-Incomplete Education:erodentina.recausente 76,78 97,72
## 2-Basic:erodentina.recausente 93,54 98,26
## 3-Medium:erodentina.recausente 94,18 99,16
## 4-High:erodentina.recausente 91,16 99,40
## 1-Incomplete Education:erodentina.recpresente 2,28 23,22
## 2-Basic:erodentina.recpresente 1,74 6,46
## 3-Medium:erodentina.recpresente 0,84 5,82
## 4-High:erodentina.recpresente 0,60 8,84
tabla8a <- svyby(~erodentina.rec, ~Nive.Educativo.de.la.Madre2.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8a[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-ELEMENTARY SCHOOL 94,15 5,85
## 2-HIGH SCHOOL 96,32 3,68
## 3-COLLEGE-UNIVERSITY 94,67 5,33
## se.erodentina.recausente
## 1-ELEMENTARY SCHOOL 1,74
## 2-HIGH SCHOOL 1,19
## 3-COLLEGE-UNIVERSITY 1,75
round(confint(tabla8a) * 1, 2)
## 2,5 % 97,5 %
## 1-ELEMENTARY SCHOOL:erodentina.recausente 0,91 0,98
## 2-HIGH SCHOOL:erodentina.recausente 0,94 0,99
## 3-COLLEGE-UNIVERSITY:erodentina.recausente 0,91 0,98
## 1-ELEMENTARY SCHOOL:erodentina.recpresente 0,02 0,09
## 2-HIGH SCHOOL:erodentina.recpresente 0,01 0,06
## 3-COLLEGE-UNIVERSITY:erodentina.recpresente 0,02 0,09
tabla9 <- svyby(~erodentina.rec, ~Tipo.de.Escuela.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla9[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Public 97,67 2,33
## 2-Private 94,53 5,47
## se.erodentina.recausente
## 1-Public 1,37
## 2-Private 1,22
round(confint(tabla9) * 100, 2)
## 2,5 % 97,5 %
## 1-Public:erodentina.recausente 94,99 100,34
## 2-Private:erodentina.recausente 92,13 96,93
## 1-Public:erodentina.recpresente -0,34 5,01
## 2-Private:erodentina.recpresente 3,07 7,87
tabla10 <- svyby(~erodentina.rec, ~FrCepDenti.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla10[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-1 vez al dia o menos 95,48 4,52
## 2- 2 veces al dia 95,90 4,10
## 3-3 o mas veces al dia 94,62 5,38
## se.erodentina.recausente
## 1-1 vez al dia o menos 1,78
## 2- 2 veces al dia 1,82
## 3-3 o mas veces al dia 1,32
round(confint(tabla10) * 100, 2)
## 2,5 % 97,5 %
## 1-1 vez al dia o menos:erodentina.recausente 92,00 98,97
## 2- 2 veces al dia:erodentina.recausente 92,32 99,47
## 3-3 o mas veces al dia:erodentina.recausente 92,05 97,20
## 1-1 vez al dia o menos:erodentina.recpresente 1,03 8,00
## 2- 2 veces al dia:erodentina.recpresente 0,53 7,68
## 3-3 o mas veces al dia:erodentina.recpresente 2,80 7,95
tabla11 <- svyby(~erodentina.rec, ~UsoDentifrico3.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla11[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Yes 95,34 4,66
## 2-No 94,29 5,71
## se.erodentina.recausente
## 1-Yes 0,99
## 2-No 4,53
round(confint(tabla11) * 100, 2)
## 2,5 % 97,5 %
## 1-Yes:erodentina.recausente 93,39 97,29
## 2-No:erodentina.recausente 85,42 103,15
## 1-Yes:erodentina.recpresente 2,71 6,61
## 2-No:erodentina.recpresente -3,15 14,58
tabla12 <- svyby(~erodentina.rec, ~Consitencia_Cepillo.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla12[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Soft 92,04 7,96
## 2-Medium 96,90 3,10
## 3-Hard 92,59 7,41
## se.erodentina.recausente
## 1-Soft 2,09
## 2-Medium 0,84
## 3-Hard 3,98
round(confint(tabla12) * 100, 2)
## 2,5 % 97,5 %
## 1-Soft:erodentina.recausente 87,93 96,14
## 2-Medium:erodentina.recausente 95,26 98,54
## 3-Hard:erodentina.recausente 84,80 100,39
## 1-Soft:erodentina.recpresente 3,86 12,07
## 2-Medium:erodentina.recpresente 1,46 4,74
## 3-Hard:erodentina.recpresente -0,39 15,20
tabla13 <- svyby(~erodentina.rec, ~IGS.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla13[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Mas de 60 99,19 0,81
## 2-De 45 a 60 95,54 4,46
## 3 -Menos o igual a 45 94,96 5,04
## se.erodentina.recausente
## 1-Mas de 60 0,86
## 2-De 45 a 60 2,52
## 3 -Menos o igual a 45 1,18
round(confint(tabla13) * 100, 2)
## 2,5 % 97,5 %
## 1-Mas de 60:erodentina.recausente 97,51 100,86
## 2-De 45 a 60:erodentina.recausente 90,61 100,47
## 3 -Menos o igual a 45:erodentina.recausente 92,65 97,26
## 1-Mas de 60:erodentina.recpresente -0,86 2,49
## 2-De 45 a 60:erodentina.recpresente -0,47 9,39
## 3 -Menos o igual a 45:erodentina.recpresente 2,74 7,35
tabla14 <- svyby(~erodentina.rec, ~RefrCola.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla14[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Nunca o raramente 97,40 2,60
## 2-Todos los dias 95,93 4,07
## 3-Mas de 3 veces al dia 93,33 6,67
## se.erodentina.recausente
## 1-Nunca o raramente 1,43
## 2-Todos los dias 1,12
## 3-Mas de 3 veces al dia 2,60
round(confint(tabla14) * 100, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente:erodentina.recausente 94,60 100,20
## 2-Todos los dias:erodentina.recausente 93,75 98,12
## 3-Mas de 3 veces al dia:erodentina.recausente 88,23 98,42
## 1-Nunca o raramente:erodentina.recpresente -0,20 5,40
## 2-Todos los dias:erodentina.recpresente 1,88 6,25
## 3-Mas de 3 veces al dia:erodentina.recpresente 1,58 11,77
tabla15 <- svyby(~erodentina.rec, ~JugFrutas.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla15[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Nunca o raramente 93,40 6,60
## 2-Todos los dias 96,15 3,85
## 3-Mas de 3 veces al dia 92,74 7,26
## se.erodentina.recausente
## 1-Nunca o raramente 4,84
## 2-Todos los dias 0,90
## 3-Mas de 3 veces al dia 2,43
round(confint(tabla15) * 100, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente:erodentina.recausente 83,91 102,88
## 2-Todos los dias:erodentina.recausente 94,39 97,91
## 3-Mas de 3 veces al dia:erodentina.recausente 87,98 97,51
## 1-Nunca o raramente:erodentina.recpresente -2,88 16,09
## 2-Todos los dias:erodentina.recpresente 2,09 5,61
## 3-Mas de 3 veces al dia:erodentina.recpresente 2,49 12,02
tabla16 <- svyby(~erodentina.rec, ~Bebidas_energizantes_b, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla16[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-No diario 95,69 4,31
## 2-Diario 95,34 4,66
## se.erodentina.recausente
## 1-No diario 1,15
## 2-Diario 1,73
round(confint(tabla16) * 100, 2)
## 2,5 % 97,5 %
## 1-No diario:erodentina.recausente 93,43 97,94
## 2-Diario:erodentina.recausente 91,95 98,74
## 1-No diario:erodentina.recpresente 2,06 6,57
## 2-Diario:erodentina.recpresente 1,26 8,05
tabla17 <- svyby(~erodentina.rec, ~Yogurt_b, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Diario 96,32 3,68
## 2-No diario 95,27 4,73
## se.erodentina.recausente
## 1-Diario 1,68
## 2-No diario 1,05
round(confint(tabla17) * 100, 2)
## 2,5 % 97,5 %
## 1-Diario:erodentina.recausente 93,03 99,61
## 2-No diario:erodentina.recausente 93,21 97,32
## 1-Diario:erodentina.recpresente 0,39 6,97
## 2-No diario:erodentina.recpresente 2,68 6,79
tabla18 <- svyby(~erodentina.rec, ~Bruxismo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla18[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-No 95,22 4,78
## 2-Yes 95,86 4,14
## se.erodentina.recausente
## 1-No 1,06
## 2-Yes 1,57
round(confint(tabla18) * 100, 2)
## 2,5 % 97,5 %
## 1-No:erodentina.recausente 93,13 97,30
## 2-Yes:erodentina.recausente 92,77 98,94
## 1-No:erodentina.recpresente 2,70 6,87
## 2-Yes:erodentina.recpresente 1,06 7,23
tabla19 <- svyby(~erodentina.rec, ~bedeportediario.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla19[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Water 95,73 4,27
## 2-Gatorade 90,28 9,72
## 3-other 98,87 1,13
## se.erodentina.recausente
## 1-Water 1,36
## 2-Gatorade 5,95
## 3-other 1,19
round(confint(tabla19) * 100, 2)
## 2,5 % 97,5 %
## 1-Water:erodentina.recausente 93,08 98,39
## 2-Gatorade:erodentina.recausente 78,63 101,94
## 3-other:erodentina.recausente 96,54 101,21
## 1-Water:erodentina.recpresente 1,61 6,92
## 2-Gatorade:erodentina.recpresente -1,94 21,37
## 3-other:erodentina.recpresente -1,21 3,46
tabla21 <- svyby(~erodentina.rec, ~BuchTragar.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla21[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-No 95,79 4,21
## 2-Yes 91,40 8,60
## se.erodentina.recausente
## 1-No 1,18
## 2-Yes 2,38
round(confint(tabla21) * 100, 2)
## 2,5 % 97,5 %
## 1-No:erodentina.recausente 93,48 98,09
## 2-Yes:erodentina.recausente 86,74 96,06
## 1-No:erodentina.recpresente 1,91 6,52
## 2-Yes:erodentina.recpresente 3,94 13,26
tabla22 <- svyby(~erodentina.rec, ~FormBeber.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla22[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-Con sorbete 95,51 4,49
## 2-Por el pico 97,10 2,90
## 3-Con vaso 94,96 5,04
## se.erodentina.recausente
## 1-Con sorbete 2,75
## 2-Por el pico 1,68
## 3-Con vaso 1,19
round(confint(tabla22) * 100, 2)
## 2,5 % 97,5 %
## 1-Con sorbete:erodentina.recausente 90,12 100,89
## 2-Por el pico:erodentina.recausente 93,80 100,39
## 3-Con vaso:erodentina.recausente 92,62 97,30
## 1-Con sorbete:erodentina.recpresente -0,89 9,88
## 2-Por el pico:erodentina.recpresente -0,39 6,20
## 3-Con vaso:erodentina.recpresente 2,70 7,38
tabla23 <- svyby(~erodentina.rec, ~Natac2vec.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla23[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-No 95,44 4,56
## 2-Yes 94,91 5,09
## se.erodentina.recausente
## 1-No 1,25
## 2-Yes 1,71
round(confint(tabla23) * 100, 2)
## 2,5 % 97,5 %
## 1-No:erodentina.recausente 93,00 97,89
## 2-Yes:erodentina.recausente 91,55 98,26
## 1-No:erodentina.recpresente 2,11 7,00
## 2-Yes:erodentina.recpresente 1,74 8,45
tabla24 <- svyby(~erodentina.rec, ~MedResp.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla24[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-No 95,03 4,97
## 2-Yes 97,00 3,00
## se.erodentina.recausente
## 1-No 1,05
## 2-Yes 1,69
round(confint(tabla24) * 100, 2)
## 2,5 % 97,5 %
## 1-No:erodentina.recausente 92,96 97,09
## 2-Yes:erodentina.recausente 93,68 100,32
## 1-No:erodentina.recpresente 2,91 7,04
## 2-Yes:erodentina.recpresente -0,32 6,32
tabla25 <- svyby(~erodentina.rec, ~AlterGastrica.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla25[, 2:4] * 100, 2)
## erodentina.recausente erodentina.recpresente
## 1-No 95,05 4,95
## 2-Yes 100,00 0,00
## se.erodentina.recausente
## 1-No 1,03
## 2-Yes 0,00
round(confint(tabla25) * 100, 2)
## 2,5 % 97,5 %
## 1-No:erodentina.recausente 93,04 97,06
## 2-Yes:erodentina.recausente 100,00 100,00
## 1-No:erodentina.recpresente 2,94 6,96
## 2-Yes:erodentina.recpresente 0,00 0,00