run_model function

Perform Bayesian pairwise or network meta-analysis