A Composite Likelihood Method for the Analysis of Multivariate Survival Data: An Application to a PBDE Study

When dealing with multivariate survival data, featuring the association structure is a key difference from the univariate survival analysis. In this paper, we explore to use the composite likelihood framework to handle multivariate survival data, where only the lower dimensional survival distributions need to be specified. The development allows us to use available modeling schemes for bivariate survival data to characterize association structures of correlated survival times. The inference procedure is based on the pseudolikelihood which is the product of the lower dimensional bivariate distributions. The proposed estimation procedure is assessed through simulation studies. As a genuine application, we apply the composite likelihood inference procedure to analyze the data from the polybrominated diphenyl ethers (PBDEs) study, where four types of PBDE congeners are available. The associations among the four PBDE congeners, and the relationships between the covariates and the PBDE congeners are of interest. The result shows that there is strong association among the concentrations of the four PBDE congeners, and statistically significant predictors on the concentrations of the four PBDE congeners are identified.