Bayesian inference for generalized linear mixed models with predictors subject to detection limits: an approach that leverages information from auxiliary variables

• Published in: Statistics in medicine Vol. 35; no. 10; pp. 1689 - 1705
• Main Authors: Yue, Yu Ryan, Wang, Xiao-Feng
• Format: Journal Article
• Language: English
• Published: England Blackwell Publishing Ltd 10.05.2016; Wiley Subscription Services, Inc
• Subjects: Algorithms; auxiliary variables; Bayes Theorem; Bayesian analysis; Bayesian modeling; censored predictors; Computer Simulation; Critical Illness; detection limits; Female; generalized linear mixed models; Hospital Mortality; Humans; Length of Stay - statistics & numerical data; Linear Models; Male; Markov analysis; Markov chain Monte Carlo; Middle Aged; Monte Carlo Method; Monte Carlo simulation; Multicenter Studies as Topic; Ohio; Retrospective Studies; Statistical methods; Variables; Vitamin D; Vitamin D Deficiency - epidemiology; auxiliary variables; detection limits; generalized linear mixed models; Bayesian modeling; Markov chain Monte Carlo; censored predictors
• ISSN: 0277-6715; 1097-0258; 1097-0258
• DOI: 10.1002/sim.6830

This paper is motivated from a retrospective study of the impact of vitamin D deficiency on the clinical outcomes for critically ill patients in multi‐center critical care units. The primary predictors of interest, vitamin D2 and D3 levels, are censored at a known detection limit. Within the context of generalized linear mixed models, we investigate statistical methods to handle multiple censored predictors in the presence of auxiliary variables. A Bayesian joint modeling approach is proposed to fit the complex heterogeneous multi‐center data, in which the data information is fully used to estimate parameters of interest. Efficient Monte Carlo Markov chain algorithms are specifically developed depending on the nature of the response. Simulation studies demonstrate the outperformance of the proposed Bayesian approach over other existing methods. An application to the data set from the vitamin D deficiency study is presented. Possible extensions of the method regarding the absence of auxiliary variables, semiparametric models, as well as the type of censoring are also discussed. Copyright © 2015 John Wiley & Sons, Ltd.