The use of restricted cubic splines to approximate complex hazard functions in the analysis of time-to-event data: a simulation study

Journal of Statistical Computation and Simulation

Volumn 85, Issue 4, 2015, Pages 777-793

  Rutherford, Mark J ^a   Crowther, Michael J ^a   Lambert, Paul C ^a,b  

^a UNIVERSITY OF LEICESTER   (United Kingdom)

^b KAROLINSKA INSTITUTE   (Sweden)

Author keywords

flexible parametric models; restricted cubic splines; simulation; survival analysis

Indexed keywords

EID: 84914665589     PISSN: 00949655     EISSN: 15635163     Source Type: Journal    
DOI: 10.1080/00949655.2013.845890     Document Type: Article

References (33)

1. Royston P, Parmar MKB. Flexible parametric proportional-hazards and proportional-odds models for censored survival data, with application to prognostic modelling and estimation of treatment effects. Stat Med. 2002;21:2175–2197. doi: 10.1002/sim.1203
2. Royston P, Lambert PC. Flexible parametric survival analysis using Stata: beyond the Cox model. College Station, TX: Stata Press; 2011.
3. Lambert PC, Holmberg L, Sandin F, Bray F, Linklater KM, Purushotham A, Robinson D, Møller H. Quantifying differences in breast cancer survival between England and Norway. Cancer Epidemiol. 2011;35:526–533. doi: 10.1016/j.canep.2011.04.003
4. Albain KS, Barlow WE, Shak S, Hortobagyi GN, Livingston RB, Yeh I-T, Ravdin P, Bugarini R, Baehner FL, Davidson NE, Sledge GW, Winer EP, Hudis C, Ingle JN, Perez EA, Pritchard KI, Shepherd L, Gralow JR, Yoshizawa C, Allred DC, Osborne CK, Hayes DF. Prognostic and predictive value of the 21-gene recurrence score assay in postmenopausal women with node-positive, oestrogen-receptor-positive breast cancer on chemotherapy: a retrospective analysis of a randomised trial. Lancet Oncol. 2010;11:55–65. doi: 10.1016/S1470-2045(09)70314-6
5. Gore ME, Griffin CL, Hancock B, Patel PM, Pyle L, Aitchison M, James N, Oliver RT, Mardiak J, Hussain T, Sylvester R, Parmar MK, Royston P, Mulders, PF. Interferon alfa-2a versus combination therapy with interferon alfa-2a, interleukin-2, and fluorouracil in patients with untreated metastatic renal cell carcinoma (MRC RE04/EORTC GU 30012): an open-label randomised trial. Lancet. 2010;375:641–648. doi: 10.1016/S0140-6736(09)61921-8
6. Nelson CP, Lambert PC, Squire IB, Jones DR. Flexible parametric models for relative survival, with application in coronary heart disease. Stat Med. 2007;26:5486–5498. doi: 10.1002/sim.3064
7. Andersson T, Dickman P, Eloranta S, Lambert P. Estimating and modelling cure in population-based cancer studies within the framework of flexible parametric survival models. BMC Med Res Methodol. 2011;11:Article 96. doi: 10.1186/1471-2288-11-96
8. Durrelman S, Simon R. Flexible regression models with cubic splines. Stat Med. 1989;8:551–561. doi: 10.1002/sim.4780080504
9. Lambert PC, Royston P. Further Development of flexible parametric models for survival analysis. Stata J. 2009;9: 265–290.
10. Lambert PC, Dickman PW, Nelson CP, Royston P. Estimating the crude probability of death due to cancer and other causes using relative survival models. Stat Med. 2010;29:885–895. doi: 10.1002/sim.3762
11. Törner A, Dickman P, Duberg AS, Kristinsson S, Landgren O, Björkholm M, Svensson A. A method to visualize and adjust for selection bias in prevalent cohort studies. Am J Epidemiol. 2011;174(8):969–976. doi: 10.1093/aje/kwr211
12. Cox DR. Regression models and life-tables. J R Stat Soc Ser B (Methodological). 1972;34:187–220.
13. McLachlan G, McGiffin D. On the role of finite mixture models in survival analysis. Stat Methods Med Res. 1994;3:211–226. doi: 10.1177/096228029400300302
14. Akaike H. Information theory as an extension of the maximum likelihood principle. In: PetrovBN, CsakiF. editors. Second international symposium on information theory. Budapest: Akademiai Kiado; 1973. p. 267–281.
15. Gelfand AE, Dey DK. Bayesian model choice: asymptotics and exact calculations. J R Stat Soc Ser B (Methodological). 1994;56:501–514.
16. Govindarajulu US, Spiegelman D, Thurston SW, Ganguli B, Eisen EA. Comparing smoothing techniques in Cox models for exposure–response relationships. Stat Med. 2007;26:3735–3752. doi: 10.1002/sim.2848
17. Govindarajulu US, Malloy EJ, Ganguli B, Spiegelman D, Eisen EA. The comparison of alternative smoothing methods for fitting non-linear exposure–response relationships with Cox models in a simulation study. Int J Biostat. 2009;5:Article 2.
18. Rosenberg PS. Hazard function estimation using B-splines. Biometrics. 1995;51:874–887. doi: 10.2307/2532989
19. Molinari N, Durand JF, Sabatier R. Bounded optimal knots for regression splines. Comput Stat Data Anal. 2004;45:159–178. doi: 10.1016/S0167-9473(02)00343-2
20. Burnham KP, Anderson DR. Multimodel inference: understanding AIC and BIC in model selection. Sociol Methods Res. 2004;33:261–304. doi: 10.1177/0049124104268644
21. Weakliem DL. A critique of the Bayesian information criterion for model selection. Sociol Methods Res. 1999;27:359–397. doi: 10.1177/0049124199027003002
22. Bantis L, Tsimikas J, Georgiou S. Survival estimation through the cumulative hazard function with monotone natural cubic splines. Lifetime Data Anal. 2012;18(3):364–396. doi: 10.1007/s10985-012-9218-4
23. Kooperberg C, Stone CJ. Logspline density estimation for censored data. J Comput Graphical Stat. 1992;1:301–328.
24. Golub GH, Loanvan CF. Matrix computations. 3rd ed. Baltimore, MD: Johns Hopkins University Press; 1996.
25. Bender R, Augustin T, Blettner M. Generating survival times to simulate Cox proportional hazards models. Stat Med. 2005;24:1713–1723. doi: 10.1002/sim.2059
26. Coleman M, Babb P, Damiecki P, Grosclaude P, Honjo S, Jones J, Knerer G, Pitard A, Quinn M, Sloggett A, De Stavola B. Cancer survival trends in England and Wales, 1971–1995: deprivation and NHS region. Series SMPS No. 61. London: The Stationery Office; 1999.
27. Silverman B. Density estimation: for statistics and data analysis. London: Chapman and Hall; 1986.
28. Crowther M, Lambert P. Simulating complex survival data. The Stata J. 2012;12(4):674–687.
29. Crowther MJ, Lambert P. STMIX: Stata module to fit two-component parametric mixture survival models; Boston College, MA: Statistical Software Components, Boston College Department of Economics 2012. Available from: http://ideas.repec.org/c/boc/bocode/s457339.html
30. Gould WW. crc46: Better numerical derivatives and integrals. Stata Tech Bull. 1997;35:3–5.
31. StataCorp. Stata Statistical Software: Release 12. College Station, TX: StataCorp LP; 2011.
32. Cox C, Chu H, Schneider MF, Muñoz A. Parametric survival analysis and taxonomy of hazard functions for the generalized gamma distribution. Stat Med. 2007;26:4352–4374. doi: 10.1002/sim.2836
33. Demiris N, Lunn D, Sharples LD. Survival extrapolation using the poly-Weibull model. Stat Methods Med Res. 2011. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21937472