Dear Heather,
You can make this test using the ordinal package. Here the function
clm fits cumulative link models where the ordinal logistic regression
model is a special case (using the logit link).
Let me illustrate how to test the parallel regression assumption for a
particular variable using
here's some code as an example hope it helps!
mod<-polr(vote~age+demsat+eusup+lrself+male+retnat+union+urban, data=dat)
summary(mod)
mod<-polr(vote~age+demsat+eusup+lrself+male+retnat+union+urban, data=dat)
levs<-levels(dat$vote)
tmpdat<-list()
for(i in 1:(nlevels(dat$vote)-1)){
tmpdat[[i]
Perhaps you should be asking whether such an algorithm exists, regardless of
whether it is already implemented in R. However, this is the wrong place to ask
such theory questions... your local statistics expert might know, or you could
ask on a statistics theory forum such as stats.stackexchange
Heather:
You are at Vanderbilt, whose statistics department under Frank Harrell
is a veritable bastion of R and statistical wisdom. I strongly
recommend that you take a stroll over there in the lovely spring
weather and seek their help. I can't imagine how you could do better
than that!
Cheers,
B
4 matches
Mail list logo
