Jointly modelling the relationship between survival and pulmonary function in cystic fibrosis patients

Modelling the relationship between pulmonary function and survival in cystic fibrosis (CF) is complicated by the fact that measures of pulmonary function commonly used such as the forced expiratory volume in one second (FEV(1)) are measured with error, and patients with the poorest lung function are increasingly censored by death, that is, data are available only for the patients who have survived to the current age. We assume a linear random effects model for FEV1 per cent predicted, where the random intercept and slope of FEV(1) per cent predicted, along with a specified transformation of the age at death follow a trivariate normal distribution. We illustrate how this model can be used to describe the relationship between age at death and parameters of the individual patient's regression of FEV(1) per cent predicted versus age, such as the slope and the intercept or true value of FEV(1) per cent predicted at a given age. We also illustrate how the model provides empirical Bayes estimates of these individual parameters. In particular, we explore how the predicted value of the age at death might be used as a prognostic or severity index. The model and methods are illustrated on a cohort of 188 cystic fibrosis patients with a common genotype (homozygous for the DeltaF508 mutation), born on or after 1965 and followed at the CF Center at the Rainbow Babies and Children's Hospital, Cleveland, OH, U.S.A.

[1]  O Siddiqui,et al.  A comparison of the random-effects pattern mixture model with last-observation-carried-forward (LOCF) analysis in longitudinal clinical trials with dropouts. , 1998, Journal of biopharmaceutical statistics.

[2]  V. De Gruttola,et al.  Modelling progression of CD4-lymphocyte count and its relationship to survival time. , 1994, Biometrics.

[3]  K. Bailey,et al.  Estimation and comparison of changes in the presence of informative right censoring: conditional linear model. , 1989, Biometrics.

[4]  N M Laird,et al.  Mixture models for the joint distribution of repeated measures and event times. , 1997, Statistics in medicine.

[5]  K. Bailey,et al.  Analysing changes in the presence of informative right censoring caused by death and withdrawal. , 1988, Statistics in medicine.

[6]  N. L. Johnson,et al.  Continuous Univariate Distributions. , 1995 .

[7]  T Park,et al.  Simple pattern-mixture models for longitudinal data with missing observations: analysis of urinary incontinence data. , 1999, Statistics in medicine.

[8]  John J. Chen,et al.  STATISTICAL ANALYSIS OF THE LONG-TERM EFFECTS OF RECOMBINANT HUMAN DEOXYRIBONUCLEASE ON PULMONARY FUNCTION IN CYSTIC FIBROSIS PATIENTS , 2000, Journal of biopharmaceutical statistics.

[9]  S. Thompson,et al.  A joint analysis of quality of life and survival using a random effect selection model. , 2000, Statistics in medicine.

[10]  Donald Hedeker,et al.  Application of random-efiects pattern-mixture models for miss-ing data in longitudinal studies , 1997 .

[11]  S J Pocock,et al.  Estimation and comparison of rates of change in longitudinal studies with informative drop-outs. , 1999, Statistics in medicine.

[12]  Yudi Pawitan,et al.  Modeling Disease Marker Processes in AIDS , 1993 .

[13]  Geert Molenberghs,et al.  Monotone missing data and pattern‐mixture models , 1998 .

[14]  M D Schluchter,et al.  Methods for the analysis of informatively censored longitudinal data. , 1992, Statistics in medicine.

[15]  M. Kenward,et al.  Informative Drop‐Out in Longitudinal Data Analysis , 1994 .

[16]  D. Rubin.  INFERENCE AND MISSING DATA , 1975 .

[17]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[18]  N M Laird,et al.  Increasing efficiency from censored survival data by using random effects to model longitudinal covariates , 1998, Statistical methods in medical research.

[19]  B. Efron.  Bootstrap Methods: Another Look at the Jackknife , 1979 .

[20]  M. Wulfsohn,et al.  A joint model for survival and longitudinal data measured with error. , 1997, Biometrics.

[21]  P. J. Byard,et al.  Effect of high-dose ibuprofen in patients with cystic fibrosis. , 1995, The New England journal of medicine.

[22]  R F Woolson,et al.  Application of empirical Bayes inference to estimation of rate of change in the presence of informative right censoring. , 1992, Statistics in medicine.

[23]  D. Thomas,et al.  Simultaneously modelling censored survival data and repeatedly measured covariates: a Gibbs sampling approach. , 1996, Statistics in medicine.

[24]  M. Lebowitz,et al.  Changes in the normal maximal expiratory flow-volume curve with growth and aging. , 1983, The American review of respiratory disease.

[25]  R Henderson,et al.  Joint modelling of longitudinal measurements and event time data. , 2000, Biostatistics.

[26]  M D Schluchter,et al.  Analysis of change in the presence of informative censoring: application to a longitudinal clinical trial of progressive renal disease , 2001, Statistics in medicine.

[27]  Raymond J. Carroll,et al.  Estimation and comparison of changes in the presence of informative right censoring by modeling the censoring process , 1988 .

[28]  R F Woolson,et al.  Slope estimation in the presence of informative right censoring: modeling the number of observations as a geometric random variable. , 1994, Biometrics.

[29]  N M Laird,et al.  Model-based approaches to analysing incomplete longitudinal and failure time data. , 1997, Statistics in medicine.

[30]  M. Corey,et al.  Longitudinal analysis of pulmonary function decline in patients with cystic fibrosis. , 1997, The Journal of pediatrics.