Extension

Aca estan las tablas de WHODFMT considerandola como variable de conteo, donde se presentan los promedios, los IC No hay que perder de vista que son variables asimetricas La codificacion de las varables independientes corresponde al mail de 15 de julio

load("~/Dropbox/odontologia/maestria_anunziatta/julio2013/datos_tana_25072013.RData")


library(survey)
## Attaching package: 'survey'
## The following object(s) are masked from 'package:graphics':
## 
## dotchart
# jueves 25 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


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


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.rec, diseniopost1, svymean, na.rm = TRUE, 
    deff = TRUE)
round(tabla9w[, 2:4] * 1, 1)
##           WHODMFT  se DEff.WHODMFT
## 1-Privada     1.1 0.1          0.8
## 2-Publica     1.9 0.1          2.7
round(confint(tabla9w) * 1, 2)
##           2.5 % 97.5 %
## 1-Privada  0.92   1.23
## 2-Publica  1.65   2.07


tabla10w <- svyby(~WHODMFT, ~AtenOdonto2.rec, diseniopost1, svymean, na.rm = TRUE, 
    deff = TRUE)
round(tabla10w[, 2:4] * 1, 1)
##                          WHODMFT  se DEff.WHODMFT
## 1-privada/convenio           1.4 0.1          2.0
## 2-publica                    2.2 0.1          1.6
## 3-nunca fue al dentista      1.9 0.2          2.7
round(confint(tabla10w) * 1, 2)
##                          2.5 % 97.5 %
## 1-privada/convenio        1.16   1.54
## 2-publica                 1.88   2.46
## 3-nunca fue al dentista   1.50   2.38

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.rec, diseniopost1, svymean, na.rm = TRUE, 
    deff = TRUE)
round(tabla13w[, 2:4] * 1, 1)
##      WHODMFT  se DEff.WHODMFT
## 1-Si     1.6 0.1          2.1
## 2-No     1.7 0.5          1.7
round(confint(tabla13w) * 1, 2)
##      2.5 % 97.5 %
## 1-Si  1.48   1.80
## 2-No  0.73   2.76

tabla14w <- svyby(~WHODMFT, ~FluorProf.rec, diseniopost1, svymean, na.rm = TRUE, 
    deff = TRUE)
round(tabla14w[, 2:4] * 1, 1)
##      WHODMFT  se DEff.WHODMFT
## 1-Si     1.7 0.1          2.0
## 2-No     1.5 0.1          1.8
round(confint(tabla14w) * 1, 2)
##      2.5 % 97.5 %
## 1-Si  1.48   1.96
## 2-No  1.31   1.72


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
round(confint(tabla16w) * 1, 2)
##                         2.5 % 97.5 %
## 1-Mas de 3 veces al dia  1.19   2.36
## 2-A veces                1.54   2.59
## 3-Nunca o raramente      1.33   1.57

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

tabla19w <- svyby(~WHODMFT, ~Masas.DulcesB.rec, diseniopost1, svymean, na.rm = TRUE, 
    deff = TRUE)
round(tabla19w[, 2:4] * 1, 1)
##                         WHODMFT  se DEff.WHODMFT
## 1-Mas de 3 veces al dia     2.1 0.2          2.4
## 2-A veces                   1.6 0.1          1.7
## 3-Nunca o raramente         1.5 0.2          1.5
round(confint(tabla19w) * 1, 2)
##                         2.5 % 97.5 %
## 1-Mas de 3 veces al dia  1.71   2.51
## 2-A veces                1.39   1.75
## 3-Nunca o raramente      1.15   1.80



tabla20w <- svyby(~WHODMFT, ~UltmVisita1.rec, diseniopost1, svymean, na.rm = TRUE, 
    deff = TRUE)
round(tabla20w[, 2:4] * 1, 1)
##                         WHODMFT  se DEff.WHODMFT
## 1 año atrás o menos         1.6 0.1          2.1
## 2- años atrás               1.7 0.4          6.7
## 3-Nunca fue al dentista     2.0 0.2          2.5
round(confint(tabla20w) * 1, 2)
##                         2.5 % 97.5 %
## 1 año atrás o menos      1.40   1.79
## 2- años atrás            0.90   2.45
## 3-Nunca fue al dentista  1.54   2.38

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
round(confint(tabla21w) * 1, 2)
##                       2.5 % 97.5 %
## 1 -Menor o igual a 20  1.27   1.84
## 2-De 20 a 60           1.53   1.87
## 3-Mas de 60            1.54   2.33


# 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-Mas de 3 veces al dia          3.4 0.4               3.2
## 2-A veces                        2.8 0.2               2.4
## 3-Nunca o raramente              2.5 0.5               4.6
round(confint(tabla19wm) * 1, 2)
##                         2.5 % 97.5 %
## 1-Mas de 3 veces al dia  2.66   4.14
## 2-A veces                2.45   3.13
## 3-Nunca o raramente      1.61   3.47



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-Mas de 3 veces al dia   6.5 0.4        1.9
## 2-A veces                 5.3 0.3        2.5
## 3-Nunca o raramente       5.5 0.5        2.8
round(confint(tabla19icdas) * 1, 2)
##                         2.5 % 97.5 %
## 1-Mas de 3 veces al dia  5.75   7.26
## 2-A veces                4.83   5.83
## 3-Nunca o raramente      4.49   6.54



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