feat: add region parameter for regional_misfit_uncertainty

mdtanker · Nov 19, 2024 · ef5239e · ef5239e
1 parent 0cb4a21
commit ef5239e
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/invert4geom/uncertainty.py b/src/invert4geom/uncertainty.py
@@ -504,6 +504,7 @@ def regional_misfit_uncertainty(
     runs: int,
     sample_gravity: bool = False,
     parameter_dict: dict[str, typing.Any] | None = None,
+    region: tuple[float, float, float, float] | None = None,
     plot: bool = True,
     plot_region: tuple[float, float, float, float] | None = None,
     true_regional: xr.DataArray | None = None,
@@ -526,6 +527,8 @@ def regional_misfit_uncertainty(
     parameter_dict : dict[str, typing.Any] | None, optional
         dictionary of parameters passes to `regional_separation` with the uncertainty
         distributions defined, by default None
+    region: tuple[float, float, float, float] | None = None,
+        region to calculate statistics within, by default None
     plot : bool, optional
         show the results, by default True
     plot_region : tuple[float, float, float, float] | None, optional
@@ -604,6 +607,7 @@ def regional_misfit_uncertainty(
     stats_ds = model_ensemble_stats(
         merged,
         weights=weights,
+        region=region,
     )
 
     if plot is True: