[R-br] Erro - hoslem.test

[Luciane Maria Pilotto]

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