Premium
Flexible parametric models for random‐effects distributions
Author(s) -
Lee Katherine J.,
Thompson Simon G.
Publication year - 2007
Publication title -
statistics in medicine
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.996
H-Index - 183
eISSN - 1097-0258
pISSN - 0277-6715
DOI - 10.1002/sim.2897
Subject(s) - random effects model , markov chain monte carlo , bivariate analysis , deviance information criterion , computer science , econometrics , markov chain , parametric statistics , mathematics , statistics , monte carlo method , meta analysis , medicine
It is commonly assumed that random effects in hierarchical models follow a normal distribution. This can be extremely restrictive in practice. We explore the use of more flexible alternatives for this assumption, namely the t distribution, and skew extensions to the normal and t distributions, implemented using Markov Chain Monte Carlo methods. Models are compared in terms of parameter estimates, deviance information criteria, and predictive distributions. These methods are applied to examples in meta‐analysis and health‐professional variation, where the distribution of the random effects is of direct interest. The results highlight the importance of allowing for potential skewing and heavy tails in random‐effects distributions, especially when estimating a predictive distribution. We describe the extension of these random‐effects models to the bivariate case, with application to a meta‐analysis examining the relationship between treatment effect and baseline response. We conclude that inferences regarding the random effects can crucially depend on the assumptions made and recommend using a distribution, such as those suggested here, which is more flexible than the normal. Copyright © 2007 John Wiley & Sons, Ltd.