outcome_binary function

Simulate from binary model given covariates