Skip to content
Prev 59615 / 398502 Next

adjacent category model in ordinal regression

Hello,

I suspect that vglm would do the job for the logit link.  Try :

     data(pneumo)
     pneumo$let = log(pneumo$exposure.time)
     vglm(cbind(normal, mild, severe) ~ let, multinomial, pneumo)

     vglm(cbind(normal, mild, severe) ~ let, acat, pneumo)

This is an adjacent category model with logit link, although the
help page speaks about log link (log of odds = logit). 
But I think the model is, in fact, the same as the previous one.
(same deviance ...). It seems that only the parametrisation is different.

Hope this help.
Claire Chabanet