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

References

• 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
   PubMed Web of Science ®Google Scholar
• Royston P, Lambert PC. Flexible parametric survival analysis using Stata: beyond the Cox model. College Station, TX: Stata Press; 2011.
   Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• Durrelman S, Simon R. Flexible regression models with cubic splines. Stat Med. 1989;8:551–561. doi: 10.1002/sim.4780080504
   PubMed Web of Science ®Google Scholar
• Lambert PC, Royston P. Further Development of flexible parametric models for survival analysis. Stata J. 2009;9: 265–290.
   Web of Science ®Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• Cox DR. Regression models and life-tables. J R Stat Soc Ser B (Methodological). 1972;34:187–220.
   Web of Science ®Google Scholar
• 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
   PubMedGoogle Scholar
• 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.
   Google Scholar
• Gelfand AE, Dey DK. Bayesian model choice: asymptotics and exact calculations. J R Stat Soc Ser B (Methodological). 1994;56:501–514.
   Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• 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.
   Google Scholar
• Rosenberg PS. Hazard function estimation using B-splines. Biometrics. 1995;51:874–887. doi: 10.2307/2532989
   PubMed Web of Science ®Google Scholar
• 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
   Web of Science ®Google Scholar
• Burnham KP, Anderson DR. Multimodel inference: understanding AIC and BIC in model selection. Sociol Methods Res. 2004;33:261–304. doi: 10.1177/0049124104268644
   Web of Science ®Google Scholar
• Weakliem DL. A critique of the Bayesian information criterion for model selection. Sociol Methods Res. 1999;27:359–397. doi: 10.1177/0049124199027003002
   Web of Science ®Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• Kooperberg C, Stone CJ. Logspline density estimation for censored data. J Comput Graphical Stat. 1992;1:301–328.
   Google Scholar
• Golub GH, Loanvan CF. Matrix computations. 3rd ed. Baltimore, MD: Johns Hopkins University Press; 1996.
   Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• 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.
   Google Scholar
• Silverman B. Density estimation: for statistics and data analysis. London: Chapman and Hall; 1986.
   Google Scholar
• Crowther M, Lambert P. Simulating complex survival data. The Stata J. 2012;12(4):674–687.
   Web of Science ®Google Scholar
• 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
   Google Scholar
• Gould WW. crc46: Better numerical derivatives and integrals. Stata Tech Bull. 1997;35:3–5.
   Google Scholar
• StataCorp. Stata Statistical Software: Release 12. College Station, TX: StataCorp LP; 2011.
   Google Scholar
• 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
   PubMed Web of Science ®Google Scholar
• 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
   Google Scholar