Regression Analysis for Multistate Models Based on a Pseudo‐value Approach, with Applications to Bone Marrow Transplantation Studies

Typically, regression analysis for multistate models has been based on regression models for the transition intensities. These models lead to highly nonlinear and very complex models for the effects of covariates on state occupation probabilities. We present a technique that models the state occupation or transition probabilities in a multistate model directly. The method is based on the pseudo-values from a jackknife statistic constructed from non-parametric estimators for the probability in question. These pseudo-values are used as outcome variables in a generalized estimating equation to obtain estimates of model parameters. We examine this approach and its properties in detail for two special multistate model probabilities, the cumulative incidence function in competing risks and the current leukaemia-free survival used in bone marrow transplants. The latter is the probability a patient is alive and in either a first or second post-transplant remission. The techniques are illustrated on a dataset of leukaemia patients given a marrow transplant. We also discuss extensions of the model that are of current research interest. Copyright 2007 Board of the Foundation of the Scandinavian Journal of Statistics..

[1]  J. Fine,et al.  Regression modeling of competing crude failure probabilities. , 2001, Biostatistics.

[2]  John P. Klein,et al.  Additive hazards Markov regression models illustrated with bone marrow transplant data , 2005 .

[3]  Robert Gray,et al.  A Proportional Hazards Model for the Subdistribution of a Competing Risk , 1999 .

[4]  David R. Cox,et al.  Regression models and life tables (with discussion , 1972 .

[5]  O. Aalen A linear regression model for the analysis of life times. , 1989, Statistics in medicine.

[6]  D. Cox Regression Models and Life-Tables , 1972 .

[7]  M. Pepe,et al.  Summarizing data on survival, relapse, and chronic graft‐versus‐host disease after bone marrow transplantation: motivation for and description of new methods , 1993, British journal of haematology.

[8]  M. Crowder Classical Competing Risks , 2001 .

[9]  Zhiliang Ying,et al.  Semiparametric analysis of the additive risk model , 1994 .

[10]  Margaret S. Pepe,et al.  Inference for Events with Dependent Risks in Multiple Endpoint Studies , 1991 .

[11]  Jerald F. Lawless,et al.  Some Simple Robust Methods for the Analysis of Recurrent Events , 1995 .

[12]  J. Klein,et al.  Generalised linear models for correlated pseudo‐observations, with applications to multi‐state models , 2003 .

[13]  J P Klein,et al.  Estimation of current leukaemia-free survival following donor lymphocyte infusion therapy for patients with leukaemia who relapse after allografting: application of a multistate model. , 2000, Statistics in medicine.

[14]  Zhiliang Ying,et al.  Semiparametric Analysis of General Additive-Multiplicative Hazard Models for Counting Processes , 1995 .

[15]  Somnath Datta,et al.  Validity of the Aalen–Johansen estimators of stage occupation probabilities and Nelson–Aalen estimators of integrated transition hazards for non-Markov models , 2001 .

[16]  Zhiliang Ying,et al.  Semiparametric regression for the mean and rate functions of recurrent events , 2000 .

[17]  R. Gray A Class of $K$-Sample Tests for Comparing the Cumulative Incidence of a Competing Risk , 1988 .

[18]  S. Zeger,et al.  Longitudinal data analysis using generalized linear models , 1986 .

[19]  M S Pepe,et al.  Kaplan-Meier, marginal or conditional probability curves in summarizing competing risks failure time data? , 1993, Statistics in medicine.

[20]  J. Klein,et al.  Estimating leukemia-free survival after allografting for chronic myeloid leukemia: a new method that takes into account patients who relapse and are restored to complete remission , 2000 .

[21]  N. Keiding,et al.  Multi-state models for event history analysis , 2002, Statistical methods in medical research.

[22]  Donald Hedeker,et al.  Longitudinal Data Analysis , 2006 .

[23]  Lee-Jen Wei,et al.  Prediction of cumulative incidence function under the proportional hazards model. , 1998, Biometrics.

[24]  John P Klein,et al.  Regression Modeling of Competing Risks Data Based on Pseudovalues of the Cumulative Incidence Function , 2005, Biometrics.

[25]  John P Klein,et al.  Modelling competing risks in cancer studies , 2006, Statistics in medicine.

[26]  Rupert G. Miller The jackknife-a review , 1974 .

[27]  J. Kalbfleisch,et al.  The Statistical Analysis of Failure Time Data: Kalbfleisch/The Statistical , 2002 .

[28]  Dario Gasbarra,et al.  Testing equality of cause-specific hazard rates corresponding to m competing risks among K groups , 2002 .

[29]  J P Klein,et al.  Summary curves for patients transplanted for chronic myeloid leukaemia salvaged by a donor lymphocyte infusion: the current leukaemia‐free survival curve , 2000, British journal of haematology.

[30]  John P. Klein,et al.  Regression Analysis of Restricted Mean Survival Time Based on Pseudo-Observations , 2004, Lifetime data analysis.

[31]  Niels Keiding,et al.  Statistical Models Based on Counting Processes , 1993 .

[32]  Ren Johansen An Empirical Transition Matrix for Non-homogeneous Markov Chains Based on Censored Observations , 1978 .

[33]  Susanne Rosthøj,et al.  Competing risks as a multi-state model , 2002, Statistical methods in medical research.

[34]  J. Klein,et al.  Statistical Models Based On Counting Process , 1994 .