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 <jrkrid...@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@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@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!
______________________________________________
R-help@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.
