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 _______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
