wip: add statistical tests between histograms

pyproject.toml

@@ -72,6 +72,7 @@ test = [
     "pytest-mpl >=0.12",
     "dask[dataframe] >=2022; python_version>='3.8'",
     "dask_histogram >=2023.1; python_version>='3.8'",
+    "numba_stats",
 ]
 dev = [
     "pytest >=6",
@@ -83,6 +84,7 @@ dev = [
     "ipykernel",
     "dask[dataframe] >=2022; python_version>='3.8'",
     "dask_histogram >=2023.1; python_version>='3.8'",
+    "numba_stats",
 ]
 docs = [
     "pytest >=6",
@@ -176,6 +178,7 @@ messages_control.disable = [
   "disallowed-name",
   "cyclic-import",
   "redefined-builtin",
+  "too-many-public-methods",
 ]
 
 [tool.ruff]

src/hist/basehist.py

@@ -630,3 +630,43 @@ def sum(self, flow: bool = False) -> float | bh.accumulators.Accumulator:
             raise AssertionError(f"Unsupported storage type {storage}")
 
         return super().sum(flow=flow)
+
+    def chisquare_1samp(
+        self,
+        distribution: str | Callable[..., Any],
+        flow: bool = False,
+        args: Any = (),
+        kwds: Any = None,
+    ) -> Any:
+        from hist import stats
+
+        return stats.chisquare_1samp(self, distribution, flow, args, kwds)
+
+    def chisquare_2samp(self, other: hist.BaseHist, flow: bool = False) -> Any:
+        from hist import stats
+
+        return stats.chisquare_2samp(self, other, flow)
+
+    def ks_1samp(
+        self,
+        distribution: str | Callable[..., Any],
+        args: Any = (),
+        kwds: Any = None,
+        alternative: str = "two-sided",
+        mode: str = "auto",
+    ) -> Any:
+        from hist import stats
+
+        return stats.ks_1samp(self, distribution, args, kwds, alternative, mode)
+
+    def ks_2samp(
+        self,
+        other: hist.BaseHist,
+        equal_bins: bool = True,
+        flow: bool = False,
+        alternative: str = "two-sided",
+        mode: str = "auto",
+    ) -> Any:
+        from hist import stats
+
+        return stats.ks_2samp(self, other, equal_bins, flow, alternative, mode)