The cumulative logistic regression model is commonly applied for ordinal outcomes in the medical literature.
On the other hand, this fascinating article describes the sequential logistic regression model:
For many ordinal variables, the assumption of a single underlying continuous variable, as in cumulative models, may not be appropriate. If the response can be understood as being the result of a sequential process, such that a higher response category is possible only after all lower categories are achieved, the sequential model proposed by Tutz (1990) is usually appropriate.
Sequential models assume that for every category k there is a latent continuous variable Yk that determines the transition between the kth and the k + 1th category.
I can imagine multiple clinical outcomes suiting these definitions, e.g.:
I didn’t read the paper but my first reaction is the sequential model may be a reinvention of the forward continuation ratio ordinal logistic model, which is a discrete hazard model. A comprehensive case study appears in RMS.
Category-specific effects requires a huge sample size as this is just the polytomous logistic regression model.
The only cohesive way to borrow information across categories is to use a Bayesian prior for the amount of borrowing, e.g., in a partial proportional odds model to specify priors for the departures from proportional odds as discussed in the first link at COVID-19.
Regarding the text below, please correct me if I misunderstood it: with a constrained partial PO model, one can not only estimate the cumulative OR as in a PO model but also a OR specific to one ordinal score (e.g., death). Would this be analogous to estimate a single category-specific effect?
From my quick experience with the Bayesian PPO model, the SD for the PPO term really determines the amount of shrinkage (e.g., SD = 5 vs. SD = 1). Can I find guidance on how to pick the exact SD value anywhere? This vignette is fascinating but doesn’t discuss this degree of detail.
I wonder if a prior predictive check approach would help in this case. rmsb::blrm() doesn’t seem to support prior predictive check, though. Please correct me if I’m mistaken.
Right it doesn’t support that directly. And yes we need some work that helps us select the amount of skepticism in the departures from PO. It’s not a very difficult problem except to put it into a clinical context. It’s easy if you stick with restrictions on ratios of odds ratios.