From e9e4381fc1a4b7d803df11315012e5a37dbf820b Mon Sep 17 00:00:00 2001
From: Levi John Wolf <levi.john.wolf@gmail.com>
Date: Thu, 4 Jul 2024 12:32:34 +0100
Subject: [PATCH] add start of testing specifications

---
 esda/tests/test_moran_local_mv.py | 57 +++++++++++++++++++++++++++++++
 1 file changed, 57 insertions(+)
 create mode 100644 esda/tests/test_moran_local_mv.py

diff --git a/esda/tests/test_moran_local_mv.py b/esda/tests/test_moran_local_mv.py
new file mode 100644
index 00000000..981f6420
--- /dev/null
+++ b/esda/tests/test_moran_local_mv.py
@@ -0,0 +1,57 @@
+import numpy
+import geodatasets
+import geopandas
+from libpysal.weights import Queen
+from sklearn.linear_model import TheilSenRegressor
+from esda.multivariate_moran import Partial_Moran_Local, Auxiliary_Moran_Local
+
+df = geopandas.read_file(geodatasets.get_path("geoda.lansing1"))
+df = df[df.FIPS.str.match("2606500[01234]...") | (df.FIPS == "26065006500")]
+
+y = df.HH_INC.values
+X = df.HSG_VAL.values
+w = Queen.from_dataframe(df) 
+
+def test_partial_runs():
+    """Check if the class computes successfully in a default configuration"""
+    m = Partial_Moran_Local(y,X,w, permutations=1)
+    # done, just check if it runs
+
+def test_partial_accuracy():
+    """Check if the class outputs expected results at a given seed"""
+    numpy.random.seed(111221)
+    m = Partial_Moran_Local(y,X,w, permutations=10)
+    # compute result by hand
+
+def test_partial_unscaled():
+    """Check if the variance scaling behaves as expected"""
+    m = Partial_Moran_Local(y,X,w, permutations=0, unit_scale=True)
+    m2 = Partial_Moran_Local(y,X,w, permutations=0, unit_scale=False)
+    # variance in the partials_ should be different
+
+def test_partial_uvquads():
+    """Check that the quadrant decisions vary correctly with the inputs"""
+    m = Partial_Moran_Local(y,X,w, permutations=0, mvquads=False)
+    ...
+
+def test_aux_runs():
+    """Check that the class completes successfully in a default configuration"""
+    m = Auxiliary_Moran_Local(y,X,w, permutations=1)
+    ...
+
+def test_aux_accuracy():
+    """Check that the class outputs expected values for a given seed"""
+    numpy.random.seed(111221)
+    m = Auxiliary_Moran_Local(y,X,w, permutations=10)
+    ...
+
+def test_aux_unscaled():
+    """Check that the variance scaling behaves as expected"""
+    m = Auxiliary_Moran_Local(y,X,w, permutations=0, unit_scale=True)
+    m2 = Auxiliary_Moran_Local(y,X,w, permutations=0, unit_scale=False)
+    ...
+
+def test_aux_transformer():
+    """Check that an alternative regressor can be used to calculate y|X"""
+    m = Auxiliary_Moran_Local(y,X,w, permutations=0, transformer=TheilSenRegressor)
+    ...
\ No newline at end of file