Aca estan las modelos de WHODFMT considerandola como variable binaria Son regresiones multivariadas Se consioderan sacar AtenOdonto2.rec y UltmVisita1.rec por problemas de multicolinealidad
#
# 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.bin1 <- svyglm(WHODMFTcateg ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
Tipo.de.Escuela.rec + AtenOdonto2.rec + FrCepDenti.4.rec + GolosinasB.rec +
Masas.DulcesB.rec + UltmVisita1.rec, diseniopost1, family = quasibinomial())
summary(modelomulti.bin1)
##
## Call:
## svyglm(formula = WHODMFTcateg ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
## Tipo.de.Escuela.rec + AtenOdonto2.rec + FrCepDenti.4.rec +
## GolosinasB.rec + Masas.DulcesB.rec + UltmVisita1.rec, diseniopost1,
## family = quasibinomial())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error
## (Intercept) 0.606 0.514
## Nivel.Socieconomico.rec2-MEDIO 0.520 0.230
## Nivel.Socieconomico.rec3-BAJO 0.603 0.337
## Nivel.Educativo.de.la.Madre1.rec2-High School 0.510 0.305
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.747 0.304
## Tipo.de.Escuela.rec2-Publica -0.224 0.257
## AtenOdonto2.rec2-publica 0.736 0.214
## AtenOdonto2.rec3-nunca fue al dentista -13.941 0.716
## FrCepDenti.4.rec2-1 vez al dias -0.622 0.318
## FrCepDenti.4.rec3-2 veces al dia -0.549 0.286
## FrCepDenti.4.rec4-3 veces al dia o mas -0.727 0.296
## GolosinasB.rec2-A veces -0.380 0.303
## GolosinasB.rec3-Nunca o raramente -0.197 0.320
## Masas.DulcesB.rec2-A veces -0.123 0.204
## Masas.DulcesB.rec3-Nunca o raramente -0.666 0.387
## UltmVisita1.rec2- años atrás -0.371 0.296
## UltmVisita1.rec3-Nunca fue al dentista 14.490 0.803
## t value Pr(>|t|)
## (Intercept) 1.18 0.2504
## Nivel.Socieconomico.rec2-MEDIO 2.26 0.0339 *
## Nivel.Socieconomico.rec3-BAJO 1.79 0.0871 .
## Nivel.Educativo.de.la.Madre1.rec2-High School 1.67 0.1086
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 2.45 0.0222 *
## Tipo.de.Escuela.rec2-Publica -0.87 0.3912
## AtenOdonto2.rec2-publica 3.44 0.0023 **
## AtenOdonto2.rec3-nunca fue al dentista -19.47 8.7e-16 ***
## FrCepDenti.4.rec2-1 vez al dias -1.96 0.0627 .
## FrCepDenti.4.rec3-2 veces al dia -1.92 0.0675 .
## FrCepDenti.4.rec4-3 veces al dia o mas -2.45 0.0221 *
## GolosinasB.rec2-A veces -1.25 0.2223
## GolosinasB.rec3-Nunca o raramente -0.62 0.5435
## Masas.DulcesB.rec2-A veces -0.61 0.5509
## Masas.DulcesB.rec3-Nunca o raramente -1.72 0.0987 .
## UltmVisita1.rec2- años atrás -1.25 0.2222
## UltmVisita1.rec3-Nunca fue al dentista 18.04 4.5e-15 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for quasibinomial family taken to be 0.8795)
##
## Number of Fisher Scoring iterations: 13
confint(modelomulti.bin1)
## 2.5 % 97.5 %
## (Intercept) -0.40114 1.612954
## Nivel.Socieconomico.rec2-MEDIO 0.06804 0.971504
## Nivel.Socieconomico.rec3-BAJO -0.05841 1.264202
## Nivel.Educativo.de.la.Madre1.rec2-High School -0.08868 1.108410
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.15013 1.343500
## Tipo.de.Escuela.rec2-Publica -0.72782 0.278900
## AtenOdonto2.rec2-publica 0.31607 1.155883
## AtenOdonto2.rec3-nunca fue al dentista -15.34514 -12.537788
## FrCepDenti.4.rec2-1 vez al dias -1.24522 0.001224
## FrCepDenti.4.rec3-2 veces al dia -1.10906 0.011724
## FrCepDenti.4.rec4-3 veces al dia o mas -1.30764 -0.146293
## GolosinasB.rec2-A veces -0.97260 0.213533
## GolosinasB.rec3-Nunca o raramente -0.82488 0.430020
## Masas.DulcesB.rec2-A veces -0.52308 0.276225
## Masas.DulcesB.rec3-Nunca o raramente -1.42545 0.092607
## UltmVisita1.rec2- años atrás -0.95174 0.208843
## UltmVisita1.rec3-Nunca fue al dentista 12.91609 16.064629
modelomulti.bin2 <- svyglm(WHODMFTcateg ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
Tipo.de.Escuela.rec + FrCepDenti.4.rec + GolosinasB.rec + Masas.DulcesB.rec +
UltmVisita1.rec, diseniopost1, family = quasibinomial())
summary(modelomulti.bin2)
##
## Call:
## svyglm(formula = WHODMFTcateg ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
## Tipo.de.Escuela.rec + FrCepDenti.4.rec + GolosinasB.rec +
## Masas.DulcesB.rec + UltmVisita1.rec, diseniopost1, family = quasibinomial())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error
## (Intercept) 0.903 0.514
## Nivel.Socieconomico.rec2-MEDIO 0.691 0.232
## Nivel.Socieconomico.rec3-BAJO 0.911 0.320
## Nivel.Educativo.de.la.Madre1.rec2-High School 0.483 0.295
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.753 0.311
## Tipo.de.Escuela.rec2-Publica -0.189 0.243
## FrCepDenti.4.rec2-1 vez al dias -0.799 0.361
## FrCepDenti.4.rec3-2 veces al dia -0.758 0.327
## FrCepDenti.4.rec4-3 veces al dia o mas -0.953 0.317
## GolosinasB.rec2-A veces -0.386 0.308
## GolosinasB.rec3-Nunca o raramente -0.254 0.286
## Masas.DulcesB.rec2-A veces -0.171 0.204
## Masas.DulcesB.rec3-Nunca o raramente -0.639 0.341
## UltmVisita1.rec2- años atrás -0.254 0.283
## UltmVisita1.rec3-Nunca fue al dentista 0.229 0.297
## t value Pr(>|t|)
## (Intercept) 1.76 0.0915 .
## Nivel.Socieconomico.rec2-MEDIO 2.98 0.0064 **
## Nivel.Socieconomico.rec3-BAJO 2.85 0.0087 **
## Nivel.Educativo.de.la.Madre1.rec2-High School 1.64 0.1138
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 2.42 0.0229 *
## Tipo.de.Escuela.rec2-Publica -0.78 0.4437
## FrCepDenti.4.rec2-1 vez al dias -2.21 0.0365 *
## FrCepDenti.4.rec3-2 veces al dia -2.32 0.0288 *
## FrCepDenti.4.rec4-3 veces al dia o mas -3.01 0.0059 **
## GolosinasB.rec2-A veces -1.25 0.2214
## GolosinasB.rec3-Nunca o raramente -0.89 0.3825
## Masas.DulcesB.rec2-A veces -0.84 0.4112
## Masas.DulcesB.rec3-Nunca o raramente -1.87 0.0727 .
## UltmVisita1.rec2- años atrás -0.90 0.3779
## UltmVisita1.rec3-Nunca fue al dentista 0.77 0.4476
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for quasibinomial family taken to be 0.9049)
##
## Number of Fisher Scoring iterations: 4
confint(modelomulti.bin2)
## 2.5 % 97.5 %
## (Intercept) -0.1054 1.91110
## Nivel.Socieconomico.rec2-MEDIO 0.2357 1.14573
## Nivel.Socieconomico.rec3-BAJO 0.2839 1.53880
## Nivel.Educativo.de.la.Madre1.rec2-High School -0.0947 1.06015
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.1440 1.36183
## Tipo.de.Escuela.rec2-Publica -0.6647 0.28683
## FrCepDenti.4.rec2-1 vez al dias -1.5070 -0.09029
## FrCepDenti.4.rec3-2 veces al dia -1.3993 -0.11765
## FrCepDenti.4.rec4-3 veces al dia o mas -1.5738 -0.33215
## GolosinasB.rec2-A veces -0.9904 0.21751
## GolosinasB.rec3-Nunca o raramente -0.8151 0.30638
## Masas.DulcesB.rec2-A veces -0.5706 0.22941
## Masas.DulcesB.rec3-Nunca o raramente -1.3080 0.02946
## UltmVisita1.rec2- años atrás -0.8099 0.30104
## UltmVisita1.rec3-Nunca fue al dentista -0.3532 0.81176
modelomulti.bin3 <- svyglm(WHODMFTcateg ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
Tipo.de.Escuela.rec + AtenOdonto2.rec + FrCepDenti.4.rec + GolosinasB.rec +
Masas.DulcesB.rec, diseniopost1, family = quasibinomial())
summary(modelomulti.bin3)
##
## Call:
## svyglm(formula = WHODMFTcateg ~ Nivel.Socieconomico.rec + Nivel.Educativo.de.la.Madre1.rec +
## Tipo.de.Escuela.rec + AtenOdonto2.rec + FrCepDenti.4.rec +
## GolosinasB.rec + Masas.DulcesB.rec, diseniopost1, family = quasibinomial())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error
## (Intercept) 0.5556 0.4819
## Nivel.Socieconomico.rec2-MEDIO 0.5203 0.2221
## Nivel.Socieconomico.rec3-BAJO 0.6826 0.3126
## Nivel.Educativo.de.la.Madre1.rec2-High School 0.4660 0.3086
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.7012 0.3011
## Tipo.de.Escuela.rec2-Publica -0.2266 0.2662
## AtenOdonto2.rec2-publica 0.7478 0.2106
## AtenOdonto2.rec3-nunca fue al dentista 0.5470 0.2711
## FrCepDenti.4.rec2-1 vez al dias -0.6070 0.3135
## FrCepDenti.4.rec3-2 veces al dia -0.5406 0.2792
## FrCepDenti.4.rec4-3 veces al dia o mas -0.7100 0.2748
## GolosinasB.rec2-A veces -0.3771 0.2998
## GolosinasB.rec3-Nunca o raramente -0.2502 0.3098
## Masas.DulcesB.rec2-A veces -0.0905 0.2059
## Masas.DulcesB.rec3-Nunca o raramente -0.6474 0.4117
## t value Pr(>|t|)
## (Intercept) 1.15 0.2599
## Nivel.Socieconomico.rec2-MEDIO 2.34 0.0274 *
## Nivel.Socieconomico.rec3-BAJO 2.18 0.0386 *
## Nivel.Educativo.de.la.Madre1.rec2-High School 1.51 0.1436
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 2.33 0.0283 *
## Tipo.de.Escuela.rec2-Publica -0.85 0.4028
## AtenOdonto2.rec2-publica 3.55 0.0016 **
## AtenOdonto2.rec3-nunca fue al dentista 2.02 0.0545 .
## FrCepDenti.4.rec2-1 vez al dias -1.94 0.0642 .
## FrCepDenti.4.rec3-2 veces al dia -1.94 0.0642 .
## FrCepDenti.4.rec4-3 veces al dia o mas -2.58 0.0160 *
## GolosinasB.rec2-A veces -1.26 0.2201
## GolosinasB.rec3-Nunca o raramente -0.81 0.4271
## Masas.DulcesB.rec2-A veces -0.44 0.6642
## Masas.DulcesB.rec3-Nunca o raramente -1.57 0.1284
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for quasibinomial family taken to be 0.8987)
##
## Number of Fisher Scoring iterations: 4
confint(modelomulti.bin3)
## 2.5 % 97.5 %
## (Intercept) -0.38897 1.500103
## Nivel.Socieconomico.rec2-MEDIO 0.08508 0.955532
## Nivel.Socieconomico.rec3-BAJO 0.06996 1.295241
## Nivel.Educativo.de.la.Madre1.rec2-High School -0.13890 1.070852
## Nivel.Educativo.de.la.Madre1.rec3-Elementary School 0.11095 1.291409
## Tipo.de.Escuela.rec2-Publica -0.74845 0.295213
## AtenOdonto2.rec2-publica 0.33491 1.160609
## AtenOdonto2.rec3-nunca fue al dentista 0.01569 1.078285
## FrCepDenti.4.rec2-1 vez al dias -1.22140 0.007365
## FrCepDenti.4.rec3-2 veces al dia -1.08776 0.006617
## FrCepDenti.4.rec4-3 veces al dia o mas -1.24866 -0.171411
## GolosinasB.rec2-A veces -0.96466 0.210512
## GolosinasB.rec3-Nunca o raramente -0.85745 0.357137
## Masas.DulcesB.rec2-A veces -0.49402 0.313101
## Masas.DulcesB.rec3-Nunca o raramente -1.45429 0.159504