py-geodetector 0.1.1__tar.gz

py_geodetector-0.1.1/.gitignore

@@ -0,0 +1,4 @@
+.vscode
+__pycache__
+tempCodeRunnerFile.py
+ref

py_geodetector-0.1.1/PKG-INFO

@@ -0,0 +1,63 @@
+Metadata-Version: 2.3
+Name: py_geodetector
+Version: 0.1.1
+Summary: A simple Python package for the geodetector
+Project-URL: Homepage, https://github.com/djw-easy/GeoDetector
+Project-URL: Bug Tracker, https://github.com/djw-easy/GeoDetector/issues
+Author-email: djw <djweasy@163.com>
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python :: 3
+Requires-Python: >=3.7
+Description-Content-Type: text/markdown
+
+# A simple Python package for the geodetector 
+
+# Install
+
+```
+pip install py-geodetector
+```
+
+# Usage
+
+A quick example of geodetector usage is given in the ./example.ipynb.
+
+```python
+from py_geodetector import load_example_data, GeoDetector
+
+# load example data
+df = load_example_data()
+
+gd = GeoDetector(df)
+# factor detect
+factor_df = gd.factor_dector()
+
+# interaction detect
+interaction_df = gd.interaction_detector()
+# or you can generate the interaction relationship as the same time
+interaction_df, interaction_relationship_df = gd.interaction_detector(relationship=True)
+
+# ecological detect
+ecological_df = gd.ecological_detector()
+
+# plot 
+# use a heatmap visualize the interaction detect result, 
+# red text means that the ecological detection results show a significant difference
+gd.plot(value_fontsize=14, tick_fontsize=16, colorbar_fontsize=14);
+```
+
+# Reference
+
+```
+@article{wang2010geographical,
+  title={Geographical detectors-based health risk assessment and its application in the neural tube defects study of the Heshun Region, China},
+  author={Wang, Jin-Feng and Li, Xin-Hu and Christakos, George and Liao, Yi-Lan and Zhang, Tin and Gu, Xue and Zheng, Xiao-Ying},
+  journal={International Journal of Geographical Information Science},
+  volume={24},
+  number={1},
+  pages={107-127},
+  year={2010},
+  publisher={Taylor \& Francis}
+}
+```

py_geodetector-0.1.1/README.md

@@ -0,0 +1,50 @@
+# A simple Python package for the geodetector 
+
+# Install
+
+```
+pip install py-geodetector
+```
+
+# Usage
+
+A quick example of geodetector usage is given in the ./example.ipynb.
+
+```python
+from py_geodetector import load_example_data, GeoDetector
+
+# load example data
+df = load_example_data()
+
+gd = GeoDetector(df)
+# factor detect
+factor_df = gd.factor_dector()
+
+# interaction detect
+interaction_df = gd.interaction_detector()
+# or you can generate the interaction relationship as the same time
+interaction_df, interaction_relationship_df = gd.interaction_detector(relationship=True)
+
+# ecological detect
+ecological_df = gd.ecological_detector()
+
+# plot 
+# use a heatmap visualize the interaction detect result, 
+# red text means that the ecological detection results show a significant difference
+gd.plot(value_fontsize=14, tick_fontsize=16, colorbar_fontsize=14);
+```
+
+# Reference
+
+```
+@article{wang2010geographical,
+  title={Geographical detectors-based health risk assessment and its application in the neural tube defects study of the Heshun Region, China},
+  author={Wang, Jin-Feng and Li, Xin-Hu and Christakos, George and Liao, Yi-Lan and Zhang, Tin and Gu, Xue and Zheng, Xiao-Ying},
+  journal={International Journal of Geographical Information Science},
+  volume={24},
+  number={1},
+  pages={107-127},
+  year={2010},
+  publisher={Taylor \& Francis}
+}
+```

py_geodetector-0.1.1/example.ipynb

@@ -0,0 +1,437 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from py_geodetector import load_example_data, GeoDetector"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>incidence</th>\n",
+       "      <th>type</th>\n",
+       "      <th>region</th>\n",
+       "      <th>level</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>5.94</td>\n",
+       "      <td>7</td>\n",
+       "      <td>5</td>\n",
+       "      <td>5</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>1</th>\n",
+       "      <td>5.87</td>\n",
+       "      <td>5</td>\n",
+       "      <td>5</td>\n",
+       "      <td>5</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>2</th>\n",
+       "      <td>5.92</td>\n",
+       "      <td>5</td>\n",
+       "      <td>5</td>\n",
+       "      <td>5</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>3</th>\n",
+       "      <td>6.32</td>\n",
+       "      <td>1</td>\n",
+       "      <td>7</td>\n",
+       "      <td>1</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>4</th>\n",
+       "      <td>6.49</td>\n",
+       "      <td>3</td>\n",
+       "      <td>2</td>\n",
+       "      <td>4</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "   incidence  type  region  level\n",
+       "0       5.94     7       5      5\n",
+       "1       5.87     5       5      5\n",
+       "2       5.92     5       5      5\n",
+       "3       6.32     1       7      1\n",
+       "4       6.49     3       2      4"
+      ]
+     },
+     "execution_count": 2,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "df = load_example_data()\n",
+    "df.head()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "gd = GeoDetector(df, 'incidence', ['type', 'region', 'level'])"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 4,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>type</th>\n",
+       "      <th>region</th>\n",
+       "      <th>level</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>q statistic</th>\n",
+       "      <td>0.385717</td>\n",
+       "      <td>0.637774</td>\n",
+       "      <td>0.606709</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>p value</th>\n",
+       "      <td>0.363236</td>\n",
+       "      <td>0.000117</td>\n",
+       "      <td>0.040804</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "                 type    region     level\n",
+       "q statistic  0.385717  0.637774  0.606709\n",
+       "p value      0.363236  0.000117  0.040804"
+      ]
+     },
+     "execution_count": 4,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "gd.factor_dector()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 5,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>type</th>\n",
+       "      <th>region</th>\n",
+       "      <th>level</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>type</th>\n",
+       "      <td>0.385717</td>\n",
+       "      <td>NaN</td>\n",
+       "      <td>NaN</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>region</th>\n",
+       "      <td>0.735681</td>\n",
+       "      <td>0.637774</td>\n",
+       "      <td>NaN</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>level</th>\n",
+       "      <td>0.663524</td>\n",
+       "      <td>0.713597</td>\n",
+       "      <td>0.606709</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "            type    region     level\n",
+       "type    0.385717       NaN       NaN\n",
+       "region  0.735681  0.637774       NaN\n",
+       "level   0.663524  0.713597  0.606709"
+      ]
+     },
+     "execution_count": 5,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "gd.interaction_detector()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 6,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>type</th>\n",
+       "      <th>region</th>\n",
+       "      <th>level</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>type</th>\n",
+       "      <td>NaN</td>\n",
+       "      <td>NaN</td>\n",
+       "      <td>NaN</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>region</th>\n",
+       "      <td>Enhance, bi-</td>\n",
+       "      <td>NaN</td>\n",
+       "      <td>NaN</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>level</th>\n",
+       "      <td>Enhance, bi-</td>\n",
+       "      <td>Enhance, bi-</td>\n",
+       "      <td>NaN</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "                type        region level\n",
+       "type             NaN           NaN   NaN\n",
+       "region  Enhance, bi-           NaN   NaN\n",
+       "level   Enhance, bi-  Enhance, bi-   NaN"
+      ]
+     },
+     "execution_count": 6,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "gd.interaction_detector(relationship=True)[1]"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 7,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>type</th>\n",
+       "      <th>region</th>\n",
+       "      <th>level</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>type</th>\n",
+       "      <td>NaN</td>\n",
+       "      <td>NaN</td>\n",
+       "      <td>NaN</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>region</th>\n",
+       "      <td>Y</td>\n",
+       "      <td>NaN</td>\n",
+       "      <td>NaN</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>level</th>\n",
+       "      <td>Y</td>\n",
+       "      <td>N</td>\n",
+       "      <td>NaN</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "       type region  level\n",
+       "type    NaN    NaN    NaN\n",
+       "region    Y    NaN    NaN\n",
+       "level     Y      N    NaN"
+      ]
+     },
+     "execution_count": 7,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "gd.ecological_detector()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW8AAAEoCAYAAACXYXDAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAy8UlEQVR4nO3dd5wU9f3H8dd777iDo0mRIogIiqBiFEsUa1QMllhjbDGgURIsUbH97BoLMYlEJaJiNFiwx2gSMRYSECuCsWAvSBFEeufq5/fHzOLecre7yN7t7O3n6WMesDPf2fnucH7uu5/5FpkZzjnn8kss1xVwzjm38Tx4O+dcHvLg7ZxzeciDt3PO5SEP3s45l4c8eDvnXB4qznUFGkPvvXa3E0ePzHU1mqQbdzs411VwbmMp1xXIhoJoea9ZvjzXVXDOuawqiODtnHNNjQdv55zLQx68nXMuD3nwds65POTB2znn8pAHb+ecy0MevJ1zLg958HbOuTzkwds55/KQB2/nnMtDHrydcy5BadfDrLTr4ZFfH7IgJqZyzrlMVS6dhlSU62qk5S1v55wLSdpdsVJASNox1/VJxVvezjkXipVtNbW41TaYVaOiFu8T4eljveXtnHOApG5WtRqVdCRW2pmaimVI2jzX9aqPt7ydcw4oat1vbqxZW6SgsV3UqjdUb/EtEW19e8vbOVfwJJXVrJtPrKzH+n1FZVtRvXYOkkpzWLV6efB2zhW84na7rS4q64H0XUhUrJhYi+40a//DdZm8h6T9JP1D0teSTNLQDM7pL2mypLXheVcr3vRPw4O3c66gSYpVr/qCopa9NzhW3GobqlZ9ToYBtRUwAzgPWJvBddsALwILgN2B3wAXAyMyqbfnvJ1zBa1Zx/2qa9bOQ0UbZkdU1AIVt6a47U41pMl9m9kEYAKApHEZXPoUoAwYYmZrgRmS+gEjJI0ys5QDhbzl7ZwraNWrPqWodZ96jxe33o7qVZ82xKX3AqaEgTvueWALoGe6kz14O+cKlqQOVlNFrFmbesvEStph1WuRdK6kaQnbsE28fBeClEmiBQnHUvK0iXOukJUoVlLrQWVdpGYYPGhmo7N8/eTUiOrZvwEP3s65AieUNgnRIF29v2HDFnan8M/kFvkGPG3inCtoEkixlFsDDdN5HdhXUvOEfYOAecBX6U724O2cK3BKG7yVQfSW1ErSzpJ2JoitPcLXPcLjIyVNTDjlYWANME7SjpKOBf4PSNvTBDx4O+cKnpCKUm4Z2g34X7i1AK4L//7b8HhXYH1ncjNbTtDS3gKYBtwB3AKMyuRinvN2zhW0eNokTam072Nmk1IVNLOhdex7H9gv7ZvXwYO3c67AKX3wjuDUVB68nXMFLn1vk0xy3o3Ng7dzrqCJTNIm0ePB2zlX2CRisXSh0FvezjkXMUqbFvG0iXPORUxGaZPoxW4P3s65QpdBb5MIRm8P3s65wpZBP+/ohW4P3jmzeuky/nH171j4+UzWLl9JWbvN2Ha/PTng7F/SvFXL9eU+fHEyr/31YZbM/pqydm3Z7WdHseepP8thzZ1rWkQMyR9YugxJMbY7YG8OOOt0yjZry5I5X/P87//MhBtHcezIqwD44tWpPHPVTRxy4dn02mt3Fn81mwk3jqK4tITdfnZ0bj+Ac01I+py3B28XKtusDQOO+8n61227dmbXn/6E18Y9un7f+8+9xLb77sWuxx8JQLvuXdlr6Em8fv/j7Hr8UWS4TqlzLqVMpoSNnvyrcRO1cuEiPvnvK/QYsNP6fdUVlRSXNKtVrllpCSu/Xcjy+Wmn+3XOZSKDKWGj2FUw8sFbKab0ynBF50h7+oob+f0+RzD6sJMoKSvjiKsvWn+s15678enk1/jyjWlYTQ2LZ83lzfFPArBq0ZJcVdm5JkUIKfUWRZFOm0gqNrMqSS2B4UAfYC7wvpn93cxMkuqa+zZcX24YQJsunZIPR8bBFwxnnzNPZcmsuUy64z5eHHUnh11+PgA7H3MYS7+ex5MXXUN1VRWlLVuy+4nHMGXsA8SKIv9717k8IWLpHlhGMIBHNnhLioWBuzUwleBbwgrgMKCZpEPMbHh9k5ab2VhgLEDX7fukndg8V1p1bE+rju3p2LMHLdq24cEzL2Cf00+mTZdOSOLAc8/kgLNOZ/XipZS1a8tXU/8HBDly59ymy2SQjqdNNoKZ1UhqBjwKzAeOMrPdga2B94BfSdonl3XMNqupAaCqsrLW/lhREa07daSoWTM+eOG/dOu/PS3bt8tFFZ1repR+JZ0oimzLO9SRYOWJPwJfhvuOAPYHLjazVySVmll5riq4MaY9/jTTHv8Hv37yPj6b8gZrl6+gS99tKSlrwaIvZzHx9rF069+P9lt2A2DNsuV89NLLbLXrTlRXVPLuP5/n44kv8/O7b8nxJ3GuKcnZAsSbJOrBe0uCPPdHZlYh6WTgIeAKM7tFUhlwhaTnzOyVnNY0A2uWrWDJrDkAFJeU8L+n/sWimbOprqykTefN6XPA3uw15MRa58yY8CL/uX0sGHTr349T7vojW+zQNxfVd67pimjrOpWoB+/5wNfA/pL6APcSBO6R4fH9gR8Ck3NUv42y37BfsN+wXwCw9Q8HsPUPB6QsX7ZZW4bcd3tjVM25giWJWCz1OpVR7HESieAtqcjMqpP3m9kcSdOB6wEDrjWzkWH3wa2BK4CFwEuNWmHnXBPiaZPvJak74AVAT+Az4AMz+xdwAvAscGBQXHsAuwCnAc2BA8KHmzEzq8nJh3DO5bWoPpRMJefBOwzcrYC3gDbAauA4oKWkm83sKkmHA48BpwDXANMJAvwvwvOLzawqRx/BOZfXFMl+3OnkLHgnpUqGEeS2jzezGZJ2JGhZXyEJM7sKOFrSFkAnYCawIhyk44HbOff9ifSdpiMY23MWvM2sOuwtcjmwAzDZzGaEx2ZIuoXgll4i6XUzm2Bm84B58fcIR1d64HbObZpY/kXvRk/0KBS+PB8YAewHfBUeLwEIA/V9wDqgX13vVd/oSuecy1i85Z1qi17sbtzgHT5UNCDeL2cU8BegHTBCUtuwP3c8gL8PLCYYqOOccw1AmFJvUdSowTvsFVICvCXpFDNbB1wE3Ab0Au5JCOCStAvQirBV7pxzWacMtgjKRc57S6AtsIekJ8JAfWlYlyHAZEl/DssdSDCL4Kgc1NM5VyhiaSJ0BAN4o+e8zewL4G8E3f7ahfsqgAsJRlDuBPwB2IZgKPzAeHfAxq6rc64AiCB4p9oiGL0bNHgrqed7PJcN3EPQn/uCMD1SHAbwS4A7gGVANfComa2T1Nx7lTjnGky64B292N0wwTshINdIKpW0M6xvYQPMIpjW9eBwf5WkZmZWSdACnxAeu01SuzA37pxz2Sel3yIoq8E73tK2QHwhhXeApyQ9I2l3SV3CKVxvAH4AnBmeUxkO3KkgGCb/BHASMLIpLHfmnIuwPHxombXgHQbuYZL+lPD6bYJBNfcC3YHngSclnQpUAJOAfSSVhd0IqxMC+CUEvVBGeX9u51yDEVCk1NvGvJ10lqSZktZJmi5p3zTlfyzpdUkrJS0KG7p90l0nmy3vEqALMFzSw8DnBF38fmpmN5rZrsDF4b5xBA8l9yJYXKF3fFKphABebmaXmNmnWayjc85tIG0/7wy//Es6gaDReRPBBHqvAc9J6lFP+a2BZ4ApYfmDgRYEqeOUsha8w7z07whGRZ5IEMxPMrOl4XJmmNm9ZvZzgqD9JvABsBlwfsLDTOqaHtY55xqE0jysTNeNsLYRwDgzu8fMPjKzcwnWJRheT/ldgWbAZWb2uZm9A4wEekvqmOpCWc15hwG8A8EDydYEy5fF89nNYP18JFOBqwgWUxgHHEQQ7J1zrvFlId8dNkB3BV5IOvQCMLCe06YBlcAZkorC54RDgLfMbFGq6zVEb5PzgMHAX4GjJN0PtR5IWnzu7TDYX0cwIOfYBqiLc86ll7bHCQBDJE1L2IYlvUtHgqk/FiTtX0CQUt6AmX0FDCKIg+XAcqA/QTo5pawHbzP7xsw+AW4E7geOTAjg1eG0rrdLii9/XgYsAnwhBedc48vogaUA7jez3RK2sfW8Y3IHC9WxLzggdSHo0PEAsDtwALASeDx5nEyyBhu1aGYLJd0Uvhwq6e/A3QTpks4Eee4i4BigPfB6Q9XFOedSSvdAMrPUySKCwYXJrexObNgajzsbWG1ml3xXFf0cmEOQaql3YfUGHWFpZt8StMDHAPsADxP8BuoXjpgsJViD8gfhsHnnnGtcWRqkE3Zxnk6QBkk0iKDXSV3KCAJ+ovjrlPG5wec2MbOFBE9P9wV+Buwff4BpZmuAv5jZhw1dD+ecq1e6+bwzN4og03CGpH6SbgO2AO4CkDRS0sSE8s8CAyRdI2lbSQMInhfOIfhFUK9GmezJzFYCH4dbfAm0yvCYD8BxzuWOyFbaBDN7TFIH4EqgKzADOMzMZoVFupKwPoGZ/UfSyQSDEi8G1gJvAIPNbHWqa+Vkpj7vx+2ciwyBpRlFaRsxPt7MxhCkius6NrSOfY8Cj2Z8gZBPs+qcK3AZ5LUjOLeJB2/nnItgcE7Hg7dzrrDFF2NIVyZiPHg751wezjrtwds5V9iUwbSvEQzuHrydc4Utk7RJBHnwds4VNAMs/2K3B2/nnPMHls45l28yGWEZQR68nXOFTYJif2DpnHP5J4LBOR0P3s65wuaDdJxzLh+FK8TnmYII3ksWVvPIXUtzXY0m6cG+/8l1FZqk2RcdmOsqFA6Rfs7uCMb2ggjezjlXLwFFaaJ3BFvmHrydc85HWDrnXB7Kv9jtwds5V+AkzHubOOdcnsloYqroRW8P3s65wiYymBK2UWqyUTx4O+cKXAZrWEaQB2/nnEub845ecPfg7ZwrbL4Yg3PO5R/LYHh8FBdr8ODtnCtsmTywjCAP3s65wpbRrILRC+4evJ1zzgfpOOdcnhGRDM7pePB2zhW4DIbHR5AHb+dcYfMpYZ1zLk9FLzanlW79COeca9IkiMVSbxvT8JZ0lqSZktZJmi5p3zTlJel8SR9LKpc0X9Lv0l3HW97OuYKXrayIpBOA24CzgFfCP5+TtL2Zza7ntFuAI4CLgfeBtkDXdNfy4O2cK2zZnZdqBDDOzO4JX58raTAwHLhsg0tL2wHnAjuZ2UcJh/6X7kKeNnHOFTQhYrHUWybBXVIJsCvwQtKhF4CB9Zx2FPAlMFjSl5K+knS/pE7prufB2zlX2MKWd6otQx2BImBB0v4FQJd6zukFbAWcCAwFTgX6Av+UlDI+e9rEOVfwUodJ4r1Rhkg6NWHvWDMbW0dpq+Ps5H1xMaAUONXMPgUIr/EJsDvwZn1V8uDtnCtoIuPW9f1mdluK44uAajZsZXdiw9Z43HygKh64Q58BVUAPUgRvT5s45wqbgqlNUm2ZxHYzqwCmA4OSDg0CXqvntFeBYkm9E/b1ImhYz0p1PQ/ezrmCJtL3896IQTyjgKGSzpDUT9JtwBbAXQCSRkqamFD+JeBt4D5Ju0jaBbiPoMU9LdWFPG3inCtsAqXJm2Qau83sMUkdgCsJ+mrPAA4zs3gruivQO6F8jaQjgNuBl4G1wIvACDOrSXUtD97OuYImMn5gmREzGwOMqefY0Dr2zQeOz/wKAQ/ezrmCF8F5p9Ly4O2cK2wZ9OWOYmz34O2cK2iZzAgbxejtwds5V9i85e2cc/lnIwbpRIoHb+dcYRPIFyB2zrn8E0u3ClrjVGOjePB2zhU0T5s451w+CucvSVcmajx458hxn07iD1PqHITF0UfexHubb1NrX8/l8/nHM5ciM/oPebAxqpj3Vk9/kZWTn6By4WxiJS1o3veHdDjp/zYoV7lwLgtu+zWY0f3GZ3NQU5dLmbS8Ixi7PXjnyr96DWRy951r7bts6oPs+u0nvNexd639zaqruO2/t/JW537s8c2HjVjL/LXyladY8Z9H2OzwYZRstT1WWU7VwrkblLOqShaPv4HSrXei/Mt3c1BTl3PKYHh8BHnwzpHy4hLKi0vWv25eVc5Bc6Yztv+RGzQDLn3rIT5pvxVvdtneg3cGatauYvmEe+k49Dqa99lt/f6Srr02KLtswj0069qL5r08eBeq+DJoqQtFr+2dh79vmqbDv3ydFpXlPNHnR7X2/2j22xw4522u2/O0HNUs/6z7dBpm1VSvXMr8P5zGvOt/xqJxV1O1eF6tcms/eoN1H71Bu6POyVFNXVSkWwYteqHbg3dknPjJS/ynxwAWlrVbv2/zNUu56dW7GbH/OawuaZHD2uWXqsXzwYwVE8ez2U+G02HIb7GaKr6960JqKtYBUL1iMUueHEX7E/+PWPOyHNfY5VT21rBsVJEP3pKKcl2Hhrbt0jns+u2nPLbdQbX2/2nSaMb3HcQ7nfrkqGZ5ymqguop2R51Di757UNqjLx1OupyaVctY9+HrACx+ZCSt9vwJpVttn+PKuihIG7wjGMAjnfOWVGRm1ZJaAr8CugH/Ad41sw2fPtU+dxgwDKCoZccGr+umOPHjl5jXssMGDzAHzp/BHt98yG/+9yQAwigy49P7TuTqgWfwaN+Dc1Db6Iu16QBAceetvtvXohWxNh2oWvYtAOWf/4/yL99lxUsPBAUMsBrmXDqIdsecR6s9j2jsarscEem7CkYwdkc7eIeBuxXBkkDtgBLgAuBxSX8ws+kpzh0LjAUo3bx3fSs351xJVQXHfP4y9+9wKJb0yHvwMX+s9frg2dM4+52nOObIm/imZfvGrGZeKe25AwBVC+dQvNnmANSUr6Vm5RKK23UGoPOFf6l1zroPXmPFxPF0+s0dFLeJ9i97l12SURxLHSKk6IWQSKZNklIlw4DZBIt4dgbOBw4BrpW0e+PX7vs79cN/8+KT59fad+hXb9C6cg1P9Dlwg/Kftu9Ra1tQ1p4aiU/b92BFaatGqnV+WPnq08z//VAAmm2+JS12GMjSZ+6gfOYMKhd8xZLHf0+s1WY077cnACVdtq61FbXpABIlXbYmVtY6h5/ENbZ4yzvVFkWRbHmHLe4yYATQB3gO+NDMDLhdUgVwA3C1pN+a2Vs5rG7G2q9bQe/ltXs8nPjJRKZ025l5rby1tylqVi+nauGc9a/bn3gZy/45hkV/vRLDKO25I52G/ZFYSfMc1tJFVR5O542CeBg9kq4FLgVWACea2X8lNTezdeHxXwPXA68BvzezV+t7r9LNe1u3o37XCLUuPFV9O+S6Ck3S7Is2/CbmsmZ9LJbUtdMP95q3x/U3pzzh9UvOZ/E7b7czs2UNXblMRbLlHRoFbAGcAVwm6XUzWyepxMwqzOwuBYmoO4FPgHqDt3PO1ccfWG6CeK+SxH1mtkLSCIL7dgwwWtK5SQH8bkkLgH/mot7OufwnoDj/BljmPnhLKjazKkmlwG5AT2AGsNjM5kq6kKCeh4fl4wG81MzKzezpxPfJzadwzuWzKPYmSSenwTtscVdJag28SBC4OwE1wMuSrjOzyZLOA24lCODVki4ws7WJ7+WB2zn3fSiDHiURbHjntqtg2KukOTAZWA2cBmwJXAl0Be6XtL+ZrQDOA54l6Dr4mxxV2TnXBMXSbFHUqC1vSTEzq0nafTBB/+0RwJQw9/07SV8Q9Ca5UtJnZjYvTKF8CfwR55zLAmHE0qRNCrrlLWlv4Jwwt52oJ9AFeCtsiZcAmNkTwAPA/sBm4b4VZjYyLNfk5zxxzjW8+APLVFsUH1g25jeCq4Ctzaw8af9cgpkljgMws4p4AAfGE3w72GD2oOTeKc45932onlGVPsIyZGaDAcKRk0cDz5rZcmA68C0wTNLHZjY1DOBFwA+ABcCcet7WOec2Wfq0SfR6ozRKy1uq9aXjFuAhYIikdmY2BzgFGAD8SdIQSW2BHxGMsPwYyIvh7865/JOvc5s0ePAOuwOu/7VlZsOB54GRBAF8MzP7L3AYsBVwD7AUGBeecoiZ1XiO2znXUNL1Noli/G7QtIkkJUwy1T4+B7eZHSrpBeAmwCTdb2aTJO1JkCrpBnwFTAzP9wE4zrkGERMZTAnbSJXZCA3W8o63uCXFgLuA2ZJ6xI+b2SHAKwQt8KFhC3yumT1rZmPN7IV4rxIP3M65hpTNtImksyTNlLRO0nRJ+2Z43raSVkpalVGdN65amQlbytWSWgDHA/F5OF+RtHW8XBjApwA3Ar8Ic921eK8S51xDEtlLm0g6AbiNIKuwC8Gsp88lNlzrOa8EeBR4OdN6Zz14h6mSqnAFnP8BpxMMd/8n0B14U1LveHkz+zFBhW8Ffpzt+jjnXGrBIJ1U20YYAYwzs3vM7CMzOxeYDwxPc97NwHvAE5leKOvBOyFV8ldgLTDMzE40s6OAnwPLgdcl9Uw451CC31ZPZbs+zjmXSib9vDPJeYet512BF5IOvQAMTHHe4cARbOS0Hw2V824G9AKmmdmseE8RM3sYuBboCExJCuAXhC32nM906JwrHBmNsAyKDpE0LWEblvRWHYEigrEpiRYQjCLf8NpSV4Iedqea2cqNqXdDBcoqoE24kdhjxMzGSzqJoGvgFEn7hAE+ZmY1/nDSOdfYMpwS9n4zuy2Dcslvpjr2xT0E3Glmb2RSgUSb3PIOUyS1hA8ZHwcODAM1Yas63le7gqDSS4Gbw7m5kyescs65BpdR2iSzt1oEVLNhK7sTG7bG4w4ErpFUJakKuBdoGb5ObtnXsknBO2xN10gqltRFUs+EeUmeAL4GLpZ0crivJnxY2QF4jKCr4A8BX67bOZcTmfQ2yYSZVRBM9zEo6dAggl4ndekP7JywXU3wrHBn0jy8/N5pk6SFFMYD2wFtgW8kXWRmLyUsonC7pOOBlcAewDoze1ZSP4J743lu51zOpJ3bJPMeJ6OAByVNJVhX99cEa/HeFbyPRgJ7mNlBAGY2o/Z1tBtQk7y/Lt87aIZ57JbAm8BiYDRBXudQYIKk881sjKTTCNagPB5YBbwBnBHOd3Io8BFBUHfOuUYnoDhN8zrTft5m9pikDny3oMwM4DAzmxUW6Qr0ru/8jbHRwTvsxx3/NXQZQUD+pZl9Gh5fS/Awcl1Y9h3gnXBJs5qwTBeCfo27Anub2epN/yjOObfxRNBFJF2ZTJnZGGBMPceGpjl3HN/N65TSRgXvhMWCW4fdWrYElgCzwuMnEXR7uczM7pPULoz1ywiftko6jKDDel9gfzP7YGPq4Jxz2RQ8sEyTFong3CYbFbwTRk5OlvR3oB8w18zKw8A9HrjczG4O+2tfCxSFKZR4F8AXCR5Qvm5ms7P2SZxz7nvKxwWIMwreSbP63U6QKrmX4Kno1ZJuA84BrjCz34Xl+ofblPi5YV/uSoKeJs45l3Px+bzzTUbBO6HF/SOCVW9eDBcEfpHgQeRZwKNmNjLsKtiXIOdTCfwW1ufKvS+3cy5SYoJmaR5YRjG4Z9ryFnANcCHB6MmjAczsXUnXEzxZPVzSOII+3F0Jctz7JUzr6rMDOuciaKMnn4qEjPqfh71LbiV4ChoD9ks49i+Clvd1BCvhLAYeBvY0s8r49LDZrbZzzmVHFkdYNqqMH1ia2deSLgfKgEskfW5mfwmPvQ28LenWxCXPfCEF51zUZdJVMIo2trfJN+GoySJgrCTiATyct8QIljWTBbzF7ZyLvCbb2ySRmS2QdE74cqykGjO7LzFQJ7a+nXMuymJAszxcw/J7DY9PCOA1wF8kLTSzf2a3as451wgyWKcygrF7k+Y2WSDpfOAz4Lms1cg55xpRk+7nXR8zm0/QTTB5II9zzuWNtC3vCAb3rE3F6oHbOZePJChKNyVsvQvh5I7Po+2cK2gxgnUqU4lgw9uDt3POFVzOO2+UFFHdo02ua9EkFc1anusqNEm97pyc6yo0WV8O37/Wa8nSp00iGNwLI3g751w9CrK3iXPONQUF1c/bOeeagkymhPW0iXPORVDa1eMbqR4bw4O3c66giQznxo4YD97OuYLnOW/nnMszQc47XVdBH2HpnHORkklXQW95O+dcBHk/b+ecyzOZPLCMYmz34O2cK2zKoB93BKO3B2/nXEHz4fHOOZenPG3inHN5xyLZFTAdD97OuYIm0resveXtnHNRk8kDywjKxyH9zjmXNQKKlHrbmOAu6SxJMyWtkzRd0r4pyh4g6RlJ8yWtkfSepNMzuY4Hb+dcwVOaLeP3kU4AbgNuAnYBXgOek9SjnlMGAu8DPwV2BO4Exko6Od21PG3inCtoIqtpkxHAODO7J3x9rqTBwHDgsuTCZnZT0q47Jf0IOA54ONWFvOXtnCt42Wh5SyoBdgVeSDr0AkELO1NtgKXpCnnwds4VvHTBOwzgQyRNS9iGJb1NR6AIWJC0fwHQJaN6SEcABwFj05X1tIlzrqApfCiZgfvN7LYMyiV3Glcd++qoh/YmSJX8xsympivvwds5V/DSDdLJMHWyCKhmw1Z2JzZsjSddX/sAE4CrzezOTC7maRPnXEHLMGWSlplVANOBQUmHBhH0Oqn7+tJ+wHPAdWZ2a6b19pa3c67gZXFWwVHAg5KmAq8Cvwa2AO4KrqORwB5mdlD4+gDgWWAMMF5SvNVebWYLU13Ig7dzrqBlcz5vM3tMUgfgSqArMAM4zMxmhUW6Ar0TThkKlAEXhVvcLKBnqmt58HbOFbxsTglrZmMIWtJ1HRtax+uhdZVNx4O3c66w5encJh68nXMFzWcVdM65POUr6TjnXJ7xlrdzzuWpWLqVdCIYvT14O+cKmz+wdM65/ONpE+ecy1P5OE+IB2/nXEHL8mIMjcaDt3Ou4ClN2zuKsd2Dd44d8+5Eznztb/ReNIc1JS3477a7M+LYi78rYMbpbzzNKdOeZculC1jeohVP7jyImwdltEZpQVr50X9YPPHPdR7revzNlHbelpqqChZPupuKhV9SuXQupV360vXY6xu5pvlr1ZsvsHziY1R+MxuVtqBshz3ZfOjlANRUlrP4kVGUz/6Mym9m0bz3jnS9IJNpsHNFSNma3aTx5FXwlhQzs5pc1yNbhr7xNGdPeYybDjmDt7v3o3lVOb0Wf12rzFXPj+XAT6dy0yFn8HGnnrQpX02nlUtyVOP80HLbvWnRY5da+5a+9gDl8z+mpNM2wQ6rQUXNaNP/UNbMepua8tU5qGl+Wv7fv7H8+fG0P+bXlG69PVZZQeW3c74rUFODiktoc8AxrJ3xBjVrV+WushmLXnBOJ/LBW1IZwQxdo8xsUVMJ4G3WruLSl/7KmSddwyu9B6zf/0nnrdf/vdeiOQx58x8MPutOPt/8u8WnP+jaqFXNO7HiUmLFpetf11SWs2bmW7QdcDQKk5uxZs3p+KNfA1CxeBYVHrwzUr1mJUufuYfOv7qBFv12W7+/pFuv9X+Plbag48kXAlDx9ReRD94K/8s3kQ/ewFnAxUB3Seeb2ZKmEMD3/eJtYlbD5quW8tKfh9F63Wre6bYdN/z4TOa0D6LzoI/fYHa7Luz/2TT+Ov5qZDW8udVO3HTIGSxutVluP0AeWf35q1hVOa36HZjrquS9tR9Ng5pqqlcsYe5vh1CzdjWlPfvS/rizaNZxi1xXbxPkX2fBfOgh8yfgduBHwGhJ7c2sRumTVJHWY+l8Ymac8/IjXP/jYQw78Wqa1VTx6LhLaV6xbn2Zbsu/5SczJnPh0RdywbGX0HvRHO59+BpUk9e/uxrVqg9epKznrhS3bJ/rquS9qkXzMDOW/ftB2h93Fp1/dT1WXcX8Wy+gJvy5zTsSUlHKLYKxO9rBW1KRmVUTtLyfBPYnwwAuaVh8lefq1csbq8oZi5lRUl3FtYcOZ/K2u/Fu9+0477hL6bh6GQd/+ub6Ms2rKrng2IuZ2rM/b221IxccezG7fP0JP5j3aY4/QX6oWDyb8m8+odX2yStTue/FDKqraH/8byjb4YeU9uxHp9OupGblUta8V+9KX5GnDP6LmkgHbzOrDgN4DXAh8AQZBnAzG2tmu5nZbkUt2zZmtTPybeugFfhZQi57ZfOWLGjdgW7Lvg3KtGpPZayImR27ry8zs0M3KmNFbLH828atcJ5a+cGLFLXqSIutdklf2KVV1LYDACVdt1q/L9aiFUVtO1K1NH9/JvMtcENEg7ekovjfw5Y3YQC/CPgbG9ECj6ppW24PQO/Fc9fvKytfS6dVS5i7WeegTI/taVZTTY8l89aX6bF0Ps1qqvm6befGrXAeqqmqYPUnk2nV78AMuoK5TDTvtSMAlQu+611Ss24N1csXU9w+X38m4wuhpdqiJ3K1SkiVIGk3SYdL6i9py3D/COAp8jCAD3nzH0wcfQYAMzt25/m+e3HNc3ex2+wP2PbbWfzx6VtY1LItE/vsAcArvXbh/a7b8Ien/8QO8z9nh/mf84en/8Tb3fvy3hbb5vKjRM6K9yYw96Fza+1b88Xr1FSsofX2B9V5TsWSOZQvnEn1uhVY5TrKF86kfOHMxqhuXlkx6SnmXncqAM06b0nZTnuz+InRrPvifSrmf8XCB2+mqPVmlPXfa/05FfO/onzOZ9SsWk5N+VrK53xG+ZzPcvURUgpGWCrlFsWkd6R6m4S9SOKBezxBgN48PPyKpNFm9rSkC4Bq4HjgT5IuNLNFual15tqtWc42i75raV9w7MVc/e+7gweQBm/12IGTh/yOdSXNAbBYjNNPuY5rJ9zJ4/ddzLpmJbzSawDXDx6GxSL/u6pRVa9dQdWy2n3kV37wIi167Exx683rPGfBP2+geuV3C3TPfyzo3tbznKcarqJ5qHrV8lot7c2HXsHiJ//MgjsvBzNKe/eny3mjiIU/twAL7riUqiUL1r+eN/JMALYeM6nR6p05pR1hGUUySzOPbQ5IugcYBFwAvAPsAPwe2BoYaGb/C1vafwCGAw8Bv66v+2Bptz7W5VejG6PqBSe2cE2uq9AkaXvvGdNQvhy+//pmtKSugw8bOO9vz/w+5TmHHXIek/87vZ2ZLWvo+mUqUi1vAEnbAnsD1wP/NrO1ktYCWxEE6U/i/bwlXQKUA3/N937fzrncSfdQMnpJkwjkvKUN5vPaEugLTA8Dd1/gA+BZ4DwzWwP8TNK2ZlZtZpebWTSTac65yMss5x09OQ3ekmRh3kZS/FvASmA10FHSdsBrwEvA6Wa2RtLBwAlAh1zU2TnXFCnNFj05TZskBO5xwHRgNPAJsBgYCfQCJgKnANWS2od/bwd8mYMqO+eaHCGK0paJmiikTVoBXYHTJHU3sxXAGUA3gjs2xsyqgF2AW4AjgbPNLH9HBDjnIiR1ykRSFGN37oO3ma0iGPq+DbB9uPtV4FcE3QHHSZoP3AvsBRxoZh/koq7OuaYqXdoketG7UdMmybMBSmpmZpVmdo+kM4GrgRfMbC3wT0k7Aj8GtgDeBj4ws6/rfHPnnPsegtDsK+mkFA/ckvqZ2UdmVpkQ0P8C/FbST83syTCwLwAeaMw6OucKTX6uH9/oaRNJY4Apkh6QtA3QMjz0L2ANcCxAGNijd8ecc02LIKZYyi2KclGrcQT5632AV4C7JQ00s3nAFcDRkg6D73qjOOdcw8lkYqrotSMbNHgnzg4Yvo6Z2VSCIL0z8CCwHcG8JQ8A/YH3gEGSmuOcc40gm7N5SzpL0kxJ6yRNl7RvmvL9JU2WtFbS15KuziTr0GA576TZAYcBvYEOkh4FpoZdAi+W1Bk4mqB3yb4Ew+BbEgR455xrYNnrTSLpBOA2guUbXwn/fE7S9mY2u47ybYAXgZeB3Qkas+MIBirekupaDRK8w5GT8cD9ZFipJUARcDpwr6S7zWxa+FDybkkvEgyLPxe4LBwG75xzDSo+PD5LRgDjzOye8PW5kgYTTKB3WR3lTwHKgCFhL7sZkvoBIySNSpU6znraJGnI+/UEgfvnwMFmthMwCvgl0CdxDm4z+9LMJpjZoWb2Trbr5ZxzdQtGWKbaMmmZSyoBdgVeSDr0AjCwntP2AqaEgTvueYLu0T1TXS8rwVtSK0lHQvCQUVKRpBYEswM+DrxpZovDGQN/DjwK/D1fFlFwzjVZy+bMmU9VVVW9BcrLK1i4cAnASfF1ccNtWFLRjgTZhQVJ+xcAXep5+y71lI8fq1e2Auc1wNOShsD6pcvaEqRBlplZRfhV4E1gMnBmOGPgJQQB3jnnGp2ZrR00aG+efvqless89tgEjjrqIMzszvi6uOE2tr63TXqtOvalK1/X/lqyFbxHEbSw/yppaLhvEbAQ6C9pc4Lk/UTgDDNbHbbCBwMDvPXtnMuVW265b+vRox+q85iZMWbMw9xww51bZPBWiwim9EhuMXdiw9Z13Df1lCfFOUCWgreZzQfOI1gc+D5JZ4STSY0hWKrsG+AZ4AQzWxkG88sIljh72hdScM7lipl91alTe6ZOfW+DYy+//BZ9+vSMx7h071NBMDvqoKRDgwimtq7L68C+SV2jBwHzgK9SXS9rLd6w18g5BAF8rKSTgbEEfbnXEnwV2FrSicDtBCMpTzGzWdmqg3POfR9PPvn8Prfeev8G+2+99X4efPCZXTfirUYBQyWdIamfpNsIHj7eBSBppKSJCeUfJhhZPk7SjpKOBf4PSNnTBLLcVdDMFkg6hyBQPwSsI+i//RVBYP8JsAL4HNjHzGZk8/rOOfc9vTZr1tfMnfsN3bsHWYwvvpjN8uUrMbO3M30TM3tMUgfgSoKprmcAhyU0UrsSjHmJl18uaRBwBzANWErQv3tUumtlvZ93GMDPJgjgTwInmdm1kkYDOxIsorDCzJZn+9rOOfd9mJk98sgo7rhjPCNHXgjA6NEP8pvf/OL7vNcYgpRxXceG1rHvfWC/jb1OgzwoDFMoZwFPAY9IOs3MFpvZZDOb44HbORc1J598YcmECZNZvXoNy5evZNKkqRxzzNmRW6Q9rsF6eYQB/GyCXij3hjlw55yLJDOrPOmkI3jwwWe4994nOe20Y+PdniOpQX+rhCmUCwhy3+805LWcc25TXXbZLR0GDNhhcXV1Ne+++3Gb8867MddVqleD968Ou9j80sw+bOhrOefcpjCzJQMH7sKgQQMxs5W5rk8qjZLPifJXD+ecSzR69IPRm7y7Dj6y0Tnn8pAHb+ecy0MevJ1zLg958HbOuTzkwds55/KQB2/nnMtDHrydcy4PefB2zrk85MHbOefykAdv55zLQx68nXMuD3nwds65PKQ0y6Q1CZIWAvm0VmZHgpWoXXb5fW04+XRvF5nZ4FxXYlMVRPDON5Kmmdluua5HU+P3teH4vW18njZxzrk85MHbOefykAfvaBqb6wo0UX5fG47f20bmOW/nnMtD3vJ2zrk85ME7yyQdLWlEruvh6ifJJF2b63pEhaRrJeX0K7ikSZIm5bIO+aZRFiAuMEcDBwOjclwPV7+9gLm5roRzm8KDt4s0SaVmVp7N9zSzN7L5fs7lgqdNskjSOGAI0C38am6S1kmqkHReHeWvlbRGUrvw9SRJr0g6StIMSeWSPpb0szrO/YGkf0haKmmtpFcl7dvgH7IBxb++S9pR0vOSVgGPSyqTdLOkmeG9nCnpCkmxpPMHSJoS3vM5ki6XdF1ySqCutImkwZJeD+/lcklPS9ouqUz83+dgSW+H/3YzJB3dQLckJyQVS7os/NkrlzRP0i2SmofHSyUtkXRLHeeeEN7fnRP27S9poqSVklaH/7Y7NuJHapI8eGfX9cAEYCHBV/P49jTwq8SCkoqAXwKPm9nShEPbALcDtwDHAp8Dj0r6UcK5A4DXgPbAmcBxwGLgJUm7NsQHa2TPAJOBI4HRwPPAGcBtwKHAX4CrgD/ET5DUEZhIcE9+AZwL/BgYmu5ikgYDzwKrgBOA4cCOwCuSuiUV7x3WYxTBv8984ElJ23yvTxpNDwFXAg8DhwMjCX5WxwOE34QeB04Of44T/RyYYWbvAEg6nODfZVV47GSgNTBF0pYN/kmaMjPzLYsbMA6Ym7TvAMCAfRP2HRnu2zNh36Q69hUBHwNTEvZNBD4CSpLKfQQ8net7sAn37trw85+XsO/UcN9+SWWvACqATuHrm8LX3RPKtAAWBD/mtc414NqE19OAz4DihH1bA5XAqKR/n0pg24R9nYBq4PJc379Nve/h3/cN788vksqcEu7fOXy9d/j6xwllNg/vzyUJ+z4HJia9VxuCeVBuTbq3k3J9L/Jp85Z3IzCzScCH1G59/wp4zzbMv85J3Gdm1cATwB6SYpJaAPuH+2rCr7jFgICXgP0a7pM0mr8n/H0wwaRir8U/a/h5XwCaAXuG5fYEXjez9Q8izWwtQYu6XpJaAgOAx8ysKuHcmcCrBPc60Wdm9llCuW+Bb4EeG/cRI2swwS/Bv9VxvyH8+TKzV4EvCH65xp1I8G1+PICkbQm+qYxPeq81wOs0jZ/VnPHg3XjuBH4qqYOkrQj+J7mrjnIL6tlXQtCyaU/Qyr6KoJWTuJ0DtEvOBeeh+Ql/7wRsxYafdWp4vEP4Z1eCIJqsrvuZqB3BL775dRz7huB+J1pSR7lyoHma6+SLTgQ/a6uofb/j97ZDQtmHgKMltQpfnwr8x8y+TngvgHvZ8N/viKT3chvJe5s0ngcIcodDCQLGWsIWSpLO9eyrIMiltwBqgDvC99yAmdVsenVzKvEB42JgJrDBQ9vQV+Gf8/kuWCSq634mWhper0sdx7qE1y8ki4F1BOmTusxL+PuDwDXAMZLeBHYneGCf+F4AlxF8K0xWsWlVLWwevLOvnCDA1mJmKySNJ0iXtAIeNrMVdZy/paQ946mT8IHQ8cDUMCivljQF+AHwdhMI1On8m+CB7Coz+zhFuTeAiyR1j6dOwhTT4ane3MxWS5oOHC/p2jBNRfjtaCDBA9NC8m/gUqCtmU1MVdDMvpD0OkGLuw+wGngqocgnBL9cdzCz3zVMdQuXB+/s+xBoL2k4wYOwdWb2fnhsDN/lvetKmUDwNf8xSdcQtLSHE/yPMTyhzAjgZeB5SfcStDo7EuRui8zs/7L4eXJtPHAaMDHsmvYuwdf63gQPfY82szUEvT+GE9yT6wh+iY4I/0w3evAqgtz4vySNIfjleh2wnKDXT8Ews0mSHiHoQTOKID1VA/QEDgMuNbNPE055gOBbYH/g72a2KuG9TNLZwDOSSgh6qCwi+DY0EJhtZj6Y7fvK9RPTprYBLYFH+O7r+FdJxz8B3qrn3EnAKwRBaQZB4PkEOKGOsv2ARwlykeUEIwb/ARyW63uwCffu2vCeFSftbx4e+zj8rEuAt8J9iT1EBoT3bx3wNUFQvg1YmvR+tXqbhPsGEzxEW0sQtJ8Btqvr36eOen8FjMv1/dvU+57wOgacR/CLcl14P94Ffk/QIk88tx3f/YI8pJ733wv4V/j/xLrwfj0K7JV0byfl+l7k0+azCjYiSX0IAtCZZnZvHccnEQSjfRq7bk1RmHJ6m2DZq4NyXR/nssnTJo1AUneCwTfXEaQ4Hs5tjZomSdcT9CueRdCT4QxgJ4Kv+841KR68G8cZwNXAp8DJFvQ/dtlnBPd5i/Dv7xHkxJ/Laa2cawCeNnHOuTyU74M5nHOuIHnwds65POTB2znn8pAHb+ecy0MevJ1zLg958HbOuTz0/1P7hrWQH5vxAAAAAElFTkSuQmCC",
+      "text/plain": [
+       "<Figure size 432x288 with 2 Axes>"
+      ]
+     },
+     "metadata": {
+      "needs_background": "light"
+     },
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "gd.plot(value_fontsize=14, tick_fontsize=16, colorbar_fontsize=14);"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "geo",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.10.12"
+  },
+  "orig_nbformat": 4
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}

py_geodetector-0.1.1/pyproject.toml

@@ -0,0 +1,29 @@
+[build-system]
+requires = [
+    "numpy",
+    "pandas",
+    "scipy",
+    "hatchling",
+    "matplotlib"
+]
+build-backend = "hatchling.build"
+
+
+[project]
+name = "py_geodetector"
+version = "0.1.1"
+authors = [
+    { name = "djw", email = "djweasy@163.com" },
+]
+description = "A simple Python package for the geodetector"
+readme = "README.md"
+requires-python = ">=3.7"
+classifiers = [
+    "Programming Language :: Python :: 3",
+    "License :: OSI Approved :: MIT License",
+    "Operating System :: OS Independent",
+]
+
+[project.urls]
+"Homepage" = "https://github.com/djw-easy/GeoDetector"
+"Bug Tracker" = "https://github.com/djw-easy/GeoDetector/issues"

py_geodetector-0.1.1/src/py_geodetector/__init__.py

@@ -0,0 +1 @@
+from .geodetector import GeoDetector, load_example_data

py_geodetector-0.1.1/src/py_geodetector/example_data/disease.csv

@@ -0,0 +1,186 @@
+incidence,type,region,level
+5.94,7,5,5
+5.87,5,5,5
+5.92,5,5,5
+6.32,1,7,1
+6.49,3,2,4
+6.46,3,2,4
+6.51,3,2,4
+6.7,3,2,4
+6.68,3,2,4
+6.65,3,2,4
+6.65,3,2,4
+6.6,3,2,4
+6.66,3,2,4
+5.86,5,5,5
+6.93,3,2,4
+6.5,3,2,4
+6.37,3,2,4
+6.64,3,2,4
+6.43,7,9,1
+7.01,2,9,1
+6.77,2,9,1
+6.36,7,9,1
+6.44,7,9,3
+6.52,7,9,1
+5.83,5,6,7
+6.32,7,9,1
+6.33,7,9,1
+6.46,7,9,1
+6.63,7,9,6
+6.53,3,9,1
+6.59,7,9,4
+6.58,2,9,1
+6.63,3,9,6
+6.61,2,9,1
+7.04,3,9,4
+5.89,5,6,5
+6.57,2,9,1
+6.53,3,9,4
+6.58,3,9,4
+6.53,3,9,1
+6.55,3,9,4
+6.57,2,9,1
+6.57,1,7,1
+6.49,2,7,4
+6.56,2,9,4
+6.57,1,7,1
+5.97,5,6,5
+6.54,2,7,4
+6.71,3,9,6
+6.58,3,9,1
+6.63,7,9,6
+7.42,3,2,6
+6.71,3,4,1
+6.82,3,2,6
+6.56,3,2,6
+6.82,3,2,6
+5.86,5,6,7
+6.8,3,2,6
+7.74,3,2,6
+6.92,3,2,6
+6.75,3,2,6
+7.06,3,2,6
+7.1,3,2,6
+6.69,3,4,6
+6.81,3,2,6
+6.99,3,2,6
+6.79,3,2,6
+5.98,5,6,5
+7.01,3,2,6
+7.2,3,2,6
+6.95,3,2,4
+6.74,3,4,6
+6.77,3,4,1
+6.73,3,4,6
+6.57,7,3,3
+7.55,7,3,3
+6.04,1,6,5
+6.28,7,3,3
+6.38,7,3,3
+6.21,7,3,3
+6.33,7,3,3
+6.19,7,3,3
+6.45,7,3,2
+6.19,7,3,2
+6.14,7,3,2
+6.2,7,3,2
+6.14,3,3,2
+5.95,1,5,5
+6.27,7,3,3
+6.14,3,3,2
+6.18,7,1,2
+6.09,7,1,2
+6.04,7,1,3
+6.01,7,5,2
+6.01,7,5,3
+6.09,7,1,3
+5.85,5,6,5
+5.88,7,5,5
+5.8,5,6,7
+5.78,5,6,7
+5.7,5,6,7
+5.82,5,6,7
+5.73,5,6,7
+5.79,5,6,7
+5.8,5,6,7
+6.45,1,8,3
+6.47,1,8,3
+6.11,3,8,1
+5.98,7,5,5
+6.12,1,8,3
+6.15,1,8,1
+6.04,3,8,1
+6.06,3,8,3
+5.95,3,8,1
+6.09,3,8,1
+6.11,3,8,1
+6.27,1,8,1
+6.28,1,8,1
+6.28,1,8,1
+5.96,1,5,5
+6.23,1,8,1
+6.53,1,8,1
+6.25,1,8,3
+6.14,3,1,3
+6.48,3,1,3
+6.02,1,1,3
+6.03,1,1,3
+6.13,3,1,3
+6.07,3,1,3
+6.09,3,1,3
+5.66,5,5,5
+6.2,3,1,1
+6.13,1,1,3
+5.99,1,1,3
+6.13,1,1,1
+6.01,7,1,3
+6.11,1,1,3
+6.17,1,1,1
+6.19,3,1,1
+6.18,1,1,3
+6.09,1,1,3
+5.74,5,5,5
+6.19,7,1,3
+6.29,1,1,3
+6.25,1,1,1
+6.24,1,1,1
+6.33,3,1,1
+6.23,3,1,3
+6.35,1,8,1
+6.32,1,1,1
+6.37,1,1,1
+6.21,1,1,1
+5.88,5,5,5
+6.2,1,1,1
+6.58,1,8,1
+6.44,1,8,1
+6.48,1,8,1
+7.4,2,8,1
+6.38,1,7,1
+6.44,1,8,1
+6.65,1,8,1
+6.86,1,8,1
+6.71,1,3,1
+6.1,5,5,5
+6.88,7,3,1
+6.72,7,3,1
+6.3,7,3,1
+6.65,1,7,1
+6.34,1,8,1
+7.22,1,3,1
+6.5,1,3,3
+6.52,1,3,1
+6.53,7,3,1
+6.43,7,3,1
+5.89,5,5,5
+6.5,7,3,3
+6.48,7,3,3
+6.5,7,3,1
+6.38,1,7,1
+6.18,3,1,3
+6.39,1,7,1
+6.45,1,8,1
+6.86,2,7,1
+6.29,2,7,1
+6.49,1,7,1

py_geodetector-0.1.1/src/py_geodetector/geodetector.py

@@ -0,0 +1,191 @@
+import warnings
+import numpy as np
+import pandas as pd
+from typing import Sequence
+from scipy.stats import f, ncf
+import matplotlib.pyplot as plt
+
+
+from pathlib import Path
+def load_example_data():
+    file_path = Path(__file__).parent / "example_data" / "disease.csv"
+    df = pd.read_csv(file_path)
+    return df
+
+
+def _plot_value(ax, interaction_df, ecological_df, value_fontsize=10):
+    length = len(interaction_df.index)
+    for i in range(length):
+        for j in range(length):
+            if not pd.isna(interaction_df.iloc[i, j]):
+                num = str(round(interaction_df.iloc[i, j], 2))
+                mark = num[-2:] if 3 == len(num) else num[-3:]
+                if 'Y'==ecological_df.iloc[i, j]:
+                    ax.text(j, i, mark, ha="center", va="center", color="r", fontsize=value_fontsize)
+                else:
+                    ax.text(j, i, mark, ha="center", va="center", color="k", fontsize=value_fontsize)
+
+
+class GeoDetector(object):
+    def __init__(self, df: pd.DataFrame, y: str, factors: Sequence[str]):
+        self.df = df
+        self.y = y
+        self.factors = factors
+        self._check_data(df, y, factors)
+        self.factor_df, self.interaction_df, self.ecological_df = None, None, None
+
+    def _check_data(self, df, y, factors):
+        for factor in factors:
+            if not factor in df.columns:
+                raise ValueError('Factor [{}] is not in data')
+            
+        for factor in factors:
+            # 检查列的数据类型
+            if df[factor].dtype not in ['int64', 'int32', 'int16', 'int8', 
+                                        'uint64', 'uint32', 'uint16', 'uint8', 
+                                        'object', 'string']:
+                # 如果数据类型不是整型或字符型，发出警告
+                warnings.warn(f"Factor '{factor}' is not of type 'int' or 'str'.")
+        
+        if y not in df.columns:
+            raise ValueError('Factor [{}] is not in data')
+            
+        for factor in factors:
+            if y==factor:
+                raise ValueError("Y variable should not in Factor variables. ")
+        
+        has_null = df.isnull().values.any()
+        if has_null:
+            raise ValueError("data hava some objects with value NULL")
+
+    @classmethod
+    def _cal_ssw(self, df: pd.DataFrame, y, factor, extra_factor=None):
+        def cal_ssw(df: pd.DataFrame, y):
+            length = df.shape[0]
+            if length==1:
+                strataVar = 0
+                lamda_1st = np.square(df[y].values[0])
+                lamda_2nd = df[y].values[0]
+            else:
+                strataVar = (length-1) * df[y].var(ddof=1)
+
+                lamda_1st = np.square(df[y].values.mean())
+                lamda_2nd = np.sqrt(length) * df[y].values.mean()
+            return strataVar, lamda_1st, lamda_2nd
+        if extra_factor==None:
+            df2 = df[[y, factor]].groupby(factor).apply(cal_ssw, y=y)
+        else:
+            df2 = df[[y]+list(set([factor, extra_factor]))].groupby([factor, extra_factor]).apply(cal_ssw, y=y)
+        df2 = df2.apply(pd.Series)
+        df2 = df2.sum()
+        strataVarSum, lamda_1st_sum, lamda_2nd_sum = df2.values
+        return strataVarSum, lamda_1st_sum, lamda_2nd_sum
+
+    @classmethod
+    def _cal_q(self, df, y, factor, extra_factor=None):
+        strataVarSum, lamda_1st_sum, lamda_2nd_sum = self._cal_ssw(df, y, factor, extra_factor)
+        TotalVar = (df.shape[0]-1)*df[y].var(ddof=1)
+        q = 1 - strataVarSum/TotalVar
+        return q, lamda_1st_sum, lamda_2nd_sum
+
+    def factor_dector(self):
+        self.factor_df = pd.DataFrame(index=["q statistic", "p value"], columns=self.factors, dtype="float32")
+        N_var = self.df[self.y].var(ddof=1)
+        N_popu = self.df.shape[0]
+        for factor in self.factors:
+            N_stra = self.df[factor].unique().shape[0]
+            q, lamda_1st_sum, lamda_2nd_sum = self._cal_q(self.df, self.y, factor)
+
+            #lamda value
+            lamda = (lamda_1st_sum - np.square(lamda_2nd_sum) / N_popu) / N_var
+            # F value
+            F_value = (N_popu - N_stra)* q / ((N_stra - 1)* (1 - q))
+            #p value
+            p_value = ncf.sf(F_value, N_stra - 1, N_popu - N_stra, nc=lamda)
+
+            self.factor_df.loc["q statistic", factor] = q
+            self.factor_df.loc["p value", factor] = p_value
+        return self.factor_df
+    
+    @classmethod
+    def _interaction_relationship(self, df):
+        out_df = pd.DataFrame(index=df.index, columns=df.columns)
+        length = len(df.index)
+        for i in range(length):
+            for j in range(i+1, length):
+                factor1, factor2 = df.index[i], df.index[j]
+                i_q = df.loc[factor2, factor1]
+                q1 = df.loc[factor1, factor1]
+                q2 = df.loc[factor2, factor2]
+
+                if (i_q <= q1 and i_q <= q2):
+                    outputRls = "Weaken, nonlinear"
+                if (i_q < max(q1, q2) and i_q > min(q1, q2)):
+                    outputRls = "Weaken, uni-"
+                if (i_q == (q1 + q2)):
+                    outputRls = "Independent"
+                if (i_q > max(q1, q2)):
+                    outputRls = "Enhance, bi-"
+                if (i_q > (q1 + q2)):
+                    outputRls = "Enhance, nonlinear"
+
+                out_df.loc[factor2, factor1] = outputRls
+        return out_df
+
+    def interaction_detector(self, relationship=False):
+        self.interaction_df = pd.DataFrame(index=self.factors, columns=self.factors, dtype="float32")
+        length = len(self.factors)
+        for i in range(0, length):
+            for j in range(0, i+1):
+                q, _, _ = self._cal_q(self.df, self.y, self.factors[i], self.factors[j])
+                self.interaction_df.loc[self.factors[i], self.factors[j]] = q
+
+        if relationship:
+            self.interaction_relationship_df = self._interaction_relationship(self.interaction_df)
+            return self.interaction_df, self.interaction_relationship_df
+        return self.interaction_df
+    
+    def ecological_detector(self):
+        self.ecological_df = pd.DataFrame(index=self.factors, columns=self.factors, dtype="float32")
+        length = len(self.factors)
+        for i in range(1, length):
+            ssw1, _, _ = self._cal_ssw(self.df, self.y, self.factors[i])
+            dfn = self.df[self.factors[i]].notna().sum()-1
+            for j in range(0, i):
+                ssw2, _, _ = self._cal_ssw(self.df, self.y, self.factors[j])
+                dfd = self.df[self.factors[j]].notna().sum()-1
+                fval = (dfn*(dfd-1)*ssw1)/(dfd*(dfn-1)*ssw2)
+                if fval<f.ppf(0.05, dfn, dfn):
+                    self.ecological_df.loc[self.factors[i], self.factors[j]] = 'Y'
+                else:
+                    self.ecological_df.loc[self.factors[i], self.factors[j]] = 'N'
+        return self.ecological_df
+
+    def plot(self, tick_fontsize=10, value_fontsize=10, colorbar_fontsize=10, show=True):
+        if isinstance(self.interaction_df, type(None)):
+            self.interaction_detector()
+        if isinstance(self.ecological_df, type(None)):
+            self.ecological_detector()
+
+        fig, ax = plt.subplots(constrained_layout=True)
+
+        im = ax.imshow(self.interaction_df.values, cmap="YlGnBu", vmin=0, vmax=1)
+        _plot_value(ax, self.interaction_df, self.ecological_df, value_fontsize=value_fontsize)
+
+        ax.set_xticks(np.arange(len(self.factors)))
+        ax.set_yticks(np.arange(len(self.factors)))
+        ax.spines['top'].set_visible(False)
+        ax.spines['right'].set_visible(False)
+
+        ax.set_xticklabels(self.factors, fontsize=tick_fontsize)
+        ax.set_yticklabels(self.factors, rotation=45, fontsize=tick_fontsize)
+        ax.tick_params(axis='y', pad=0.1)
+
+        colorbar = fig.colorbar(im, ax=ax, shrink=0.9, pad=0.01, aspect=25, extend="both")
+        colorbar.ax.tick_params(labelsize=colorbar_fontsize)
+
+        if show:
+            plt.show()
+            return ax
+        else:
+            return ax