[R] help - hoslem.test

[Luciane Maria Pilotto]

load("id3.rda")
attach(id3)

#transformando q13 em binária
q131<-ifelse(q13==1,1,ifelse(q13==2,2,ifelse(q13==3,2, ifelse(q13==4,2,ifelse(q13==5,2,NA)))))
id3<-cbind(id3,q131)
id3$q131 <- as.factor(id3$q131)

tp1 <- glm(q131 ~ q11 + q10+q12+edcat + q08+q06+ q14, family = binomial(link = "logit"), data=id3) 
tp1

library(ResourceSelection)
hoslem.test(tp1$q131, fitted(tp1), g=10)

dataframe: https://www.dropbox.com/s/9qrdf4mhd6tzypi/id3.rda?dl=0


__________________________________________________
Luciane Maria Pilotto
Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS 
NDE Odontologia - UNIVATES
Telefone: (51) 84512344

--------------------------------------------
Em qui, 30/4/15, John Kane <jrkrideau em inbox.com> escreveu:

 Assunto: RE: [R] help - hoslem.test

.org
 Data: Quinta-feira, 30 de Abril de 2015, 7:52

 http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example

 John Kane
 Kingston ON Canada


 > -----Original Message-----

 > Sent: Wed, 29 Apr 2015 18:45:26 -0700
 > To: r-help em r-project.org
 > Subject: [R] help - hoslem.test
 > 
 > Hello,
 > 
 > I'm working with
 ordinal logistic regression model (polr) and would like
 > to test the proportional odds assumption.
 For this, I ran the binary
 > logistic
 regressions with varying cutpoints on the dependent
 variable, as
 > described in the following
 commands. When running the test of Hosmer and
 > Lemeshow (hoslem.test) for residuals gives
 error.
 > 
 > Thanks,
 > Luciane
 > 
 >
 ______________________________________________
 > R-help em r-project.org
 mailing list -- To UNSUBSCRIBE and more, see
 > https://stat.ethz.ch/mailman/listinfo/r-help
 > PLEASE do read the posting guide
 > http://www.R-project.org/posting-guide.html
 > and provide commented, minimal,
 self-contained, reproducible code.

 ____________________________________________________________
 FREE ONLINE PHOTOSHARING - Share your photos
[[elided Yahoo spam]]
 Visit
 http://www.inbox.com/photosharing to
 find out more!