spacr 0.3.61__tar.gz → 0.3.64__tar.gz

{spacr-0.3.61/spacr.egg-info → spacr-0.3.64}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: spacr
-Version: 0.3.61
+Version: 0.3.64
 Summary: Spatial phenotype analysis of crisp screens (SpaCr)
 Home-page: https://github.com/EinarOlafsson/spacr
 Author: Einar Birnir Olafsson

spacr-0.3.64/notebooks/old/.ipynb_checkpoints/Untitled-checkpoint.ipynb

@@ -0,0 +1,337 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "0a6d929d",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from spacr.plot import graph_importance\n",
+    "%matplotlib inline\n",
+    "\n",
+    "csvs = ['/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240810_132824/plate1/results/feature_importance_compartment.csv',\n",
+    "       '/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240824_072829/plate2/results/feature_importance_compartment.csv',\n",
+    "       '/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240825_094106/plate3/results/feature_importance_compartment.csv',\n",
+    "       '/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240826_140251/plate4/results/feature_importance_compartment.csv']\n",
+    "\n",
+    "settings = {'csvs':csvs,\n",
+    "           'grouping_column':'compartment',\n",
+    "           'data_column':'compartment_importance_sum',\n",
+    "           'graph_type':'jitter_bar',\n",
+    "            'save':False}\n",
+    "\n",
+    "graph_importance(settings)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "858a77a9",
+   "metadata": {
+    "scrolled": false
+   },
+   "outputs": [],
+   "source": [
+    "import pandas as pd\n",
+    "from spacr.plot import spacrGraph\n",
+    "import matplotlib.pyplot as plt\n",
+    "%matplotlib inline\n",
+    "\n",
+    "srcs = ['/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240810_132824/plate1/results/feature_importance_compartment.csv',\n",
+    "       '/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240824_072829/plate2/results/feature_importance_compartment.csv',\n",
+    "       '/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240825_094106/plate3/results/feature_importance_compartment.csv',\n",
+    "       '/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240826_140251/plate4/results/feature_importance_compartment.csv']\n",
+    "\n",
+    "dfs = []\n",
+    "for path in srcs:\n",
+    "    dft = pd.read_csv(path)\n",
+    "    dfs.append(dft)\n",
+    "\n",
+    "df = pd.concat(dfs)\n",
+    "display(df)\n",
+    "\n",
+    "spacr_graph = spacrGraph(\n",
+    "    df=df,                                     \n",
+    "    grouping_column='compartment',                         \n",
+    "    data_column='compartment_importance_sum',         \n",
+    "    graph_type='jitter_bar',          \n",
+    "    graph_name='compartment',          \n",
+    "    summary_func='mean',                         \n",
+    "    colors=None,                                \n",
+    "    output_dir='/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240810_132824/plate1/results',                              \n",
+    "    save=True,                       \n",
+    "    y_lim=None,                     \n",
+    "    error_bar_type='std',                       \n",
+    "    representation='object',\n",
+    "    theme='pastel',                    \n",
+    ")\n",
+    "\n",
+    "\n",
+    "# Create the plot\n",
+    "spacr_graph.create_plot()\n",
+    "\n",
+    "# Get the figure object if needed\n",
+    "fig = spacr_graph.get_figure()\n",
+    "plt.show()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "0f756080",
+   "metadata": {
+    "scrolled": false
+   },
+   "outputs": [],
+   "source": [
+    "from spacr.submodules import analyze_endodyogeny\n",
+    "%matplotlib inline\n",
+    "\n",
+    "settings = {\n",
+    "    'src': '/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240810_132824/plate1/',\n",
+    "    'tables': ['cell', 'nucleus', 'pathogen', 'cytoplasm'],\n",
+    "    'cell_types': ['Hela'],\n",
+    "    'cell_plate_metadata': None,\n",
+    "    'pathogen_types': ['nc', 'pc'],\n",
+    "    'pathogen_plate_metadata': [['c1'], ['c2']],\n",
+    "    'treatments': None,\n",
+    "    'treatment_plate_metadata': None,\n",
+    "    'min_area_bin': 500,\n",
+    "    'group_column': 'pathogen',\n",
+    "    'compartment': 'pathogen',\n",
+    "    'pathogen_limit': 1,\n",
+    "    'nuclei_limit': 1000,\n",
+    "    'level': 'well',\n",
+    "    'um_per_px':0.1,\n",
+    "    'max_bins':5,\n",
+    "    'save':False,\n",
+    "    'verbose': True}\n",
+    "\n",
+    "# Run the analysis\n",
+    "chi2, p = analyze_endodyogeny(settings)\n",
+    "print(f\"Chi-squared statistic: {chi2}, p-value: {p}\")\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 13,
+   "id": "0d73aa58",
+   "metadata": {
+    "scrolled": false
+   },
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "cell: 60816\n",
+      "nucleus: 80889\n",
+      "pathogen: 111317\n",
+      "cytoplasm: 60816\n",
+      "png_list: 60816\n",
+      "cells: 60816, cells grouped: 60816\n",
+      "cytoplasms: 60816, cytoplasms grouped: 60816\n",
+      "nucleus: 80889, nucleus grouped: 60816\n",
+      "pathogens: 37541, pathogens grouped: 37541\n",
+      "png_list: 60816, png_list grouped: 60816\n",
+      "Added png_list columns: Index(['test', 'test_2'], dtype='object'), Index(['png_path'], dtype='object')\n",
+      "Generated dataframe with: 794 columns and 37541 rows\n",
+      "Chi-squared test statistic (raw data): 0.0000\n",
+      "p-value (raw data): 1.0000e+00\n",
+      "Chi-squared results saved to /nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240810_132824/plate1/results/class_chi_squared_results.csv\n",
+      "Annotated data saved to /nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240810_132824/plate1/results/class_chi_squared_data.csv\n"
+     ]
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAABLMAAALACAYAAABo5DA3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABj3klEQVR4nO3dd5xU1f0//teCsDSBKE0UBXtHI18R0cREFEuwxFiwoNgSu6JGjQU7Ro2axEI0sSUSjT3Ggog1lhh7NFYsGCMgIqAgIHB/f/hjPqws1ZXlmufz8djHgzlz7p33nbkzw772nHOriqIoAgAAAAAl0KC+CwAAAACABSXMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAJYQJ07d85+++1X32XM05tvvpmtt946rVq1SlVVVe64447F9thleH7qwn777ZfOnTsv9sds0aLFYn3Mb6unn346jRs3znvvvVffpcA3pqqqKqeffnrl9rXXXpuqqqq8++67C7yPjz/+OM2bN88999xT9wUCwNckzAIWyaz/GM/6adKkSVZfffUcfvjhGT16dH2Xt8ieeOKJnH766Rk/fnx9l7JI9t133/zrX//KOeeckz/+8Y/p1q3bHH0uuuiiVFVV5YEHHpjrfq666qpUVVXlr3/96zdZ7hJj9nO5qqoqzZs3z9prr52zzz47kydPru/yFquZM2fm+uuvz1ZbbZU2bdqkUaNGadeuXbbeeutceeWVmTp1an2X+LWcfPLJ6du3b1ZaaaVK2xZbbJGqqqqsttpqtW4zbNiwyrlxyy23LK5S691HH32Uo446KmuuuWaaNm2adu3aZeONN84JJ5yQzz77rNJvv/32q/H+adGiRVZeeeX85Cc/ya233pqZM2fW41GwqJZddtkceOCBOfXUU+u7FACYw1L1XQBQbmeeeWa6dOmSKVOm5O9//3uuuOKK3HPPPXn55ZfTrFmz+i5voT3xxBM544wzst9++6V169Y17nv99dfToMGS+zeAzz//PE8++WROPvnkHH744XPtt8cee+T444/PkCFD0qtXr1r7DBkyJMsuu2y23Xbbb6rcJc5WW22Vfv36JUk+++yzPPbYYzn11FPz4osv5uabb670u+qqq761v5x//vnn2XnnnTN06NBsuummOe6449K+ffuMGzcujzzySA499ND84x//yB/+8If6LnWRvPDCC3nggQfyxBNPzHFfkyZN8tZbb+Xpp5/OxhtvXOO+G264IU2aNMmUKVMWV6n1bty4cenWrVsmTpyY/fffP2uuuWY+/vjjvPTSS7niiityyCGH1BgtWF1dnd///vdJvjyP3nvvvdx11135yU9+ki222CJ33nlnWrZsWV+HwyL62c9+lt/85jd58MEH88Mf/rC+ywGACmEW8LVsu+22ldE/Bx54YJZddtlcdNFFufPOO9O3b99at5k0aVKaN2++OMucrwWpqbq6ejFVs2g++uijJJkjhPuqjh075gc/+EFuu+22XHHFFXMc1wcffJBHH300Bx98cBo1avRNlbvEWX311bP33ntXbv/sZz/LtGnTctttt2XKlClp0qRJknyrn5NjjjkmQ4cOzSWXXJKjjjqqxn3HHnts3nzzzQwbNmye+5g+fXpmzpyZxo0bf5OlLpJrrrkmK664YjbZZJM57ltllVUyffr0/PnPf64RZk2ZMiW33357tt9++9x6662Ls9w6sd9+++Xdd9/Nww8/vFDb/eEPf8jIkSPz+OOPZ9NNN61x38SJE+d4fZdaaqka758kOfvss3PeeeflpJNOykEHHZSbbrppkY6hriyJ3z1LurXWWivrrrturr32WmEWAEuUJXeIAVBKs/6z+8477yT5v7V+RowYke222y5LL7109tprryRf/mJx7LHHplOnTqmurs4aa6yRCy+8MEVR1NhnVVVVDj/88Nxwww1ZY4010qRJk2y00UZ59NFH53j8559/Pttuu21atmyZFi1aZMstt8xTTz1Vo8+sKZKzRpq0a9cuK6ywQk4//fQcf/zxSZIuXbpUpszMWmOktjWh3n777ey6665ZZpll0qxZs2yyySa5++67a/R5+OGHU1VVlb/85S8555xzssIKK6RJkybZcsst89Zbby3Q8zq/4zr99NMr06aOP/74VFVVzXNdp7333jsTJkyYo9YkufHGGzNz5syFfp2+6vTTT09VVdUc7bWt3dK5c+f86Ec/ysMPP5xu3bqladOmWW+99Sq/gN92221Zb731Kq/9888/P8d+X3vttfzkJz/JMssskyZNmqRbt25fe5pkhw4dUlVVlaWW+r+//Xx1zax33303VVVVufDCC3PllVdmlVVWSXV1df7f//t/+ec//1ljf6NGjUr//v2zwgorpLq6Osstt1x23HHHBV7H5u23307v3r3TvHnzdOzYMWeeeWbldSiKIp07d86OO+44x3ZTpkxJq1at8tOf/nSu+37//ffz+9//Pttss80cQdYsq622Wg499NBaj/2SSy6pHPu///3vJMmDDz6YzTffPM2bN0/r1q2z44475tVXX62xz7mtQVbb+bMwnwW1ueOOO/LDH/6w1vMySfr27Zubbrqpxsi7u+66K5MnT85uu+1W6zYffPBB9t9//7Rv3z7V1dVZZ511cvXVV9foM23atJx22mnZaKON0qpVqzRv3jybb755HnrooRr9FuZc+qaNGDEiDRs2rDX4a9myZSXcnZ8TTzwxW2+9dW6++ea88cYb8+z70ksvZb/99svKK6+cJk2apEOHDtl///3z8ccfz9H3gw8+yAEHHJCOHTumuro6Xbp0ySGHHJJp06Ylmfvn/CyXX3551llnnVRXV6djx4457LDD5phe/uabb2aXXXZJhw4d0qRJk6ywwgrZY489MmHChEqfYcOGZbPNNkvr1q3TokWLrLHGGvnFL34xz+P88Y9/nO9+97s12vr06TPH1O5//OMfqaqqyr333ltpGz9+fI4++ujK5/Gqq66aX/7yl4s0WvSZZ55J796906ZNmzRt2jRdunTJ/vvvP0e/rbbaKnfdddd8P/MBYHEyMguoUyNGjEjy5Vobs0yfPj29e/fOZpttlgsvvDDNmjVLURTZYYcd8tBDD+WAAw7IBhtskKFDh+b444/PBx98kIsvvrjGfh955JHcdNNNOfLII1NdXZ3LL78822yzTZ5++umsu+66SZJXXnklm2++eVq2bJmf//znadSoUX73u99liy22yCOPPJLu3bvX2Oehhx6atm3b5rTTTsukSZOy7bbb5o033sif//znXHzxxWnTpk2SpG3btrUe6+jRo7Pppptm8uTJOfLII7Psssvmuuuuyw477JBbbrklO++8c43+5513Xho0aJDjjjsuEyZMyPnnn5+99tor//jHP+b5nC7Icf34xz9O69atc8wxx6Rv377Zbrvt5rlg+I9//OMccsghGTJkSH784x/XuG/IkCFZaaWV0rNnz4V+nb6Ot956K3vuuWd++tOfZu+9986FF16YPn36ZPDgwfnFL35RCVEGDRqU3Xbbrca0z1deeSU9e/bM8ssvnxNPPDHNmzfPX/7yl+y000659dZb53gtajNlypSMHTs2yZcB3uOPP57rrrsue+65Z40wa26GDBmSTz/9ND/96U9TVVWV888/Pz/+8Y/z9ttvV0Zz7bLLLnnllVdyxBFHpHPnzhkzZkyGDRuWkSNHzndR+RkzZmSbbbbJJptskvPPPz/33XdfBg4cmOnTp+fMM89MVVVV9t5775x//vkZN25clllmmcq2d911VyZOnDjHyJnZ3XvvvZkxY8Y8+8zNNddckylTpuTggw9OdXV1lllmmTzwwAPZdttts/LKK+f000/P559/nt/+9rfp2bNnnnvuuUVeRH9BPgtq88EHH2TkyJFzhAiz23PPPXP66afn4YcfrgTzQ4YMyZZbbpl27drN0X/06NHZZJNNKiFb27Ztc++99+aAAw7IxIkTc/TRRyf5ciTT73//+/Tt2zcHHXRQPv300/zhD39I79698/TTT2eDDTaosd8FOZdqM3PmzIwbN65G29SpU/PFF19Uzu1ZWrVqNc99rbTSSpkxY0b++Mc/Zt99951rvwWxzz775P7778+wYcOy+uqrz7XfsGHD8vbbb6d///7p0KFDXnnllVx55ZV55ZVX8tRTT1VCyP/+97/ZeOONM378+Bx88MFZc80188EHH+SWW27J5MmTa4wa++rnfPJlUHrGGWekV69eOeSQQ/L666/niiuuyD//+c88/vjjadSoUaZNm5bevXtn6tSpOeKII9KhQ4d88MEH+dvf/pbx48enVatWeeWVV/KjH/0o66+/fs4888xUV1fnrbfeyuOPPz7P52PzzTfPnXfemYkTJ6Zly5YpiiKPP/54GjRokMceeyw77LBDkuSxxx5LgwYN0rNnzyTJ5MmT8/3vfz8ffPBBfvrTn2bFFVfME088kZNOOikffvhhLrnkkgV+TcaMGZOtt946bdu2zYknnpjWrVvn3XffzW233TZH34022igXX3xxXnnllXm+xwBgsSoAFsE111xTJCkeeOCB4qOPPiref//94sYbbyyWXXbZomnTpsV//vOfoiiKYt999y2SFCeeeGKN7e+4444iSXH22WfXaP/JT35SVFVVFW+99ValLUmRpHjmmWcqbe+9917RpEmTYuedd6607bTTTkXjxo2LESNGVNr++9//FksvvXTxve99b47aN9tss2L69Ok1Hv+CCy4okhTvvPPOHMe80korFfvuu2/l9tFHH10kKR577LFK26efflp06dKl6Ny5czFjxoyiKIrioYceKpIUa621VjF16tRK31//+tdFkuJf//rXnE/wbBb0uN55550iSXHBBRfMc3+z7LrrrkWTJk2KCRMmVNpee+21Iklx0kknFUWxcK/TV5+fgQMHFrV9zcx6/md/jldaaaUiSfHEE09U2oYOHVokKZo2bVq89957lfbf/e53RZLioYceqrRtueWWxXrrrVdMmTKl0jZz5sxi0003LVZbbbX5PhezzrGv/uy000419lkUX57TK620UuX2rOd92WWXLcaNG1dpv/POO4skxV133VUURVF88sknC/X6fPUxkxRHHHFEjePbfvvti8aNGxcfffRRURRF8frrrxdJiiuuuKLG9jvssEPRuXPnYubMmXN9jGOOOaZIUrzwwgs12qdOnVp89NFHlZ+xY8fOcewtW7YsxowZU2O7DTbYoGjXrl3x8ccfV9pefPHFokGDBkW/fv1qHNvsz+cstZ0/C/pZUJsHHnigxusxu+9///vFOuusUxRFUXTr1q044IADiqL48jVr3Lhxcd1111XexzfffHNluwMOOKBYbrnlajwnRVEUe+yxR9GqVati8uTJRVEUxfTp02u892ftu3379sX+++9faVvQc2luZm2/ID+zv39qM2rUqKJt27ZFkmLNNdcsfvaznxVDhgwpxo8fP0fffffdt2jevPlc9/X8888XSYpjjjlmno856/ma3Z///OciSfHoo49W2vr161c0aNCg+Oc//zlH/1nn+Nw+58eMGVM0bty42HrrrSuf0UVRFJdeemmRpLj66qtr1Dz76/1VF198cZGk8v5bUP/85z+LJMU999xTFEVRvPTSS0WSYtdddy26d+9e6bfDDjsUG264YeX2WWedVTRv3rx44403auzvxBNPLBo2bFiMHDmy0pakGDhwYOX2Vz93b7/99iJJrc/hVz3xxBNFkuKmm25aqOMEgG+SaYbA19KrV6+0bds2nTp1yh577JEWLVrk9ttvz/LLL1+j3yGHHFLj9j333JOGDRvmyCOPrNF+7LHHpiiKGtMqkqRHjx7ZaKONKrdXXHHF7Ljjjhk6dGhmzJiRGTNm5P77789OO+2UlVdeudJvueWWy5577pm///3vmThxYo19HnTQQWnYsOEiH/s999yTjTfeOJtttlmlrUWLFjn44IPz7rvvVqZazdK/f/8aIwY233zzJF9OHZubRTmuBbX33ntnypQpNf4SP2TIkCSpTDFc2Nfp61h77bXTo0ePyu1ZI+l++MMfZsUVV5yjfdbzNm7cuDz44IPZbbfd8umnn2bs2LEZO3ZsPv744/Tu3TtvvvlmPvjgg/k+/o477phhw4Zl2LBhufPOO3PSSSflvvvuy5577rlA02t23333fOc736nc/urr27Rp0zRu3DgPP/xwPvnkk/nurzazL+w/azTQtGnTKlemXH311dO9e/fccMMNlX7jxo3Lvffem7322muu0+uSVM6jr47ou+eee9K2bdvKz+xXAZxll112qTGC8cMPP8wLL7yQ/fbbr8YIsfXXXz9bbbVV7rnnnoU88v8zv8+CuZk1VW3216g2e+65Z2677bZMmzYtt9xySxo2bFjryL6iKHLrrbemT58+KYqict6NHTs2vXv3zoQJE/Lcc88lSRo2bFh5788aPTV9+vR069at0md28zuX5qZDhw6Vc3jWz9Zbb531119/jvauXbvOc1/t27fPiy++mJ/97Gf55JNPMnjw4Oy5555p165dzjrrrIWacjbrnPr000/n2a9p06aVf88aKTlrmuOs52nmzJm544470qdPn1qv1vrVc/yrn/MPPPBApk2blqOPPrrGBT0OOuigtGzZsjL1ulWrVkmSoUOHzvWKprPWJ7zzzjsXaprfhhtumBYtWlSmxz722GNZYYUV0q9fvzz33HOZPHlyiqLI3//+98prnyQ333xzNt9883znO9+pcb716tUrM2bMWODptrPX/re//S1ffPHFPPvOOhe/OroPAOqTMAv4Wi677LIMGzYsDz30UP79739X1vSZ3VJLLVVjrZIkee+999KxY8csvfTSNdrXWmutyv2zW2211eZ47NVXXz2TJ0/ORx99lI8++iiTJ0/OGmusMUe/tdZaKzNnzsz7779fo71Lly4LfqC1eO+99+b6eLPun93sgUzyf78gzCvYWJTjWlDbbrttlllmmUqAlSR//vOf07Vr16yzzjqVY1iY1+nr+OrzM+uXyU6dOtXaPut5e+utt1IURU499dQaoUvbtm0zcODAJF9OqZmfFVZYIb169UqvXr2yww475Nxzz83ZZ5+d2267LX/7298Wuv6vvr7V1dX55S9/mXvvvTft27fP9773vZx//vkZNWrUfPedJA0aNKgRaCapTNmafc2tfv365fHHH6+8NjfffHO++OKL7LPPPvPc/6zX+LPPPqvR3rNnzxrBSG2++l6a9dhzO2/Hjh1bmfK1sOb3WTA/8wthZq2JdO+99+aGG27Ij370oznO/+TL9+b48eNz5ZVXznHe9e/fP0nN8+66667L+uuvnyZNmmTZZZdN27Ztc/fdd9dYf2mWRfmsSL68IuOsc3jWz3LLLZfvfOc7c7TPL9RLvgzNr7jiinz44Yd5/fXX85vf/KYyZW9hrmg565yq7Xmc3bhx43LUUUelffv2adq0adq2bVs5t2Y9Tx999FEmTpy4wNPdFvTcbNy4cVZeeeXK/V26dMmAAQPy+9//Pm3atEnv3r1z2WWX1Xi9dt999/Ts2TMHHnhg2rdvnz322CN/+ctf5htsNWzYMD169Mhjjz2W5Mswa/PNN89mm22WGTNm5Kmnnsq///3vjBs3rkaY9eabb+a+++6b43ybdVXaBfmcm+X73/9+dtlll5xxxhlp06ZNdtxxx1xzzTWZOnXqHH1nvWfmFYYDwOJmzSzga9l4441r/ev47Kqrq2v8BXxJMfsogMVhbqPAFmaEQ11q1KhRdtttt1x11VUZPXp0Ro4cmTfffDPnn39+nex/br/4zG30zNyen/k9b7N+cTzuuOPmCFJnWXXVVedZ69xsueWWSZJHH300ffr0mWffBXl9jz766PTp0yd33HFHhg4dmlNPPTWDBg3Kgw8+mA033HCRavyqPfbYI8ccc0xuuOGG/OIXv8if/vSndOvWrdZgaXZrrrlmkuTll1+uMWpn9l+W//SnP9W67dd5Ly3sebKoZq3jN79AaLnllssWW2yRX/3qV3n88cfnegXDWefd3nvvPdc1pdZff/0kXz5v++23X3baaaccf/zxadeuXRo2bJhBgwZV1hmc3ZL2WVFVVZXVV189q6++erbffvusttpqueGGG3LggQcu0PYvv/xykvm/D3fbbbc88cQTOf7447PBBhukRYsWmTlzZrbZZptFWuA8+Xrn5q9+9avst99+ufPOO3P//ffnyCOPzKBBg/LUU09lhRVWSNOmTfPoo4/moYceyt1335377rsvN910U374wx/m/vvvn+fI38022yznnHNOpkyZksceeywnn3xyWrdunXXXXTePPfZY2rdvnyQ1wqyZM2dmq622ys9//vNa9zmv9ci+qqqqKrfcckueeuqp3HXXXRk6dGj233///OpXv8pTTz1VY4TmrPfMrHUkAWBJIMwC6sVKK62UBx54IJ9++mmNv9a/9tprlftn9+abb86xjzfeeCPNmjWrTG9q1qxZXn/99Tn6vfbaa2nQoMEcI3xqszB/eV5ppZXm+niz7v+62rZtWyfHNTd77bVXBg8enJtuuinvvPNOqqqq0rdv38r9C/s6zW7WyI/x48dXprQkdTuaK0lltFKjRo0qoUtdmT59epI5Ryt9HausskqOPfbYHHvssXnzzTezwQYb5Fe/+tVcg6JZZs6cmbfffrvGL6yzrg43+2LqyyyzTLbffvvccMMN2WuvvfL4448v0MLQ2267bRo2bFjZ7uuYdV7M7bxt06ZNmjdvnuTL8+SrV5FL5n6eLMhnQW1mhXWzrrQ6L3vuuWcOPPDAtG7dOtttt12tfdq2bZull146M2bMmO95d8stt2TllVfObbfdVuMzZtbIwTJZeeWV853vfCcffvjhAm/zxz/+MVVVVdlqq63m2ueTTz7J8OHDc8YZZ+S0006rtH/19W7btm1atmxZCcgW1uzn5uwjHadNm5Z33nlnjtdyvfXWy3rrrZdTTjklTzzxRHr27JnBgwfn7LPPTvLliMktt9wyW265ZS666KKce+65Ofnkk/PQQw/N87zYfPPNM23atPz5z3/OBx98UAmtvve971XCrNVXX70SaiVffnZ89tlndfo5t8kmm2STTTbJOeeckyFDhmSvvfbKjTfeWCOonPWemTUiFwCWBEveUAngf8J2222XGTNm5NJLL63RfvHFF6eqqirbbrttjfYnn3yyxtoy77//fu68885svfXWadiwYRo2bJitt946d955Z40pV6NHj86QIUOy2WabpWXLlvOta9Yv2LX9cl3bMTz99NN58sknK22TJk3KlVdemc6dO2fttdee7z7mp66Oa2569uyZzp07509/+lNuuummfP/7368xJXRhX6fZrbLKKklSYx2XSZMm5brrrlvkemvTrl27bLHFFvnd735X6y/YCzL1bG7uuuuuJJnv+kILYvLkyZkyZUqNtlVWWSVLL710rVN7ajP761AURS699NI0atSoMoJsln322Sf//ve/c/zxx6dhw4bZY4895rvvFVdcMfvvv3/uvffeOV7v2R9zQSy33HLZYIMNct1119V4L7388su5//77awREq6yySiZMmJCXXnqp0vbhhx/m9ttvr3Xf8/ssmJvll18+nTp1yjPPPDPf+n/yk59k4MCBufzyy2uscze7hg0bZpdddsmtt95aa7Ay+3k3q67Zn79//OMfNT47vinXXnttHn744YXe7h//+EetU0GffvrpfPzxx/Md6TfLeeedl/vvvz+77757rVNEZ6ntOUoyRxDboEGD7LTTTrnrrrtqfS3nd4726tUrjRs3zm9+85saff/whz9kwoQJ2X777ZN8uYbcrDB7lvXWWy8NGjSovF+/euXIJJUrU87vPd29e/c0atQov/zlL7PMMstUpnZvvvnmeeqpp/LII4/UGJWVfDly7cknn8zQoUPn2N/48ePnqHdePvnkkzmeq7nV/uyzz6ZVq1aVGgFgSWBkFlAv+vTpkx/84Ac5+eST8+6776Zr1665//77c+edd+boo4+uBCGzrLvuuundu3eOPPLIVFdX5/LLL0+SnHHGGZU+Z599doYNG5bNNtsshx56aJZaaqn87ne/y9SpUxd46tyshaVPPvnk7LHHHmnUqFH69OlTCblmd+KJJ+bPf/5ztt122xx55JFZZpllct111+Wdd97JrbfeWmdTK+viuOamqqoqe+65Z84999wkyZlnnlnj/oV9nWa39dZbZ8UVV8wBBxxQCVWuvvrqtG3bNiNHjvxadX/VZZddls022yzrrbdeDjrooKy88soZPXp0nnzyyfznP//Jiy++ON99vPHGG5XRUZMnT85TTz2V6667Lquuuup815taEG+88Ua23HLL7Lbbbll77bWz1FJL5fbbb8/o0aMXKGxq0qRJ7rvvvuy7777p3r177r333tx99935xS9+MceIpO233z7LLrtsbr755my77bZp167dAtV4ySWX5J133skRRxyRG2+8MX369Em7du0yduzYPP7447nrrrsWOMS44IILsu2226ZHjx454IAD8vnnn+e3v/1tWrVqldNPP73Sb4899sgJJ5yQnXfeOUceeWQmT56cK664Iquvvnqti6MvyGfB3Oy44465/fbbUxTFPEdhfrXGuTnvvPPy0EMPpXv37jnooIOy9tprZ9y4cXnuuefywAMPVMKOH/3oR7ntttuy8847Z/vtt88777yTwYMHZ+21167TUX+TJk2aawj4VVtttVWNUT9f9cc//jE33HBDdt5552y00UZp3LhxXn311Vx99dVp0qRJfvGLX9ToP3369Mr7Z8qUKXnvvffy17/+NS+99FJ+8IMf5Morr5xnPS1btqysI/fFF19k+eWXz/3331/rSLpzzz03999/f77//e/n4IMPzlprrZUPP/wwN998c/7+97/XGAn6VW3bts1JJ52UM844I9tss0122GGHvP7667n88svz//7f/8vee++dJHnwwQdz+OGHZ9ddd83qq6+e6dOn549//GMlxEy+/Lx89NFHs/3222ellVbKmDFjcvnll2eFFVaocWGQ2jRr1iwbbbRRnnrqqfTp06dyPn7ve9/LpEmTMmnSpDnCrOOPPz5//etf86Mf/Sj77bdfNtpoo0yaNCn/+te/csstt+Tdd99d4KmA1113XS6//PLsvPPOWWWVVfLpp5/mqquuSsuWLecYjThs2LAaNQLAEmFxXjoR+PaYdZnv+V3We16XbP/000+LY445pujYsWPRqFGjYrXVVisuuOCCyqXVZ0lSHHbYYcWf/vSnYrXVViuqq6uLDTfcsNZLyz/33HNF7969ixYtWhTNmjUrfvCDHxRPPPHEQtV+1llnFcsvv3zRoEGDGpcyX2mllYp99923Rt8RI0YUP/nJT4rWrVsXTZo0KTbeeOPib3/7W40+Dz30UK2XeH/nnXeKJMU111xTax0Le1yz9nfBBRfMd3+ze+WVV4okRXV1dfHJJ5/Mcf+Cvk61PT/PPvts0b1796Jx48bFiiuuWFx00UVzXCJ+1rbbb7/9HI8967VfkOMcMWJE0a9fv6JDhw5Fo0aNiuWXX7740Y9+VNxyyy3zfQ6S1Php2LBhscIKKxQHH3xwMXr06Bp9991332KllVaabz2z9jtw4MCiKIpi7NixxWGHHVasueaaRfPmzYtWrVoV3bt3L/7yl7/Mt75Z76MRI0YUW2+9ddGsWbOiffv2xcCBA4sZM2bUus2hhx5aJCmGDBky3/3Pbvr06cU111xT/PCHPyyWWWaZYqmlliratGlTbLnllsXgwYOLzz//fIGOvSiK4oEHHih69uxZNG3atGjZsmXRp0+f4t///vcc/e6///5i3XXXLRo3blysscYaxZ/+9Kdi4MCBxVf/m7IwnwW1ee6554okxWOPPVaj/fvf/36xzjrrzHPbub2PR48eXRx22GFFp06dikaNGhUdOnQottxyy+LKK6+s9Jk5c2Zx7rnnFiuttFKl5r/97W+LfC7NzaztF+Rnfs/ZSy+9VBx//PHFd7/73cp5sNxyyxW77rpr8dxzz9Xou++++9bYd7NmzYrOnTsXu+yyS3HLLbfM9Rz9qv/85z/FzjvvXLRu3bpo1apVseuuuxb//e9/az329957r+jXr1/Rtm3borq6ulh55ZWLww47rJg6dWpRFPP/nL/00kuLNddcs2jUqFHRvn374pBDDqnx+ff2228X+++/f7HKKqsUTZo0KZZZZpniBz/4QfHAAw9U+gwfPrzYcccdi44dOxaNGzcuOnbsWPTt27d44403Fuh4jz/++CJJ8ctf/rJG+6qrrlokKUaMGDHHNp9++mlx0kknFauuumrRuHHjok2bNsWmm25aXHjhhcW0adMq/b76nH31c/e5554r+vbtW6y44opFdXV10a5du+JHP/pR8cwzz9R4vFdffbVIUuO4AWBJUFUU9bSaKMACqqqqymGHHTbXqU/AnI455pj84Q9/yKhRo9KsWbP6LqdO1MVnwZZbbpmOHTvmj3/8Yx1WBt9ORx99dB599NE8++yzRmYBsESxZhYAfMtMmTIlf/rTn7LLLrt8a4KsunLuuefmpptuqvMLEcC3zccff5zf//73OfvsswVZACxxrJkFAN8SY8aMyQMPPJBbbrklH3/8cY466qj6LmmJ071790ybNq2+y4Al3rLLLluna7oBQF0SZgHAt8S///3v7LXXXmnXrl1+85vfVK5OBgAA3yb1umbWo48+mgsuuCDPPvts5TLcO+200zy3efjhhzNgwIC88sor6dSpU0455ZTst99+i6VeAAAAAOpXva6ZNWnSpHTt2jWXXXbZAvV/5513sv322+cHP/hBXnjhhRx99NE58MADM3To0G+4UgAAAACWBEvM1QyrqqrmOzLrhBNOyN13352XX3650rbHHntk/Pjxue+++xZDlQAAAADUp1KtmfXkk0+mV69eNdp69+6do48+eq7bTJ06NVOnTq3cnjlzZsaNG5dll13WlVkAAAD4xhVFkU8//TQdO3ZMgwb1OkEKvhVKFWaNGjUq7du3r9HWvn37TJw4MZ9//nmaNm06xzaDBg3KGWecsbhKBAAAgFq9//77WWGFFeq7DCi9UoVZi+Kkk07KgAEDKrcnTJiQFVdcMe+//35atmxZj5VRl3Zs1a++S2AhfZrxmZRPF3q75lk6S6d13RfEN+bOCdfXdwnAPPgOLR/fof87fId+e0ycODGdOnXK0ksvXd+lwLdCqcKsDh06ZPTo0TXaRo8enZYtW9Y6KitJqqurU11dPUd7y5YthVnfIktVNarvElhI30nbfCdt67sMFgOftbBk8x1aPr5D/3f4Dv32sdQN1I1STdbt0aNHhg8fXqNt2LBh6dGjRz1VBAAAAMDiVK9h1meffZYXXnghL7zwQpLknXfeyQsvvJCRI0cm+XKKYL9+/zf0/Wc/+1nefvvt/PznP89rr72Wyy+/PH/5y19yzDHH1Ef5AAAAACxm9RpmPfPMM9lwww2z4YYbJkkGDBiQDTfcMKeddlqS5MMPP6wEW0nSpUuX3H333Rk2bFi6du2aX/3qV/n973+f3r1710v9AAAAACxe9bpm1hZbbJGiKOZ6/7XXXlvrNs8///w3WBUAAACwMGbMmJEvvviivsugxBo1apSGDRsuUN9SLQAPAAAALDmKosioUaMyfvz4+i6Fb4HWrVunQ4cO871YgjALAAAAWCSzgqx27dqlWbNmrtjIIimKIpMnT86YMWOSJMstt9w8+wuzAAAAgIU2Y8aMSpC17LLL1nc5lFzTpk2TJGPGjEm7du3mOeWwXheABwAAAMpp1hpZzZo1q+dK+LaYdS7Nb/01YRYAAACwyEwtpK4s6LkkzAIAAACgNIRZAAAAwBKnc+fOueSSS+q7jG+1ffbZJ+eee+43/jhjx45Nu3bt8p///KdO9ifMAgAAAOpMnz59ss0229R632OPPZaqqqq89NJLi7mqurXFFlukqqoqVVVVadKkSVZfffUMGjQoRVHUd2kL7MUXX8w999yTI488ss72OWnSpOyzzz7Zdttt06tXr9xyyy1JkjZt2qRfv34ZOHBgnTyOMAsAAACoMwcccECGDRtW6yica665Jt26dcv6669fD5XVrYMOOigffvhhXn/99Zx00kk57bTTMnjw4Poua4H99re/za677poWLVrU2T6bN2+es88+O1OnTs1zzz2XO++8s3Jf//79c8MNN2TcuHFf+3GEWQAAAECd+dGPfpS2bdvm2muvrdH+2Wef5eabb84BBxyQJLn11luzzjrrpLq6Op07d86vfvWrue7z3XffTVVVVV544YVK2/jx41NVVZWHH344SfLwww+nqqoqQ4cOzYYbbpimTZvmhz/8YcaMGZN77703a621Vlq2bJk999wzkydPruxn5syZGTRoULp06ZKmTZuma9eulRFF89KsWbN06NAhK620Uvr375/1118/w4YNq9w/YsSI7Ljjjmnfvn1atGiR//f//l8eeOCBGvvo3Llzzj333Oy///5Zeumls+KKK+bKK6+s0eeJJ57IBhtskCZNmqRbt26544475nguXn755Wy77bZp0aJF2rdvn3322Sdjx46da+0zZszILbfckj59+sxRz1lnnZW+ffumefPmWX755XPZZZdV7n/44YfTuHHjPPbYY5W2888/P+3atcvo0aOTJCuttFIefPDBDBkyJLvvvnul3zrrrJOOHTvm9ttvn+9zOz/CLAAAAKDOLLXUUunXr1+uvfbaGtPubr755syYMSN9+/bNs88+m9122y177LFH/vWvf+X000/PqaeeOkcAtihOP/30XHrppXniiSfy/vvvZ7fddssll1ySIUOG5O67787999+f3/72t5X+gwYNyvXXX5/BgwfnlVdeyTHHHJO99947jzzyyAI9XlEUeeyxx/Laa6+lcePGlfbPPvss2223XYYPH57nn38+22yzTfr06ZORI0fW2P5Xv/pVunXrlueffz6HHnpoDjnkkLz++utJkokTJ6ZPnz5Zb7318txzz+Wss87KCSecUGP78ePH54c//GE23HDDPPPMM7nvvvsyevTo7LbbbnOt+aWXXsqECRPSrVu3Oe674IIL0rVr1zz//PM58cQTc9RRR1VCui222CJHH3109tlnn0yYMCHPP/98Tj311Pz+979P+/btM2XKlMp+GjRokFNOOaXGvjfeeOMaQdiiWupr7wEAAABgNvvvv38uuOCCPPLII9liiy2SfDnFcJdddkmrVq1y0UUXZcstt8ypp56aJFl99dXz73//OxdccEH222+/r/XYZ599dnr27JnkyymPJ510UkaMGJGVV145SfKTn/wkDz30UE444YRMnTo15557bh544IH06NEjSbLyyivn73//e373u9/l+9///lwf5/LLL8/vf//7TJs2LV988UWaNGlSY/2prl27pmvXrpXbZ511Vm6//fb89a9/zeGHH15p32677XLooYcmSU444YRcfPHFeeihh7LGGmtkyJAhqaqqylVXXZUmTZpk7bXXzgcffJCDDjqosv2ll16aDTfcsMZC7ldffXU6deqUN954I6uvvvoctb/33ntp2LBh2rVrN8d9PXv2zIknnpjky9fl8ccfz8UXX5ytttqq8vwOGzYsBx98cF5++eXsu+++2WGHHZJ8GZIdddRRadSoUcaNGzfHGlkdO3bM888/P9fndEEJswAAAIA6teaaa2bTTTfN1VdfnS222CJvvfVWHnvssZx55plJkldffTU77rhjjW169uyZSy65JDNmzEjDhg0X+bFnX4+rffv2adasWSXImtX29NNPJ0neeuutTJ48uRLUzDJt2rRsuOGG83ycvfbaKyeffHI++eSTDBw4MJtuumk23XTTyv2fffZZTj/99Nx999358MMPM3369Hz++edzjMyavd6qqqp06NAhY8aMSZK8/vrrWX/99dOkSZNKn4033rjG9i+++GIeeuihWte+GjFiRK1h1ueff57q6upUVVXNcd+sUG/227NfVbJx48a54YYbsv7662ellVbKxRdfXKO2J598co59ztK0adMaUzwXlTALAAAAqHMHHHBAjjjiiFx22WW55pprssoqq8xzpNO8NGjw5SpJs09b/OKLL2rt26hRo8q/q6qqatye1TZz5swkXwZOSXL33Xdn+eWXr9Gvurp6njW1atUqq666apLkL3/5S1ZdddVssskm6dWrV5LkuOOOy7Bhw3LhhRdm1VVXTdOmTfOTn/wk06ZNm2u9X61vQXz22Wfp06dPfvnLX85x33LLLVfrNm3atMnkyZMzbdq0GlMjF9QTTzyRJBk3blzGjRuX5s2bL9B248aNS9u2bRf68b7KmlkAAABAndttt93SoEGDDBkyJNdff33233//ykigtdZaK48//niN/o8//nhWX331WkdlzQpAPvzww0rb7AugL6q111471dXVGTlyZFZdddUaP506dVrg/bRo0SJHHXVUjjvuuErg9vjjj2e//fbLzjvvnPXWWy8dOnTIu+++u1D1rbHGGvnXv/6VqVOnVtr++c9/1ujz3e9+N6+88ko6d+48xzHMLWTaYIMNkiT//ve/57jvqaeemuP2WmutVbk9YsSIHHPMMbnqqqvSvXv37Lvvvgscvr388svzHfG2IIRZAAAAQJ1r0aJFdt9995x00kn58MMPa6yFdeyxx2b48OE566yz8sYbb+S6667LpZdemuOOO67WfTVt2jSbbLJJzjvvvLz66qt55JFH5lhcfFEsvfTSOe6443LMMcfkuuuuy4gRI/Lcc8/lt7/9ba677rqF2tdPf/rTvPHGG7n11luTJKuttlpuu+22vPDCC3nxxRez5557LtSIqySVbQ4++OC8+uqrGTp0aC688MIkqQSDhx12WMaNG5e+ffvmn//8Z0aMGJGhQ4emf//+mTFjRq37bdu2bb773e/m73//+xz3Pf744zn//PPzxhtv5LLLLsvNN9+co446KsmXV0Hce++907t37/Tv3z/XXHNNXnrppXleiXKWyZMn59lnn83WW2+9UM9BbYRZAAAAwDfigAMOyCeffJLevXunY8eOlfbvfve7+ctf/pIbb7wx6667bk477bSceeaZ81z8/eqrr8706dOz0UYb5eijj87ZZ59dJzWeddZZOfXUUzNo0KCstdZa2WabbXL33XenS5cuC7WfZZZZJv369cvpp5+emTNn5qKLLsp3vvOdbLrppunTp0969+6d7373uwu1z5YtW+auu+7KCy+8kA022CAnn3xyTjvttCSprKPVsWPHPP7445kxY0a23nrrrLfeejn66KPTunXryvTM2hx44IG54YYb5mg/9thj88wzz2TDDTfM2WefnYsuuii9e/dOkpxzzjl577338rvf/S7Jl9MYr7zyypxyyil58cUX53ksd955Z1ZcccVsvvnmC/Uc1KaqmH3C6f+AiRMnplWrVpkwYUJatmxZ3+VQR7ZqsGt9lwDMxbCZN9d3CcA8+A6FJZfv0G+Pb+vvoVOmTMk777yTLl261FignG/WDTfckP79+2fChAlp2rTpIu/n888/zxprrJGbbrqpsuh7586dc/TRR+foo4+uo2r/zyabbJIjjzwye+6551z7LOg5ZQF4AAAAgCXU9ddfn5VXXjnLL798XnzxxZxwwgnZbbfdvlaQlXw5dfP666/P2LFj66jSuRs7dmx+/OMfp2/fvnWyP2EWAAAAwBJq1KhROe200zJq1Kgst9xy2XXXXXPOOefUyb632GKLOtnP/LRp0yY///nP62x/wiwAAACAJdTPf/7zOg2C5mVhr7ZYXywADwAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGkvVdwEAAAAAX9dWDXZdrI83bObNC9X/0UcfzQUXXJBnn302H374YW6//fbstNNOC/24RVFk4MCBueqqqzJ+/Pj07NkzV1xxRVZbbbV5bnfZZZflggsuyKhRo9K1a9f89re/zcYbb7zQj78kMDILAAAA4Bs2adKkdO3aNZdddtnX2s/555+f3/zmNxk8eHD+8Y9/pHnz5undu3emTJky121uuummDBgwIAMHDsxzzz2Xrl27pnfv3hkzZszXqqW+CLMAAAAAvmHbbrttzj777Oy8886LvI+iKHLJJZfklFNOyY477pj1118/119/ff773//mjjvumOt2F110UQ466KD0798/a6+9dgYPHpxmzZrl6quvXuRa6pMwCwAAAGAJcPrpp6dz585zvf+dd97JqFGj0qtXr0pbq1at0r179zz55JO1bjNt2rQ8++yzNbZp0KBBevXqNddtlnTCLAAAAIAlQJs2bbLKKqvM9f5Ro0YlSdq3b1+jvX379pX7vmrs2LGZMWPGQm2zpBNmAQAAACwBDj/88AwfPry+y1jiCbMAAAAASqBDhw5JktGjR9doHz16dOW+r2rTpk0aNmy4UNss6YRZAAAAACXQpUuXdOjQocborYkTJ+Yf//hHevToUes2jRs3zkYbbVRjm5kzZ2b48OFz3WZJJ8wCAAAA+IZ99tlneeGFF/LCCy8k+XIx9xdeeCEjR46s9Ln00kuz5ZZbznUfVVVVOfroo3P22Wfnr3/9a/71r3+lX79+6dixY3baaadKvy233DKXXnpp5faAAQNy1VVX5brrrsurr76aQw45JJMmTUr//v3r/DgXh6XquwAAAACAb7tnnnkmP/jBDyq3BwwYkCTZd999c+211yb5crH2ESNGzHM/P//5zzNp0qQcfPDBGT9+fDbbbLPcd999adKkSaXPiBEjMnbs2Mrt3XffPR999FFOO+20jBo1KhtssEHuu+++ORaFL4uqoiiK+i5icZo4cWJatWqVCRMmpGXLlvVdDnVkqwa71ncJwFwMm3lzfZcAzIPvUFhy+Q799vi2/h46ZcqUvPPOO+nSpUuNIAUW1YKeU6YZAgAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAALDI/seuK8c3aEHPJWEWAAAAsNAaNWqUJJk8eXI9V8K3xaxzada5NTdLLY5iAAAAgG+Xhg0bpnXr1hkzZkySpFmzZqmqqqrnqiijoigyefLkjBkzJq1bt07Dhg3n2V+YBQAAACySDh06JEkl0IKvo3Xr1pVzal6EWQAAAMAiqaqqynLLLZd27drliy++qO9yKLFGjRrNd0TWLMIsAAAA4Gtp2LDhAgcR8HVZAB4AAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA06j3Muuyyy9K5c+c0adIk3bt3z9NPPz3P/pdccknWWGONNG3aNJ06dcoxxxyTKVOmLKZqAQAAAKhP9Rpm3XTTTRkwYEAGDhyY5557Ll27dk3v3r0zZsyYWvsPGTIkJ554YgYOHJhXX301f/jDH3LTTTflF7/4xWKuHAAAAID6UK9h1kUXXZSDDjoo/fv3z9prr53BgwenWbNmufrqq2vt/8QTT6Rnz57Zc88907lz52y99dbp27fvfEdzAQAAAPDtUG9h1rRp0/Lss8+mV69e/1dMgwbp1atXnnzyyVq32XTTTfPss89Wwqu3334799xzT7bbbru5Ps7UqVMzceLEGj8AAAAAlNNS9fXAY8eOzYwZM9K+ffsa7e3bt89rr71W6zZ77rlnxo4dm8022yxFUWT69On52c9+Ns9phoMGDcoZZ5xRp7UDAAAAUD/qfQH4hfHwww/n3HPPzeWXX57nnnsut912W+6+++6cddZZc93mpJNOyoQJEyo/77///mKsGAAAAIC6VG8js9q0aZOGDRtm9OjRNdpHjx6dDh061LrNqaeemn322ScHHnhgkmS99dbLpEmTcvDBB+fkk09OgwZzZnPV1dWprq6u+wMAAAAAYLGrt5FZjRs3zkYbbZThw4dX2mbOnJnhw4enR48etW4zefLkOQKrhg0bJkmKovjmigUAAABgiVBvI7OSZMCAAdl3333TrVu3bLzxxrnkkksyadKk9O/fP0nSr1+/LL/88hk0aFCSpE+fPrnooouy4YYbpnv37nnrrbdy6qmnpk+fPpVQCwAAAIBvr3oNs3bfffd89NFHOe200zJq1KhssMEGue+++yqLwo8cObLGSKxTTjklVVVVOeWUU/LBBx+kbdu26dOnT84555z6OgQAAAAAFqOq4n9sft7EiRPTqlWrTJgwIS1btqzvcqgjWzXYtb5LAOZi2Myb67sEYB58h8KSy3fot4ffQ6FulepqhgAAAAD8bxNmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlEa9h1mXXXZZOnfunCZNmqR79+55+umn59l//PjxOeyww7Lccsuluro6q6++eu65557FVC0AAAAA9Wmp+nzwm266KQMGDMjgwYPTvXv3XHLJJendu3def/31tGvXbo7+06ZNy1ZbbZV27drllltuyfLLL5/33nsvrVu3XvzFAwAAALDY1WuYddFFF+Wggw5K//79kySDBw/O3XffnauvvjonnnjiHP2vvvrqjBs3Lk888UQaNWqUJOncufPiLBkAAACAelRv0wynTZuWZ599Nr169fq/Yho0SK9evfLkk0/Wus1f//rX9OjRI4cddljat2+fddddN+eee25mzJixuMoGAAAAoB7V28issWPHZsaMGWnfvn2N9vbt2+e1116rdZu33347Dz74YPbaa6/cc889eeutt3LooYfmiy++yMCBA2vdZurUqZk6dWrl9sSJE+vuIAAAAABYrOp9AfiFMXPmzLRr1y5XXnllNtpoo+y+++45+eSTM3jw4LluM2jQoLRq1ary06lTp8VYMQAAAAB1qd7CrDZt2qRhw4YZPXp0jfbRo0enQ4cOtW6z3HLLZfXVV0/Dhg0rbWuttVZGjRqVadOm1brNSSedlAkTJlR+3n///bo7CAAAAAAWq3oLsxo3bpyNNtoow4cPr7TNnDkzw4cPT48ePWrdpmfPnnnrrbcyc+bMStsbb7yR5ZZbLo0bN651m+rq6rRs2bLGDwAAAADlVK/TDAcMGJCrrroq1113XV599dUccsghmTRpUuXqhv369ctJJ51U6X/IIYdk3LhxOeqoo/LGG2/k7rvvzrnnnpvDDjusvg4BAAAAgMWo3haAT5Ldd989H330UU477bSMGjUqG2ywQe67777KovAjR45Mgwb/l7d16tQpQ4cOzTHHHJP1118/yy+/fI466qiccMIJ9XUIAAAAACxGVUVRFPVdxOI0ceLEtGrVKhMmTDDl8Ftkqwa71ncJwFwMm3lzfZcAzIPvUFhy+Q799vB7KNStUl3NEAAAAID/bcIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUxlKLstGMGTNy7bXXZvjw4RkzZkxmzpxZ4/4HH3ywTooDAAAAgNktUph11FFH5dprr83222+fddddN1VVVXVdFwAAAADMYZHCrBtvvDF/+ctfst1229V1PQAAAAAwV4u0Zlbjxo2z6qqr1nUtAAAAADBPixRmHXvssfn1r3+doijquh4AAAAAmKtFmmb497//PQ899FDuvfferLPOOmnUqFGN+2+77bY6KQ4AAAAAZrdIYVbr1q2z884713UtAAAAADBPixRmXXPNNXVdBwAAAADM1yKFWbN89NFHef3115Mka6yxRtq2bVsnRQEAAABAbRZpAfhJkyZl//33z3LLLZfvfe97+d73vpeOHTvmgAMOyOTJk+u6RgAAAABIsohh1oABA/LII4/krrvuyvjx4zN+/PjceeedeeSRR3LsscfWdY0AAAAAkGQRpxneeuutueWWW7LFFltU2rbbbrs0bdo0u+22W6644oq6qg8AAAAAKhZpZNbkyZPTvn37OdrbtWtnmiEAAAAA35hFCrN69OiRgQMHZsqUKZW2zz//PGeccUZ69OhRZ8UBAAAAwOwWaZrhr3/96/Tu3TsrrLBCunbtmiR58cUX06RJkwwdOrROCwQAAACAWRYpzFp33XXz5ptv5oYbbshrr72WJOnbt2/22muvNG3atE4LBAAAAIBZFinMSpJmzZrloIMOqstaAAAAAGCeFjjM+utf/5ptt902jRo1yl//+td59t1hhx2+dmEAAAAA8FULHGbttNNOGTVqVNq1a5eddtpprv2qqqoyY8aMuqgNAAAAAGpY4DBr5syZtf4bAAAAABaXBouy0fXXX5+pU6fO0T5t2rRcf/31X7soAAAAAKjNIoVZ/fv3z4QJE+Zo//TTT9O/f/+vXRQAAAAA1GaRwqyiKFJVVTVH+3/+85+0atXqaxcFAAAAALVZ4DWzkmTDDTdMVVVVqqqqsuWWW2appf5v8xkzZuSdd97JNttsU+dFAgAAAECykGHWrKsYvvDCC+ndu3datGhRua9x48bp3LlzdtlllzotEAAAAABmWagwa+DAgZkxY0Y6d+6crbfeOsstt9w3VRcAAAAAzGGh18xq2LBhfvrTn2bKlCnfRD0AAAAAMFeLtAD8uuuum7fffruuawEAAACAeVqkMOvss8/Occcdl7/97W/58MMPM3HixBo/AAAAAPBNWKg1s2bZbrvtkiQ77LBDqqqqKu1FUaSqqiozZsyom+oAAAAAYDaLFGY99NBDdV0HAAAAAMzXIoVZ3//+9+u6DgAAAACYr0UKs5Jk/Pjx+cMf/pBXX301SbLOOutk//33T6tWreqsOAAAAACY3SItAP/MM89klVVWycUXX5xx48Zl3Lhxueiii7LKKqvkueeeq+saAQAAACDJIo7MOuaYY7LDDjvkqquuylJLfbmL6dOn58ADD8zRRx+dRx99tE6LBAAAAIBkEcOsZ555pkaQlSRLLbVUfv7zn6dbt251VhwAAAAAzG6Rphm2bNkyI0eOnKP9/fffz9JLL/21iwIAAACA2ixSmLX77rvngAMOyE033ZT3338/77//fm688cYceOCB6du3b13XCAAAAABJFnGa4YUXXpiqqqr069cv06dPT5I0atQohxxySM4777w6LRAAAAAAZlmkMKtx48b59a9/nUGDBmXEiBFJklVWWSXNmjWr0+IAAAAAYHaLFGbN0qxZs7Ru3brybwAAAAD4Ji3SmlnTp0/PqaeemlatWqVz587p3LlzWrVqlVNOOSVffPFFXdcIAAAAAEkWcWTWEUcckdtuuy3nn39+evTokSR58sknc/rpp+fjjz/OFVdcUadFAgAAAECyiGHWkCFDcuONN2bbbbettK2//vrp1KlT+vbtK8wCAAAA4BuxSNMMq6ur07lz5znau3TpksaNG3/dmgAAAACgVosUZh1++OE566yzMnXq1Erb1KlTc8455+Twww+vs+IAAAAAYHaLNM3w+eefz/Dhw7PCCiuka9euSZIXX3wx06ZNy5Zbbpkf//jHlb633XZb3VQKAAAAwP+8RQqzWrdunV122aVGW6dOneqkIAAAAACYm0UKs6655pq6rgMAAAAA5muRwqxZPvroo7z++utJkjXWWCNt27atk6IAAAAAoDaLtAD8pEmTsv/++2e55ZbL9773vXzve99Lx44dc8ABB2Ty5Ml1XSMAAAAAJFnEMGvAgAF55JFHctddd2X8+PEZP3587rzzzjzyyCM59thj67pGAAAAAEiyiNMMb7311txyyy3ZYostKm3bbbddmjZtmt122y1XXHFFXdUHAAAAABWLNDJr8uTJad++/Rzt7dq1M80QAAAAgG/MIoVZPXr0yMCBAzNlypRK2+eff54zzjgjPXr0qLPiAAAAAGB2izTN8JJLLsk222yTFVZYIV27dk2SvPjii2nSpEmGDh1apwUCAAAAwCyLFGatt956efPNN3PDDTfktddeS5L07ds3e+21V5o2bVqnBQIAAADALAsdZn3xxRdZc80187e//S0HHXTQN1ETAAAAANRqodfMatSoUY21sgAAAABgcVmkBeAPO+yw/PKXv8z06dPruh4AAAAAmKtFWjPrn//8Z4YPH577778/6623Xpo3b17j/ttuu61OigMAAACA2S1SmNW6devssssudV0LAAAAAMzTQoVZM2fOzAUXXJA33ngj06ZNyw9/+MOcfvrprmAIAAAAwGKxUGtmnXPOOfnFL36RFi1aZPnll89vfvObHHbYYd9UbQAAAABQw0KFWddff30uv/zyDB06NHfccUfuuuuu3HDDDZk5c+Y3VR8AAAAAVCxUmDVy5Mhst912ldu9evVKVVVV/vvf/9Z5YQAAAADwVQsVZk2fPj1NmjSp0daoUaN88cUXdVoUAAAAANRmoRaAL4oi++23X6qrqyttU6ZMyc9+9rM0b9680nbbbbfVXYUAAAAA8P9bqDBr3333naNt7733rrNiAAAAAGBeFirMuuaaa76pOgAAAABgvhZqzSwAAAAAqE/CLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaSwRYdZll12Wzp07p0mTJunevXuefvrpBdruxhtvTFVVVXbaaadvtkAAAAAAlgj1HmbddNNNGTBgQAYOHJjnnnsuXbt2Te/evTNmzJh5bvfuu+/muOOOy+abb76YKgUAAACgvtV7mHXRRRfloIMOSv/+/bP22mtn8ODBadasWa6++uq5bjNjxozstddeOeOMM7LyyisvxmoBAAAAqE/1GmZNmzYtzz77bHr16lVpa9CgQXr16pUnn3xyrtudeeaZadeuXQ444IDFUSYAAAAAS4il6vPBx44dmxkzZqR9+/Y12tu3b5/XXnut1m3+/ve/5w9/+ENeeOGFBXqMqVOnZurUqZXbEydOXOR6AQAAAKhf9T7NcGF8+umn2WeffXLVVVelTZs2C7TNoEGD0qpVq8pPp06dvuEqAQAAAPim1OvIrDZt2qRhw4YZPXp0jfbRo0enQ4cOc/QfMWJE3n333fTp06fSNnPmzCTJUkstlddffz2rrLJKjW1OOumkDBgwoHJ74sSJAi0AAACAkqrXMKtx48bZaKONMnz48Oy0005Jvgynhg8fnsMPP3yO/muuuWb+9a9/1Wg75ZRT8umnn+bXv/51rSFVdXV1qqurv5H6AQAAAFi86jXMSpIBAwZk3333Tbdu3bLxxhvnkksuyaRJk9K/f/8kSb9+/bL88stn0KBBadKkSdZdd90a27du3TpJ5mgHAAAA4Nun3sOs3XffPR999FFOO+20jBo1KhtssEHuu+++yqLwI0eOTIMGpVraCwAAAIBvSFVRFEV9F7E4TZw4Ma1atcqECRPSsmXL+i6HOrJVg13ruwRgLobNvLm+SwDmwXcoLLl8h357+D0U6pYhTwAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJTGEhFmXXbZZencuXOaNGmS7t275+mnn55r36uuuiqbb755vvOd7+Q73/lOevXqNc/+AAAAAHx71HuYddNNN2XAgAEZOHBgnnvuuXTt2jW9e/fOmDFjau3/8MMPp2/fvnnooYfy5JNPplOnTtl6663zwQcfLObKAQAAAFjc6j3Muuiii3LQQQelf//+WXvttTN48OA0a9YsV199da39b7jhhhx66KHZYIMNsuaaa+b3v/99Zs6cmeHDhy/mygEAAABY3Oo1zJo2bVqeffbZ9OrVq9LWoEGD9OrVK08++eQC7WPy5Mn54osvsswyy9R6/9SpUzNx4sQaPwAAAACUU72GWWPHjs2MGTPSvn37Gu3t27fPqFGjFmgfJ5xwQjp27FgjEJvdoEGD0qpVq8pPp06dvnbdAAAAANSPep9m+HWcd955ufHGG3P77benSZMmtfY56aSTMmHChMrP+++/v5irBAAAAKCuLFWfD96mTZs0bNgwo0ePrtE+evTodOjQYZ7bXnjhhTnvvPPywAMPZP31159rv+rq6lRXV9dJvQAAAADUr3odmdW4ceNstNFGNRZvn7WYe48ePea63fnnn5+zzjor9913X7p167Y4SgUAAABgCVCvI7OSZMCAAdl3333TrVu3bLzxxrnkkksyadKk9O/fP0nSr1+/LL/88hk0aFCS5Je//GVOO+20DBkyJJ07d66srdWiRYu0aNGi3o4DAAAAgG9evYdZu+++ez766KOcdtppGTVqVDbYYIPcd999lUXhR44cmQYN/m8A2RVXXJFp06blJz/5SY39DBw4MKeffvriLB0AAACAxazew6wkOfzww3P44YfXet/DDz9c4/a77777zRcEAAAAwBKp1FczBAAAAOB/izALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAQAAAFAawiwAAAAASkOYBQAAAEBpCLMAAAAAKA1hFgAAAAClIcwCAAAAoDSEWQAAAACUhjALAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGkIswAAAAAoDWEWAAAAAKUhzAIAAACgNJaIMOuyyy5L586d06RJk3Tv3j1PP/30PPvffPPNWXPNNdOkSZOst956ueeeexZTpQAAAADUp3oPs2666aYMGDAgAwcOzHPPPZeuXbumd+/eGTNmTK39n3jiifTt2zcHHHBAnn/++ey0007Zaaed8vLLLy/mygEAAABY3KqKoijqs4Du3bvn//2//5dLL700STJz5sx06tQpRxxxRE488cQ5+u++++6ZNGlS/va3v1XaNtlkk2ywwQYZPHjwfB9v4sSJadWqVSZMmJCWLVvW3YFQr7ZqsGt9lwDMxbCZN9d3CcA8+A6FJZfv0G8Pv4dC3VqqPh982rRpefbZZ3PSSSdV2ho0aJBevXrlySefrHWbJ598MgMGDKjR1rt379xxxx219p86dWqmTp1auT1hwoQkX36Y8O0xvfiivksA5sLnLSzZfIfCkst36LfHrNeynseSwLdGvYZZY8eOzYwZM9K+ffsa7e3bt89rr71W6zajRo2qtf+oUaNq7T9o0KCcccYZc7R36tRpEasGYGG0atWqvksAgFLyHfrt8+mnn3pdoQ7Ua5i1OJx00kk1RnLNnDkz48aNy7LLLpuqqqp6rAyozcSJE9OpU6e8//77hmADwELwHQpLrqIo8umnn6Zjx471XQp8K9RrmNWmTZs0bNgwo0ePrtE+evTodOjQodZtOnTosFD9q6urU11dXaOtdevWi140sFi0bNnSf8QBYBH4DoUlkxFZUHfq9WqGjRs3zkYbbZThw4dX2mbOnJnhw4enR48etW7To0ePGv2TZNiwYXPtDwAAAMC3R71PMxwwYED23XffdOvWLRtvvHEuueSSTJo0Kf3790+S9OvXL8svv3wGDRqUJDnqqKPy/e9/P7/61a+y/fbb58Ybb8wzzzyTK6+8sj4PAwAAAIDFoN7DrN133z0fffRRTjvttIwaNSobbLBB7rvvvsoi7yNHjkyDBv83gGzTTTfNkCFDcsopp+QXv/hFVltttdxxxx1Zd9116+sQgDpUXV2dgQMHzjE9GACYN9+hAPyvqCpcGxQAAACAkqjXNbMAAAAAYGEIswAAAAAoDWEWAAAAAKUhzAIAAACgNIRZAAAAAJSGMAsAAACA0hBmAfXunnvuydChQ+doHzp0aO699956qAgAAIAllTALqHcnnnhiZsyYMUd7URQ58cQT66EiACiHa665JjfffPMc7TfffHOuu+66eqgIAL55wiyg3r355ptZe+2152hfc80189Zbb9VDRQBQDoMGDUqbNm3maG/Xrl3OPffceqgIAL55wiyg3rVq1Spvv/32HO1vvfVWmjdvXg8VAUA5jBw5Ml26dJmjfaWVVsrIkSProSIA+OYJs4B6t+OOO+boo4/OiBEjKm1vvfVWjj322Oywww71WBkALNnatWuXl156aY72F198Mcsuu2w9VAQA3zxhFlDvzj///DRv3jxrrrlmunTpki5dumTNNdfMsssumwsvvLC+ywOAJVbfvn1z5JFH5qGHHsqMGTMyY8aMPPjggznqqKOyxx571Hd5APCNqCqKoqjvIgCKosgDDzyQF154IU2bNs3666+f733ve/VdFgAs0aZNm5Z99tknN998c5ZaaqkkyYwZM7Lvvvtm8ODBady4cT1XCAB1T5gFLBGGDx+e4cOHZ8yYMZk5c2aN+66++up6qgoAyuHNN9/M888/X/mD0EorrVTfJQHAN2ap+i4A4IwzzsiZZ56Zbt26ZbnllktVVVV9lwQApfGHP/whF198cd58880kyWqrrZajjz46Bx54YD1XBgDfDGEWUO8GDx6ca6+9Nvvss099lwIApXLaaafloosuyhFHHJEePXokSZ588skcc8wxGTlyZM4888x6rhAA6p5phkC9W3bZZfP0009nlVVWqe9SAKBU2rZtm9/85jfp27dvjfY///nPOeKIIzJ27Nh6qgwAvjmuZgjUuwMPPDBDhgyp7zIAoHS++OKLdOvWbY72jTbaKNOnT6+HigDgm2eaIVDvpkyZkiuvvDIPPPBA1l9//TRq1KjG/RdddFE9VQYAS7Z99tknV1xxxRzflVdeeWX22muveqoKAL5ZphkC9e4HP/jBXO+rqqrKgw8+uBirAYDyOOKII3L99denU6dO2WSTTZIk//jHPzJy5Mj069evxh+I/HEIgG8LYRYAAJTUvP4gNDt/HALg20SYBQAAAEBpWAAeAAAAgNIQZgEAAABQGsIsAAAAAEpDmAUAAABAaQizAOB/wKhRo3LUUUdl1VVXTZMmTdK+ffv07NkzV1xxRSZPnlzf5QEAwAJbqr4LAAC+WW+//XZ69uyZ1q1b59xzz816662X6urq/Otf/8qVV16Z5ZdfPjvssMMc233xxRdp1KhRPVQMAABzZ2QWAHzLHXrooVlqqaXyzDPPZLfddstaa62VlVdeOTvuuGPuvvvu9OnTJ0lSVVWVK664IjvssEOaN2+ec845J0lyxRVXZJVVVknjxo2zxhpr5I9//GNl3++++26qqqrywgsvVNrGjx+fqqqqPPzww0mShx9+OFVVVbn77ruz/vrrp0mTJtlkk03y8ssvL7bnAACAbw9hFgB8i3388ce5//77c9hhh6V58+a19qmqqqr8+/TTT8/OO++cf/3rX9l///1z++2356ijjsqxxx6bl19+OT/96U/Tv3//PPTQQwtdy/HHH59f/epX+ec//5m2bdumT58++eKLLxb52AAA+N8kzAKAb7G33norRVFkjTXWqNHepk2btGjRIi1atMgJJ5xQad9zzz3Tv3//rLzyyllxxRVz4YUXZr/99suhhx6a1VdfPQMGDMiPf/zjXHjhhQtdy8CBA7PVVltlvfXWy3XXXZfRo0fn9ttv/9rHCADA/xZhFgD8D3r66afzwgsvZJ111snUqVMr7d26davR79VXX03Pnj1rtPXs2TOvvvrqQj9mjx49Kv9eZpllssYaayzSfgAA+N9mAXgA+BZbddVVU1VVlddff71G+8orr5wkadq0aY32uU1FnJsGDb78u1hRFJU2UwcBAPgmGZkFAN9iyy67bLbaaqtceumlmTRp0kJvv9Zaa+Xxxx+v0fb4449n7bXXTpK0bds2SfLhhx9W7p99MfjZPfXUU5V/f/LJJ3njjTey1lprLXRNAAD8bzMyCwC+5S6//PL07Nkz3bp1y+mnn571118/DRo0yD//+c+89tpr2Wijjea67fHHH5/ddtstG264YXr16pW77rort912Wx544IEkX47s2mSTTXLeeeelS5cuGTNmTE455ZRa93XmmWdm2WWXTfv27XPyySenTZs22Wmnnb6JQwYA4FtMmAUA33KrrLJKnn/++Zx77rk56aST8p///CfV1dVZe+21c9xxx+XQQw+d67Y77bRTfv3rX+fCCy/MUUcdlS5duuSaa67JFltsUelz9dVX54ADDshGG22UNdZYI+eff3623nrrOfZ13nnn5aijjsqbb76ZDTbYIHfddVcaN278TRwyAADfYlXF7ItcAADUsYcffjg/+MEP8sknn6R169b1XQ4AACVnzSwAAAAASkOYBQAAAEBpmGYIAAAAQGkYmQUAAABAaQizAAAAACgNYRYAAAAApSHMAgAAAKA0hFkAAAAAlIYwCwAAAIDSEGYBAAAAUBrCLAAAAABKQ5gFAAAAQGn8f92/rLjpCkFDAAAAAElFTkSuQmCC",
+      "text/plain": [
+       "<Figure size 1200x800 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "import spacr\n",
+    "import os\n",
+    "from scipy.stats import chi2_contingency\n",
+    "import pandas as pd\n",
+    "import matplotlib.pyplot as plt\n",
+    "%matplotlib inline\n",
+    "\n",
+    "def analyze_class_proportion(settings):\n",
+    "    \n",
+    "    from spacr.utils import annotate_conditions, save_settings\n",
+    "    from spacr.io import _read_and_merge_data\n",
+    "    #from spacr.settings import set_analyze_class_proportion_defaults\n",
+    "    \n",
+    "    def _plot_proportion_stacked_bars(settings, df, group_column, bin_column, prc_column='prc', level='object'):\n",
+    "        # Always calculate chi-squared on raw data\n",
+    "        raw_counts = df.groupby([group_column, bin_column]).size().unstack(fill_value=0)\n",
+    "        chi2, p, dof, expected = chi2_contingency(raw_counts)\n",
+    "        print(f\"Chi-squared test statistic (raw data): {chi2:.4f}\")\n",
+    "        print(f\"p-value (raw data): {p:.4e}\")\n",
+    "\n",
+    "        # Plot based on level setting\n",
+    "        if level == 'well':\n",
+    "            # Aggregate by well for mean ± SD visualization\n",
+    "            well_proportions = (\n",
+    "                df.groupby([group_column, prc_column, bin_column])\n",
+    "                .size()\n",
+    "                .groupby(level=[0, 1])\n",
+    "                .apply(lambda x: x / x.sum())\n",
+    "                .unstack(fill_value=0)\n",
+    "            )\n",
+    "            mean_proportions = well_proportions.groupby(group_column).mean()\n",
+    "            std_proportions = well_proportions.groupby(group_column).std()\n",
+    "\n",
+    "            ax = mean_proportions.plot(\n",
+    "                kind='bar', stacked=True, yerr=std_proportions, capsize=5, colormap='viridis', figsize=(12, 8)\n",
+    "            )\n",
+    "            plt.title('Proportion of Volume Bins by Group (Mean ± SD across wells)')\n",
+    "        else:\n",
+    "            # Object-level plotting without aggregation\n",
+    "            group_counts = df.groupby([group_column, bin_column]).size()\n",
+    "            group_totals = group_counts.groupby(level=0).sum()\n",
+    "            proportions = group_counts / group_totals\n",
+    "            proportion_df = proportions.unstack(fill_value=0)\n",
+    "\n",
+    "            ax = proportion_df.plot(kind='bar', stacked=True, colormap='viridis', figsize=(12, 8))\n",
+    "            plt.title('Proportion of Volume Bins by Group')\n",
+    "\n",
+    "        plt.xlabel('Group')\n",
+    "        plt.ylabel('Proportion')\n",
+    "\n",
+    "        # Update legend with formatted labels, maintaining correct order\n",
+    "        plt.legend(title=f'Classes', bbox_to_anchor=(1.05, 1), loc='upper left')\n",
+    "        plt.ylim(0, 1)\n",
+    "        fig = plt.gcf() \n",
+    "        return chi2, p, dof, expected, raw_counts, fig\n",
+    "    \n",
+    "    #settings = set_analyze_endodyogeny_defaults(settings)\n",
+    "    #save_settings(settings, name='analyze_endodyogeny', show=True)\n",
+    "    output = {}\n",
+    "\n",
+    "    # Process data\n",
+    "    if not isinstance(settings['src'], list):\n",
+    "        settings['src'] = [settings['src']]\n",
+    "    \n",
+    "    locs = []\n",
+    "    for s in settings['src']:\n",
+    "        loc = os.path.join(s, 'measurements/measurements.db')\n",
+    "        locs.append(loc)\n",
+    "        \n",
+    "    if 'png_list' not in settings['tables']:\n",
+    "        settings['tables'] = settings['tables'] + ['png_list']\n",
+    "            \n",
+    "    df, _ = _read_and_merge_data(\n",
+    "        locs, \n",
+    "        tables=settings['tables'], \n",
+    "        verbose=settings['verbose'], \n",
+    "        nuclei_limit=settings['nuclei_limit'], \n",
+    "        pathogen_limit=settings['pathogen_limit']\n",
+    "    )\n",
+    "        \n",
+    "    df = annotate_conditions(\n",
+    "        df=df, \n",
+    "        cells=settings['cell_types'], \n",
+    "        cell_loc=settings['cell_plate_metadata'], \n",
+    "        pathogens=settings['pathogen_types'],\n",
+    "        pathogen_loc=settings['pathogen_plate_metadata'],\n",
+    "        treatments=settings['treatments'], \n",
+    "        treatment_loc=settings['treatment_plate_metadata']\n",
+    "    )\n",
+    "    \n",
+    "    if settings['group_column'] not in df.columns:\n",
+    "        print(f\"{settings['group_column']} not found in DataFrame, please choose from:\")\n",
+    "        for col in df.columns:\n",
+    "            print(col)\n",
+    "    \n",
+    "    df[settings['class_column']] = df[settings['class_column']].fillna(0)\n",
+    "    output['data'] = df\n",
+    "    \n",
+    "    # Perform chi-squared test and plot\n",
+    "    chi2, p, dof, expected, raw_counts, fig = _plot_proportion_stacked_bars(settings, df, settings['group_column'], bin_column=settings['class_column'], level=settings['level'])\n",
+    "\n",
+    "    # Create a DataFrame with chi-squared test results and raw counts\n",
+    "    results_df = pd.DataFrame({\n",
+    "        'chi_squared_stat': [chi2],\n",
+    "        'p_value': [p],\n",
+    "        'degrees_of_freedom': [dof]\n",
+    "    })\n",
+    "    \n",
+    "    output['chi_squared'] = results_df\n",
+    "    \n",
+    "    if settings['save']:\n",
+    "        output_dir = os.path.join(settings['src'][0], 'results')\n",
+    "        os.makedirs(output_dir, exist_ok=True)\n",
+    "        output_path_chi = os.path.join(output_dir, 'class_chi_squared_results.csv')\n",
+    "        output_path_data = os.path.join(output_dir, 'class_chi_squared_data.csv')\n",
+    "        output_path_fig = os.path.join(output_dir, 'class_chi_squared.pdf')\n",
+    "        fig.savefig(output_path_fig, dpi=300, bbox_inches='tight')\n",
+    "        results_df.to_csv(output_path_chi, index=False)\n",
+    "        df.to_csv(output_path_data, index=False)\n",
+    "        print(f\"Chi-squared results saved to {output_path_chi}\")\n",
+    "        print(f\"Annotated data saved to {output_path_data}\")\n",
+    "\n",
+    "    plt.show()     \n",
+    "\n",
+    "    return output\n",
+    "\n",
+    "%matplotlib inline\n",
+    "\n",
+    "settings = {\n",
+    "    'src': '/nas_mnt/carruthers/Einar/tsg101_screen/TSG101SCREEN_20240810_132824/plate1/',\n",
+    "    'tables': ['cell', 'nucleus', 'pathogen', 'cytoplasm'],\n",
+    "    'cell_types': ['Hela'],\n",
+    "    'cell_plate_metadata': None,\n",
+    "    'pathogen_types': ['nc', 'pc'],\n",
+    "    'pathogen_plate_metadata': [['c1'], ['c2']],\n",
+    "    'treatments': None,\n",
+    "    'treatment_plate_metadata': None,\n",
+    "    'group_column': 'pathogen',\n",
+    "    'class_column':'test',\n",
+    "    'pathogen_limit': 1,\n",
+    "    'nuclei_limit': 1000,\n",
+    "    'level': 'well',\n",
+    "    'save':True,\n",
+    "    'verbose': True}\n",
+    "\n",
+    "# Run the analysis\n",
+    "output = analyze_class_proportion(settings)\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "adc1df67",
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "spacr",
+   "language": "python",
+   "name": "spacr"
+  },
+  "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.9.19"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}