Frailty Survival Model Analysis of the National Deceased Donor Kidney Transplant Dataset Using Poisson Variance Structures

In a recent study of transplant outcomes, donor age, cerebrovascular accident as the cause of death (CVA), renal insufficiency (serum creatinine >1.5 mg/dL), and history of hypertension have been identified as donor factors associated with elevated risk of kidney transplant failure. It is of great interest to know whether there remain other unmeasured donor factors associated with elevated risk of graft failure. In this article we study a sample of 6,024 deceased donor kidney transplants performed in 194 centers from 1995 to 2000. In addition to variation among transplant recipients, there are two other random effects: unmeasured donor and unrecorded center factors (data not available at the physician level). These two random effects are crossed, because the two kidneys from the same donor can be transplanted in different centers. Multivariate frailty models are applied to analyze the data. The likelihood functions of both parametric (e.g., with piecewise constant baseline hazard) and semiparametric multivariate frailty models are shown to be proportional to the likelihood functions of a class of mixed Poisson regression models. The penalized quasi-likelihood method is used as the numerical procedure for these mixed Poisson regression models. Thus we are able to estimate and model crossed random-effects structures for survival analysis. Although about 30% of recipient graft survival rate variation due to donor factors is explained by the measured donor characteristics, the remaining variation among donors in graft survival rate is still statistically significant, suggesting that there may be other unmeasured donor factors associated with a reduced graft survival rate. We also find significant variation of graft failure rates among transplant centers due to unrecorded center factors. Therefore, this study suggests that practice patterns at transplant centers and identification of other donor factors may merit further investigation.

[1]  Norman Breslow,et al.  Discussion of Professor Cox''s paper , 1974 .

[2]  Nan M. Laird,et al.  Covariance Analysis of Censored Survival Data Using Log-Linear Analysis Techniques , 1981 .

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

[4]  N. Breslow,et al.  Approximate inference in generalized linear mixed models , 1993 .

[5]  N. Sastry A nested frailty model for survival data, with an application to the study of child survival in northeast Brazil. , 1997, Journal of the American Statistical Association.

[6]  R Ecochard,et al.  Multi-level modelling of conception in artificial insemination by donor. , 1998, Statistics in medicine.

[7]  J. F. Lawless,et al.  Analysis of interval‐grouped recurrent‐event data using piecewise constant rate functions , 1998 .

[8]  Erik T. Parner,et al.  Asymptotic theory for the correlated gamma-frailty model , 1998 .

[9]  T. Louis,et al.  Triple‐goal estimates in two‐stage hierarchical models , 1998 .

[10]  R. Henderson,et al.  Effect of frailty on marginal regression estimates in survival analysis , 1999 .

[11]  R. Wolfe,et al.  Comparison of mortality in all patients on dialysis, patients on dialysis awaiting transplantation, and recipients of a first cadaveric transplant. , 1999, The New England journal of medicine.

[12]  F. Vaida,et al.  Proportional hazards model with random effects. , 2000, Statistics in medicine.

[13]  J. Palmgren,et al.  Estimation of Multivariate Frailty Models Using Penalized Partial Likelihood , 2000, Biometrics.

[14]  R. Wolfe,et al.  Survival in recipients of marginal cadaveric donor kidneys compared with other recipients and wait-listed transplant candidates. , 2001, Journal of the American Society of Nephrology : JASN.

[15]  Robert M. Merion,et al.  Donor characteristics associated with reduced graft survival: an approach to expanding the pool of kidney donors1 , 2002, Transplantation.

[16]  Tze Leung Lai,et al.  A hybrid estimator in nonlinear and generalised linear mixed effects models , 2003 .

[17]  R. Merion Preface: SRTR Report on the State of Transplantation , 2003 .

[18]  Daniel Krewski,et al.  Random effects Cox models: A Poisson modelling approach , 2003 .