A model for overdispersed hierarchical ordinal data

• Published in: Statistical modelling Vol. 14; no. 5; pp. 399 - 415
• Main Authors: Ivanova, Anna, Molenberghs, Geert, Verbeke, Geert
• Format: Journal Article
• Language: English
• Published: New Delhi, India SAGE Publications 01.10.2014; Sage Publications Ltd
• Subjects: Case studies; Cholesterol; Diabetes; Low density lipoprotein; Maximum likelihood method; Normal distribution; Parameter estimation; Patients; Quality improvement; Statistical analysis; Studies; generalized linear mixed model; proportional odds model; overdispersion; beta distribution; maximum likelihood
• ISSN: 1471-082X; 1477-0342
• DOI: 10.1177/1471082X14522910

Non-Gaussian outcomes are frequently modelled using members of the exponential family. In particular, the Bernoulli model for binary data and the Poisson model for count data are well-known. Two reasons for extending this family are (1) the occurrence of overdispersion, implying that the variability in the data is not adequately described by the models, and (2) the incorporation of hierarchical structure in the data. These issues are routinely addressed separately, the first one through overdispersion models, the second one, for example, by means of random effects within the generalized linear mixed models framework. Molenberghs et al. (2007, 2010) introduced a so-called ‘combined model’ that simultaneously addresses both. In these and subsequent papers, a lot of attention was given to binary outcomes, counts, and time-to-event responses. While common in practice, ordinal data have not been studied from this angle. In this article, a model for ordinal repeated measures, subject to overdispersion, is formulated. It can be fitted without difficulty using standard statistical software. The model is exemplified using data from an epidemiological study in diabetic patients and using data from a clinical trial in psychiatric patients.