A Unified Framework for Fitting Bayesian Semiparametric Models to Arbitrarily Censored Survival Data, Including Spatially Referenced Data. Issue 522 (3rd April 2018)

• Record Type: Journal Article
• Authors: Zhou, Haiming; Hanson, Timothy
• Is Part Of: Journal of the American Statistical Association. Volume 113:Issue 522(2018)
• Page Start: 571
• Page End: 581
• Publication Date: 2018-04-03
• Subjects: Bernstein polynomial -- Interval censored data -- Spatial frailty models -- Variable selection; Statistics -- Periodicals; Statistics -- Periodicals; Statistiques -- Périodiques; États-Unis -- Statistiques -- Périodiques; 519.5
• DOI: 10.1080/01621459.2017.1356316
• Languages: English
• ISSNs: 0162-1459

Abstract

ABSTRACT: A comprehensive, unified approach to modeling arbitrarily censored spatial survival data is presented for the three most commonly used semiparametric models: proportional hazards, proportional odds, and accelerated failure time. Unlike many other approaches, all manner of censored survival times are simultaneously accommodated including uncensored, interval censored, current-status, left and right censored, and mixtures of these. Left-truncated data are also accommodated leading to models for time-dependent covariates. Both georeferenced (location exactly observed) and areally observed (location known up to a geographic unit such as a county) spatial locations are handled; formal variable selection makes model selection especially easy. Model fit is assessed with conditional Cox–Snell residual plots, and model choice is carried out via log pseudo marginal likelihood (LPML) and deviance information criterion (DIC). Baseline survival is modeled with a novel transformed Bernstein polynomial prior. All models are fit via a new function which calls efficient compiled C++ in the R packagespBayesSurv . The methodology is broadly illustrated with simulations and real data applications. An important finding is that proportional odds and accelerated failure time models often fit significantly better than the commonly used proportional hazards model. Supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplement.