Identifying counterfactual probabilities using bivariate distributions and uplift modeling