Competing risks regression for clustered data with covariate-dependent censoring

Competing risks data in clinical trials or observational studies often suffer from cluster effects such as center effects and matched pairs design. The proportional subdistribution hazards (PSH) model is one of the most widely used methods for competing risks data analyses. However, the current literature on the PSH model for clustered competing risks data is limited to covariate-independent censoring and the unstratified model. In practice, competing risks data often face covariate-dependent censoring and have a non PSH structure. Thus, we propose a marginal stratified PSH model with a covariate-adjusted censoring weight for clustered competing risks data. We use a marginal stratified proportional hazards model to estimate the survival probability of censoring by taking clusters and non proportional hazards structures into account. Our simulation results show that, in the presence of covariate-dependent censoring, the parameter estimates of the proposed method are unbiased with approximate 95% coverage rates. We apply the proposed method to stem cell transplant data of leukemia patients to evaluate the clinical implications of donor-recipient HLA matching on chronic graft-versus-host disease.