fix: negbinom consistency plots now have the correct boundaries (#245)

SCECcode · Nov 16, 2023 · d764e5f · d764e5f
1 parent 65451fb
commit d764e5f
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