diff --git a/cara/montecarlo.py b/cara/montecarlo.py
index b993ea02..2d4aaef3 100644
--- a/cara/montecarlo.py
+++ b/cara/montecarlo.py
@@ -833,8 +833,13 @@ exposure_models = [MCExposureModel(
 ) for qid in (100, 60)]
 
 
+rs = [model.expected_new_cases() for model in large_population_baselines]
 
-compare_infection_probabilities_vs_viral_loads(*exposure_models)
+print(f"R0 - original variant:\t{np.mean(rs[0])}")
+print(f"R0 - english variant:\t{np.mean(rs[1])}")
+print(f"Ratio between R0's:\t\t{np.mean(rs[1]) / np.mean(rs[0])}")
+
+# compare_infection_probabilities_vs_viral_loads(*exposure_models)
 
 
 # present_model(exposure_models[0].concentration_model)