SciELO - Scientific Electronic Library Online

 
vol.43 issue1Generalized Poisson Hidden Markov Model for Overdispersed or Underdispersed Count DataNested and Repeated Cross Validation for Classification Model With High-dimensional Data author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Revista Colombiana de Estadística

Print version ISSN 0120-1751

Rev.Colomb.Estad. vol.43 no.1 Bogotá Jan./June 2020  Epub June 05, 2020

https://doi.org/10.15446/rce.v43n1.77851 

ARTÍCULOS ORIGINALES DE INVESTIGACIÓN

A Birnbaum-Saunders Model for Joint Survival and Longitudinal Analysis of Congestive Heart Failure Data

Un modelo Birnbaum-Saunders para el análisis conjunto de datos de supervivencia y longitudinales de insuficiencia cardíaca congestive

Diana e. Franco-Soto1 

Antonio e. Pedroso-de-Lima2 

Julio M. Singer3 

1 Departamento de Estadística, Facultad de Ciencias, Universidad Nacional de Colombia, Bogotá D.C., Colombia. PhD. E-mail: dcfrancos@unal.edu.co

2 Departamento de Estatística, Instituto de Matemática e Estatística, Universidade de São Paulo, São Paulo, Brazil. PhD. E-mail: acarlos@ime.usp.br

3 Departamento de Estatística, Instituto de Matemática e Estatística, Universidade de São Paulo, São Paulo, Brazil. PhD. E-mail: jmsinger@ime.usp.br


Abstract

We consider a parametric joint modelling of longitudinal measurements and survival times, motivated by a study conducted at the Heart Institute (Incor), São Paulo, Brazil, with the objective of evaluating the impact of B-type Natriuretic Peptide (BNP) collected at different instants on the survival of patients with Congestive Heart Failure (CHF). We employ a linear mixed model for the longitudinal response and a Birnbaum-Saunders model for the survival times, allowing the inclusion of subjects without longitudinal observations. We derive maximum likelihood estimators of the joint model parameters and conduct a simulation study to compare the true survival probabilities with dynamic predictions obtained from the fit of the proposed joint model and to evaluate the performance of the method for estimating the model parameters. The proposed joint model is applied to the cohort of 1609 patients with CHF, of which 1080 have no BNP measurements. The parameter estimates and their standard errors obtained via: i) the traditional approach, where only individuals with at least one measurement of the longitudinal response are included and ii) the proposed approach, which includes survival information from all individuals, are compared with those obtained via marginal (longitudinal and survival) models.

Key words: Linear mixed model; Parametric survival model; Repeated measures

Resumen

Consideramos una modelación conjunta paramétrica de mediciones longitudinales y tiempos de supervivencia, motivados por un estudio realizado en el Instituto do Coração (Incor), São Paulo, Brasil, con el objetivo de evaluar el impacto del Péptido Natriurético tipo B (BNP) recolectado en diferentes instantes, sobre la supervivencia de pacientes con Insuficiencia Cardíaca Congestiva (ICC). Empleamos un modelo lineal de efectos mixtos para la respuesta longitudinal y un modelo Birnbaum-Saunders para los tiempos de supervivencia, permitiendo la inclusión de sujetos sin observaciones longitudinales. Obtenemos los estimadores de máxima verosimilitud de los parámetros del modelo conjunto y realizamos un estudio de simulación para comparar las probabilidades de supervivencia verdaderas con las predicciones dinámicas obtenidas al ajustar el modelo conjunto propuesto y para evaluar el desempeño del método para estimar los parámetros del modelo. El modelo conjunto propuesto se aplica a la cohorte de 1609 pacientes con ICC, de los cuales 1080 no tienen mediciones de BNP. Las estimaciones de los parámetros y sus errores estándar obtenidos por medio de: i) el enfoque tradicional, donde únicamente se incluyen individuos con al menos una medición de la respuesta longitudinal y ii) el enfoque propuesto, que incluye la información de supervivencia de todos los individuos; se comparan con los obtenidos por medio de los modelos marginales (longitudinal y de supervivencia).

Palabras clave: Medidas repetidas; Modelo de supervivencia paramétrico; Modelo lineal de efectos mixtos

Full text available only in PDF format.

References

Albert, P. S. & Shih, J. H. (2010), 'On estimating the relationship between longitudinal measurements and time-to-event data using a simple two-stage procedure', Biometrics 66, 983-991. [ Links ]

Balakrishnan, N., Leiva, V., Sanhueza, A. & Vilca, F. (2009), 'Estimation in the Birnbaum-Saunders distribution based on scale-mixture of normals and the EM algorithm', Statistics and Operations Research Transactions 33, 171-192. [ Links ]

Barros, M., Paula, G. A. & Leiva, V. (2008), 'A new class of survival regression models with heavy-tailed errors: robustness and diagnostics', Lifetime Data Analysis 14, 316-332. [ Links ]

Crowther, M. J., Abrams, K. R. & Lambert, P. C. (2012), 'Flexible parametric joint modelling of longitudinal and survival data', Statistics in Medicine 31, 4456-4471. [ Links ]

DeGruttola, V. & Tu, X. M. (1994), 'Modelling progression of CD-4 lymphocyte count and its relationship to survival time', Biometrics 50, 1003-1014. [ Links ]

Diggle, P., Farewell, D. & Henderson, R. (2007), 'Analysis of longitudinal data with drop-out: objectives, assumptions and a proposal', Journal of the Royal Statistical Society: Series C (Applied Statistics) 56, 499-550. [ Links ]

Diggle, P. J., Sousa, I. & Chetwynd, A. G. (2008), 'Joint modelling of repeated measurements and time-to-event outcomes: the fourth armitage lecture', Statistics in Medicine 27, 2981-2998. [ Links ]

Franco-Soto, D. C. (2014), Modelos Birnbaum-Saunders para sobrevivência com dados longitudinais, PhD dissertation, Departamento de Estatística, Instituto de Matemática e Estatística, Universidade de São Paulo, Brasil. *https://www.ime.usp.br/Links ]

Galea, M., Leiva, V. & Paula, G. (2004), 'Influence diagnostics in log-Birnbaum-Saunders regression models', Journal Applied Statistics 31, 1049-1064. [ Links ]

Gay, D. M. (1990), 'Usage summary for selected optimization routines', Computing science technical report 153, 1-21. [ Links ]

Gould, A. L., Boye, M. E., Crowther, M. J., Ibrahim, J. G., Quartey, G., Micallef, S. & Bois, F. Y. (2015), 'Joint modeling of survival and longitudinal non-survival data: current methods and issues. Report of the DIA Bayesian joint modeling working group', Statistics in Medicine 34(14), 2181-2195. [ Links ]

Greene, W. F. & Cai, J. (2004), 'Measurement error in covariates in the marginal hazards model for multivariate failure time data', Biometrics 60, 987-996. [ Links ]

Henderson, R., Diggle, P. & Dobson, A. (2000), 'Joint modelling of longitudinal measurements and event time data', Biostatistics 1, 465-480. [ Links ]

Hogan, J. W. & Laird, N. M. (1997a), 'Mixture models for the joint distributions of repeated measures and event times', Statistics in Medicine 16, 239-257. [ Links ]

Hogan, J. W. & Laird, N. M. (1997b), 'Model-based approaches to analysing incomplete longitudinal and failure time data', Statistics in Medicine 16, 259-272. [ Links ]

Hsieh, F., Tseng, Y. K. & Wang, J. L. (2006), 'Joint modelling of survival and longitudinal data: likelihood approach revisited', Biometrics 62, 1037-1043. [ Links ]

Hu, P., Tsiatis, A. A. & Davidian, M. (1998), 'Estimating the parameters in the Cox model when covariate variables are measured with error', Biometrics 54, 1407-1419. [ Links ]

Kalbfleisch, J. D. & Prentice, R. L. (2002), The Statistical Analysis of Failure Time Data, 2 nd edn, John Wiley & Sons, New Jersey. [ Links ]

Klein, J. P. & Moeschberger, M. L. (2003), Survival Analysis: Techniques for Censored and Truncated Data, 2 nd edn, Springer-Verlag, New York. [ Links ]

Leiva, V., Athayde, E., Azevedo, C. & Marchant, C. (2011), 'Modelling wind energy flux by a Birnbaum-Saunders distribution with unknown shift parameter', Journal of Applied Statistics 54, 2819-2838. [ Links ]

Leiva, V., Barros, M. K., Paula, G. A. & Galea, M. (2007), 'Influence diagnostics in log-Birnbaum-Saunders regression models with censored data', Computational Statistics and Data Analysis 51, 5694-5707. [ Links ]

Pawitan, Y. & Self, S. (1993), 'Modelling disease marker processes in AIDS', Journal of the American Statistical Association 83, 719-726. [ Links ]

R Development Core Team (2013), R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria. *https://www.R-project.orgLinks ]

Rizopoulos, D. (2010), 'JM: An R package for the joint modelling of longitudinal and time-to-event data', Journal of Statistical Software 35, 1-33. [ Links ]

Rizopoulos, D. (2011), 'Dynamic predictions and prospective accuracy in joint models for longitudinal and time-to-event data', Biometrics 67, 819-829. [ Links ]

Rizopoulos, D. (2012a), JM: Shared parameter models for the joint modelling of longitudinal and time-to-event data, R package version 1.2-0. *https://cran.r-project.org/web/packages/JM/index.htmlLinks ]

Rizopoulos, D. (2016), Joint Models for Longitudinal and Time-to-Event Data With Applications in R, 1 st edn, Chapman-Hall/CRC Biostatistics Series, Boca Raton. [ Links ]

Rizopoulos, D., Verbeke, G. & Lesaffre, E. (2009), 'Fully exponential Laplace approximations for the joint modelling of survival and longitudinal data', Journal of the Royal Statistical Society: Series B (Methodological) 71, 637-654. [ Links ]

Rizopoulos, D., Verbeke, G. & Molenberghs, G. (2010), 'Multiple-imputation-based residuals and diagnostic plots for joint models of longitudinal and survival outcomes', Biometrics 66, 20-29. [ Links ]

Schluchter, M. D. (1992), 'Methods for the analysis of informatively censored longitudinal data', Statistics in Medicine 11, 1861-1870. [ Links ]

Tsiatis, A. A. & Davidian, M. (2004), 'Joint modelling of longitudinal and time-to-event data: an overview', Statistica Sinica 14, 809-834. [ Links ]

Tsiatis, A. A., DeGruttola, V. & Wulfsohn, M. S. (1995), 'Modelling the relationship of survival to longitudinal data measured with error: applications to survival and CD4 counts in patients with AIDS', Journal of the American Statistical Association 90, 27-37. [ Links ]

Wu, L., Liu, W., Yi, G. Y. & Huang, Y. (2012), 'Analysis of longitudinal and survival data: joint modeling, inference methods, and issues', Journal of Probability and Statistics 2012, 1-17. [ Links ]

Wulfsohn, M. S. & Tsiatis, A. A. (1997), 'A joint model for survival and longitudinal data measured with error', Biometrics 53, 330-339. [ Links ]

Ye, W., Lin, X. & Taylor, J. M. G. (2008), 'Semiparametric modeling of longitudinal measurements and time-to-event data: a two-stage regression calibration approach', Biometrics 64, 1238-1246. [ Links ]

Yu, M., Law, N. J., Taylor, J. M. G. & Sandler, H. M. (2004), 'Joint longitudinal-survival-cure models and their application to prostate cancer', Statistica Sinica 14, 835-862. [ Links ]

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License