Aca estan las modelos de WHODFMT considerandola como variable de conteo Son regresiones multiivariadas se hacen 3 escenarios cambiando atencion odontologica y visita al dentisita
#
# load('~/Dropbox/odontologia/maestria_anunziatta/julio2013/datos_tana_25072013.RData')
load("C:/Users/usuario/Dropbox/odontologia/maestria_anunziatta/julio2013/datos_tana_25072013.RData")
library(survey)
## Attaching package: 'survey'
## The following object(s) are masked from 'package:graphics':
##
## dotchart
# 20 de agosto 2013
# WHO
modelomulti.poi1 <- svyglm(WHODMFT ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
Tipo.de.Escuela.rec + AtenOdonto2.rec + FrCepDenti.4.rec + RefrColaB.rec +
GolosinasB.rec + IGS.rec1 + Masas.DulcesB.rec + UltmVisita1.rec, diseniopost1,
family = quasipoisson())
summary(modelomulti.poi1)
##
## Call:
## svyglm(formula = WHODMFT ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
## Tipo.de.Escuela.rec + AtenOdonto2.rec + FrCepDenti.4.rec +
## RefrColaB.rec + GolosinasB.rec + IGS.rec1 + Masas.DulcesB.rec +
## UltmVisita1.rec, diseniopost1, family = quasipoisson())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error
## (Intercept) 0.0324 0.2629
## Nivel.Socieconomico.rec2-MEDIO 0.4981 0.1465
## Nivel.Socieconomico.rec3-BAJO 0.3985 0.1716
## Nivel.Educativo.de.la.Madre1.rec2-High School 0.4037 0.2141
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.4677 0.2071
## Tipo.de.Escuela.rec2-Publica 0.1385 0.1222
## AtenOdonto2.rec2-publica 0.1905 0.1103
## AtenOdonto2.rec3-nunca fue al dentista -0.8908 0.2317
## FrCepDenti.4.rec2-1 vez al dias -0.2044 0.1715
## FrCepDenti.4.rec3-2 veces al dia -0.1674 0.1755
## FrCepDenti.4.rec4-3 veces al dia o mas -0.2992 0.1044
## RefrColaB.rec2-A veces -0.0818 0.0919
## RefrColaB.rec3-Nunca o raramente -0.1776 0.1469
## GolosinasB.rec2-A veces -0.1294 0.1345
## GolosinasB.rec3-Nunca o raramente -0.1782 0.1859
## IGS.rec12-De 20 a 60 0.0412 0.1186
## IGS.rec13-Mas de 60 0.2113 0.1058
## Masas.DulcesB.rec2-A veces -0.0639 0.0835
## Masas.DulcesB.rec3-Nunca o raramente -0.3298 0.1951
## UltmVisita1.rec2- años atrás -0.1151 0.2258
## UltmVisita1.rec3-Nunca fue al dentista 0.8993 0.2542
## t value Pr(>|t|)
## (Intercept) 0.12 0.9032
## Nivel.Socieconomico.rec2-MEDIO 3.40 0.0030 **
## Nivel.Socieconomico.rec3-BAJO 2.32 0.0314 *
## Nivel.Educativo.de.la.Madre1.rec2-High School 1.89 0.0747 .
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 2.26 0.0358 *
## Tipo.de.Escuela.rec2-Publica 1.13 0.2714
## AtenOdonto2.rec2-publica 1.73 0.1003
## AtenOdonto2.rec3-nunca fue al dentista -3.84 0.0011 **
## FrCepDenti.4.rec2-1 vez al dias -1.19 0.2482
## FrCepDenti.4.rec3-2 veces al dia -0.95 0.3522
## FrCepDenti.4.rec4-3 veces al dia o mas -2.87 0.0099 **
## RefrColaB.rec2-A veces -0.89 0.3850
## RefrColaB.rec3-Nunca o raramente -1.21 0.2415
## GolosinasB.rec2-A veces -0.96 0.3481
## GolosinasB.rec3-Nunca o raramente -0.96 0.3497
## IGS.rec12-De 20 a 60 0.35 0.7322
## IGS.rec13-Mas de 60 2.00 0.0604 .
## Masas.DulcesB.rec2-A veces -0.77 0.4532
## Masas.DulcesB.rec3-Nunca o raramente -1.69 0.1072
## UltmVisita1.rec2- años atrás -0.51 0.6162
## UltmVisita1.rec3-Nunca fue al dentista 3.54 0.0022 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for quasipoisson family taken to be 1.771)
##
## Number of Fisher Scoring iterations: 6
confint(modelomulti.poi1)
## 2.5 % 97.5 %
## (Intercept) -0.482859 0.54768
## Nivel.Socieconomico.rec2-MEDIO 0.210977 0.78520
## Nivel.Socieconomico.rec3-BAJO 0.062244 0.73475
## Nivel.Educativo.de.la.Madre1.rec2-High School -0.015910 0.82330
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.061871 0.87359
## Tipo.de.Escuela.rec2-Publica -0.101113 0.37804
## AtenOdonto2.rec2-publica -0.025648 0.40656
## AtenOdonto2.rec3-nunca fue al dentista -1.344971 -0.43668
## FrCepDenti.4.rec2-1 vez al dias -0.540525 0.13182
## FrCepDenti.4.rec3-2 veces al dia -0.511406 0.17661
## FrCepDenti.4.rec4-3 veces al dia o mas -0.503734 -0.09460
## RefrColaB.rec2-A veces -0.261936 0.09843
## RefrColaB.rec3-Nunca o raramente -0.465533 0.11032
## GolosinasB.rec2-A veces -0.393066 0.13422
## GolosinasB.rec3-Nunca o raramente -0.542481 0.18609
## IGS.rec12-De 20 a 60 -0.191196 0.27353
## IGS.rec13-Mas de 60 0.003897 0.41871
## Masas.DulcesB.rec2-A veces -0.227510 0.09968
## Masas.DulcesB.rec3-Nunca o raramente -0.712168 0.05252
## UltmVisita1.rec2- años atrás -0.557560 0.32745
## UltmVisita1.rec3-Nunca fue al dentista 0.401057 1.39750
modelomulti.poi2 <- svyglm(WHODMFT ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
Tipo.de.Escuela.rec + FrCepDenti.4.rec + RefrColaB.rec + GolosinasB.rec +
IGS.rec1 + Masas.DulcesB.rec + UltmVisita1.rec, diseniopost1, family = quasipoisson())
summary(modelomulti.poi2)
##
## Call:
## svyglm(formula = WHODMFT ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
## Tipo.de.Escuela.rec + FrCepDenti.4.rec + RefrColaB.rec +
## GolosinasB.rec + IGS.rec1 + Masas.DulcesB.rec + UltmVisita1.rec,
## diseniopost1, family = quasipoisson())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error
## (Intercept) 0.1750 0.2529
## Nivel.Socieconomico.rec2-MEDIO 0.5628 0.1488
## Nivel.Socieconomico.rec3-BAJO 0.5051 0.1682
## Nivel.Educativo.de.la.Madre1.rec2-High School 0.3356 0.2148
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.4101 0.2097
## Tipo.de.Escuela.rec2-Publica 0.1354 0.1115
## FrCepDenti.4.rec2-1 vez al dias -0.2655 0.1731
## FrCepDenti.4.rec3-2 veces al dia -0.2399 0.1797
## FrCepDenti.4.rec4-3 veces al dia o mas -0.3726 0.1048
## RefrColaB.rec2-A veces -0.0717 0.0946
## RefrColaB.rec3-Nunca o raramente -0.1237 0.1364
## GolosinasB.rec2-A veces -0.1492 0.1322
## GolosinasB.rec3-Nunca o raramente -0.2045 0.1782
## IGS.rec12-De 20 a 60 0.0537 0.1112
## IGS.rec13-Mas de 60 0.2028 0.0926
## Masas.DulcesB.rec2-A veces -0.0830 0.0887
## Masas.DulcesB.rec3-Nunca o raramente -0.2870 0.1788
## UltmVisita1.rec2- años atrás -0.0791 0.2188
## UltmVisita1.rec3-Nunca fue al dentista -0.0841 0.1276
## t value Pr(>|t|)
## (Intercept) 0.69 0.4966
## Nivel.Socieconomico.rec2-MEDIO 3.78 0.0011 **
## Nivel.Socieconomico.rec3-BAJO 3.00 0.0068 **
## Nivel.Educativo.de.la.Madre1.rec2-High School 1.56 0.1331
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 1.96 0.0639 .
## Tipo.de.Escuela.rec2-Publica 1.21 0.2382
## FrCepDenti.4.rec2-1 vez al dias -1.53 0.1400
## FrCepDenti.4.rec3-2 veces al dia -1.33 0.1962
## FrCepDenti.4.rec4-3 veces al dia o mas -3.56 0.0019 **
## RefrColaB.rec2-A veces -0.76 0.4569
## RefrColaB.rec3-Nunca o raramente -0.91 0.3745
## GolosinasB.rec2-A veces -1.13 0.2718
## GolosinasB.rec3-Nunca o raramente -1.15 0.2639
## IGS.rec12-De 20 a 60 0.48 0.6343
## IGS.rec13-Mas de 60 2.19 0.0400 *
## Masas.DulcesB.rec2-A veces -0.94 0.3597
## Masas.DulcesB.rec3-Nunca o raramente -1.61 0.1233
## UltmVisita1.rec2- años atrás -0.36 0.7212
## UltmVisita1.rec3-Nunca fue al dentista -0.66 0.5173
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for quasipoisson family taken to be 1.772)
##
## Number of Fisher Scoring iterations: 6
confint(modelomulti.poi2)
## 2.5 % 97.5 %
## (Intercept) -0.3207689 0.67077
## Nivel.Socieconomico.rec2-MEDIO 0.2711638 0.85450
## Nivel.Socieconomico.rec3-BAJO 0.1754758 0.83469
## Nivel.Educativo.de.la.Madre1.rec2-High School -0.0853647 0.75661
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School -0.0009104 0.82110
## Tipo.de.Escuela.rec2-Publica -0.0831912 0.35394
## FrCepDenti.4.rec2-1 vez al dias -0.6047395 0.07375
## FrCepDenti.4.rec3-2 veces al dia -0.5920686 0.11231
## FrCepDenti.4.rec4-3 veces al dia o mas -0.5779809 -0.16718
## RefrColaB.rec2-A veces -0.2572550 0.11376
## RefrColaB.rec3-Nunca o raramente -0.3910183 0.14353
## GolosinasB.rec2-A veces -0.4084057 0.10994
## GolosinasB.rec3-Nunca o raramente -0.5537869 0.14469
## IGS.rec12-De 20 a 60 -0.1642252 0.27153
## IGS.rec13-Mas de 60 0.0212095 0.38436
## Masas.DulcesB.rec2-A veces -0.2568311 0.09075
## Masas.DulcesB.rec3-Nunca o raramente -0.6374301 0.06334
## UltmVisita1.rec2- años atrás -0.5079478 0.34970
## UltmVisita1.rec3-Nunca fue al dentista -0.3341953 0.16607
modelomulti.poi3 <- svyglm(WHODMFT ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
Tipo.de.Escuela.rec + AtenOdonto2.rec + FrCepDenti.4.rec + RefrColaB.rec +
GolosinasB.rec + IGS.rec1 + Masas.DulcesB.rec, diseniopost1, family = quasipoisson())
summary(modelomulti.poi3)
##
## Call:
## svyglm(formula = WHODMFT ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
## Tipo.de.Escuela.rec + AtenOdonto2.rec + FrCepDenti.4.rec +
## RefrColaB.rec + GolosinasB.rec + IGS.rec1 + Masas.DulcesB.rec,
## diseniopost1, family = quasipoisson())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error
## (Intercept) 0.0191 0.2721
## Nivel.Socieconomico.rec2-MEDIO 0.4957 0.1432
## Nivel.Socieconomico.rec3-BAJO 0.3993 0.1607
## Nivel.Educativo.de.la.Madre1.rec2-High School 0.3933 0.2110
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.4551 0.2065
## Tipo.de.Escuela.rec2-Publica 0.1353 0.1280
## AtenOdonto2.rec2-publica 0.1922 0.1085
## AtenOdonto2.rec3-nunca fue al dentista 0.0307 0.1185
## FrCepDenti.4.rec2-1 vez al dias -0.2176 0.1904
## FrCepDenti.4.rec3-2 veces al dia -0.1642 0.1901
## FrCepDenti.4.rec4-3 veces al dia o mas -0.2997 0.1105
## RefrColaB.rec2-A veces -0.0794 0.0909
## RefrColaB.rec3-Nunca o raramente -0.1534 0.1381
## GolosinasB.rec2-A veces -0.1382 0.1364
## GolosinasB.rec3-Nunca o raramente -0.1919 0.1783
## IGS.rec12-De 20 a 60 0.0288 0.1154
## IGS.rec13-Mas de 60 0.1856 0.1069
## Masas.DulcesB.rec2-A veces -0.0454 0.0887
## Masas.DulcesB.rec3-Nunca o raramente -0.2919 0.2046
## t value Pr(>|t|)
## (Intercept) 0.07 0.9446
## Nivel.Socieconomico.rec2-MEDIO 3.46 0.0023 **
## Nivel.Socieconomico.rec3-BAJO 2.49 0.0215 *
## Nivel.Educativo.de.la.Madre1.rec2-High School 1.86 0.0764 .
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 2.20 0.0388 *
## Tipo.de.Escuela.rec2-Publica 1.06 0.3025
## AtenOdonto2.rec2-publica 1.77 0.0908 .
## AtenOdonto2.rec3-nunca fue al dentista 0.26 0.7981
## FrCepDenti.4.rec2-1 vez al dias -1.14 0.2660
## FrCepDenti.4.rec3-2 veces al dia -0.86 0.3976
## FrCepDenti.4.rec4-3 veces al dia o mas -2.71 0.0130 *
## RefrColaB.rec2-A veces -0.87 0.3924
## RefrColaB.rec3-Nunca o raramente -1.11 0.2792
## GolosinasB.rec2-A veces -1.01 0.3226
## GolosinasB.rec3-Nunca o raramente -1.08 0.2940
## IGS.rec12-De 20 a 60 0.25 0.8055
## IGS.rec13-Mas de 60 1.74 0.0971 .
## Masas.DulcesB.rec2-A veces -0.51 0.6137
## Masas.DulcesB.rec3-Nunca o raramente -1.43 0.1685
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for quasipoisson family taken to be 1.803)
##
## Number of Fisher Scoring iterations: 6
confint(modelomulti.poi3)
## 2.5 % 97.5 %
## (Intercept) -0.51414 0.55241
## Nivel.Socieconomico.rec2-MEDIO 0.21506 0.77630
## Nivel.Socieconomico.rec3-BAJO 0.08439 0.71419
## Nivel.Educativo.de.la.Madre1.rec2-High School -0.02029 0.80697
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.05032 0.85980
## Tipo.de.Escuela.rec2-Publica -0.11562 0.38628
## AtenOdonto2.rec2-publica -0.02034 0.40481
## AtenOdonto2.rec3-nunca fue al dentista -0.20148 0.26288
## FrCepDenti.4.rec2-1 vez al dias -0.59068 0.15556
## FrCepDenti.4.rec3-2 veces al dia -0.53683 0.20846
## FrCepDenti.4.rec4-3 veces al dia o mas -0.51626 -0.08314
## RefrColaB.rec2-A veces -0.25747 0.09877
## RefrColaB.rec3-Nunca o raramente -0.42413 0.11727
## GolosinasB.rec2-A veces -0.40554 0.12918
## GolosinasB.rec3-Nunca o raramente -0.54136 0.15754
## IGS.rec12-De 20 a 60 -0.19746 0.25501
## IGS.rec13-Mas de 60 -0.02386 0.39514
## Masas.DulcesB.rec2-A veces -0.21927 0.12838
## Masas.DulcesB.rec3-Nunca o raramente -0.69293 0.10920