diff --git a/src/bayesvalidrox/bayes_inference/bayes_inference.py b/src/bayesvalidrox/bayes_inference/bayes_inference.py
index 1898a8ae619597d92bc355ac4249f57019f0aed7..00971b92af1bf574d1c61a470cfdde7506cffecb 100644
--- a/src/bayesvalidrox/bayes_inference/bayes_inference.py
+++ b/src/bayesvalidrox/bayes_inference/bayes_inference.py
@@ -446,6 +446,9 @@ class BayesInference:
             # BME = BME + BMECorrFactor
             if self.emulator:
                 self.log_BME += self.log_BME_corr_factor
+                
+        if self.bmc: 
+            return self
 
         # ---------------- Parameter Bayesian inference ----------------
         if self.inference_method.lower() == 'mcmc':
diff --git a/src/bayesvalidrox/bayes_inference/bayes_model_comparison.py b/src/bayesvalidrox/bayes_inference/bayes_model_comparison.py
index 828613556e90ec0c529b91f2592eec148c98136b..fd01689d70e59031434cea0696e82079d03b83d2 100644
--- a/src/bayesvalidrox/bayes_inference/bayes_model_comparison.py
+++ b/src/bayesvalidrox/bayes_inference/bayes_model_comparison.py
@@ -169,6 +169,7 @@ class BayesModelComparison:
             # Pass justifiability data as perturbed data
             BayesOpts.perturbed_data = justData
             BayesOpts.just_analysis = justifiability
+            BayesOpts.bmc = True
 
             bayes_dict[model] = BayesOpts.create_inference()
             print("-"*20)