Bivariate survival models with Clayton aging functions