Olá,
estou realizando os ajustes de um modelo de regressão logistica ordinal e
testando se há violação do pressuposto de ordinalidade. Para isto, dicotomizei
o desfecho e rodei uma regressão logística normal, como descrito nos comandos
abaixo. Ao rodar o teste do Hosmer e Lemeshow (hoslem.test) para os resíduos dá
erro.
Agradeço quem puder ajudar!
########################################################################################
tp1 <- glm(q131[getario == "15 a 19 anos"] ~ q11[getario == "15 a 19 anos"] +
q10[getario == "15 a 19 anos"]
+rendacat[getario == "15 a 19 anos"]+q12[getario == "15 a 19
anos"] + q08[getario == "15 a 19 anos"]
+q06[getario == "15 a 19 anos"] + q14[getario == "15 a 19 anos"],
family = binomial(link = "logit"), data=id, na.action=na.exclude)
summary(tp1)
hoslem.test(tp1$q131, fitted(tp1), g=10)
Error in quantile.default(yhat, probs = seq(0, 1, 1/g)) :
missing values and NaN's not allowed if 'na.rm' is FALSE
###################################################
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.
