fit_model function

Fit probit model to choice data