Adding covariates to Twin Ace model