fit_models_RA_5bm function

Model fitting for reference analysis using 5 benchmarks: Posterior inference for benchmark and actual heterogeneity priors