[R-br] ajuste Regressão logística com o uso do cluster

[Luciane Maria Pilotto]

Bom dia,

gostaria de ajuda para realizar o ajuste de um modelo de regressão logística usando cluster. Como estimar a variância robusta ou o odds proporcional?


#Utilizei os seguintes comandos: 

dstrat<-svydesign(ids=~1, probs=cod_mun, strata = NULL, variable = NULL, fpc=NULL,
                  data = id15, nest = FALSE, weights=NULL,pps=FALSE,method="logistic")
m15<-svyolr(q13~q11,design=dstrat, method="logistic") #modelo bruto
#m156<-svyolr(q13~q11+q06+q10+q08+q12+q14,design=dstrat, method="logistic") #modelo ajustado
(ctable <- coef(summary(m15)))
p <- pnorm(abs(ctable[, "t value"]), lower.tail = FALSE) * 2
(ctable <- cbind(ctable, "p value" = p))
(ci <- confint(m15)) # default method gives profiled CIs
exp(coef(m15))
exp(cbind(OR = coef(m15), ci))

regTermTest(m15, ~q11, method="Wald")

LogRegR2(m15) #erro
Error in if (!(model$family$family == "binomial" && (model$family$link ==  : 
  missing value where TRUE/FALSE needed

Atenciosamente,
Luciane Pilotto