fit_models_RA function

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