Modelling spatially correlated survival data for individuals with multiple cancers

Abstract

Epidemiologists and biostatisticians investigating spatial variation in diseases are often interested in estimating spatial effects in survival data, where patients are monitored until their time to failure (for example, death, relapse). Spatial variation in survival patterns often reveals underlying lurking factors that could assist public health professionals in their decision–making process to identify regions requiring attention. The Surveillance Epidemiology and End Results (SEER) database of the National Cancer Institute provides a fairly sophisticated platform for exploring novel approaches in modelling cancer survival, particularly with models accounting for spatial clustering and variation. Modelling survival data for patients with multiple cancers poses unique challenges in itself and in capturing the spatial associations of the different cancers. This paper develops the Bayesian hierarchical survival models for capturing spatial patterns within the framework of proportional hazard. Spatial variation is introduced in the form of county–cancer level frailties. The baseline hazard function is modelled semiparametrically using mixtures of beta distributions. We illustrate with data from the SEER database, perform model checking and comparison among competing models, and discuss implementation issues.

Keywords

Bayesian hierarchical models frailty models Markov Chain Monte Carlo (MCMC)mixture of beta functions spatial association survival modelling

Get full access to this article

View all access options for this article.

References

Agarwal D and Gelfand AE (2005) Slice Gibbs sampling for simulation based fitting of spatial models . Statistics and Computing, 15, 61–69 .

Banerjee S and Carlin BP (2003) Semi–parametric spatiotemporal analysis . Environmetrics, 14, 523–35 .

Banerjee S , Carlin BP and Gelfand AE (2004). Hierarchical Modeling and Analysis for Spatial Data. Boca Raton: Chapman and Hall/CRC .

Banerjee S and Dey DK (2005) Semiparametric proportional odds model for spatially correlated survival data . Lifetime Data Analysis, 11, 175–91 .

Banerjee S , Wall M and Carlin BP (2003) Frailty modelling for spatially correlated survival data with application to infant mortality in Minnesota . Biostatistics, 4, 123–42 .

Biggeri A , Marchi M , Lagazio C , Martuzzi M and Böhning D (2000) Nonparametric maximum likelihood estimators for disease mapping . Statistics in Medicine, 19, 2539–554 .

Carlin BP and Banerjee S (2003) Hierarchical multivariate CAR models for spatio–temporally correlated survival data. In Bernardo JM , Bayarri MJ , Berger JO , Dawid AP et al . eds. Bayesian Statistics 7. Oxford: Oxford University Press , 45–64.

Celeux G , Forbes F , Robert CP and Titterington DM (2006) Deviance information criteria for missing data models (with discussion). Bayesian Analysis, 1, 651–706 .

Chen MH , Shao QM and Ibrahim JG (2000) Monte Carlo Methods in Bayesian Computation. New York: Springer–Verlag .

10.

Clayton DG (1995) Some approaches to the analysis of recurrent event data . Statistics in Medical Research, 3, 244–62 .

11.

Gelfand AE and Dey DK (1994) Bayesian model choice: assymptotics and exact calculations . Journal of Royal Statistical Society, Series B, 56, 501–14 .

12.

Gelfand AE and Mallick BK (1995) Bayesian analysis of proportional hazards models built from monotone functions . Biometrics, 51, 843–52 .

13.

Gelfand AE and Vounatsou P (2002) Proper multivariate conditional autoregressive models for spatial data analysis . Biostatistics, 4, 11–25 .

14.

Gilks WR and Wild P (1992) Adaptive rejection sampling for Gibbs sampling . Journal of Royal Statistical Society, Series C, 41, 337–48 .

15.

Harville DA (1999) Matrix Algebra from a Statistician's Perspective. New York: Springer–Verlag .

16.

Heinävaara S (2003) Modelling Survival of Patients with Multiple Cancers. Ph.D. thesis, University of Helsinki. The Finnish Statistical Society , Statistical Research Reports No. 18.

17.

Held L and Best NG (2001) A shared component model for detecting joint and selective clustering of two diseases . Journal of the Royal Statistical Society, Series A, 164, 73–85 .

18.

Held L , Natario I , Fenton S , Rue H and Becker N (2005) Towards joint disease mapping . Statistical Methods in Medical Research, 14, 61–82 .

19.

Ibrahim JG , Chen MH and Sinha D (2001) Bayesian Survival Analysis. New York: Springer–Verlag .

20.

Jin X , Carlin BP and Banerjee S (2005) Generalized hierarchical multivariate CAR models for areal data . Biometrics, 61, 950–61 .

21.

Kelderman H (1984) Loglinear Rasch model tests . Psychometrika, 49, 223–45 .

22.

Li Y and Ryan L (2002) Modeling spatial survival data using semiparametric frailty models . Biometrics, 58, 287–97 .

23.

Lichstein JW , Simons TR , Shriner SA and Franzreb KE (2002) Spatial autocorrelation and autoregressive models in ecology . Ecological Monographs, 72, 3, 445–63 .

24.

Mardia KV (1998) Multidimensional multivariate Gaussian Markov random fields with application to image processing . Journal of Multivariate Analysis, 24, 45–64 .

25.

Ramsay T , Burnett R and Krewski D (2003) Exploring bias in a generalized additive model for spatial air pollution data . Environmental Health Perspectives, 111, 1283–88 .

26.

Ries LAG , Eisner MP , Kosary CL , Hankey BF et al. , (eds) (2005) SEER Cancer Statistics Review, 1975–2002, National Cancer Institute. Bethesda, MD , http://seer.cancer.gov/csr/1975_2002/, based on November 2004 SEER data submission, posted to the SEER website 2005.

27.

Sankila R and Hakulinen T (1998) Survival of patients with colorectal carcinoma: effect of prior breast cancer . Journal of the National Cancer Institute, 90, 6365.

28.

Smith BA (2005) Bayesian Output Analysis Program (BOA) Version 1.1 User's Manual. http://www.public–health.uiowa.edu/boa/.

29.

Spiegelhalter DJ , Best N , Carlin BP and van der Linde A (2002) Bayesian measures of model complexity and fit (with discussion). Journal of the Royal Statistical Society, Series B, 64, 583–639 .

30.

Surveillance, Epidemiology, and End Results (SEER) Program (www.seer.cancer.gov) SEER*Stat Database: Incidence – SEER 11 Regs+AK Public–Use, Nov 2003 Sub (1973–2001 varying), National Cancer Institute, DCCPS, Surveillance Research Program, Cancer Statistics Branch, released April 2004, based on the November 2003 submission.

31.

Turechek WW and Madden LV (2002) A generalized linear modeling approach for characterizing disease incidence in spatial hierarchy . Phytopathology, 93, 458–66 .