Compute IQR more manual

shakyground2/workflows.py

@@ -258,7 +258,8 @@ def shakemaps_from_quakeml(
                 "maximum": np.exp(np.max(mean_shakemaps[imt])),
                 "median": np.exp(np.mean(mean_shakemaps[imt])),
                 "minimum": np.exp(np.min(mean_shakemaps[imt])),
-                "IQR": stats.iqr(mean_shakemaps[imt]),
+                "IQR": np.exp(np.quantile(mean_shakemaps[imt], 0.75))
+                - np.exp(np.quantile(mean_shakemaps[imt], 0.25)),
             }
         if export_folder:
             filestem = "{:s}_{:s}".format(results.earthquake.id, imt)