![[PDF]](/pb-assets/Icons/adobe-pdf-icon-1498649896243.png){kind=icon}
Abstract
Abstract:
The evaluation of peritoneal dialysis (PD) programmes requires the use of statistical methods that suit the complexity of such programmes. Multi-state regression models taking competing risks into account are a good example of suitable approaches. In this work, multi-state structured additive regression (STAR) models combined with penalized splines (P-splines) are proposed to evaluate peritoneal dialysis programmes. These models are very flexible since they may consider smooth estimates of baseline transition intensities and the inclusion of time-varying and smooth covariate effects at each transition. A key issue in survival analysis is the quantification of the time-dependent predictive accuracy of a given regression model, which is typically assessed using receiver operating characteristic (ROC)’based methodologies. The main objective of the present study is to adapt the concept of time-dependent ROC curve, and their corresponding area under the curve (AUC), to a multi-state competing risks framework. All statistical methodologies discussed in this work were applied to PD survival data. Using a multi-state competing risks framework, this study explored the effects of major clinical covariates on survival such as age, sex, diabetes and previous renal replacement therapy. Such multi-state model was composed of one transient state (peritonitis) and several absorbing states (death, transfer to haemodialysis and renal transplantation). The application of STAR models combined with time-dependent ROC curves revealed important conclusions not previously reported in the nephrology literature when using standard statistical methodologies. For practical application, all the statistical methods proposed in this article were implemented in R and we wrote and made available a script named as NestedCompRisks.
Keywords Key words:time-dependent ROC curve, competing risks, multi-state models, peritoneal dialysis, STAR model, survival analysis
References
| Abraham, G, Kumar, V, Nayak, KS, Ravichandran, R, Srinivasan, G, Krishnamurthy, M, Prasath, AK, Kumar, S, Thiagarajan, T, Mathew, M, Lesley, N (2010) Predictors of long-term survival on peritoneal dialysis in south india: A multicenter study. Peritoneal Dialysis International, 30, 29–34. doi: 10.3747/pdi.2008.00028 Google Scholar | SAGE Journals | ISI | |
| Akritas, MG (1994) Nearest-neighbor estimation of a bivariate distribution under random censoring. Annals of Statistics, 22, 1299–327. doi: DOI 10.1214/aos/1176325630 Google Scholar | Crossref | ISI | |
| Andersen, PK, Keiding, N (2012) Interpretability and importance of functionals in competing risks and multistate models. Statistics in Medicine, 31, 1074–88. doi: 10.1002/sim.4385 Google Scholar | Crossref | Medline | ISI | |
| Belitz, C, Brezger, A, Kneib, T, Lang, S, Umlauf, N (2012) Bayesx ’ software for bayesian inference in structured additive regression models. Version 3.0. URL: http://www.bayesx.org Google Scholar | |
| Beyersmann, J, Schumacher, M, Allignol, A (2012) Competing risks and multistate models with R. Use R! New York: Springer. Google Scholar | Crossref | |
| Blanche, P, Dartigues, JF, Jacqmin-Gadda, H (2013) Estimating and comparing time-dependent areas under receiver operating characteristic curves for censored event times with competing risks. Statistics in Medicine, 32, 5381–97. Google Scholar | Crossref | Medline | ISI | |
| Brezger, A, Lang, S (2006) Generalized structured additive regression based on bayesian P-splines. Computational Statistics & Data Analysis, 50, 967–91. doi: DOI 10.1016/j.csda.2004.10.011 Google Scholar | Crossref | ISI | |
| Cadarso-Suárez, C, Meira-Machado, L, Kneib, T, Gude, F (2010) Flexible hazard ratio curves for continuous predictors in multi-state models: An application to breast cancer data. Statistical Modelling, 10, 291–314. doi: 10.1177/1471082x0801000303 Google Scholar | SAGE Journals | ISI | |
| Chen, HC, Kodell, R, Cheng, K, Chen, J (2012) Assessment of performance of survival prediction models for cancer prognosis. BMC Medical Research Methodology, 12, 102. Google Scholar | |
| Davenport, A (2009) Peritonitis remains the major clinical complication of peritoneal dialysis: The London, UK, peritonitis audit 2002–2003. Peritoneal Dialysis International, 29, 297–302. Google Scholar | Abstract | ISI | |
| Davison, A, Hinkley, D (1997) Bootstrap Methods and Their Application. UK: Cambridge University Press. ISBN 9780521574716 Google Scholar | Crossref | |
| de Wreede, LC, Fiocco, M, Putter, H (2011) Mstate: An r package for the analysis of competing risks and multi-state models. Journal of Statistical Software, 38, 1–30. Google Scholar | Crossref | ISI | |
| Efron, B (1981) Censored data and the bootstrap. Journal of the American Statistical Association, 76, 312–19. doi: 10.1080/01621459.1981.10477650 Google Scholar | Crossref | ISI | |
| Eilers, PHC, Marx, BD (1996) Flexible smoothing with b-splines and penalties. Statistical Science, 11, 89–102. doi: DOI 10.1214/ss/1038425655 Google Scholar | Crossref | ISI | |
| Eilers, PHC, Marx, BD (2002) Generalized linear additive smooth structures. Journal of Computational and Graphical Statistics, 11, 758–83. Google Scholar | Crossref | ISI | |
| Etzioni, R, Pepe, M, Longton, G, Hu, C, Goodman, G (1999) Incorporating the time dimension in receiver operating characteristic curves: A case study of prostate cancer. Medical Decision Making, 19, 242–51. Google Scholar | SAGE Journals | ISI | |
| Fahrmeir, L, Kneib, T (2011) Bayesian Smoothing and Regression for Longitudinal, Spatial and Event History Data. Oxford, UK: Oxford University Press. ISBN 9780199533022. Google Scholar | Crossref | |
| Fine, JP, Gray, RJ (1999) A proportional hazards model for the subdistribution of a competing risk. Journal of the American Statistical Association, 94, 496–509. Google Scholar | Crossref | ISI | |
| Gooley, TA, Leisenring, W, Crowley, J, Storer, BE (1999) Estimation of failure probabilities in the presence of competing risks: New repre- sentations of old estimators. Statistics in Medicine, 18, 695–706. Google Scholar | Crossref | Medline | ISI | |
| Hastie, T, Tibshirani, R (1990a) Exploring the nature of covariate effects in the proportional hazards model. Biometrics, 46, 1005–16. Google Scholar | Crossref | Medline | ISI | |
| Hastie, T, Tibshirani, R (1990b) Generalized additive models (Vol. 43). Boca Raton, FL: CRC Press. Google Scholar | |
| Heagerty, PJ, Lumley, T, Pepe, MS (2000) Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics, 56, 337–44. Google Scholar | Crossref | Medline | ISI | |
| Heagerty, PJ, Zheng, Y (2005) Survival model predictive accuracy and ROC curves. Biometrics, 61, 92–105. Google Scholar | Crossref | Medline | ISI | |
| Hennerfeind, A, Brezger, A, Fahrmeir, L (2006) Geoadditive survival models. Journal of the American Statistical Association, 101, 1065–75. Google Scholar | Crossref | ISI | |
| Hougaard, P (1999) Multi-state models: A review. Lifetime Data Analysis, 5, 239–64 Google Scholar | Crossref | Medline | ISI | |
| Kneib, T, Hennerfeind, A (2008) Bayesian semiparametric multi-state models. Statistical Modelling, 8, 169–98. Google Scholar | SAGE Journals | ISI | |
| Kneib, T, Fahrmeir, L (2007) A mixed model approach for geoadditive hazard regression. Scandinavian Journal of Statistics, 34, 207–28. Google Scholar | Crossref | ISI | |
| Kotsanas, D, Polkinghorne, KR, Korman, TM (2007) Risk factors for peritoneal dialysis-related peritonitis: Can we reduce the incidence and improve patient selection? Nephrology, 12, 239–45. Google Scholar | Crossref | Medline | ISI | |
| Levey, AS, Andreoli, SP, DuBose, T, Provenzano, R, Collins, AJ (2007) Chronic kidney disease: Common, harmful, and treatable—World Kidney Day 2007. Journal of the American Society of Nephrology, 18, 374–78. Google Scholar | Crossref | Medline | ISI | |
| Martins, M, Rodrigues, A, Pedrosa, JM, Carvalho, MJ, Cabrita, A, Oliveira, R (2013) Update on the challenging role of biofilms in peritoneal dialysis. Biofouling, 29, 1015–27. Google Scholar | Crossref | Medline | ISI | |
| Martinussen, T, Scheike, T (2010) Dynamic regression models for survival data. USA: Springer. ISBN 9781441919045. Google Scholar | |
| Meira-Machado, L, de Una-Alvarez, J, rez, C, Andersen, PK (2009) Multi- state models for the analysis of time-to- event data. Statistical Methods in Medical Research, 18, 195–222. Google Scholar | SAGE Journals | ISI | |
| Meira-Machado, L, Cadarso-Suárez, C, Gude, F, Arajo, A (2013) Smooth HR: An R package for pointwise nonparametric estimation of hazard ratio curves of continuous predictors. Computational and Mathematical Methods in Medicine, 28, 1–11. Google Scholar | Crossref | ISI | |
| Pepe, (2004) The statistical evaluation of medical tests for classification and prediction. USA: Oxford University Press. ISBN 9780198565826. Google Scholar | |
| Pepe, M, Zheng, Y, Jin, Y, Huang, Y, Parikh, C, Levy, W (2008) Evaluating the ROC performance of markers for future events. Lifetime Data Analysis, 14, 86–114. Google Scholar | Crossref | Medline | ISI | |
| Perl, J, Wald, R, Bargman, JM, Na, Y, Jassal, SV, Jain, AK, Moist, L, Nessim, SJ (2012) Changes in patient and technique survival over time among incident peritoneal dialysis patients in Canada. Clinical Journal of the American Society of Nephrology, 7, 1145–54. Google Scholar | Crossref | Medline | ISI | |
| Prentice, RL, Kalbfleisch, JD, Peterson, AV, Flournoy, N, Farewell, VT, Breslow, NE (1978) Analysis of failure times in presence of competing risks. Biometrics, 34, 541–54. Google Scholar | Crossref | Medline | ISI | |
| Putter, H, Fiocco, M, Geskus, RB (2007) Tutorial in biostatistics: Competing risks and multi-state models. Statistics in Medicine, 26, 2389–430. Google Scholar | Crossref | Medline | ISI | |
| R Development CoreTeam (2013) R: A language and environment for statistical computing. Vienna, Austria: The R Foundation for Statistical Computing. ISBN: 3-900051-07-0. URL: http://www.R-project.org/ Google Scholar | |
| Rocha, A, Rodrigues, A, Teixeira, L, Carvalho, MJ, Mendonça, D, Cabrita, A (2012) Temporal trends in peritonitis rates, microbiology and outcomes: The major clinical complication of peritoneal dialysis. Blood Purification, 33, 284–91. Google Scholar | Crossref | Medline | ISI | |
| Rodriguez-Girondo, M, Kneib, T, Cadarso-Suárez, C, Abu-Assi, E (2013) Model building in nonproportional hazard regression. Statistics in Medicine, 32, 5301–14. Google Scholar | Crossref | Medline | ISI | |
| Saha, P, Heagerty, PJ (2010) Time-dependent predictive accuracy in the presence of competing risks. Biometrics, 66, 999–1011. Google Scholar | Crossref | Medline | ISI | |
| Schemper, M (2003) Predictive accuracy and explained variation. Statistics in Medicine, 22, 2299–2308. Google Scholar | Crossref | Medline | ISI | |
| Swets, J, Pickett, R (1982) Evaluation of diagnostic systems: methods from signal detection theory. USA: Elsevier Science & Technology Books. ISBN 9780126790801. Google Scholar | |
| Teixeira, L, Rodrigues, A, Carvalho, MJ, Cabrita, A, Mendonça, D (2013) Modelling competing risks in nephrology research: An example in peritoneal dialysis. BMC Nephrology, 14, 110. Google Scholar | Medline | ISI | |
| Yang, X, Yi, C, Liu, X, Guo, Q, Yang, R, Cao, P, Lin, J, Mao, H, Yu, X (2013) Clinical outcome and risk factors for mortality in chinese patients with diabetes on peritoneal dialysis: A 5-year clinical cohort study. Diabetes Research and Clinical Practice, 100, 354–61. Google Scholar | Crossref | Medline | ISI | |
| Zheng, YY, Cai, TX, Jin, YY, Feng, ZD (2012) Evaluating prognostic accuracy of biomarkers under competing risk. Biometrics, 68, 388–96. Google Scholar | Crossref | Medline | ISI | |
| Zhou, X, Obuchowski, N, McClish, D (2011) Statistical Methods in Diagnostic Medicine. USA: Wiley. ISBN 9780470906507. Google Scholar | Crossref |
