outcome_count function

Simulate from count model given covariates