Frailty modeling via the empirical Bayes-Hastings sampler

详细信息	查看全文 | 推荐本文 |

• 作者：Richard A. Levine^a ; ^{ralevine@sciences.sdsu.edu} ; Juanjuan Fan^a ; ^{jjfan@sciences.sdsu.edu} ; Pamela Ohman Strickland^b ; ^{ohmanpa@umdnj.edu} ; Shaban Demirel^c ; ^{sdemirel@deverseye.org}
• 关键词：Multivariate survival analysis ; Nonparametric Pó ; lya tree prior ; Gibbs sampler ; Metropolis&ndash ; Hastings sampler ; Goodness of fit ; Glaucoma and ophthalmology data
• 刊名：Computational Statistics and Data Analysis
• 出版年：2012
• 期刊代码：19_01679473
• 类别：cp
• 出版时间：June, 2012
• 卷：56
• 期：6
• 页码：1303-1318
• 文件大小：440 K

摘要

Studies of ocular disease and analyses of time to disease onset are complicated by the correlation expected between the two eyes from a single patient. We overcome these statistical modeling challenges through a nonparametric Bayesian frailty model. While this model suggests itself as a natural one for such complex data structures, model fitting routines become overwhelmingly complicated and computationally intensive given the nonparametric form assumed for the frailty distribution and baseline hazard function. We consider empirical Bayesian methods to alleviate these difficulties through a routine that iterates between frequentist, data-driven estimation of the cumulative baseline hazard and Markov chain Monte Carlo estimation of the frailty and regression coefficients. We show both in theory and through simulation that this approach yields consistent estimators of the parameters of interest. We then apply the method to the short-wave automated perimetry (SWAP) data set to study risk factors of glaucomatous visual field deficits.