Meta-analysis of individual patient data with semi-competing risks under the Weibull joint frailty–copula model

In meta-analysis of individual patient data with semi-competing risks, the joint frailty–copula model has been proposed, where frailty terms account for the between-study heterogeneity and copulas account for dependence between terminal and nonterminal event times. In the previous works, the baseline hazard functions in the joint frailty–copula model are estimated by the nonparametric model or the penalized spline model, which requires complex maximization schemes and resampling-based interval estimation. In this article, we propose the Weibull distribution for the baseline hazard functions under the joint frailty–copula model. We show that the Weibull model constitutes a conjugate model for the gamma frailty, leading to explicit expressions for the moments, survival functions, hazard functions, quantiles, and mean residual lifetimes. These results facilitate the parameter interpretation of prognostic inference. We propose a maximum likelihood estimation method and make our computer programs available in the R package, joint.Cox. We also show that the delta method is feasible to calculate interval estimates, which is a useful alternative to the resampling-based method. We conduct simulation studies to examine the accuracy of the proposed methods. Finally, we use the data on ovarian cancer patients to illustrate the proposed method.