fix hpd_plot interpolation (#1062)

arviz-devs · Feb 10, 2020 · c2a2589 · c2a2589
1 parent ce0cf87
commit c2a2589
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/arviz/plots/hpdplot.py b/arviz/plots/hpdplot.py
@@ -102,8 +102,8 @@ def plot_hpd(
             smooth_kwargs = {}
         smooth_kwargs.setdefault("window_length", 55)
         smooth_kwargs.setdefault("polyorder", 2)
-        eps = np.finfo(float).eps
-        x_data = np.linspace(x.min() + eps, x.max() - eps, 200)
+        x_data = np.linspace(x.min(), x.max(), 200)
+        x_data[0] = (x_data[0] + x_data[1]) / 2
         hpd_interp = griddata(x, hpd_, x_data)
         y_data = savgol_filter(hpd_interp, axis=0, **smooth_kwargs)
     else: