RMS Semiparametric Ordinal Longitudinal Model

Regression Modeling Strategies: Semiparametric Ordinal Longitudinal Model

This is the 22nd of several connected topics organized around chapters in Regression Modeling Strategies. This topic is for a chapter that is not in the book but is in the course notes. The purposes of these topics are to introduce key concepts in the chapter and to provide a place for questions, answers, and discussion around the chapter’s topics.

Overview | Course Notes

A common cause of disappointment (e.g., uninformative nulls) is pursuing low-information (insensitive) outcomes. Thoughtful effort given to understanding and choosing high-resolution high-information Y will likely improve PTS.

A high-resolution high-information Y can flexibly accommodate the timing and severity of a variety of outcomes (terminal events, non-terminal events, and recurrent events); and the more levels of Y the better (fharrell.com/post/ordinal-info). The longitudinal ordinal model is a general and flexible way to capture severity and timing of outcomes.

The proportional odds longitudinal ordinal logistic model with covariate adjustment is recommended (the Markov model better still). With this ordinal model there is no assumption about Y distribution, and random effects (intercepts) handle intra-patient correlation.

The proportional odds ordinal logistic model can estimate the probability that Y=y or worse as a function of time and treatment. This modeling approach provides estimates of efficacy for individual patients by addressing the fundamental clinical question: ‘If I compared two patients who have the same baseline variables but were given different treatments, by how much better should I expect the outcome to be with treatment B instead of treatment A?

With this ordinal longitudinal model one can obtain a variety of estimates: such as time until a condition, and, expected time in state. The ordinal model does assume proportional odds but the partial proportional odds model relaxes this.

The model provides a correct basis for analysis of heterogeneity of treatment effect.

Bayesian partial proportional odds model, moreover, can compute more complex probabilities of special interest, such as the probability that the treatment affects mortality differently than it affects nonfatal outcomes.

Additional links


Q&A From May 2021 Course

  1. Where do you think is the best place to start learning about Bayes, coming from a frequentist perspective? I know McElreath’s course is highly recommended, but I don’t think it ‘tells you’ the parallel approaches from the frequentist world - is there an (intro level) course or resource that does? dgl-I am working my way through the new Gelman book (Regression and Other Stories, Gelman, et al, (2021).), and am impressed with it. I recommend it. It has a balanced approach. fh-Kruschke does a lot of side-by-side Bayesian/frequentist analyses.
  2. Is there any limit to the number of states that the Bayesian Markov model can handle? Do you know if the models fitted using the Bayesian approach are comparable with the non-Bayesian Markov models fitted in the msm() package? Great questions. msm does not handle ordinal states, so every state is its own category and needs its own large sample size. With ordinal states there is no limit to the number of states as long as the proportional odds assumption is reasonably satisfied. You just need an good overall sample size.