Abstract
In clinical and epidemiological studies, when the time-to-event(s) and the longitudinal outcomes are associated, modelling them separately may give biased estimates. A joint modelling approach is required to obtain unbiased results and to evaluate their association. In the joint model, a subject may be exposed to more than one type of failure event (competing risks). Considering the competing event as an independent censoring of the time-to-event process may underestimate the true survival probability and give biased results. Within the joint model, longitudinal outcomes may have nonlinear (irregular) trajectories over time and exhibit skewness with heavy tails. Accordingly, fully parametric mixed-effect models may not be flexible enough to model this type of complex longitudinal data. In addition, assuming a Gaussian distribution for model errors may be too restrictive to adequately represent within-individual variations and may lack robustness against deviation from distributional assumptions. To simultaneously overcome these issues, in this paper, we presented semiparametric joint models for competing risks failure time and skewed-longitudinal data by using a smoothing spline approach and a multivariate skew-t distribution. We also considered different parameterization approaches in the formulation of joint models and used a Bayesian approach to make the statistical inference. We illustrated the proposed methods by analyzing real data on a chronic kidney disease. To evaluate the performance of the methods, we also carried out simulation studies. The results of both the application and simulation studies revealed that the joint modelling approach proposed in this study performed well when the semiparametric, random-effects parameterization, and skew-t distribution specifications were taken into account.