Problems expanding univariate liability threshold script to bivariate and multivariate script