diff --git a/benchmarks/scripts/analyze.py b/benchmarks/scripts/analyze.py
index 8f50b3003..8e0772cd1 100755
--- a/benchmarks/scripts/analyze.py
+++ b/benchmarks/scripts/analyze.py
@@ -234,7 +234,7 @@ def qrde_hd(samples, precision=0.01):
     quantiles = np.linspace(precision, 1 - precision, num_quantiles - 1)
     hd_quantiles = [min_sample] + list(hdquantiles(samples, quantiles)) + [max_sample]
     width = [hd_quantiles[idx + 1] - hd_quantiles[idx] for idx in range(num_quantiles)]
-    p = 1.0 / num_quantiles
+    p = 1.0 / precision
     height = [1.0 / (p * w) for w in width]
     return width, height