fix: negbinom consistency plots now have the correct boundaries

SCECcode · Nov 15, 2023 · 3101a2f · 3101a2f
1 parent 1e2badb
commit 3101a2f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/csep/utils/plots.py b/csep/utils/plots.py
@@ -1856,9 +1856,9 @@ def plot_consistency_test(eval_results, normalize=False, axes=None,
             mean = res.test_distribution[1]
             upsilon = 1.0 - ((var - mean) / var)
             tau = (mean ** 2 / (var - mean))
-            phigh = scipy.stats.nbinom.ppf((1 - percentile / 100.) / 2., tau,
+            plow = scipy.stats.nbinom.ppf((1 - percentile / 100.) / 2., tau,
                                            upsilon)
-            plow = scipy.stats.nbinom.ppf(1 - (1 - percentile / 100.) / 2.,
+            phigh = scipy.stats.nbinom.ppf(1 - (1 - percentile / 100.) / 2.,
                                           tau, upsilon)
 
         # empirical distributions