Estan las recodificaciones planteadas en el documento de 28 de julio Estan tratadas como variables de conteo por lo cual se hacen modelos de regresion de poisson
domingo 11 de AGOSTO de 2013
# domingo 11 de AGOSTO de 2013
library(survey)
## Attaching package: 'survey'
## The following object(s) are masked from 'package:graphics':
##
## dotchart
options(OutDec = ",")
library(survey)
load("~/Dropbox/odontologia/maestria licet/julio_2013/datos_licet_25072013.RData")
library(car)
## Loading required package: MASS
## Loading required package: nnet
model_multi6.poi <- svyglm(bewes ~ Sexo.rec + Nivel.Socieconomico.rec + FrCepDenti.rec +
UsoDentifrico3.rec + JugFrutas.rec + Bebidas_energizantes.rec + BuchTragar.rec +
FormBeber.rec, diseniopost1, family = quasipoisson())
names(model_multi6.poi)
## [1] "coefficients" "residuals" "fitted.values"
## [4] "effects" "R" "rank"
## [7] "qr" "family" "linear.predictors"
## [10] "deviance" "aic" "null.deviance"
## [13] "iter" "weights" "prior.weights"
## [16] "df.residual" "df.null" "y"
## [19] "converged" "boundary" "model"
## [22] "na.action" "call" "formula"
## [25] "terms" "data" "offset"
## [28] "control" "method" "contrasts"
## [31] "xlevels" "naive.cov" "cov.unscaled"
## [34] "survey.design"
summary(model_multi6.poi)
##
## Call:
## svyglm(formula = bewes ~ Sexo.rec + Nivel.Socieconomico.rec +
## FrCepDenti.rec + UsoDentifrico3.rec + JugFrutas.rec + Bebidas_energizantes.rec +
## BuchTragar.rec + FormBeber.rec, diseniopost1, family = quasipoisson())
##
## Survey design:
## postStratify(disenio1, ~categor.rec + Sexo, tabla.pob)
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 1,9906 0,1863 10,69
## Sexo.rec2-M 0,0507 0,0779 0,65
## Nivel.Socieconomico.rec2-MEDIO 0,0477 0,0974 0,49
## Nivel.Socieconomico.rec3-ALTO -0,3251 0,1601 -2,03
## FrCepDenti.rec2- 2 veces al dia 0,1242 0,1900 0,65
## FrCepDenti.rec3-3 o mas veces al dia 0,2926 0,2308 1,27
## UsoDentifrico3.rec2-No -0,9143 0,5696 -1,61
## JugFrutas.rec3-Nunca o raramente -0,3254 0,3279 -0,99
## JugFrutas.recMas de 3 veces al dia 0,2336 0,0710 3,29
## Bebidas_energizantes.rec2-Todos los dias -0,0303 0,1176 -0,26
## Bebidas_energizantes.recMas de 3 veces al dia -1,0133 0,6536 -1,55
## BuchTragar.rec2-Yes 0,2828 0,1177 2,40
## FormBeber.rec2-con sorbete 0,2892 0,1096 2,64
## FormBeber.rec3-botella -0,1436 0,2178 -0,66
## Pr(>|t|)
## (Intercept) 6e-10 ***
## Sexo.rec2-M 0,5225
## Nivel.Socieconomico.rec2-MEDIO 0,6296
## Nivel.Socieconomico.rec3-ALTO 0,0551 .
## FrCepDenti.rec2- 2 veces al dia 0,5206
## FrCepDenti.rec3-3 o mas veces al dia 0,2187
## UsoDentifrico3.rec2-No 0,1234
## JugFrutas.rec3-Nunca o raramente 0,3323
## JugFrutas.recMas de 3 veces al dia 0,0035 **
## Bebidas_energizantes.rec2-Todos los dias 0,7993
## Bebidas_energizantes.recMas de 3 veces al dia 0,1360
## BuchTragar.rec2-Yes 0,0256 *
## FormBeber.rec2-con sorbete 0,0154 *
## FormBeber.rec3-botella 0,5168
## ---
## Signif. codes: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1
##
## (Dispersion parameter for quasipoisson family taken to be 9,044)
##
## Number of Fisher Scoring iterations: 6
confint(model_multi6.poi)
## 2,5 % 97,5 %
## (Intercept) 1,62545 2,35570
## Sexo.rec2-M -0,10199 0,20330
## Nivel.Socieconomico.rec2-MEDIO -0,14325 0,23861
## Nivel.Socieconomico.rec3-ALTO -0,63883 -0,01138
## FrCepDenti.rec2- 2 veces al dia -0,24831 0,49664
## FrCepDenti.rec3-3 o mas veces al dia -0,15974 0,74493
## UsoDentifrico3.rec2-No -2,03059 0,20205
## JugFrutas.rec3-Nunca o raramente -0,96801 0,31720
## JugFrutas.recMas de 3 veces al dia 0,09442 0,37286
## Bebidas_energizantes.rec2-Todos los dias -0,26076 0,20020
## Bebidas_energizantes.recMas de 3 veces al dia -2,29432 0,26778
## BuchTragar.rec2-Yes 0,05214 0,51336
## FormBeber.rec2-con sorbete 0,07441 0,50408
## FormBeber.rec3-botella -0,57057 0,28331