This is an R Markdown document. Markdown is a simple formatting syntax for authoring web pages (click the MD toolbar button for help on Markdown).
When you click the Knit HTML button a web page will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
load("~/Dropbox/odontologia/maestria_anunziatta/julio2013/datos_tana_20072013.RData")
library(survey)
## Attaching package: 'survey'
## The following object(s) are masked from 'package:graphics':
##
## dotchart
# sabado 20 de julio 2013
# extension
tabla6w <- svyby(~WHODMFT, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE, deff = TRUE)
tabla6w[, 2:4] * 1
## WHODMFT se DEff.WHODMFT
## 1-M 1.616 0.1404 3.291
## 2-F 1.680 0.1131 2.014
round(confint(tabla6w) * 1, 2)
## 2.5 % 97.5 %
## 1-M 1.34 1.89
## 2-F 1.46 1.90
svychisq(~WHODMFT + Sexo.rec, diseniopost1, statistic = "Chisq")
## Warning: Chi-squared approximation may be incorrect
##
## Pearson's X^2: Rao & Scott adjustment
##
## data: svychisq(~WHODMFT + Sexo.rec, diseniopost1, statistic = "Chisq")
## X-squared = 12.5, df = 13, p-value = 0.796
tabla7w <- svyby(~WHODMFT, ~Nivel.Socieconomico.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
tabla7w[, 2:4] * 1
## WHODMFT se DEff.WHODMFT
## 1-ALTO 0.9142 0.1005 1.666
## 2-MEDIO 1.9631 0.1404 2.718
## 3-BAJO 2.0723 0.1340 1.315
round(confint(tabla7w) * 1, 2)
## 2.5 % 97.5 %
## 1-ALTO 0.72 1.11
## 2-MEDIO 1.69 2.24
## 3-BAJO 1.81 2.33
svychisq(~WHODMFT + Nivel.Socieconomico, diseniopost1, statistic = "Chisq")
## Error: sistema es computacionalmente singular: número de condición
## recíproco = 8.60475e-18
tabla8w <- svyby(~WHODMFT, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1-College-University 0.8 0.1 1.0
## 2-High School 1.5 0.2 2.4
## 3-Elementary School 2.0 0.1 2.3
round(confint(tabla8w) * 1, 2)
## 2.5 % 97.5 %
## 1-College-University 0.59 0.99
## 2-High School 1.07 1.86
## 3-Elementary School 1.75 2.18
tabla9w <- svyby(~WHODMFT, ~Tipo.de.Escuela, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla9w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 2-Privada 1.1 0.1 0.8
## 1-Publica 1.9 0.1 2.7
tabla10w <- svyby(~WHODMFT, ~AtenOdonto.1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla10w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1 1.1 0.1 1.7
## 2 1.6 0.1 2.0
## 3 2.2 0.1 1.6
## 4 1.9 0.2 2.7
## 99 2.2 0.2 0.9
tabla11w <- svyby(~WHODMFT, ~AtenOdonto2, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla11w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1 1.4 0.1 2.0
## 2 2.2 0.1 1.6
## 3 1.9 0.2 2.7
## 99 2.2 0.2 0.9
tabla12w <- svyby(~WHODMFT, ~FrCepDenti.4.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla12w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1-algun dia de la semana 2.3 0.3 2.5
## 2-1 vez al dias 1.8 0.2 1.5
## 3-2 veces al dia 1.7 0.2 2.7
## 4-3 veces al dia o mas 1.4 0.1 2.1
round(confint(tabla12w) * 1, 2)
## 2.5 % 97.5 %
## 1-algun dia de la semana 1.68 2.95
## 2-1 vez al dias 1.40 2.12
## 3-2 veces al dia 1.40 2.06
## 4-3 veces al dia o mas 1.23 1.67
tabla13w <- svyby(~WHODMFT, ~UsoDentifrico3, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla13w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1 1.6 0.1 2.1
## 2 1.7 0.5 1.7
## 99 2.8 0.7 1.2
tabla14w <- svyby(~WHODMFT, ~FluorProf., diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla14w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1 1.7 0.1 2.0
## 2 1.5 0.1 1.8
## 99 1.9 0.2 0.9
tabla15w <- svyby(~WHODMFT, ~RefrColaB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla15w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1-Mas de 3 veces al dia 2.0 0.2 2.0
## 2-A veces 1.5 0.1 1.4
## 3-Nunca o raramente 1.5 0.2 1.3
round(confint(tabla15w) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 1.64 2.37
## 2-A veces 1.39 1.70
## 3-Nunca o raramente 1.20 1.80
tabla16w <- svyby(~WHODMFT, ~MateDulceB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla16w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1-Mas de 3 veces al dia 1.8 0.3 3.5
## 2-A veces 2.1 0.3 5.7
## 3-Nunca o raramente 1.4 0.1 0.8
tabla17w <- svyby(~WHODMFT, ~GolosinasB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1-Mas de 3 veces al dia 2.0 0.2 2.6
## 2-A veces 1.6 0.1 1.6
## 3-Nunca o raramente 1.5 0.1 1.5
round(confint(tabla17w) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 1.60 2.39
## 2-A veces 1.37 1.74
## 3-Nunca o raramente 1.19 1.72
tabla18w <- svyby(~WHODMFT, ~GolosinasB, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla18w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1 1.5 0.1 1.7
## 2 2.0 0.2 2.6
tabla19w <- svyby(~WHODMFT, ~Masas.DulcesB, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla19w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1 1.6 0.1 1.5
## 2 2.1 0.2 2.4
## 99 1.4 0.4 1.4
tabla20w <- svyby(~WHODMFT, ~UltmVisita1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla20w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1 1.6 0.1 2.1
## 2 1.7 0.4 6.7
## 3 2.0 0.2 2.5
## 99 1.6 0.3 1.4
tabla21w <- svyby(~WHODMFT, ~IGS.rec1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla21w[, 2:4] * 1, 1)
## WHODMFT se DEff.WHODMFT
## 1 -Menor o igual a 20 1.6 0.1 3.6
## 2-De 20 a 60 1.7 0.1 1.0
## 3-Mas de 60 1.9 0.2 1.1
# extension
tabla6wm <- svyby(~WHO.modified, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
tabla6wm[, 2:4] * 1
## WHO.modified se DEff.WHO.modified
## 1-M 2.881 0.3082 5.255
## 2-F 2.805 0.1860 2.452
round(confint(tabla6wm) * 1, 2)
## 2.5 % 97.5 %
## 1-M 2.28 3.48
## 2-F 2.44 3.17
tabla7wm <- svyby(~WHO.modified, ~Nivel.Socieconomico.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
tabla7wm[, 2:4] * 1
## WHO.modified se DEff.WHO.modified
## 1-ALTO 2.024 0.3390 4.692
## 2-MEDIO 3.238 0.2389 3.305
## 3-BAJO 3.254 0.3084 2.666
round(confint(tabla7wm) * 1, 2)
## 2.5 % 97.5 %
## 1-ALTO 1.36 2.69
## 2-MEDIO 2.77 3.71
## 3-BAJO 2.65 3.86
tabla8wm <- svyby(~WHO.modified, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-College-University 1.7 0.4 3.4
## 2-High School 2.7 0.3 2.9
## 3-Elementary School 3.2 0.2 3.6
round(confint(tabla8wm) * 1, 2)
## 2.5 % 97.5 %
## 1-College-University 0.89 2.46
## 2-High School 2.04 3.38
## 3-Elementary School 2.81 3.65
tabla9wm <- svyby(~WHO.modified, ~Tipo.de.Escuela.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla9wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-Privada 2.0 0.2 1.1
## 2-Publica 3.2 0.3 7.0
round(confint(tabla9wm) * 1, 2)
## 2.5 % 97.5 %
## 1-Privada 1.64 2.28
## 2-Publica 2.62 3.71
tabla11wm <- svyby(~WHO.modified, ~AtenOdonto2.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla11wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-privada/convenio 2.5 0.3 6.2
## 2-publica 3.3 0.3 3.0
## 3-nunca fue al dentista 3.6 0.4 3.0
round(confint(tabla11wm) * 1, 2)
## 2.5 % 97.5 %
## 1-privada/convenio 1.94 3.05
## 2-publica 2.66 3.89
## 3-nunca fue al dentista 2.71 4.43
tabla12wm <- svyby(~WHO.modified, ~FrCepDenti.4.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla12wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-algun dia de la semana 3.8 0.6 2.8
## 2-1 vez al dias 3.1 0.4 2.1
## 3-2 veces al dia 2.8 0.3 3.1
## 4-3 veces al dia o mas 2.6 0.3 3.9
round(confint(tabla12wm) * 1, 2)
## 2.5 % 97.5 %
## 1-algun dia de la semana 2.58 4.97
## 2-1 vez al dias 2.37 3.87
## 3-2 veces al dia 2.31 3.38
## 4-3 veces al dia o mas 2.13 3.12
tabla13wm <- svyby(~WHO.modified, ~UsoDentifrico3.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla13wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-Si 2.8 0.2 5.1
## 2-No 2.6 0.8 1.9
round(confint(tabla13wm) * 1, 2)
## 2.5 % 97.5 %
## 1-Si 2.44 3.25
## 2-No 0.93 4.24
tabla14wm <- svyby(~WHO.modified, ~FluorProf.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla14wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-Si 2.9 0.3 4.8
## 2-No 2.8 0.2 3.3
round(confint(tabla14wm) * 1, 2)
## 2.5 % 97.5 %
## 1-Si 2.31 3.47
## 2-No 2.30 3.23
tabla15wm <- svyby(~WHO.modified, ~RefrColaB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla15wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-Mas de 3 veces al dia 3.3 0.3 2.5
## 2-A veces 2.7 0.2 3.7
## 3-Nunca o raramente 2.6 0.3 2.2
round(confint(tabla15wm) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 2.70 3.99
## 2-A veces 2.31 3.14
## 3-Nunca o raramente 1.96 3.17
tabla16wm <- svyby(~WHO.modified, ~MateDulceB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla16wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-Mas de 3 veces al dia 3.1 0.5 3.5
## 2-A veces 3.5 0.5 6.1
## 3-Nunca o raramente 2.5 0.2 3.4
round(confint(tabla16wm) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 2.14 4.03
## 2-A veces 2.62 4.40
## 3-Nunca o raramente 2.13 2.94
tabla17wm <- svyby(~WHO.modified, ~GolosinasB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-Mas de 3 veces al dia 3.4 0.4 5.1
## 2-A veces 2.8 0.2 2.2
## 3-Nunca o raramente 2.3 0.2 1.7
round(confint(tabla17wm) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 2.50 4.24
## 2-A veces 2.43 3.14
## 3-Nunca o raramente 1.88 2.78
tabla19wm <- svyby(~WHO.modified, ~Masas.DulcesB.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla19wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1-No diario 2.7 0.2 3.8
## 2-Diario 3.4 0.4 3.2
round(confint(tabla19wm) * 1, 2)
## 2.5 % 97.5 %
## 1-No diario 2.36 3.13
## 2-Diario 2.66 4.14
tabla20wm <- svyby(~WHO.modified, ~UltmVisita1.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla20wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1 año atrás o menos 2.7 0.3 6.4
## 2- años atrás 3.1 0.5 3.6
## 3-Nunca fue al dentista 3.6 0.4 2.9
round(confint(tabla20wm) * 1, 2)
## 2.5 % 97.5 %
## 1 año atrás o menos 2.17 3.22
## 2- años atrás 2.20 3.99
## 3-Nunca fue al dentista 2.74 4.40
tabla21wm <- svyby(~WHO.modified, ~IGS.rec1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla21wm[, 2:4] * 1, 1)
## WHO.modified se DEff.WHO.modified
## 1 -Menor o igual a 20 2.6 0.3 6.6
## 2-De 20 a 60 3.0 0.2 2.4
## 3-Mas de 60 3.7 0.3 1.0
round(confint(tabla21wm) * 1, 2)
## 2.5 % 97.5 %
## 1 -Menor o igual a 20 2.01 3.19
## 2-De 20 a 60 2.52 3.41
## 3-Mas de 60 3.08 4.37
# extension ICDAS
tabla6icdas <- svyby(~ICDAS, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
tabla6icdas[, 2:4] * 1
## ICDAS se DEff.ICDAS
## 1-M 5.376 0.3198 3.203
## 2-F 5.790 0.3510 3.817
round(confint(tabla6icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-M 4.75 6.00
## 2-F 5.10 6.48
tabla7icdas <- svyby(~ICDAS, ~Nivel.Socieconomico.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
tabla7icdas[, 2:4] * 1
## ICDAS se DEff.ICDAS
## 1-ALTO 3.971 0.1754 0.5891
## 2-MEDIO 6.126 0.2062 1.3095
## 3-BAJO 6.866 0.3277 1.6458
round(confint(tabla7icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-ALTO 3.63 4.32
## 2-MEDIO 5.72 6.53
## 3-BAJO 6.22 7.51
tabla8icdas <- svyby(~ICDAS, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(tabla8icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-College-University 3.6 0.3 0.8
## 2-High School 4.8 0.3 1.1
## 3-Elementary School 6.4 0.2 2.4
round(confint(tabla8icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-College-University 2.99 4.22
## 2-High School 4.23 5.33
## 3-Elementary School 5.97 6.92
tabla9icdas <- svyby(~ICDAS, ~Tipo.de.Escuela.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla9icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-Privada 4.1 0.4 1.9
## 2-Publica 6.1 0.2 2.7
round(confint(tabla9icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-Privada 3.45 4.84
## 2-Publica 5.64 6.55
tabla11icdas <- svyby(~ICDAS, ~AtenOdonto2.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla11icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-privada/convenio 4.8 0.2 1.4
## 2-publica 6.4 0.4 2.2
## 3-nunca fue al dentista 7.3 0.5 2.7
round(confint(tabla11icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-privada/convenio 4.41 5.19
## 2-publica 5.73 7.13
## 3-nunca fue al dentista 6.31 8.25
tabla12icdas <- svyby(~ICDAS, ~FrCepDenti.4.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla12icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-algun dia de la semana 6.6 0.5 1.6
## 2-1 vez al dias 5.5 0.5 2.1
## 3-2 veces al dia 5.7 0.3 1.5
## 4-3 veces al dia o mas 5.3 0.2 1.5
round(confint(tabla12icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-algun dia de la semana 5.53 7.67
## 2-1 vez al dias 4.60 6.47
## 3-2 veces al dia 5.11 6.21
## 4-3 veces al dia o mas 4.82 5.74
tabla13icdas <- svyby(~ICDAS, ~UsoDentifrico3.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla13icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-Si 5.6 0.2 2.2
## 2-No 4.8 1.1 1.7
round(confint(tabla13icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-Si 5.20 5.95
## 2-No 2.74 6.95
tabla14icdas <- svyby(~ICDAS, ~FluorProf.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla14icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-Si 5.5 0.3 1.8
## 2-No 5.5 0.3 2.6
round(confint(tabla14icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-Si 5.00 6.01
## 2-No 4.96 6.08
tabla15icdas <- svyby(~ICDAS, ~RefrColaB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla15icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-Mas de 3 veces al dia 6.0 0.4 2.5
## 2-A veces 5.3 0.2 1.7
## 3-Nunca o raramente 6.1 0.5 1.8
round(confint(tabla15icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 5.21 6.87
## 2-A veces 4.86 5.66
## 3-Nunca o raramente 5.19 7.01
tabla16icdas <- svyby(~ICDAS, ~MateDulceB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla16icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-Mas de 3 veces al dia 7.0 0.5 2.1
## 2-A veces 6.4 0.4 3.4
## 3-Nunca o raramente 4.9 0.2 1.0
round(confint(tabla16icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 6.05 7.92
## 2-A veces 5.55 7.21
## 3-Nunca o raramente 4.56 5.21
tabla17icdas <- svyby(~ICDAS, ~GolosinasB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla17icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-Mas de 3 veces al dia 6.1 0.4 2.0
## 2-A veces 5.5 0.2 1.7
## 3-Nunca o raramente 5.0 0.4 2.2
round(confint(tabla17icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-Mas de 3 veces al dia 5.39 6.87
## 2-A veces 5.05 5.92
## 3-Nunca o raramente 4.20 5.84
tabla19icdas <- svyby(~ICDAS, ~Masas.DulcesB.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla19icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1-No diario 5.4 0.2 2.1
## 2-Diario 6.5 0.4 1.9
round(confint(tabla19icdas) * 1, 2)
## 2.5 % 97.5 %
## 1-No diario 4.96 5.77
## 2-Diario 5.75 7.26
tabla20icdas <- svyby(~ICDAS, ~UltmVisita1.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla20icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1 año atrás o menos 5.2 0.2 1.8
## 2- años atrás 6.3 0.5 2.6
## 3-Nunca fue al dentista 7.2 0.5 2.6
round(confint(tabla20icdas) * 1, 2)
## 2.5 % 97.5 %
## 1 año atrás o menos 4.76 5.57
## 2- años atrás 5.41 7.28
## 3-Nunca fue al dentista 6.28 8.20
tabla21icdas <- svyby(~ICDAS, ~IGS.rec1, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla21icdas[, 2:4] * 1, 1)
## ICDAS se DEff.ICDAS
## 1 -Menor o igual a 20 5.5 0.3 3.4
## 2-De 20 a 60 5.5 0.3 2.7
## 3-Mas de 60 6.3 0.6 1.6
round(confint(tabla21icdas) * 1, 2)
## 2.5 % 97.5 %
## 1 -Menor o igual a 20 4.87 6.18
## 2-De 20 a 60 4.89 6.13
## 3-Mas de 60 5.20 7.42