Bayesian model choice in cumulative link ordinal regression models

Change log
McKinley, TJ 
Morters, M 
Wood, JLN 

The use of the proportional odds (PO) model for ordinal regression is ubiquitous in the literature. If the assumption of parallel lines does not hold for the data, then an alternative is to specify a non-proportional odds (NPO) model, where the regression parameters are allowed to vary depending on the level of the response. However, it is often difficult to fit these models, and challenges regarding model choice and fitting are further compounded if there are a large number of explanatory variables. We make two contributions towards tackling these issues: firstly, we develop a Bayesian method for fitting these models, that ensures the stochastic ordering conditions hold for an arbitrary finite range of the explanatory variables, allowing NPO models to be fitted to any observed data set. Secondly, we use reversible-jump Markov chain Monte Carlo to allow the model to choose between PO and NPO structures for each explanatory variable, and show how variable selection can be incorporated. These methods can be adapted for any monotonic increasing link functions. We illustrate the utility of these approaches on novel data from a longitudinal study of individual-level risk factors affecting body condition score in a dog population in Zenzele, South Africa.

Bayesian inference, ordinal regression, Markov chain Monte Carlo, reversible-jump, Bayesian model choice
Journal Title
Bayesian Analysis
Conference Name
Journal ISSN
Volume Title
Institute of Mathematical Statistics
Biotechnology and Biological Sciences Research Council (BB/I012192/1)
International Fund for Animal Welfare (IFAW) (unknown)
International Fund for Animal Welfare (IFAW) (unknown)
TJM is supported by Biotechnology and Biological Sciences Research Council grant number BB/I012192/1. MM is supported by a grant from the International Fund for Animal Welfare (IFAW) and the World Society for the Protection of Animals (WSPA), with additional support from the Charles Slater Fund and the Jowett Fund. JW is supported by the Alborada Trust and the RAPIDD program of the Science and Technology Directorate, Department of Homeland Security and the Fogarty International Centre.