update_model function

Fit causal model using 'stan'