compare_methods_with_gibbs function

Compare estimation methods with and without Gibbs sampling