Premium
Variational methods for fitting complex Bayesian mixed effects models to health data
Author(s) -
Lee Cathy Yuen Yi,
Wand Matt P.
Publication year - 2015
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.6737
Subject(s) - computer science , markov chain monte carlo , bayesian inference , inference , bayes' theorem , bayesian probability , context (archaeology) , algorithm , econometrics , statistics , mathematics , artificial intelligence , paleontology , biology
We consider approximate inference methods for Bayesian inference to longitudinal and multilevel data within the context of health science studies. The complexity of these grouped data often necessitates the use of sophisticated statistical models. However, the large size of these data can pose significant challenges for model fitting in terms of computational speed and memory storage. Our methodology is motivated by a study that examines trends in cesarean section rates in the largest state of Australia, New South Wales, between 1994 and 2010. We propose a group‐specific curve model that encapsulates the complex nonlinear features of the overall and hospital‐specific trends in cesarean section rates while taking into account hospital variability over time. We use penalized spline‐based smooth functions that represent trends and implement a fully mean field variational Bayes approach to model fitting. Our mean field variational Bayes algorithms allow a fast (up to the order of thousands) and streamlined analytical approximate inference for complex mixed effects models, with minor degradation in accuracy compared with the standard Markov chain Monte Carlo methods. Copyright © 2015 John Wiley & Sons, Ltd.