2
$\begingroup$

I have done an ordinal regression using lrm from the rms package.

load(url("http://www.produnis.de/R/NEXT-Q0-OR.RData"))
require(rms)
mydata <- NEXT.or
fit <- lrm(QuitJob~QD,data=mydata);fit

I would like to predict the rank in "QuitJob" (1=won't quit ... 5=will quit very soon) by giving a value for "QD" (how much work to do). For example, if "QD" is 78, what rank of "QuitJob" is likely?

If I type in

predict.lrm(fit, data.frame(QD=78))

I get

        1 
0.3191933 

Do I understand it right, that the probability of being in rank 1 is 31,92% ?! If so, how can I get the probabilities for ranks 2-5 ?

$\endgroup$
  • 4
    $\begingroup$ predict(fit, data.frame(QD=78), type="fitted.ind") gives predicted class probabilities, what you got is the value of the linear predictor, see this answer. $\endgroup$ – caracal Jul 24 '13 at 8:39
  • $\begingroup$ thx, caracal... If you post your commend as an answer, I would mark it as "accepted" $\endgroup$ – Produnis Jul 24 '13 at 8:45

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.