Estimating and testing for center effects in competing risks.

The problems of fitting Gaussian frailties proportional hazards models for the subdistribution of a competing risk and of testing for center effects are considered. In the analysis of competing risks data, Fine and Gray proposed a proportional hazards model for the subdistribution to directly assess the effects of covariates on the marginal failure probabilities of a given failure cause. Katsahianbiet al. extended their model to clustered time to event data, by including random center effects or frailties in the subdistribution hazard. We first introduce an alternate estimation procedure to the one proposed by Katsahian et al. This alternate estimation method is based on the penalized partial likelihood approach often used in fitting Gaussian frailty proportional hazards models in the standard survival analysis context, and has the advantage of using standard survival analysis software. Second, four hypothesis tests for the presence of center effects are given and compared via Monte-Carlo simulations. Statistical and numerical considerations lead us to formulate pragmatic guidelines as to which of the four tests is preferable. We also illustrate the proposed methodology with registry data from bone marrow transplantation for acute myeloid leukemia (AML).