septum-mec/actions/stimulus-response/data/20_stimulus-spike-response....

766 lines
597 KiB
Plaintext
Raw Permalink Normal View History

2019-10-17 17:50:31 +00:00
{
"cells": [
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": null,
2019-10-17 17:50:31 +00:00
"metadata": {},
"outputs": [],
"source": [
"%load_ext autoreload\n",
"%autoreload 2"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
2021-03-10 12:30:50 +00:00
"outputs": [],
2019-10-17 17:50:31 +00:00
"source": [
"import os\n",
"import expipe\n",
"import pathlib\n",
"import numpy as np\n",
"import spatial_maps.stats as stats\n",
"import septum_mec\n",
"import septum_mec.analysis.data_processing as dp\n",
"import septum_mec.analysis.registration\n",
"import head_direction.head as head\n",
"import spatial_maps as sp\n",
"import speed_cells.speed as spd\n",
"import re\n",
"import joblib\n",
"import multiprocessing\n",
"import shutil\n",
"import psutil\n",
"import pandas as pd\n",
"import matplotlib.pyplot as plt\n",
"import matplotlib\n",
2021-03-10 12:30:50 +00:00
"import seaborn as sns\n",
2019-10-17 17:50:31 +00:00
"from distutils.dir_util import copy_tree\n",
"from neo import SpikeTrain\n",
"import scipy\n",
"\n",
2019-12-13 10:43:57 +00:00
"from tqdm.notebook import tqdm_notebook as tqdm\n",
"tqdm.pandas()\n",
2019-10-17 17:50:31 +00:00
"\n",
2021-03-10 12:30:50 +00:00
"from septum_mec.analysis.statistics import load_data_frames, make_paired_tables, make_statistics_table\n",
2019-10-17 17:50:31 +00:00
"\n",
"from spikewaveform.core import calculate_waveform_features_from_template, cluster_waveform_features\n",
"\n",
"from septum_mec.analysis.plotting import violinplot, despine"
2019-10-17 17:50:31 +00:00
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"plt.rc('axes', titlesize=12)\n",
"plt.rcParams.update({\n",
" 'font.size': 12, \n",
" 'figure.figsize': (6, 4), \n",
" 'figure.dpi': 150\n",
"})\n",
"\n",
"output_path = pathlib.Path(\"output\") / \"stimulus-response\"\n",
"(output_path / \"statistics\").mkdir(exist_ok=True, parents=True)\n",
"(output_path / \"figures\").mkdir(exist_ok=True, parents=True)\n",
"output_path.mkdir(exist_ok=True)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
2021-03-10 12:30:50 +00:00
"project_path = dp.project_path()\n",
"project = expipe.get_project(project_path)\n",
"actions = project.actions"
2019-10-17 17:50:31 +00:00
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"colors = ['#d95f02','#e7298a']\n",
"labels = ['11 Hz', '30 HZ']\n",
2021-03-10 12:30:50 +00:00
"queries = ['frequency==11 and stim_location==\"ms\"', 'frequency==30 and stim_location==\"ms\"']"
2019-10-17 17:50:31 +00:00
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 6,
"metadata": {},
2019-10-17 17:50:31 +00:00
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Number of sessions above threshold 194\n",
2021-03-10 12:30:50 +00:00
"Number of animals 4\n",
"Number of individual gridcells 139\n",
"Number of gridcell recordings 230\n"
2019-10-17 17:50:31 +00:00
]
}
],
"source": [
2021-03-10 12:30:50 +00:00
"data, labels, colors, queries = load_data_frames(queries, labels, colors)"
2019-10-17 17:50:31 +00:00
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 7,
2019-10-17 17:50:31 +00:00
"metadata": {},
"outputs": [],
"source": [
2021-03-10 12:30:50 +00:00
"keys = [\n",
" 't_i_peak',\n",
" 't_e_peak',\n",
" 'p_i_peak',\n",
" 'p_e_peak',\n",
"]"
2019-10-17 17:50:31 +00:00
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 8,
2019-12-13 10:43:57 +00:00
"metadata": {},
"outputs": [],
"source": [
2021-03-10 12:30:50 +00:00
"results, labels = make_paired_tables(data, keys, queries=queries, labels=labels)"
2019-12-13 10:43:57 +00:00
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 9,
2019-12-13 10:43:57 +00:00
"metadata": {},
2019-10-17 17:50:31 +00:00
"outputs": [
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5yU1bnA8d+U7Wyhd+k8AjYsoIAFe4vdmCgabNGo8cZ7bYkm6r0xWJKoMSYYExsGI/ZK7CKgKLGLehBpItKXZfvszLz3j/PO7OwyOzu7O7Oz5fl+PnzeOe973jNndoedZ071OI6DUkoppVRTvJmugFJKKaU6Ng0WlFJKKZWQBgtKKaWUSkiDBaWUUkolpMGCUkoppRLSYEEppZRSCWmwoJRSSqmENFhQSimlVEIaLCillFIqIQ0WlFJKKZWQBgtKKaWUSkiDBaWUUkolpMGCUkoppRLSYEEppZRSCfkzXQGlRGQPY8ynaSo7sgf7y8aYo9PxHJ2diMwEHnCTPzbG/CuD1WmTtv6+G/0sphtj3kpR1dqsqdcmIocAb7rJXxpjbmnvuqmuT4MFlTEiUgz8L3Ap+l5USqkOS/9Aq0z6I3BepiuhlFIqMQ0WVCb50v0ExhhPup9DdRz6+1YqPXSAo1JKKaUS0mBBKaWUUglpN4RqdyJyI3BDo3ORkd4LjDGHpPC50jYbQkQeBH4CGGPMriKyN/ArYBrQC9gAvAH83hjzeTNlZQEzgdOBPdz7twOfAE8ADxhjAs2UMQ04y33+QUARUA6sAxYA9xhjvmrF6+wNvAXs5p76tTHmt60oZzhwJXAUsAtQCXwMzDbGzBOR2cBFwBpjzPBG961ykycD32PHu+wNVANfAJcYYz5t7vctIh7gOOBnwO5AX+Bb7M/41iRfh88t41zs72qwW49lwDzgb8aYmgT3H+XeewDQ3713BfAScLcxZksy9VCqPWnLglIpICKnA+8Cp2I/ALKAodhg4iMRmZHg3tHYoOBvwBEx9/cFDgdmA5+IiDRxf56IPAEsBC7Gfqj3wn4Z6In9ULwM+ExEzmrh6yoCXqbtgcLR2A/TS4HRQLZbt+nAYyLyKMl9eZmInSY4Bch1y9gTWJlEHfKAp4DngWOxv59cYAzwS+AzYGQzZQzABoDPAie5+XOAEmAqcBfwoYjsEufeAhF5Cvg3cAY2YIrcuy/wG+AbEflBc69FqfamwYLKhNnYP/rPx5yb6P67ICM1apsBwEPYD/h/YD/wDwR+B9RgPwQfFpFjG9/ofvgsBMYBAeAe7LfWScCJbrkhYFfgTREZGOf5Z2ODFID/AOcDhwAHYVsr3nGv+YF73ZaCZolIAfbb7j7uqdYGCvsDzwH52Nf4B2yQcCC2hakc+BFwdhLFXQ94sB/uU4FzgBuNMRVJ3DsH+wEP8Ck2kNsf+7N7FRs8XJfgdeRiA5WD3FNvAzPcMk7BBlVgf5fPua1FkXu92J/Bye6pF7ABw2RsQPg7oAzbGvS0iByWxOtRqt1oN4Rqd8aYDcAGEdkWc+7jDFaprYrd4wxjzD9jzi8SkRex30RzgLtE5FVjTF1MntnYYGMHcLgxZmmjsp9zWw2eAwYCd2A/WAEQkTHUf8i+BxxsjKmNuX8h8JCIPA6cBhQAxwCPJHpB7gfjs9gPZGh9oOAF/owNpALAEcaYt2OyLBKRedgP3r5JFOkFLjfG3Oem30mUOaYeR1AfUL0OHNfo5/SUiNwFXJ6gmOuxQRvY7oLGeZ+O6ZraE/t7muNeuxw41H18sTHm3kb3vi4i/wAWY98PD4jIqEbvFaUyRlsWlEqNhxsFCgAYY94BbneTo7HfIgEQkbHACW7y5jiBQqSMF7AtDACni8igmMu7AV9jWzB+1+gDMFZs3QYneiHuN+IngMi321YFCq5p1LdM/KFRoACAO47iv5Msr5r6n0VLXOQeQ8B5TfycrsSOHdiJG/REWr2+Af6niee5Eoh8wB8Vc2/k9c2PEygAYIxZCVzrJodiWyuU6hA0WFAqNe5JcO3+mMex/dHHYpvUwTaDJ/KSe/RiuxgAMMY8bYwRbBP/83Hui9gQ8zgnQT4fMBfbFQJtCxSg4et9oMlc8Bh2QGdzPmxuoGdjIuLHdg0BLDbGrI2Xz/0W31SLy37YsSQAjzb1jd8dnLg30McYExmnsgf2wx+a/z3Pj3msXRGqw9BuCKXargb4oKmLxphVIlKKHYy3a8yliTGPP2xi/GI8Ow3CM8ZEZgEgIv3dPKOB8djBc1Nisif6kvBbYHhMemuylWrCXu6xzBjzdVOZjDF1IvIRdixDIt+2og6RmSFgZ18k8n4T58fGPP4wUQFxZr7E/p7/KCJ/bKYOEQkHWyrVnjRYUKrtNhljQs3k2YINFgbEnOvTyufr2fiEOx3vEuBg6sdQxAonWfZw91iLbYG4VUReMMa05kMa6r+NJxN0bEoiz4421CGZemxMQRmNpez3rFSmaLCgVNsFk8gTWdo6tgk99v/fJOr7upuzOfLAXTfgb+w8i+Rb4CvsqP93sE38rydZ/v9i++UfAgqBe7FdJq2R7R6T6fJ0ms+SVJ5E9zS3HHRTv4O2/K2MvfcS7BTbZFS34TmVSikNFpRqu15J5Il8u/w+5ty2mMffGWPWt+K5L6Y+UPgaO1f/FWNMbNmRbYyTcacx5gb3np9gR/AfIyJnG2PmJL41ri2AAH1ExBPbXRJHa7+BNyf2Z97cjIumfpexP8+kpp42cW9ZJ5/5o7opDRaUarsSERlmjFkT76K7mFKkz/yTmEuxfdv7YxcMiktEJmMHNq7GDtJb51661D2GgKOMMat2vhuwCwAl472YxxdjWyZygTtE5GVjTDJdBbE+xk6/7IFd/Gh5vEzujIG94l1LgQ3YD+xe2IGKiTRVhy9iHu8JPN1UASIyFzvI8SvsugqNf89zE9zbF9v6sBr42BjzSVN5lWpPOhtCZVKy/eidQaIFhWK34Y4NCF6OefyzZsq/HbgF+BcNB76Ndo+lCQIFsMtARyT1JcEdkBiZCdEbu15CS8XO0Ei0euRxpKllwR1P8oybnCQiu8fL53bpNLXS5lLqx0uc4QY38crIxXbZCNDbbUlZSn3rwpkiEm9MScTPgRuBB7GLcinVIWiwoDIpOtddRHpksiIp8EsR2afxSRGZDlzhJhcbY6Kj7Y0x/8EuRgRwuIj8Kl7BIvI/2NUOwX5TXxhzObKPQB8R2TfOvR4RuQk4MuZ0oqmTjd1G/bfq00XkpESZ43gl5v6rRGRSnDoOAe5uYbktdTe29QXsapolcfL8EjtzZCfuugx/c5O7Yj/Q4/kD9QNM73PvDVAfaPUG5ojITr8DEZkKXOUmqyP3K9URaDdEBonIcwDGmBOay9tFxfYl/05EHgZCxpiPMlWhNsgHFrjT4l7B/t86HvtNMQs7vfKiOPddgF2iuQi4WUQOxi4ZvQY75e8s6hfnCQA/bdTvP4/6YORFEbnNLc8BJmA3LGrc9J7om20D7pTGi7BBjQf4i4i8ZYxJZk0EjDGOe/+bQB7wloj8Cbs/QgA7pfNqGo4laM0gxubq8bH7s/kltqvhYxG5BfgI26IxE7vCZSV2lct4bsIuojUW+LW7cdgD2I26hmOX2Y6s5/A2DddsmIV9P+yNXXviYxG5E9stVYJdrOsSbJcPwLXGmNj/H0pllLYsZNao0aNH/wD7x7Hb/XvmmWdu8vkikwT4ObB00KBBH6b4eQCYNm3aUamu/8knn/yTSPkXXnghHo+nAPg19pv/m9hV/rJ79+7NY489luvOv29QhjFm+TPPPFM0eHB0UcUjsQsULcF2WZwKeIqLi7nvvvuy3ZaJ6P0ffPDBFbvtFtnjiX7A77E7RC4A/gLsl5WVxbXXXkv//nb238SJEy+MLWPWrFnRxZL++Mc/PhqnjgvPOOOMyCyCgaeeemppS35OxpiFd9xxhz83NxdswHCN+/NZjO1e6XvmmWcybNgwAEaMGDE89v7XX399VaR+Z5xxxkWt/X1/9dVXvzz//PMjWYYBf3V/zi8Ap5WUlDBr1qxooPDwww+
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5xcVfn48c/M7GxN75UkpDwEQgm9ExABCUiXrqiggBW/iGD5Cd+viogFRRQrYAQFaSJFlEhvQiAkBHjSSCe9sLvZ3dmZub8/zp2ym2m7O7P1eb9e+7r3zD33zJnJZueZUwOe52GMMcYYk02wqytgjDHGmO7NggVjjDHG5GTBgjHGGGNysmDBGGOMMTlZsGCMMcaYnCxYMMYYY0xOFiwYY4wxJicLFowxxhiTkwULxhhjjMnJggVjjDHG5GTBgjHGGGNysmDBGGOMMTlZsGCMMcaYnCxYMMYYY0xOZV1dAdO3iMg+qrqgq+thchORS4A7/OT5qvrXLqxOh4iI558+qaonteP+S0i9F8eq6jNFqlqHZXttIjILeNpPXqeqP+zsupnexYIF0ylEZCDwv8AXsN87Y4zpUeyPtuksPwU+09WVMMYY03YWLJjOEurqCpi+SVUDXV0HY3o6G+BojDHGmJwsWDDGGGNMTtYNYUpKRK4HvtvqscQI7mdVdVYJnrMauAI4HdgDGABsAV4H7gHuU9V4sZ/Xf+47gU8Bqqp7iMj+wDeBI4EhwHrgP8CPVfXtPGWFgUuAc4B9/Pu3A28B9wN3qGokTxlHAhf6zz8G917UAmuAZ4HbVPW9drzOocAzwAz/oe+o6vfaUc5E4GrgRGA3oB6YD9yuqveJyO3A54GVqjqx1X3v+8kzgA9w42L2BxqAd4ArVXVBvtkQIhIAZuN+Z/YGhgOrce/xTQW+jpBfxqdx/1Zj/XosAu4DfquqjTnuP9G/9zBgpH/vUuBx4FZV3VxIPYwpFWtZML2KiBwELAZ+jPuAHAaUA6OBU4G/AM+LyOhOqMs5wMvAWbgPgDAwHhdMvCkiF+W4dwouKPgt8NG0+4cDxwO3A2+JiGS5v0pE7geeBy7HfagPwX1BGIz7UPwisFBELmzj6xoAPEnHA4WTcB+mXwCm4P6dBgPHAveKyF8o7AvNTNw0wcOBSr+MfYHlBdShCngQ+AdwMu7fpxKYClwHLAR2z1PGKFwA+HdcgLo7UAEMAo4Afg68ISK7Zbi3RkQeBP4JnIsLmBL3Hgj8P2CZiJya77UYU0oWLJhSux33x/wfaY/N9H8uLeYTicheuA+NsbhvzzcBJwCHAJ8AHvazHg78S0Rqivn8rYwC7sJ9wP8B94F/FPADoBH3IfgnETk5w+sYhfuQnw5EgNtw31oPBk7zy43hWk2ezhL43I4LUsC1qHwWmAUcjWuteMm/Vgb8xm8pyMt/zx4HDvAfam+gcCjwCFCNe40/wQUJR+FaomqB84CLCyju20AA9+F+BPBJ4HpVrSvg3jm4D3iABbhA7lDce/dvXPDwrRyvoxL3O3e0/9BzwEV+GWfigipw/5aP+K1FiXuDuPfgDP+hR3EBwyG4gPAHwA5ca9BDIvKRAl6PMSVh3RCmpFR1PbBeRLamPTa/RE/3Z6AG18Q+S1WXpV37L/A3EbkS9+E7A/gOcG2J6jLQP16kqnenPf6CiDyG+yZaAfxcRP6tqs1peW7HBRsfAser6mutyn7EbzV4BNdi8jPcBysAIjKV1Ifsq8AxqtqUdv/zwF0i8jfgbNx79jHc+5eV/8H4d9wHMrQ/UAgCv8QFUhHgo6r6XFqWF0TkPtwH7/ACigwCX1bV3/npl3JlTqvHR0kFVHOB2a3epwdF5OfAl3MU821c0Aauu6B13ofSuqb2xf07zfGvfRk4zj+/XFV/0+reuSLyB+BF3O/DHSIyudXvijGdwloWTK/g/+Hfz0/+T6tAIUlVf4Xrawe4Iv2bXgn8qVWgkKjDS8DNfnIK7lskACIyDfi4n/x+hkAhUcajuBYGgHNEZEza5RnAElwLxg9afQCmS6/b2FwvxH+f7gcS327bFSj4jiTVMvGTVoECAP44iq8VWF4DqfeiLT7vH2PAZ7K8T1fjxg7swg96Eq1jy4D/yfI8VwOJD/gT0+5NvL4nMgQKAKjqclIB7Xhca4Uxnc6CBdNbzE47/3eevI/7xwHAQaWpDuBaMLL5Y9p5en/0ybgmdSj8dQRxXQwAqOpDqiq4Jv5/ZLgvYX3aeUWOfCHcwNDEe9yRQAFavt47suaCe3EDOvN5I99Az9ZEpAzXNQTwoqquypTP/xafrcXlINxYEoC/ZPvG7w9O3B8YpqqJcSr74D78If+/8xNp59YVYbqEdUOY3mJm2vnWLOP+MtmdAput26gRmJftoqq+LyLbcIPx9ki7lP463mjj62j9HIlZAIjISD/PFGBP3OC5w9Oy5/ri8D1gYlp6S6GVyiLRArRDVZdky6SqzSLyJm4sQy6r21GHxMwQcLMvcvlvlsenpZ2/kauADDNf0v+dfyoiP81Th4Scgy2NKRULFkxvMayd9w0uai1SNqpqLE+ezf7zj0p7rGivw5+OdyVwDKkxFOkKnT460T824VogbhKRR1W1PR/SkPo2XkjQsbGAPB92oA6F1GNDEcporbv9vhqTkwULprdI/C5vJtW8XIi1JagLQLSAPIklsNOb0NP/Tx5Mqq87n02JE3/dgN+y62yT1cB7uFH/L+Ga+OcWWP7/4vrl7wL6A7/BdZm0R7l/LKQb1MufpaA8ue7Jtxx0tn+Djvz9TL/3StwU20I0dOA5jWk3CxZMb5GYbdEPWFCqRZfaYEgBeRLfLj9Ie2xr2vlaVV3Xjue+nFSgsAQ3V/9fqppedmIb40Lcoqrf9e/5FG4E/8dE5GJVnZP71ow2AwIME5FAendJBu39Bp5P+nueb8ZFtn/L9PezoKmnWe7dUcIZQsYUhQULprd4Gze3vRLXJ561D9lfDGgvYAVuFclSrI43SEQmqOrKLHUQUn3mb6VdSu/bPhS3YFBGInIIbmDjCtwgvTX+pS/4xxhwoqq+v+vdgFsAqBCvpp1fjmuZqAR+JiJPqmohXQXp5uOmX/bDLX60OFMmf8bAfpmuFcF63Af2EPIPcs1Wh3fSzvcFHspWgIjcgxvk+B5uXYXW/8735Lh3OK71YQUwX1XfypbXmFKx2RCms5T6m/6TaedXZsvkL8v7K9wKj/fj1hgolVwLCqVv150eEKS/jivylH8z8EPgr7Qc+DbFP27LESiAWwY6oaAvDv6AxMRMiKG49RLaKn2GRq7VI2dTopYFfzxJYpGug0Vk70z5/C6dbCttvkZqvMS5fnCTqYxKXJeNAEP9lpTXSLUuXCAimcaUJHwJuB64E7colzGdzoIF01mSc9hFpF8Jyn+Y1PK+n8mxhPHPgEn++d+zffMvkutE5IDWD4rIscBVfvJFVU2OtlfV13GLEQEcLyLfzFSwiPwPbrVDcN/Un0+7nGgpGSYiB2a4NyAiN+BWt0zINXWytR+R+lZ9joicnitzBv9Ku//rInJwhjqOA25tY7ltdSuu9QXcapqDMuS5DjdzZBf+ugy/9ZN74D7QM/kJqQGmv/PvjZAKtIYCc0Rkl38DETkC+LqfbEjcb0xns26ILiQijwCo6sfz5e0F0vuIfyAifwJiqvpmMQpX1aiIfBK39G4Y98c3sRfEetyI/s+RWjFvG/DVYjx3DtXAs/60uH/h/r+dgvumGMZNr/x8hvsuxS3RPAD4vogcg1syeiVuyt+FpBbniQCfa9Xvfx+pYOQxEfmRX56H6375NLs2vef6ZtuCP6Xx87igJgD8SkSeUdVC1kRAVT3//qeBKuAZEfkFbn+ECG5K5zW0HEvQnkGM+eox339vrsN1NcwXkR8Cb+JaNC7BrXBZT/YWqBtwi2hNA77jbxx2B24V0Ym4ZbYTA26fo+WaDTfifh/2x609MV9EbsF1Sw3CLdZ1Ja7LB+BaVU3/f2RMp7GWha41ecqUKafi/hD26p+HH374hlAoMfifLwGvjRkz5o1iPoeqvnD
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd3Rc1bXA4d/MqMuWJfcCNtjGm9CLccCY3kOvIXQILfTewgskDwKE5OGEEiCEEBwIJdTQewdDaKGYbVNsg3uRZFuWVWbu++PcKZKl0UiaIo32t5bWnXPLuWc0tu6eUwOe52GMMcYY055grgtgjDHGmJ7NggVjjDHGJGXBgjHGGGOSsmDBGGOMMUlZsGCMMcaYpCxYMMYYY0xSFiwYY4wxJikLFowxxhiTlAULxhhjjEnKggVjjDHGJGXBgjHGGGOSsmDBGGOMMUlZsGCMMcaYpCxYMMYYY0xSBbkugDHZIiL3AMf7yRGqujCHxemxRCS6bv3zqrp3TgvTDen4vBN+F39X1RPSVLRuS/beRGQ2MAZQVd0w64UzeclqFowxxhiTlAULxhhjjEnKmiFMn+FXI5+Q42KYLLHP25j0sZoFY4wxxiRlwYIxxhhjkrJmCJMVCb3KLwduAE4GTgE2xP07/Ap4GPiTqtZlqAz3kKHRECKyM/Cqn9wO+AK4DDgMGA2sAT4G/grcr6peG9kk5rctcCqwEzASaAJmAy/gfkdzO7h+GO73uxsgwECgGVgGfAA8ADzSUTnayfu3uM8R4D/AHqpa08k8ivzyHQVsChTi3t+jwO+AHwHv+qfvoqqvJVz7Gu738gRwDPB/wKFAKfA9cLuq3pTK5y0i44DzgF2B9YCVwJvAjar6forvZUP/veyBG4VQCMwFXgL+qKqzkly7LnA2sJd//yJggV+G21X13fauNSabLFgw2RYEHsI9RBNt6f8cLyJ7dvQw7OGGAn/DBUJRJcAu/s9hIvIzVV3T+kIRKQD+BPyi1aES3EN1U+AsETlHVe9s6+YicjzwZ9zDM1ExUI4LXg4F/i0iB6tqONU3JiK/pPuBwiDgGWBSq0MbAlcAxwFXppBVoZ/PDgn7JgDLUyzHsbjgrTBhdxnu3+bBInJWB9cH/PL+hrVracX/OVlETlHVaW1c/3PgFtxnm2h9/+c4EbkTOEtVm1J5T8ZkijVDmGw7G/fHeDHuG91k4EDgWf+4AK+JSFluipcWd+AefB8DRwPb4jrazfCPHwT8vZ1r7yIeKLyF+2a8LbAz7sG0EPfQv8MPCloQkV2Be3CBwnLcg2wfP49DgVuBRv/0/XE1PCkRkfOAa/xkVwOFAuBF4oHCK7h/D9vigoSPgHWANgOhVvbGBQpPArvjvp3/GVdD1VE59sd9BoVAHfBrP69dcDUbzbjfVTJX4X4fQdy/5ytwNR67Ab8F6nGf1d9EZLtW9z8B91mXAN8BF/r3n4yrpfjEP/VUUvtdGJNRVrNgsm04rrp5B1X9IWH/kyLyJ1wwsT5wKe6PcW80HPeN92BVjT6Yp4vIv3BV09sCR4jI7aoabbpARA4iXm1+nape0Srf10XkLuA1YCPgVhF5SlWXJZzzG3/bDOylqv9plcejIvIc8G8/fTguuElKRE4FbvKTXQoUfGfiapAA/qyqZyQcmy4iDwAPAgenkFcQF1AdrKoRf98LHV0kIsXAVCAArML9W/wk4ZTXROQJ4GXW/tYfzeNHuOAAYBawk6ouSDjlFRF5FXgOCOGCit38a0cSD0ReBg5s1fT2rt+Eci/wM+AEEXlQVZ/r6L0ZkylWs2By4aRWgULUhcC3/utT/Wre3qgaODYhUADAfyAcB0QfbKe3uu5if/sF8Mu2MlbVJUD0AVtOQs2AXxtTjKtR+HcbgUI0j6eA6IN+VEdvRkSOwX1jh+4FCgDn+tvZuJql1mVrAk7E9a1IxR0JgUKq9gDG+q+vbxUoRMvxDq6GoT0nEm++OK1VoBDN4yVcrQfAtiJS4b8+E9fc0Qwc11YfHVVtxtUw1fq7zm19jjHZZMGCybYZid+mE/kPin/4yeHANlkrVXo9oKpttpv7nd3e8pN7+9XyiEgVrmMkwMsddDx8C9cRD/xvq37eq1V1G1UdxNp9QlqLdvYrTnaSiByKa9YIAh/SjUBBRDbG1RqB6+TZ2NZ5qlqLq11IxXtdKMpPEl7fn+S8u5Mc29/f/tDev2ff2bjApL+qrvD37etvv1TV+e1d6P8e3vaTO4pIYXvnGpNp1gxhsu3tDo5/lPB6QyClHuk9TCrvcUegAjfSYS6wBa5aHOAcETknxXuNbWtn9Nu2iJTjHtDjcL/PzYEpwLr+qcm+MGyKe5iG/HQjsKL90zu0RcLrNms9EkwnXoOSzPddKEe042mtqn7X3kmqOkdEFuM6rMb4NV7j/eRHa13YMo95ra4twP1eATZLGCXUkTJgGNBWjZwxGWfBgsm2dr9J+ZYkvB6eyYJkUGff41xgcBfvVdV6h4isg2vSOYB2gglcU0hHNYsj/W0DrgZiO9w35T92qaTuYRe1tINzF6eQX6OqNnSjHKmMmlhEq2ABGET8b2eqzSVRVXS9RrcKCxZMjliwYLKtuYPjoYTXbVZT9wJdeY+J/xevAR5J8V4thj2KyN7Av3D9GaJW4kZifIGbY+El4DFg4xTy/wT4Ka62ZDBwrYg8mewbeRJFCa87emCm8o2703NEtLoulT4xbQ1Z7M7fzcRrnyM+DDUVX3fjvsZ0iwULJtsGdnB8SMLrtTqN9RJdeY+J33Lr2+p01xERGY6bbKkc95C7Dvgnbqlir9W5/VLI8ltgN1VdLiIX4HrnlwN/wQ1V7KzE2oQh7Z7ldLWmJRULcIHSYBEJdNA/pK3Psjrh9aBO3jvxcw515XM2Jhesg6PJti06OJ7YqfHTTBYkg1J9jwuIV7d/nnB822QXi0iRiFwlIj9vNX7/GGCA//oaVb1KVb9qI1AopmWTQHtmRTtq+pMKvezv382fUKizEh+MW3Vw7tZdyD9VX/jbfrScOKsFEanEzcjYgt/08Y2f3DzZjURkPxGZJyJvisj2/rXRGoKtox1ck1x/uoicKSL7dnSuMZlkwYLJth1EZK0/wBB7iB3lJ2eo6lfZK1Za/UxEQm0dEJGNgB/7ycejD3K/I1z0Iba3iKzf1vW+Y4CrcZP6nJKwf3zC6w+TXH8o8fkDOvMAOh03bTXA70VkRCeuBTdJVbTD30/be/iJSAmu6SNTEpt4kgU9x9J+U8VL/naMiExOksd+uL4fU4iPQHne3w4EjmzvQhEZj5uP4RbgZn84pTE5YcGCybYQcI//QIjxe5jfjJuKGNx8/72VANeutdONs78X9wBqwj0EEv3B3xYCD/jfbFvnMYH4+H8P9zuLSqzm36fNgolManVN0qGTiVT1a+IzOFYSn3sh1evDCfceC1zfRvmi/w5Gtj6WRm8TH41xrois1aTiD/P8Tev9CW4l3vfhz/7Q19Z5TCYejLyqqtHaiD8S72syVUQ2bePaMtww4ujf6D8lKYsxGWfVWjkkIk8CqOoBuS5Llu0MfCAiN+IWkFoXOMvfD/A6bs7+3uxS/4FzJ66pYTPcwlKJkwF92eqae3Df+vfFTYf8mYjchBtGWIL7dnoe7kENMFVVP064/l+4WQUDwBn+A+cRXBAxCjet9s9ouRZCRQrt9ol+5+exMXCgiPxUVVOdEwHcLJBH4X4fF4rIJrjf0Q+4hZTOxA0rTdTVjoxtUtWI34zyPi5YekZE/gw8jutwuitwEdAfN2Vz6zU2UNXPROR6XAfFzYBPROQPuA6kFbhpo8/F/Y2tJ2FSJVWdJSJX4BZUG4SbufJWXIfHNbihlRcAG/iXvE/HU08bk1FWs5Bb48aPH78/7o9hvv8AsPXWW7P11lsDbIKbm3867iG3M8BOO+3ERx99tJM/T0Bay3DwwQfH1lJ46623FqQz73vvvTc2Mc/BBx/MkCFDwFVBP4mbOOhO/EDhzDPPRFX/p3Ueqhr59NNP991vv/2iWa2Dq214C1ftfTV+oHD00UczY8aM81td//H
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5hU5fXA8e/M9l1YeldAQA6KDQuKvcTYxW5ij1GjiTXFFhNNfho1ppjYjSYaIrF3RWPvYkFQAQ8iAiq9LmXb7NzfH++dssvs7Ozu3J0t5/M8PHfeufe+98zssnPmvW8JeZ6HMcYYY0xjwrkOwBhjjDHtmyULxhhjjEnLkgVjjDHGpGXJgjHGGGPSsmTBGGOMMWlZsmCMMcaYtCxZMMYYY0xaliwYY4wxJi1LFowxxhiTliULxhhjjEnLkgVjjDHGpGXJgjHGGGPSsmTBGGOMMWlZsmCMMcaYtPJzHYAxpmMQkdh69i+q6sE5DaYVROQ+4HS/OEhVl7Sgjth7cb+qnpGl0Fot3WsTkfnAMEBVdUybB2c6NGtZMMYYY0xaliwYY4wxJi27DWGM6VL82wZn5DgMYzoUa1kwxhhjTFqWLBhjjDEmLbsNYQKR1Fv8CuBG4CzgbGAM7vfuC+AR4O+quqEN4jkI+BEwARgAVAJzgeeBW1R1RUDX3Rd4zS9OAGYClwPHAUOBKuAT4F5gsqp6KapJrm834BxgH2AwUAvMB/6Hey8XNnH+ANzP4QBAgN5ABFgJfAg8CDzWVByN1P0H3M8b4CPgQFVd08w6Cv34TgK2BQpwr+9x4I/AVsB7/uH7qerrSee+jntfngJOAf4CHAuUAN8Ad6rqXzMZDSEiI4GLgf2B4cA64C3gJlX9IMPXMsZ/LQfiRiEUAAuBl4G/qeqXac7dHLgAOMi/fiGw2I/hTlV9r7FzjQmCJQsmaGHgYdyHY7Jx/r/TReT7TX3ItZSIlAGTgKMb7CoCdvb/XSwip6jqM0HEkKQ/8C9cwhRTDOzn/ztORH6oqlUNTxSRfODvwHkNdhXjPlS3Bc4XkQtV9e5UFxeR04E7cB+eyYqAMlzycizwjIgcrap1mb4wEfk1rU8U+uCSt/ENdo0BrgROA67KoKoCv569kp4bDazKMI5TcclbQdLTpbjf4aNF5Pwmzg/58f6eTVtvxf93loicraqTUpz/Y+BW3M822Rb+v9NE5G7gfFWtzeQ1GdNadhvCBO0C3B/ZZbhvarsDE4Ep/n4BXheR0mxfWETCwNMkEoVngROBXYHvAX8A1gLlwBMickC2Y2jgLtwH3yfAycBuuI52s/39RwH3N3LuPSQShbdx34x3A/bFfTAtwX3o3+UnBfWIyP7AfbhEYRXug+wQv45jgduAGv/wI3AtQRkRkYuBa/1iSxOFfOAlEonCq7jfm91wScI0YDMgZSLUwMG4ROFp3M/5IFyS9EgGcRyB+xkUABuA3/l17Ydr2Yjg3qt0rsa9H2Hc7/2VuBaPA3C/c5W4n9W/RGRCg+ufgftZFwNfA7/wr787rpViun/oOWT2XhiTFdayYII2ENeMvJeqfpv0/NMi8ndcMrEFcBnuj2w2XYhrRgY4V1XvarD/FRG5F3jHj/NfIjIywG9rA3HfeI9W1dgH81QReRTXNL0bcIKI3KmqsVsXiMhRJJrNr1fVKxvU+4aI3AO8DmwN3CYiz6rqyqRjfu9vI8BBqvpRgzoeF5EXgFjryvG45CYtETkH+KtfbFGi4PsZrqUJ4A5V/WnSvqki8iDwEJu2EKUSxiVUR6tq1H/uf02dJCJFwM1ACFiP+52dnnTI6yLyFPAKm37rj9WxFS45APgS2EdVFycd8qqIvAa8AOThkooD/HMHk0hEXgEmNrhF955/C+XfwA+BM0TkIVV9oanXZkxrWcuCaQtnNkgUYn4BzPMfn+M332aF36rwc784JUWiAICqzsP1IQDYHDgmWzGksBo4NSlRiMWwAfftOfbBdm6D837lb2cCv05VsaouB2IfsGUktQz4rTZFuBaFZ1IkCrE6ngViH/RDmnoxInIK7hs7tC5RALjI387HtUA1jK0W1+dkZcN9jbgrKVHI1IHACP/xDQ0ShVgc7+JaGBrzIxK3L37SIFGI1fEyrtUDYDcRKfcf/wx3uyMCnJaqL4+qRnAtTGv9py5qeIwxQbBkwQRtdvK35GT+B8B//OJAYJcsXnc73Ic/uObtdKYkPQ7yVsSDqpryvrnf2e1tv3iw3yyPiPTCdYwEeKWJjodv4zriQdLrUNWNqrqLqvZh074jDcU6+xWlO0hEjsXd1ggDH9OKREFExuJal8B18qxJdZyqrsW1LmTi/RaEcmjS48lpjvtnmn1H+NtvG/u9912AS0y6q2qF/9xh/naWqi5q7ET/fXjHL+4tIgWNHWtMtthtCBO0d5rYPy3p8Rggo57mGRiX9PgvIvKXDM8b0fQhLZbJe7E3rg/FYFzP+R1wzeIAF4rIhRleK+XriH3b9jt+bgGMxL3v2wN7kkiw0n2R2Bb3YZrnl2uAisYPb9IOSY9TtnokmUqiBSWdb1oQR6zj6VpV/bqxg1R1gYgsw3VYjfNbxkb5xWmbnFi/ju8anJuPe18BtksaTdSUUtzonlQtd8ZkjSULJmiNfkPyLU96PDCL1+3bwvN6ZTGGhpr7Xiwki69DRDbD3fo5ksaToihNtzgO9rfVuBaICbhvyn9rUaTuwy6mqSGsyzKor0ZVq1sRRyajJpbSIFkA+pD4m5rp7ZKYXrS8pbcXliyYgFmyYIIWaWJ/XtLjlM3PLZT8u/1TEmPzm1KZxRgaasl7kfw6rgUey/Ba9YY9isjBwKO4/gwx63AjMWbi5lh4GXgCGJtB/dNxI0vewSU014nI0+m+kadRmPS4qQ/MTL5xN3uOiAbnZdJ3JlUn2Nb8PU0+9wUSw1AzMbcV1zUmI5YsmKD1bmJ/v6THm3QGa4Xkb4drU3VWy4GWvBfJr6OyJa9DRAbiJlsqw33IXQ/8F7dUsdfg2G4ZVDkPOEBVV4nIz3G988uAf+CGKjZXcmtCv0aPclra0pKJxbhEqa+IhJroH5LqZ7k66XGfZl47+eec105+X42Js2TBBG2HJvYnd2qckcXrfp70eDfSdFgTkX641of5wHRVzWYcyXbAzS7YmNh7sZhEc3vD19Eof/bDK3BN0rOSZvk7BejhP75WVX/fyPlF1L8l0JgvYx01VXWSP6/DAcABIvJjVb03gzqSJX8w7ohrAWnMTs2suzlm4pKdbrj+C7NTHSQiPXEzMtajqtUi8hWuH8j26S4kIofjhqbOAy5X1XdEZC6uz8NOIpLvj3xo7PxzcS1R84EX0x1rTDbYaAgTtL1EZJM/rBD/cDrJL85W1S+yeN0PSXxbO0lEeqQ59gLgGlzv/olZjKGhH4pIXqodIrI1brIogCdj32r9jnAz/ecPFpEtUp3vOwX3Ou7BTeATMyrp8cdpzj+WxPwBzfkicS5u2mqAP4nIoGacC26SqliHvxNjI0EaEpFi3K2PoCTf4vlxmuNOpfFbFS/722EisnuaOg7H9f3Yk8QIlBf9bW/gB42dKCKjcPMx3IqbqtwSBRM4SxZM0PKA+/w/9HF+z/FbcFMMg5vHP2v84Xe3+sU+wCQ/OalHRPYgMY9BJa4pPSgCXJcihnJcU34Id5vg1gaH/NnfFgAP+t9sG9YxmsT4fw/33sYkN/MfkjIwkfENzkk7dDKZqs4lMYNjTxJzL2R6fl3StUcAN6SIL/b7Mrjhvix6h8RojItEZJNbKv4wz5QtM77bSPR9uMMf+tqwjt1JJCOvqepX/uO/kehrcrOIbJvi3FLccOPY3+6/p4nFmKyx2xA5JCJPA6jqkbmOJWD7Ah+KyE24BaQ2B873nwd4AzcXf7Zdj/sGtyNu/Pt0EbkZd7ujJ67J+ackvk1fnmoSnSy7zP/AuRt3q2E73KRQyZMBzWpwzn24b/2H4aZD/kxE/oobRliM+3Z6Me41Adysqp8knf8oblbBEPBT/wPnMVwSMQTXmvJD6q+FUJ7Bfftkf/TrGAtMFJE
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5hb5ZX48a+mN/eOOzY+2GBqMM0GU03AVFNCSyBAQnYJG5KQkLIb+O0SQmAJGyC0hBB6C8UUU0KzaQZMB3Nsg20wxrjb0zWS7u+P96rMWKPRzKhMOZ/nmUf36t776tXMSDp6y3kDnudhjDHGGNOagnxXwBhjjDFdmwULxhhjjEnJggVjjDHGpGTBgjHGGGNSsmDBGGOMMSlZsGCMMcaYlCxYMMYYY0xKFiwYY4wxJiULFowxxhiTkgULxhhjjEnJggVjjDHGpGTBgjHGGGNSsmDBGGOMMSlZsGCMMcaYlIryXQGTPhHZRVU/yFLZ0bXKn1HVI7LxGN2diJwF/N3fPVVV78tw+SuAsYCq6o6ZLLs3yubrJVdEZByw3N+9WVXP70AZlwK/83fHq+qKjFSuk1I9t2y/1kz7WbDQDYhIP+D/Af+O/c2MScleL8Zknr2QuodrgO/nuxLGdBP2ejEmwyxY6B4Ks/0AqhrI9mMYkyNZf73kit9lYK9Nk3c2wNEYY4wxKVmwYIwxxpiUrBuiC2sxijl6X3TWwsuqOjODj5W12RAicjvwPfxR/iKyB/BrYDowEFgDvABcraoftVFWMXAWcBKwi3/9ZuB94CHg76oabKOM6cDp/uNvB/QFqoFVwMvADar6aQee5yDgJWBn/67/VNX/aW85CeUdC1wI7AEUA58DjwHXquqGFNd9C/ghcCAwGogA64DXgfuAuarqtXZ9J+p7Ke7/dYuq9heR4cBPgaOBMUATsBi4H7hRVRtTlFUKfBc4EdgNGID7O38MPArcoqr1rTx+4n1Zeb0klL8LcBFwMDDCr+NC4M+q+pyIPA3Mavn4IjITeNHf3R33f3g5MBn3v/gBcCZQQhuzIUSkCDgZOBfYEegPfAbcDfwpzedRivtdn+nXYTiwFfe6ugu4Q1UjrVwb8B//NOBbwGD/OSzG/b/epKo16dTDdF3WsmBySkROwn1ozQGG4T4ER+OCiXdF5IwU107EvXndAhyWcP0Q4FDgJuB9EZFWri8XkYeABcD5uA/1gbigeQAwFbgA+FBETm/n8+oLPEPmAoU/4T4UD8a9+Vf69fstsFREZrRy3S+AN3EfHDsAZUAFbkrmd/wynxWRyo7WLc36T8d9sF+M+wCrAPoB++A+wN4UkYGtXLsL8b/z4cBQ4n/nmcC1gPpBUd6IyNnAO7jgdQzxOs7G/Y6vSrOoI4HHcUFRKe7DdhzwdRp1GIgLcO8GDsIFLOW4/8MrcP8Lg9ooYxLwNi4omOU/lxK/HofgpjDOF5H+Sa4dins93Qccgwt6SvzHnA5chftb7dvWczFdmwULXdtNuG8djyfct7v/c25eatQ5w4F/4N5U/4b7wJ8B/B5owH1o3yEiR7a80P+WugD3rScI3AAcBUwDjvXLDeM+mF4UkRFJHv8mXJAC7s3xHNyHzwG4N/zX/GNFwM1+S0Gb/A/ep4A9/bs6FSgAk4CfAJuAS4D9gSOAO/3jA4Cn/HnqifU4APgDbkDc+8DZwH64N+1zgU/8Uw/FTS3MlnJgLq7F5ka/7vvhArHV/jm74L5JNyMiOwD/AqIB3xO4VqS9geOBB/z7RwMvicjOCZfn7PUiInOA23CDKbfiWjRm4IK7/8X9j/4c9wHelv8GtuB+P/vjAtn/bqv1R0QKgHm43y3AfNw3/H2AM4C3cL/nf0tRxhDgFeJB7lxcC8M+wKm4YAO/Xne3uLYS1zqyP+Dhgo3jca/JbwPX4V7X2+GCp51SPR/TtVk3RBemqmuANSKyMeG+9/JYpc7q59+eoaqJbzyviMiTuK6IUuD/ROQ5VW1KOOcm4k2jh6rqWy3Knuu3GszFfbv6E+6bNBD7EDrT310IHNiiGXwB8A8ReRD3ZlmJe8O7K9UTEpEyXFPr/v5dnQ0UwH3YrwGmq+pnCfc/IyILgeuBKty3tpMSjp/tX7sB9/y2JBx71X9u7wATgHNE5BeqGu5kXZMpwX2IzlbVpxPuf11EngI+wrU0nCYiF7b4O/8F9+0c4Geqek2Lsh/1y7gd9ze6W0R2U1UvV68X/0PyWn93I+7vtDjhlBdF5FHgOVzLTlsKgFNU9Tl//7VUJyc4B/fBDC6Q/F5CgLFQRO7HBVfHpyjjT8R/3z9X1f9NOLbQf03NwwWYR4rIDFVd4B+/HJgChIDjVfWJFmU/LSJ34Fo+qnBfEPZJ87mZLsZaFkyu3dEiUABAVV/DffgBTMS9OQGxZtJj/N3LkwQK0TKewLUwAJwkItslHN4ZWIr7pvP7FP3liXUbmeqJ+OMnHsI11UJmAoWoH7cIFABQ1Rtw4yIAjvNbXKKi22tbBArRa7fiujGuAX5Deh9kHfVwi0AhWofluA9RcC0P46PHRGQ34n/3J5MECtEy/oELFsB9c851xtE5wCh/+1ctAgUAVPUVkrSctOKzhEChPX7o324C/q1lS4SqhnD5Jrb5X4BY8qqT/d35LQKFxDJ+lnDXLP/a/sB5/n23JgkUote/DfzR391bRPZu60mZrsmCBZNrN6Q4dlvC9tEJ20cSn2ve1pvqU/5tAa6LAQBVfURVBfeN9vEk10WtSdguTXFeIXAPrisEMhsorAceSXE8GhAV4Zq9o6KDMieLyK0iMrblhap6n6r+TFVvUNXazFQ3qWdTHEsMgvokbM9K2L65jfL/krD97XQrlSHR/80mWjTNt3BLmuW90d4KiMgw3MBXgMdbG0Coqptp/X/pCFyXIMAdrT2WnzJ7d6Cfqv7Wv3sm7rUE6b8mIR5Ym27GuiFMLjUAi1o7qKrLRWQTrk8+cW2E3RO232ll/GIy2yd5jNi3L/8Nd3tcS8YU3Eju/RJOTxVM/w9uEFpUq7MTOuDtNroH3k3YnpqwfQOuaboPro/+XBH5GPdm/i/gpSwHCIlWpDiW+MGW+B6UOP6grQ/Qd3HjAkpo/jvIhd3828Wpfp+qulZElpPQetKKLztQByEeQLfV1fImbkxOS5MStt9JVUCS7pzE1+TDnXlNmu7BWhZMLq1No498vX+b2Lw+uIOPN6DlHSIyS0QeE5HNuFaE13Dfqi7BNYGn2zQ/zr+NdmdcKSKjO1jPlta2cXx9wnZsRoGqLsN9O0+c9rkTbrDkE8AGEXlcRI4i+1JNlUtsLk/MThgdUBqh+XPcht88Hg3Qks6qyKJh/m3KOvra+luCG4fT0TpA24HqNxkoo6WMvSZN92AtCyaXQmmcE03Vm5grIfH/dBqu+Tcd66Ib/lzwW9h2VPyXuA/XD3CBw2bg+TTL/3+4JvV/4L7N34zrMsm2xA/YZmMvVPV1f9T5TOAEXFPzBP9wKW5a32wRuQs3IC7p3Pk8aW9a4+j/Sq6fQ4l/m86XrXTyWXQk50VrAVcyrb1eOvP+n3jtccDKNK9LOn7CdH0WLJhcSucbYPQbS+Ic840J21+p6mra73zigcJS4L+AZ1U1sexospx0XKuqv/Ov+R5u7MC3ReRMVb0z9aVtauv3NDRhe5tvhH4A8IL/E10K+BBcX/tRuNf9GbhxBZ2tayZF/xYFuP+Dda2d6A8ujX5L3djaeVmyHjf4dUhbJ9Lxb+BtSXx9tFWP1v6fEn9vg4Av2vH4ideu6+aztEwarBvC5FL/ZIPuovxkSn393fcTDiVmdUw59UpE9haRX4rIKSIyKuHQv/u3YWCWP9Av2YfMmFTlJ1iYsH0+bjwGwJ/8RDWd0dZ89GkJ27G+ZhHpIyJ7iUizfmFVXaGqf1PV44jnmQDXytCVfJCw3dao+T2JD85rd7bNTop+ME5KldzKnzGQrT76T3H/ywB7tXHubq3c/0nC9q6pChCRBSLykYjc5N/VntfkJBH5jYic7k9hNt2QBQvdQ1dqKu6sM1McS1xW+OGE7WcStn/
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd3hb5dn48a8kyzN7D0ImuUMCgbBngQ5ogTJKKbRAgZa20EHLW2jppu9bSikd/NpCobSUsimFskfZO0ADhBVuSCCBLLKHE9uypfP74zmSjh1JlmzJsuX7c12+zjk6Q49sS+fWM+4n5HkexhhjjDHZhMtdAGOMMcb0bhYsGGOMMSYnCxaMMcYYk5MFC8YYY4zJyYIFY4wxxuRkwYIxxhhjcrJgwRhjjDE5WbBgjDHGmJwsWDDGGGNMThYsGGOMMSYnCxaMMcYYk5MFC8YYY4zJyYIFY4wxxuRkwYIxxhhjcqoqdwFMZiIyW1VfLXc5TG4ichrwd3/z86p6c5GvvxiYCKiqzijmtfujSnhficgk4D1/80pVPbML17gA+Jm/OVlVFxelcN2U67WV+r1mcrNgoZcRkcHA/wLfwP4+xhSFva+M6R570/Q+vwO+VO5CGFNh7H1lTDdYsND7RMpdAGMqUMW8r/wmg1C5y2H6F+vgaIwxxpicLFgwxhhjTE7WDNFLdOidnHzM81efUNWDS/Cc9cBZwDHADGAQsBb4L3Aj8E9VTRT7ef3nvgY4Fb+Xv4jsBvwQOAAYBqwEHgV+o6qvd3KtKHAacDww2z9/AzAf+Bfwd1WNdXKNA4CT/Ocfh/tdbAaWAk8Al6nqW114ncOBx4Gd/Id+oqq/KPQ6gesdDZwN7AZEgXeBO4FLVXVtjvP2AL4GHARMABLAauA54GbgLlX1sp3fjfJegPu/3qiqQ0RkDPA/wKeB7YFWYAFwC/BnVW3Jca0a4IvAZ4FdgaG4v/MbwB3AX1S1KcvzBx8r9ftqNnAO8FFgrF/G54E/qOpDIvIAcFjH5xeRg4HH/M05uP/DC4Edcf+LrwKnANV0MhpCRKqAzwFn4N7bQ4BFwA3A7/N8HTW43/UpfhnGAJtw76vrgWuzfT6ISMh//i8AewAj/NewAPf/eoWqNuZTDtM7WM1CPyUiewJvA7/B3SBH4D6ExuI+yG8CnhKRsT1QluNxN63jgNG4m+AEXDDxsoicnOPcabgPr78AnwicPxL4OHAFMF9EJMv5dSLyL+Ap4EzcTX0YLpAeCuwMfBN4TUROKvB1DQIepHiBwu9xN8WP4j78G/zy/Rh4R0QOzHLe94AXcDeOHYBaoB43JPNE/5r/EZGGrpYtz/IfgLuxn4e7gdUDg4F9cDewF0RkWJZzZ5P+Ox8KjCL9dz4YuBRQPygqGxE5HXgJF7xuT7qMR+J+x5fkeanDgbtxQVEN7v05CViRRxmG4QLcG4BDcO/pOtz/4UW4/4XhnVxjOu5Lw/W4wGZ73OfDCOBjuCGMT4rIkAznjsK9n24GjsIFPdX+cx4AXIL7W+3b2WsxvYcFC73HFbhvE3cHHpvj/5xRzCcSkVm4bzDjcdH+xbgP4L1x3wbu8A/dj9LfRMYA/8B9qP4Nd8M/EPgl0Iy7aV8rIodneB1jcB9KOwIx4DLgCGAv4Gj/unHcjemxLIHPFbggBdyH45dxN5+P4D7wn/X3VQFX+jUFnfJ/Z/cBu/sPdStQAKYD3wHWA+cD+wOfBK7z9w8F7vPHqQfL8RHgV7gOcfOB03F/1wNw/1dv+od+HDe0sFTqgLtwNTZ/9su+Hy4QW+4fMxv3TbodEdkBeBhIBnz34GqR9gaOBf7pPz4BeFxEdgqc3pPvq+OAq3GdKTfhajQOxAV3v8X9j56Lu4F35v+Ajbjfz/64QPb/Oqv9EZEwcD/udwvwJO49vQ9wMvAi7vf89RzXGAk8TTrIvQtXw7AP8HlcsIFfrhs6nNuA+2zZH/BwwcaxuPfkp4A/4t7X43CfLbNyvR7Te1gzRC+hqiuBlSKyLvDYKyV6uutx30qXAger6qLAvheAW0Xk67ib707AT3A3qFIY7C9PVtXgB8/TInIvrimiBvh/IvKQqrYGjrmCdNXox1X1xQ7XvsuvNbgL9+3q97hv0kDqJnSKv/k8cFCHavCngH+IyK24D8sG3Afe9blekIjU4qpa9/cf6m6gAO5mvxI4oMPf60EReR74EzAA963t+MD+0/1z1+Je38bAvmf81/YSMBX4soh8T1Xj3SxrJtW4m+iRqvpA4PHnROQ+4HVcTcMXROTsDn/ny3HfzgG+q6q/63DtO/xrXIP7G90gIruqqtdT7yv/Jnmpv7kO93daEDjkMRG5A3gIV7PTmTBwgqo+5G8/m+vggC/jbszgAslTAwHG8yJyCy64OjbHNX5P+vd9rqr+NrDvef89dT8uwDxcRA5U1af8/RcCM4E24FhVvafDtR8QkWtxNR8DcF8Q9snztZkyspqFfkZEPoGr2gT3wbso03GqejmurR3gLL9fQKlc2yFQSJbhWdzND2Aa7sMJSFWTHuVvXpghUEhe4x5cDQPA8SIyLrB7J+Ad3DedX+ZoLw+WbXyuF+L/nv6Fq6qF4gQKSd/K9PdS1ctI/62O8WtckpLrqzoECslzN+GaMX4H/Ij8bmRddXuHQCFZhvdwN1FwNQ+Tk/tEZFfSf/d7MwQKyWv8AxcsgPvm/MkilTlfxwHb+es/6BAoAKCqT5Oh5iSLRYFAoRBf85frga93rIlQ1TZcvolt/hcglbzqc/7mkx0CheA1vht46DD/3CHAV/zHrsoQKCTP/y/wa39zbxHZu7MXZcrPgoX+54jAemcfRvf5y0HAnqUpDuBqMLK5OrD+6cD64aTHmuf7OsK4JgYAVPXfqiq4b7R3ZzgvaWVgvSbHcRFcx9Dk77iYgcIa4N859icDoipctXdSslPmjiJylYhM7Hiiqt6sqt9V1ctUdUtxipvRf3LsCwZBAwPrhwXWr+zk+pcH1j+Vb6GKJPm/2UqHqvkO/pLn9eYWWgARGY3r+Apwd7YOhKq6gez/S5/ENQkCXJvtufyU2XOAwar6Y//hg3HvJcj/PQnpwNr0YtYM0f/MCayvy9LvL5Mp5F8VWohmYF62nar6noisx7XJB+dGCL6Olwp8HR2fI/Xty//AnYKryZiJ68m9X+DwXAH2L3Cd0JKyjk7ogv920jzwcmB958D6Zbiq6YG4NvozROQN3If5w8DjJQ4Qghbn2Be8sQU/l4L9Dzq7gb6M6xdQTfvfQU9I1tYtyPX7VNVVIvIegdqTLD7oQhmEdADdWVPLC7g+OR1ND6y/lOsCGZpzgu/J27vznjS9j9Us9D8junje0KKWIm1VHm3ka/xlsHq9aK9DRA4TkTtFZAOuFuFZ3Leq83FV4PlWzU/yl8nmjItFZEIXy9nRqk72rwmsp0YUqOpC3Lfz4LDPWbjOkvcAa0XkbhEJ1jiVSq6hcsHq8mB2wmSH0gTtX+M2/OrxZICWcVRFCY32lznL6OvsbwmuH05XywCdB6ofFuEaHfW2zxZTRFaz0P8k/+ZrcCMP8rWsBGUB1xGqM8lUvcFcCcH/3b1w1b/5WJ1c8ceC/4Vte8V/gLu5vooLHDYAj+R5/f/FVan/A/dt/kpck0mpBW+w7fpeqOpzfq/zg4HP4Kqap/q7a3DD+o4UketxHeJKklujiwpNa5z8X+np11DtL/P5ApZPPouu5LzIFnBlku390p17QvDcY4AleZ6Xsf+E6V0sWOh/kr3CBwCv9oIbQz7fAJPfWIJjzNcF1pep6nIKdybpQOEd4KfAf1Q1eO1kspx8XKqqP/PPORXXd+BTInKKql6X+9ROdfZ7GhVY3+Ybof93ftT/SU4F/DFcW/sRuM+Ck3H9Crpb1mJK/i3CuP+D1dkO9DuXJr+lrst2XImswXV+HdnZgXT9G3hngu+PzsqR7f8p+HsbDrxfwPMHz11dwtFcpgwsWOh/XscNVarFtbNmbZcUkU/iqqwX47LN5VPFWqghIjJRVTN+C/GTKQ3yN+cHdgWzOu4D3J7tCfze1gf
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5hU1fnA8e+07bssvUrXV8HeCwqK0VgJIZbYNRb82Y1GjSZqEmOMKSbGgrF3Y0cRjWJHJYiigHCwAUqRtsD23Sm/P86d2dlldnZ2d+7MlvfzPPvMnFvP7JT73lM9kUgEpZRSSqnmeLOdAaWUUkp1bBosKKWUUiopDRaUUkoplZQGC0oppZRKSoMFpZRSSiWlwYJSSimlktJgQSmllFJJabCglFJKqaQ0WFBKKaVUUhosKKWUUiopDRaUUkoplZQGC0oppZRKSoMFpZRSSiWlwYJSSimlkvJnOwOqYxKRB4HTneRAY8yaLGanwxKR6BzvrxljfpzmY98AXO8k9zPGfJTO43c3IlIE9DPGfJPtvLSHiLwNjAdqjTF5bdh/OPCtk7zRGHND2jLXTslem5vfNdUyLVlQSnV5InICsAQ4KNt5Uaoz0pIFpVSXJiIHAk9mOx9KdWYaLKiEjDFnAGdkORtKpYMv2xlIJ2PMhGznQXU/Wg2hlFJKqaQ0WFBKKaVUUloN0QHFtfq9BrgFOBs4B9ge+54tAZ4G/mmMqXQpDw/iUm8IEZkAvOUk9wMWAVcDPwOGAjXAp8B9wOPGmEiCw8Qfb1/gXGwr6kFAPbAM+C/2f7Sihf37Y/+/EwEBegFBYAMwF1vf/WxL+Wjm2H/Evo8AHwM/MsZsau1xnGMNA34D/BjoA6wD3gHuNMZ8kGS/nsD5wNHAGKAAKMN+jl4BprU1Ty3kdzgNre4nA9OxVVunAWOBYmAl8BrwV2PM1y0cbz/s+3wgMBj7Hq0AZgH/MsYsTXL+qAdE5AHn+QhjzLLWv7KkeSwGLsJ+lgWIAF8CTwD/wP4fnkh0fhFZBgxztrsZuB37Xnuwn+ebjTGPp9IbQkR2AS7Dfr+GYj/LrwN/AmpTfC17Yb8XBzj5AvgaeBW4zRizOsm+OwAXYr9TQ5zX8B32e3+7MeaLVPKgOg4tWejYvMB/gHuAvbA/rvnAbsAfgXkiMjR72UuLfsD/gF8D2wF5QClwMPAo8JyINPeD6BeRO4EPgTOBkc7+xcBOwC+BpSJybnMnF5HTsReU3wMTgIFALlCI/ZGdgg3MXhSRVtV9i8i1pClQwF4gFwC/wF4oc7E/wicDs0Xk1mbysDM2GLsJe+HoAQSw//eDsBePr50LsZsKgDewAeB4bLCTi33PzgcWiciRzbyGXCd4/QAbbIzCvs9F2ODnImf/X7v7EpITkRHYIPcm7He0APs52hUb9H8M9E/hUCXAu8Bx2M9yEbAjNrBKJR9XO/k4nYbv1GDs/24+cGwL+wdE5Hbs9/Ic7P+40PnbGfgVsFhEDm1m/99gP6v/hw2YCrH/CwGmAgtE5AYR8aTyelTHoMFCxxa9Q1kLXArsD0wCZjrrBXhbRAqyk720mIYtMfkUe+HbF/ujtthZ/xPgoWb2vRd7oQF4H/vjuC/2ov9rYA32gjTNCQoaEZFDgAexAdhG4HfAEc4xpgB3AHXO5sdgS3hSIiKXAn9wku0NFMBebIqBx5w87g9c6eQb4AoR+WWTPPiAZ7ABUCVwA/ZOb2/s5+gxZ9NewNMikt+O/LXkr9gA8CPgVGAf7Hv7urM+F3jQGQsh/jV4gcdpKOVaDlwCjMO+z9cDm7Albjc5Y1NErcJetM+JW3a9s2w3Z31aiEgP7AV+lLPoWexnZj/gPOwd+Y7An1M4XPQiHw2sjsXejb+TQj4uxJZKeID1wBXYz8ph2O9aDvD3Fg7zb2ypANhA+mJs6cKPsd+JIDbofF5ERjY5/w3Y75EP+BwbHOyPfb8uwf4fvNj34XpUp6HVEB3bAGzx44HGmO/jlk8XkX9ig4kRwFV03i/eAGxR+GRjTPTCPEdEnsHeie4LHC8idxtjolUXiMhPaLiA3GyMaXpX+Y6I3Au8jb0zukNEXjbGbIjb5nfOYxA43BjzcZNjPCcirwIvOenjsD+4STklGdEf5HQECmB//C8wxtwZt+xDEXkemI29Y/29iDxkjFnvrB8HbOs8P88Y8xiNTReRVdigYzBwJPYi54YBwCPAGcaYcHShiEwHXnbO3Rc4Cngqbr+fAz91ns8GjjDGlMetf0dEHsK+z8OB3zrv88fO52m+iJTGbb/CGDM/ra/MugFb0gNwtTHmlrh1H4nIE9jqllRKcLzY6rf44PSl5jaOEpG+NASoq4D9jTHL4zZ5XURex5aUNXeMQ2n4Xn2I/V7E/79fE5G52CC7CLgWW9qFiOyOrSYD+16fZYwJxu07W0Tuw77fE7Dv1X+0SqJz0JKFju+sJoFC1C+B6Eh053biIr0y4NS4QAEApy3GaUD0wjK1yX5XOo+LsD9YWzHGrMMWhYItCo39+DqlMbnYO/OXEgQK0WO8jL1zBXtBTUpETgHucpLpChQAZjQJFKL5+xobLIItITkjbvWAuOdfNnPcf2DvJK8Bvmp/NptVA1waHygAOO1A/h23aJcm+0Xf51rgxCYXrugxltPwuj1x+2SEU00Wrer6oEmgAICT75OwgWkq7mp5k62cgL3jB/hVk0Ahmo9nsRfy5kRfRxg4rZn/90PYkkCAH8X99vwSe03ZAExtEihE960EzsK25fBgb3hUJ6DBQse2OP5uOp4xph5bpw/2orBXxnKVXk8aYzYmWmGM+RJbvQDwYxHxQ6zBXvQObVYLDQ/fB6I/eBPjjl1ljNnLGNMbW9WTTLRxZ26yjURkCvaOywvMI32BAth2K815ioZGa4fFLV8S9/wBETnEKdaPMcasNMaca4z5kzHmszTlNZF5zb3P2KLpqOLoExEZQEPwMKOZoBkAp4g+eod6WNPX6bKJ2Dp5sO9/Qk5jxv+mcLwgNtBsrWibj1qSlxDdn2ih8/2KDqP8gTEmWfA4BdgGGGaMiTgBwxHOutnGmKrmdjTGfEtDNePE5rZTHYtWQ3Rss1tY/0nc8+2xDZI6m1Re40HYRl+DsK3fd8XelQBcLCIXp3iukYkWRu92RaQQW60zCvv/3AVblL+Ns2myC9BO2Lr1aCPIOmBLivlKxZzmVhhjakRkMfb/slPc8s9EZCb2R3wMttfABhGZha3i+W+iu0+XLEuyriLuefxv0ti456nMi/ER9nWWYt+zTL22XeOet3SRn0PDRb05640xNW3Ix/bOo2lh/49puLOPN5CGYO0TknAu+PGGAz2d58fG9ehqyYgUt1NZpsFCx9ZSA6x1cc8HNLtVx9ba17gC25K+LXo2XSAiQ7DFp8fSTDCBLZJt6U51kPNYiy2B2A9bxPqPNuV0a+taWB9tp9CryfITsY3STsZeHHoDxzt/iMgCbJDzL2NMBe5Jduz4C0v8Bax33PO1KZzjh7jnvchcsBDfw2F9s1tZqbyOtgaZ0XxsSLaRMaZSRCqIK8Vpsn+Lx0igrd9Jv4gUJ6ruUB2LBgsdW0v1m/Fd+eqa3apja8trjP/c/oHUG+WF4hMi8mNsb4HCuMXl2CLSRdgxFt4AnqfxXW5z5mPrjWdjfzxvEpHpCe7C2qKlO7XoRbZRH3pjzBbgVBH5LbaB5tHYRqMBZ5OdsK3n/09EDm5prIMMa207nPjPSrjZrdIvJ+55S0FlKnfcrR7Po8l+qfzf6hMsa8/1IH7f+7FjRKSq2SoL1XFosNCxNb1LbKpv3PNmB0jp4NryGuPrvqvb0rrdqQ9/Ehso1GMvmE9gi3AjTbYt2voIW/kGmGiM2SgilwMPO8f+N5CwP3or9SL53V4/5zHhNk7A8mfgz87rORBbP30C9o5yG2y7iI5Uhxz/PvdrdqsG8XfGzbWPcEN8aUJfkpdotPUOPBWrgdE0/s5sxelSW5JgVfz/rHeC9cnE7xtyqceJyiINFjq2XVtYH9+o0c3GaW7aFXgxyfroa1xNQxHuwrj1+yY
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5hTVfrA8W/KJNMYGHoTEJBXwY6iIIi9rYpl7Y3dn73Xta5lV9e6q67ddRUb9oYFC7g2FAt2xIOKgNLbUKZlJsnvj3OTyQyZTGYmdzLl/TzPPMm5uffmZGaS++aU93ii0ShKKaWUUvXxZrsCSimllGrdNFhQSimlVEoaLCillFIqJQ0WlFJKKZWSBgtKKaWUSkmDBaWUUkqlpMGCUkoppVLSYEEppZRSKWmwoJRSSqmUNFhQSimlVEoaLCillFIqJQ0WlFJKKZWSBgtKKaWUSkmDBaWUUkql5M92BZRqb0Qktu77W8aY/TJ87muBa5ziaGPMzEyev6MRkUKgpzFmXrbr0hwi8h4wHqg0xuQ24fhBwK9O8TpjzLUZq1wzpXptbr7XVG3asqCU6pBE5CjgR2DXbNdFqdZOWxaUUh2OiIwDns52PZRqKzRYUEp1RL5sVyCTjDG7ZbsOqn3TbgillFJKpaTBglJKKaVS0m6ILEsYzXs5cDNwMnAKsDn27/Mj8Bzwb2NMaQvUZ1/gT8BooBdQDvwMvAHcZYxZ6dLz7gb8zymOBmYDlwF/BAYAFcBXwH+BycaYaJLTJJ5vZ+BU7CjqvkAVMB94G/u7XNjA8b2wf4c9AQG6AtXAKuBzbH/3Cw3Vo55z/wP79wb4AtjbGFPS2PM45xoI/BXYD+gOrADeB+41xnyc4rhi4AzgQGA4kA+swf6/vQE80NQ6NVDfQdSMuj8UmAJMBE4ERgCdgEXAW8A/jTG/NHC+0di/8zigH/ZvtBCYDtxtjJmb4vljHhGRR5z7mxpj5jf+laWsYyfgHOz/sgBR4CfgKeBO7O/hqWTPLyLzgYHOfjcCd2H/1h7s//ONxpjJ6cyGEJFtgAuw768B2P/ld4CbgMo0X8uO2PfFLk69AH4B3gTuMMYsSXHsFsDZ2PdUf+c1/IZ9399ljPkhnTqo7NCWhdbDCzwLPAjsiP3QzAO2A/4BzBKRAW49uYgUiMiL2Df9UdgPkyDQBdgBuBr4RUQOcqsOCXoCnwFXAMOAXKceuwNPAC+KSH0fiH4RuRf4BBv0DHaO7wRsBVwEzBWRU+t7chE5CXtB+TuwG9AH+7sowP5eDscGcK+ISKP6vkXkSjIUKGAvkN8B/4e9UAaxH8LHATNE5NZ66rA1Nhi7AXvh6AzkYH/vu2IvHr84F2I35QPTsAHgeGywE8T+zc4AZovIAfW8hqCITAI+xgYbQ7B/50Js8HOOc/wV7r6E1ERkU2yQewP2vZyP/T/aFvvl4AtsUN6QIuAD4Ajs/3IhsCU2sEqnHpc59TiJmvdUP+zv7mvg4AaOzxGRu7Dvy1Owv+MC52dr4C/AHBHZq57j/4r9Xz0TGzAVYH8XApwOfCci14qIJ53Xo1qeBgutR+ybx3LgfGAMMAGY6jwuwHsikp/pJxYRL/Yb3qHOptewAcNOwF7YYGUt9gPrJRHZM9N1qOMBbMvKV9gL387YD7U5zuOHAI/Wc+xD2AsNwEfYD8edsRf9K4Cl2AvSA05QUIuI7AFMwgZqq4G/Afs75zgcuAcIObsfhG0JSouInA9c7xSbGyiAvdh0Ap506jgGuMSpN8DFInJRnTr4gOexAVApcC32m94o7P/bk86uXYHnRCSvGfVryD+xAeBM4ATs/9sh2G+7YP9Ok5xcCImvwQtMxv5tARYA5wFjsX/na4ASbMvcDU5uipjF2Iv2KQnbrnG2bec8nhEi0hl7gR/ibHoB+z8zGjgN+418S+CWNE4Xu8jHAquDsd/G30+jHmdjWyU8wErgYuz/yj7Y91oAuL2B0/wH2yoANpA+F9u6sB/2PVGNDTpfEpHBdZ7/Wuz7yAd8iw0OxmD/Xudhfw9e7N/hGlSrpN0QrUdvbLPiOGPM7wnbp4jIv7HBxKbApWT+DXUusIdz/3RjzAN1Hp8uIv8FZjj1fEREhhhjqjJcj5je2KbwQ40xsQvzpyLyPPab6M7AkSJyvzEm1nWBiBxCzQXkRmNM3W+V74vIQ8B72G9G94jIa8aYVQn7/M25rQb2NcZ8UeccL4rIm8CrTvkI7AduSk5LRuwDOROBAtgP/7OMMfcmbPtERF7C/q16AX8XkUcTuo/GAps5908zxjxJbVNEZDE26OgHHIC9yLmhN/A4MNEYE4ltFJEp2ID1AKAH8AfgmYTjjgEOc+7PAPY3xqxPePx9EXkU+3ceBFzt/J2/cP6fvhaRLgn7LzTGfJ3RV2Zdi23pAbjMGHNzwmMzReQpbHdLOi04Xmz3W2Jw+mp9O8eISA9qAtTFwBhjzIKEXd4RkXewLWX1nWMvat5Xn2DfF4m/77dE5HNskF0IXIlt7UJEtsd2k4H9W//ZGFOdcOwM57PlNWygd7WIPKtdEq2Ptiy0Ln+uEyjEXATEMsydmsmmOudb2oVOcWqSQAEAJ8PdZU5xE2o+rN2wBjghIVCI1aEU27cdu7CcXue4S5zb2dgPrI0YY1Zgm0LBNoXGP3ydVpsg9pv5q0kChdg5XsN+cwV7QU1JRI4H7nOKmQoUAF6vEyjE6vcLNqgE20IyMeHh3gn3f6rnvHdiv0lejh2v4pYK4PzEQAHAGQfyn4RN29Q5LvZ3rgSOrnPhip1jATWv25NwTItwusliXV0f1wkUAHDqfSw2ME3HfQ3vspGjsN/4Af5SJ1CI1eMF7IW8PrHXEQFOrOf3/Si2JRBg74TPqIuw15lV2C8iG71W5339Z+xYDg/2i5FqZTRYaD3mJH5LTuR8g3/CKfbGjmnIlK2xF3+oaf6tz9SE+252RTxtjFmd7AFjzE/Y7gWA/UTED/EBe7FvaNMbGHj4ERD7wIu/DmNMmTFmR2NMN2yXUCpLndtgqp1E5HDsNy4vMIvMBQpgx7fU5xlqBq3tk7D9x4T7j4jIHk7AGGeMWWSMOdUYc5Mx5psM1TWZWfX9nbFN0zGdYndEpDc1wcPr9QTXADhN9LFvqPvUfZ0u2xPbJw/275+UM5jx7TTOV40NNBsrNuajktQtRA8n2+i8v2JplD82xqQKHg/HfpYMNMZEnYBhf+exGcaYsvoONMb8Sk03o9vdnKoJtBui9ZjRwONfJtzfHDvQKBO2S7j/LxH5V5rHDW54lyZL53exK3YMRV/s6Pdtsd9KAM4VkXPTfK6kryP2bVdECrDdP0Owv/dtsE35sQAr1QVoK2zfemwQZAhYl2a90vFpfQ8YYypEZA7297JVwvZvRGQq9kN8OHbWwCoRmY7t4nk72bdPl8xP8diGhPuJn1MjEu6nsy7GTOzr7IL9m7XUa9s24X5DF/lPqbmo12elMaaiCfXY3Lk1DRz/BTXf7BP1oSZY+5IUnAt+okFAsXP/4ISZXw3ZNM39VAvSYKH1aGhg1YqE+73r3avxujfxuOKGd2myxv4uFpLB1yEi/bHNpwdTf1AUoeGWub7ObSW2BWI0ton1zibVdGMrGng8Nk6ha53tR2MHpR2HvTh0A450fhCR77BBzt3GmA24J9W5Ey8siRewbgn3l6fxHMsS7nel5YKFxBkODU03Tud1NDXIjNVjVaqdjDGlIrKBhFacOsc3eI4kmvqe9ItIp2TdHSp7NFhoPRrqt0ycoheqd6/GS/wfOBM7gCkd5RmsQ11N+V0kvo7rSX9QXjixICL7YWcLFCRsXo9tIp2NzbEwDXiJ2t9y6/M1tt94BvbD8wYRmZLkW1hTNPRNLXaRrTWH3hizDjhBRK7GDtA8EDtoNMfZZSvs6PkzRWT3hnIdtLDGjtdJ/F+J1LtX5gUS7jcUVKbzjbvR+TzqHJfO7y3ZgOXmXCMSj30YmyMiXfV2Wajs0GCh9aj77a+uHgn360180gSJfcZrXRoV3lhN+V0kvo7yprwOpz/8aWygUIW9YD6FbcK
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
2019-10-17 17:50:31 +00:00
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdeZxcVZ3//1d19Z49ZCUECAl8CHtYgggIKq4oLoAbIODgiLuMgrgC81URmZ/D6IiooyIoIyAIDJvsu+z79kkCJAFCyEr23qrq98e5tXSnq7q6u6qrq/r9fDz6ce+te8+5p7qquz511lgqlUJEREQkn7pKF0BERESGNwULIiIiUpCCBRERESlIwYKIiIgUpGBBREREClKwICIiIgUpWBAREZGCFCyIiIhIQQoWREREpCAFCyIiIlKQggUREREpSMGCiIiIFKRgQURERApSsCAiIiIF1Ve6AFJZZraXuz9dprzT65//w93fX457VDszOwn4Y3T4aXf/awWLM+yY2dnAWdHhQe7+YInzH9R71MwOB+6MDr/j7j/NOXc2gyy7md0FHAa0u3tznmvK9jdcrJzf45/c/aRKlkXKQ8HCCGVm44B/B76M3gciVUd/wzKU9AYbuX4OfK7ShRCRAdPfsAwZBQsjV7zcN3D3WLnvIbXN3c8Gzq5wMfJy97uAsr3P3f3wAqfL/jcskqYOjiIiIlKQggUREREpSM0QI0yPHtrpx9I9me/uo9qzv/cq22gIM7sYOBFwd9/VzPYFvgscAkwElgN3AP/h7s/2kVcDcBJwLLBXlP4t4Cngb8Af3b2jjzwOAY6L7r8tMBbYALwG3A38yt1fHMDz3Aa4C9gjeugH7v6jfqRfDOwA/Je7f8PM3gV8CTgImASsAe4DLnT3Owvksz/wBULP/JlAElgJ/BP4K3Cdu6fypR+oQiMKckYKXOvuHzWz/YCvAYcDU4H1wCPA/7j734u831GE53kA4TVcTvj9/5e7P9HL9YeTZzREL9fuAPwAeD/hd7+S8N640N0fyJMm/RwzoyEG8jdsZu8DTia87lOBLcAi4Ebgl+6+Kl+5o/SzgW8A7wJ2JLy37wXOd/eHC6WV2qCaBal6ZnYs4UPraMI/wgbCB9qJwBNmdnyBtHMIQcFvgffkpJ8MHAFcBDxlZpYnfYuZ/Y3wj/NUwof6REIgPgHYE/gK8IyZHdfP5zUW+AcDDBR6ye884HbC72lboBGYBhwD3GFm/y9PujOAh4FTgJ2BZqCVEIR8CrgGuMXMRg20bINlZl8FHgQ+C2wPNBFeww8CV5vZJWZWqG9Bg5n9Bbg2SjM5ymMHwvvoITM7cRBFPBR4BvgXYEaU93aEAPN+Mzt/EHnnZWajzOxq4Gbgk2R/N+OB/YEfAi+Z2YcL5HEC8ALhfbwb4bWfSnjfPGBmp5aj7DK8KFgYeS4C5gH/l/PYvOjnlIqUaHCmAX8ifMD/nvCBfyjwE6CN8KF9iZl9sGdCM5tG+JCfC3QAvwKOBOYDH4nyTQC7Anea2fRe7n8R4cMX4FHCh8HhwDsItRXpb4z1wG+imoI+RR+8NwL7RQ8NKlAgfFCcAbwKnAa8nRAMXQikv5V+38ze1qMc7wB+SujE9xTh2+nbCTUopwDPR5ceQRjGVwlvA/6LUJPwA8Lv/jDCeyBdI3QC8IkCebwL+AzwMqF24lDgY4SaJQjvr1+b2Y4DLON5wBjgL8AHCL/D0wk1OwDfMrNvFplXUX/DZlYHXEd4HgDXE94HBxJer58A6wg1KH83s3f3vFEURKT/vjYB5xB+N+8EfgZ0Ef5upMapGWKEcfflwHIzW5Pz2JMVLNJgjYu2x7v7X3Iev8/MbiA0RTQB/2Vmt7p7Z841FxGCjfXAEe7+SI+8r4tqDa4DpgP/SfgmDYCZ7Uz4EAJ4CDjM3dtz0t8L/MnMriR8CxtF+KD4c6EnZGbNhG+4B0cPDTZQgPA8nwMOdfe1OY/fbmZLCB9mEL5F504edDIhUFhNeH7rcs7dHz23x4HZwL+Y2RnunhhkWftrKvAG8DZ3X5rz+D1m9hRweXR8Us5+bx4APuDu63Meu8bMLiI0TbQQmqoGUgsQA77s7hfmPPZPM/s7cH/0HP6fmf2pryaBfvwNf40QBAGc6u6/6XH+djP7fXT/acAfzWx2+m/EzJqAC6KybyS8d3Lvc5eZXUuorep1wiipHapZkFpwSY9AAYCoHTj9j30O4dsUAGa2C3BUdPjjXgKFdB7XE75ZARxrZtvmnN4DWEiowfhJj0AhV27ZZhR6IlH/ib8B6W95pQgU0s7oESik/Y5s7cLePc5Ni7YregQKAEQfrN8njPn/HpX70PhRj0Ah7Uog/Zx7PreevtAjUEjL/ea8Ry/ni3FDj0ABAHd/Cfh2dNhCCGgGLapV+Lfo8KZeAoX0/V8GzowOZwIfzzn9HmCnaP+nvQUk0d/Yz0pRZhneFCxILShUDfqHnP3cdtkPkh0ff2sf+d8YbesITQwAuPvf3d0Ibbj/10u6tOU5+00FrosDlxGaQqC0gUIXoZZlK1EAkf5AHdPjdLpT5lwz+13USa9n+r+6+zfd/VfuvqlE5e2vW3p7MOp0+Up02PO55Xq1QEfYhTn7EwdQNgh9YvK5HEgHmu8dYP497UX48Ie+39835eznNkXkNt1dViD9HwqckxqhZgipdm3AY/lOuvsrZraW0Nlw15xT83L2H8/Tf7E3O/V8IHcUgJlNja6ZQ+gMtj+hfTqtUID+I0JP87TVxRaqCG+6e1uB8xvJdszM9StCP4wxhPbwU8zsOcIH0G3AXRUMEHItLnBuY7Qt9P/utXwn3L0t5/0x0P+ZD/WR/wvAPoQOsaWQ+/7+uZn9vMh0ue/v9N/LOnd/pbeLAdx9iZmtAKb0s4xSRRQsSLVbUUQb+SpCsDAt57FJA7zfhJ4PRMPSvkToVDduqxRhmGExdoy27YQaiPPM7Hp3f3UA5expYx/n0wFPtxED7r4oen5/IPvhsXv08w2g3cxuBS5y9xtKUM6BaHf3rgLne31uPWwo8l4Dna1xZR/n0/0UBlpz0VMp3t9To+2a3i7s4U0ULNQ0BQtS7Qp9SKSlp8XNnSsh970/H8jt+FhI5p9+NBTvt2w9iuRVQvX904ROc28ROoEV49+Blwj9JMYAv6F7dfCQc/d/mtnuhCaYjxPmCZgdnW4CPgR8yMz+DJzo7sUGR8NJyeeI6Gf+6SAkX7+X/sp9f3+JMLS4GFty9osJstKK/fuRKqVgQapdMd/E0t+y3sh5LPfb0uvuvmwA9z6VbKCwkDBm/RZ37/ZNLJq4pxgXuPtZUZoTCT3ZP2BmJ7j7pQMoX8lEAcAd0Q/REMJ3E/qBHEn4X3I8oe9ARcs6TE2kcLNS+lt5qZqect+D6wY44ukNQg3SJDOL9THpVqlqRGSYUrAg1W68me3g7kt6OxlNpjQ2Onwq51RuZ7a3AVfnu4GZHUj4Vr0YuN/d0+3bX462CeB9Bdp1ty/0BHLktmufSqiZaAb+08z+4e4risynZMxsDKH5YXXUcx4Ad19MmNfi99Gsh9dGpz6EgoXe7A7c09uJ6Hc8Nzp8vET36/n+zttB0cwmE2ofFgNPunv67+Q5wgii0YT3wAt50o8nTF4lNUyjIUauaqwqzueEAudyl/DNDQj+kbP/xT7yP58wMdFf6d4BbE60XVuoAxhhlr60ogJ0d19I6PAIsA3w38WkKyUzm0mYg+JhCq/8+A+y7yeNt+/dJwucO5Hs++LmfuRZ6G/4EbK1C58xs9760qR9lfD6XkyYjCztqpz9fymQ/gTKuPKmDA8KFkauTNuomY2uZEFK4DvRugDdmNk7CbMVQqgRyMxh7+6Pkv2md4SZfbe3jKNZ9Q6NDp8kTLSUlu6UNilaO6Fn2piZnUP34XCFhk729DOyMyQea2Yf7UfaQYs6Vqa/ZX7CzA7Kc+mnyP4v6XW+CuELZrbV+ihmthfZoHAFfUzY1UPev+FoLZN0gLkNcGk0yVLP+x9MmEkSQn+F3+Wcvp8wKynA183sCHqI+rJUauZOGUJqhqg
2019-10-17 17:50:31 +00:00
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
2019-10-17 17:50:31 +00:00
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd3yb1fX48Y+8V+xMkhAgQEIOGYxQIGWHvTcp/BgFWmgpUFpaoJRSRlugQDeFwpdZoLSsJFBC2DtsKAFCOAmEBEjIThzvIT2/P+6j5diyLUuWZJ/36+WX9KyrK1uWju44N+B5HsYYY4wxHcnLdAWMMcYYk90sWDDGGGNMQhYsGGOMMSYhCxaMMcYYk5AFC8YYY4xJyIIFY4wxxiRkwYIxxhhjErJgwRhjjDEJWbBgjDHGmIQsWDDGGGNMQhYsGGOMMSYhCxaMMcYYk5AFC8YYY4xJyIIFY4wxxiRkwUI/IiLbZ7oOpnMicoaIeP7PSZmuT7YRkatifj/fTkP54bKfSvL6qTFlXNrmWI/rLiIv+dc3Jjgn4//rMc/znkzXxfRcQaYrYNJPRKqA3wDnYX9zY/os+1836WIvpv7hT8D3Ml0JY0za2f+6SQsLFvqH/ExXwJhUUdWrgKsyXI0OqepLQCCN5U9NcNj+101a2JgFY4wxxiRkwYIxxhhjErJuiD5MRK4Crmyzz/PvvtxJc2ayj1kG/Ag4BtgWqATWAO8CDwAPqWoo1Y/rP/Y9wOmAquq2IrITcBmwJzAYWA68APxBVT/upKxC4AxgGrC9f/16YC7wCHC3qjZ3UsaewCn+42+K+13UAF8DLwM3q+qnSTzPIcBLwCR/169V9XfduH4xMBr4q6r+VET2A84FdgOGAmuB14BbVPXFBOXsDPwQ2AfYHAgBq4A3gP8Aj6uq19H1yWrzut5NVd+MOfaSX5/HVPUYEfkWcAEwFRgObADeAe5Q1RldfLyjcM9zF9zfcDnu9/9XVf1fO+dPBcK/t1+q6u8TlD0a+DVwCO53vwr32rhFVV/v4Jrwc2xS1RJ/31V0839dRA4GzsT93YcDDcBnwJPATaq6uqN6+9ePAX4K7AdsiXttvwrcqKpvJ7rW5B5rWTApIyK7AAuAP+A+IIcCRcBI4Ejg38CrIjKyF+oyDfehdTzujbAQ94F2OvA/ETk1wbVjcUHB/wEHxlw/DDgAuBWYKyLSwfWlIvII7o3zHNyH+mBccD4I2A44H/hIRE7p5vOqBJ4myUChnfKuB57H/Z42xf29RgAnAC+IyG87uO4S4G3gLGAboAQowwUhJwEzgWdEpDzZuvWUiPwYeBP4LrAFUIz7Gx4GTBeRe0Uk0diCQhH5F/CYf80wv4zRuNfRWyJyeg+quBfwEfB9YJRf9ma4AHOOiNzYg7I7JCLlIjIdeAo4kejvZiCwM3AF8LmIHJmgjNOA+bjX8QTc33447nXzuoick466m8yxYKFvuxWYDPw3Zt9k/+esVD6QiEzEfZsahfuGcT1wEDAF+A7uwwNgd9L/ITIC+CfuA/5O3Af+XsC1QCPuQ/teETmsnecxAvchPx5oBm4GDgd2BY72yw3iWk1e7CDwuRX34QuuReX7uG+2e+NaK8LfGAuA2/yWgk75v7MngW/5u3oUKOA+KC4BvgIuxP1tDgBuAcLfSi9vmw9ARPYGfo8bxDcX9+10d1yAeBbwiX/qAbhpfJnwbeCvuJaEX+N+9/vgXgPhFqHTcK/NjuwHnAwswrVO7AUci2tZAvf6+oeIbJlkHa8HBgD/Ag7F/Q4vxrXsAFwkIj/vYlld+l8XkTzgcdzzAHgC9zqYgvt7XQtU41pQZojI/m0fyA8iwv9fdcDVuN/NvsANQCvu/8b0IdYN0Yep6nJguYisjdn3QZoe7n6gHNfEPlVVP4859jbwsIici3sTmYR7A790o1JSo8q/PVVV/xWz/zURmYXriigG/ioiz6pqS8w5t+KCjQ3AAar6TpuyH/dbDR7HtZj8GfdNGgAR2Qb3IQTwFrCPqjbFXP8q8E8ReRj3Lawc90Fxf6InJCIluG+4e/i7ehoogHue84C9VHVdzP7nRWQJ7sMM3LfoN2OOn4kLFNbgnl91zLE5/nN7HxgDfF9ELlHVYA/r2l3DgW+Ab6vqlzH7XxGRucCD/vYZMffb8zpwqKpuiNk3U0RuxXVNlOK6qpJpBQgA56nqLTH73hCRGcAc/zn8VkT+2VmXQDf+1y/ABUEA56jqbW2OPy8id/qPPwK4W0TGhP9HRKQY+Itf91rcayf2cV4SkcdwrVUliepscou1LJgeE5EDgR39zZ+3CRQi/DfFl/zNH/njAtLl3jaBQrgOrxN9Yx+L+zYFgIiMA47yN69pJ1AIl/EE7psVwDQR2TTm8CRgIa4F49o2gUKs2LqNSvRE/N/TI0D4W14qAoWwS9oECmG3E21d2KHNsRH+7co2gQIA/gfr5bg5/78icx8av2sTKIQ9DISfc9vn1tYP2wQKYbHfnCe1c7wrZrUJFADw/39+4W+W4gKaHvNbFX7mb85uJ1AIP/4iooH85sBxMYcPBLb27/++vYDE/x+7IRV1NtnDggWTCofH3H+2k3Of9G8rcQPG0iVRM+hdMfdj+2UPIzo/vqvPIw/XxQCAqs5QVcH14f63nevClsfcL05wXj5uYGj4d5zKQKEV18qyET+ACH+gDmhzODwoc7yI3O4P0mt7/X9U9eeqerOq1qWovt31THs7/UGXX/ibbZ9brK8SDIRdGHN/cBJ1AzcmpiMPAuFA86Aky29re9yHP3T++p4dcz+2KyK26+6BBNffleCYyUHWDWFSYXLM/bUdjPtrz9ZE++9TqRF4r6ODqvqFiKzDDTbcNuZQ7PN4v5vPo+1jRGYBiMhw/5yxuMFgO+P6p8MSBe2/w400D1vT1Up1wQpV7XB9AVwzc3hgZqybceMwBuD6w88SkXm4D6DngJcyGCDEWpzgWK1/m+g98OuODqhqY8zrI9n30bc6KX8+rsVuuyTLbyv29f0nEflTF6+LfX2H/1+qVfWL9k4GUNUlIrIS2KSbdTRZyoIFkwpDk7xuUEprEbWyC33kq/3HHxGzL2XPw5+Wdi5uUF3VRle4aYZdsaV/24RrgbheRJ5Q1a+SqGdbtZ0cDwc8cTMGVPUz//ndRfTDY6L/81OgSUSeBW5V1VkpqGcymlS1NcHxdp9bGzVdfKxkszWu6uR4eJxCsi0XbaXi9T3cv13b3oltrMCChT7DggWTCuHX0Wpcn2ZXLU1DXcA1r3cmnBY3NldC7P/DrkDswMdEIm/6/lS8/2Pj2SZf4ZrvP8S1pqzHDQLrit8An+PGSQwAbiO+ObjXqeob/gyYqbg+7UNwAxrBBTVHAEeIyP3A6enKrZFmKc8R0c3yw0FIR+Neuiv29X0ubmpxVzTE3O9KkBXW1f8fkwMsWDCpEP6WUQF8mAUfDF35Jhb+lvVNzL7Yb0tLVXVZEo99DtFAYSFuzvozqhr3TcxP3NMVf1HVK/1rTseNZD9URE5T1fuSqF/K+H/nF/wf/CmE++PGgRyOe385FTd2IKN1zVKDSdytFP5Wnqqup9jXYHWSM6O+wbUgDRWRQCdJt1LVImKygAULJhU+xs1rL8H1sb7f0YkicgjuzWYxLrNcwilhSRooIqNVdUkHdRDcAEtweQLCYgezfRuY3tEDiMgU3LfqxcAcVQ33b5/n3waBgxP0626R6AnEiO3XPgfXMlEC/FlEnlbVlV0sJ2VEZACu+2GNP3IeAFVdjMtrcaef9fAx/9ARWLDQnonAK+0d8H/H4/3NDv+fuqnt67vDAYoiMgzX+rAY+EBVw/8n83AziCpwr4H5HVw/EJe8yvQRNhuif0j3N/2nY+6f29FJIpKPS/jzB9xUwHQmZjotwbHYJXxjA4LY5/GjTsq/EZeY6D/EDwAb69+uSzQADJelL6xLQbuqLsQNeAQYAvy9K9elkohsjstB8TaJV358mujrzubbt+/EBMdOJ/q6eKobZSb6X3+HaOvCySLS3liasB/j/r734JKRhT0ac//7Ca4/jTS
2019-10-17 17:50:31 +00:00
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
2019-10-17 17:50:31 +00:00
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5yU1dXA8d/MttldytKbVIEjYscSFQVBxYIi1tg1sUaNJbHG1/ImdpOYWKK+9hZ7wW7EiKKIDRQpF5TeYRfYPjvlef+4z8zOLrOzbWZny/l+Pnxm71PvLLMzZ2451+M4DkoppZRSdfGmuwJKKaWUat00WFBKKaVUQhosKKWUUiohDRaUUkoplZAGC0oppZRKSIMFpZRSSiWkwYJSSimlEtJgQSmllFIJabCglFJKqYQ0WFBKKaVUQhosKKWUUiohDRaUUkoplZAGC0oppZRKSIMFpZRSSiWUme4KqI5DRJ4CznaL/Ywx69NYnVZLRCLrxn9ojDkirZVpZURkCLDMLT5ijLkoydd/ima+RkVkOTAYMMaYnWK2D6GZdReRc4An3eKpxpgX4xzTCehtjFna2OsnS+zv0RjjSVc9VPJoy4JSSrUTInIKsAg4ON11Ue2LtiwopVQ7ICIHAdu1NCiVDBosqBZjjDkHOCfN1VBtmDFmOdCqm7WNMUNSeO2ngKfq2J2Rqvsqpd0QSimllEpIgwWllFJKJaTdEB1EzAj764G7gPOA84GdsK+DRcArwD+NMWUpqsNTpGg2hIiMB/7rFvcH5gPXAScCg4BKYA7wOPCCMcaJc5nY6/0KuAAYB/QHAsBy4CPs72hlPef3wf5+JwICdAeCQCHwDbZv+bX66lHHtW/H/j8CfAscZozZ2sBzbwFuBrYZYwpEpC9wFXAM9vcUABYCLwH/Msb467hON+BiYDKwM5AHbMG+jt7DjvZvUJ0aI9GMglozBboBYeAK4HhgGLb74mfgdeAfxpjiBt7vauAIYABQDCzAdgU8He//r67ZEHGOywQuw/5NjMS+Rn8C/g08ZowJxDkn9jmeaox5sdbvJOJJEYkcN9Ttvom9zkD33pOAIUA2sA74HHjYGDOrrnq75+diuxTPBEYAPuz//ePAo4nOVW2Ttix0PF7gZewf9D5AZyAX2BO4HfhORAalr3pJ0Rv4GrgB+ybsAwqAQ4DngNdFxBfvRBHJFJGHgFnAudgPGR/297Qr8AdgsYhcUNfNReRs7Jv3n4HxQD8gB8jHfiCfgA3M3hKRRvUzi8ifaGKgEOdaY7FB1dXYoDEP6Ar8Cvg78LWIdI9z3m7uebdhA7OuQBb2934wcCfwi4js35R6JYkA84Bbgd2x/3+dgD2A/wV+cj9kEzkBGzj9Dvs6yAF6YQPIJ4F33Q/8pugEfAr8za1fLjbAOQh4CPt32K+J105IRH4LLMb+v+8GdMG+xocCZwFfisgjIpJVx/mDge/deu4P9HSfz97Av7ABddy/L9V2abDQ8VyG/ba9Efut6wBgCvC+u1+AT0UkLz3VS4pHsB9+c4DTsR9+52Df+AGOA56u49zHsN+YAWZiv/X9CvuhfwOwHvuh8YgbFNQgIhOw3zpzgSLsB9OR7jVOAB4EqtzDj8G28DSIiFwB/MUtNitQcOs3DftB8S/sN+cDgEuBte4xu2EDgtg6ZACvYgOgMuAWbOvJvtjX0fPuod2BV9xvoOnwJjYw+zdwLPb3fy72QxJgIPZ5J/IAtjXiQezv5xDgT9jnDfb/9aom1u804EBsK8W5bv1OBr5w9+9Kw4ORtdhg//yYbTe72/ak+v8z0jLxGPbDfBk2+D0I+39/PjDXPfQC4rQQiEg+MAP79wXwFtW/3/Owv9+JwEkNqLdqQ7QbouPpi21OP8gYszpm+zQR+Sc2mBgKXIt9w2mL+mKbwqcaYyIfzLNF5FXgY9w3ZhF52BgT6bpARI6jupvkDmPMDbWuO0NEHsN+I9wZeFBE3jHGFMYc87/uYxCYZIz5ttY1XheRD4C33fJJ2OAmIbcl4+9usbmBAthm5wxgsjHmg5jts0TkPWxzeB5wmoj8PqZJfCy22RngQmPM89Q0TUTWYr+1DgCOAl5rRj2bqi9wvjHmsZhts0XkTewHdD9gkoj0TdAd5gcmGGO+jNn2qYh8gm15Avt6ubsJ9fNgP3SPNMZUxNTvNWyAczL2g/5CbLBSJ/c1PldECmI2rzTGzI09TkT6x1xrOjClVpfjLLer8BngVOAcEXmp1uvjBmwXC8BtxpgbY/bNFpEXgf9gWxxUO6ItCx3Tb2oFChF/ACJZ3y4QkVY9RS2BLcCZMYECAO4b41nYvmyA2hn0rnYf52O/QW7HGLMJ2ywNtlsh2jLgtsbkYFsU3o4TKESu8Q4Q+aAfUN+TEZEzqP4WnIxAIeL1Wh8Ekfotw77hg215GBqzu2/Mz0vquO4/gP/Ddpf8nIR6NsU3tQIFANzf28tu0YNtPanLA7UChcg1vsI2wwPs1MSuCD9wekygELl2GPsNf5u76eLaJzbDJdgAMAicFW9skjEm6N4zcv/LI/tExEt164UhzpeJOH9jqp3QYKHjWRj7bTqW++3xObfYFzumoS160RhTFG+HMWYJtnsB4IjIG707YC/ybWh6PQMPZwIl7s8TY65dbozZxxjTA9vVk0jk22xOooNE5ARst4YX+I7kBQpg+5br8kvMz51jfl4U8/OTIjLB/RCJMsasMcZcYIy50xjzQzIq2gRNeW61vZdgXyRQ8mLHbDTWe8aYNfF2uAMv33SLo90WgWQ42n1cYIxZW9dBxphtVHeHHBwzdmFv7JgNgJeMMaE6zv8Z+CwJ9VWtiHZDdDxf1LP/+5ifd8IOFGxrGvIcD8Z+a+4PrMQOfIu0pPxeRH7fwHsNi7fR/YYY6eMdCuyI/X3ujm3KH+gemihg3xV4gepkO1XY0fjJsjzBvtKYn6PvE8aYH0TkfWx//c7Y5uxCEZmO7eL5yBizIol1bKrlCfbFfW5xxGt9i6hs4DXqMrue/XOo7hLblZhxB03hBsW7usXdYmZH1ScP6IP9XcTO7Jgb//Cor7HjfFQ7ocFCx1Pfm86mmJ/71nlU69bY57gSO6K7KbrV3iAiO2C7dI6ljmAC20xbX8te5BulH9sCsT92TMk/mlTT7ZUm2Bf7YVK7O+rX2L7v0919PbB97CcDiMg8bJDzgDEm0T1SqanPLVZJgn2xmtJdt7Ge/Ztjft5uRkoTdKPpLcndsMFCn5hthXUcG7GhifdSrZQGCx1PsJ79sVP5quo8qnVrynOM/Vv4Cw0flFejKVZEjsDOFsiP2VyCnYkxH5tj4WPgDWB0A64/FzgF21rSE7hNRKa54wrSwm0mP1NEbsIO0JyMHTQaaa7eFbgD+J2IHGKM+SX+lVq9RufASKLYACRurotGin19f0D19NuGiIw7aWiQBTZfh2pHNFjoeOr7ltIr5ud1qaxICjXlOcaOcaioPZK8IdwERy9iA4UA9gPz39jkPE6tYzs14JJLgYnGmCIRuQo7Sj0fO3jw0MbWL9ncgOVu4G73+RyEnWJ4CvZb6EDs9LuJdV6k46rvNdo75uf6vsU3ROzrO6Mpr29qvh/0qvMoKxmtIaoV0QGOHc8e9eyPHdSYrsFpzdXQ57iO6ubgn2L2/yrRySKSLSI3i8hvayUeOoPqwW5/McbcbIxZFCdQyKFmk25dlkQGahpjnsWODwCY6CbWaXFu0qqRInJg7HZjTKkx5n1jzOXAKKpn1UxIY66F1qy+VqV93ccw9Y8PqJebiTPSQjCmvhkcInKRiFwiIkfHHDs/5pD6Bj/X9zeo2hgNFjqeg9wMbNtxP8ROc4sLjTGL4h3XBpxaV2ZEEdkZ2M8tvhn5IHdHpkfeDI8QkaHxznedgU1G9Bg1E+EMj/n5uwTnn0B1hrvGtO5dRPXAuntTleGvHh9gp8197A7e3I4xZgvVeQhAs/nFM7muxGci0hubsAvgK3d2QkPUN13xQ/exO3bcSVwiMhw7JuUB4H53OiVua0QkCDwtQRbUXtgWJtWOaLDQ8WQAT9X+Q3d
2019-10-17 17:50:31 +00:00
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
2019-10-17 17:50:31 +00:00
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd3zTdf7A8VeSTsoqs1A2yNvDLSoOUNw4ce99DjzHieP0PH+n3um5PT2358AtDlSOoZ56bnFPxA8iArJXodCVJvn+/vh8821a0nQlTVPez8ejj+TzXXknbZN3PtPnOA5KKaWUUvXxpzsApZRSSrVtmiwopZRSKiFNFpRSSimVkCYLSimllEpIkwWllFJKJaTJglJKKaUS0mRBKaWUUglpsqCUUkqphDRZUEoppVRCmiwopZRSKiFNFpRSSimVkCYLSimllEpIkwWllFJKJaTJglJKKaUSykp3AEqpholIdC35N4wx49IaTBsjIoOAX93iQ8aYCUm+/iTgdLfYxxizvBnXWAAMBIwxZsuY7YNoYewicgbwuFs80RjzfJxjOgK9jDHzm3r9ZIl9HY0xvnTFoZpHaxaUUqodE5HjgZ+APdMdi8pcWrOglFLtlIiMATapaVCqqTRZUEplNGPMAqBNV2sbYwal8NqTgEn17A6k6nHV5kWbIZRSSimVkCYLSimllEpImyHaoZie838GbgHOBs4BtsT+zn8CXgT+ZYwpa4V4DgTOBHYDegMVwDxgBnCPMWZ1ih53LPA/t7gbMBu4CjgGGABUAl8DjwLPGmOcOJeJvd6uwLnAXkBfoBpYALyJfS0XNXB+b+zvYV9AgG5ACFgDfI5tW365oTjqufY/sL9vgC+A/Y0x6xp57nXAtcB6Y0xXESkCLgUOw75O1cAcYDLwgDGmqp7rFALnA4cCI4AOQAn2720Gtrd/o2JqikQjCuqMFCgEIsAlwFHAEGzzxTxgCnC3Maa0kY93BTAOKAZKgR+xTQFPxPv91TcaIs5xWcBF2FEDw7F/oz8AzwGPGGOq45wT+xxPNMY8X+c1iXpcRKLHDXabb2Kv09997AOBQUAOsAz4AHjQGPNJfXG75+cDZwCnAlsAedjf/aPAw4nOVW2f1iy0b37gBew/6s5AJyAf2AH4B/CliAxI1YOLSIGITAFeB47HfvDkAl2BnYC/Ar+IyGGpiiFGL+Az4Grsm3CeG8fewNPAFBHJi3eiiGSJyP3AJ9ikZ4h7fidgG+AyYK6InFvfg4vI6dg3778DY4E+2NeiAPu6HI1N4F4TkSa1M4vIX2hmohDnWqOxSdUV2OSyA9AF2BX4J/CZiHSLc9627nk3YhOzLkA29nXfE7gZ+7verTlxJYkA3wPXA9thf38dge2BvwE/uB+yiRyNTZz+gP07yAV6YhPIx4Hp7gd+c3QE3gXudOPLxyY4Y4D7sf+vfZp57YRE5PfAXOzvfVugM/ZvfDBwGvCxiDwkItn1nD8Q+MqNczegh/t8dgIewCbUcf+/VGbQZKF9uwj7LXol9tvU7sB4YKa7X4B3RaRDsh9YRPzAVOBId9M0bMIwCtgPm6ysx74pvSIi+yY7hjoewn74fQ2cjP3wOwP7xg9wBPBEPec+gv3GDPAh9lvfrtgP/auB5dgPjYfcpKAWEdkH+60zH1iL/WA6yL3G0cB9QNA9/DBsTVCjiMglwA1usUWJghvfVOzv5AHsN+fdgQuBpe4x22ITgtgYAsBL2ASoDLgOW3uyC/bv7Rn30G7Ai+430HR4FZuYPQccjn39z8R+SAL0xz7vRO7F1kbch3199gb+gn3eYH+vlzYzvpOAPbC1FGe68R0HfOTu34bGJyNLsV8KzonZdq27bQdqfp/RmolHsB/mv2KT3zHY3/05wDfuoecSp4ZARAqA97D/XwCvUfP6no19ffcFjm1E3KqN0maI9q0IW00+xhizOGb7VBH5FzaZGAxciX0jSaaLgX3c+xOMMQ/V2f+2iDyKfSMswlaRDo1XzZokRdiq8CONMdEP5k9F5CXgLdw3ZhF50BgTbbpARI6gZkKem4wxV9e57nsi8gj2G+EI4D4RmWaMWRNzzN/c2xBwoDHmizrXmCIirwP/ccvHYpObhNyajH+6xZYmCmCrnQPAocaY12O2fyIiM7DV4R2Ak0Tk4pjf1WhstTPAecaYZ6htqogsxX5rLQYOBl5uQZzNVQScY4x5JGbbpyLyKvYDug9woIgUJZh4qQrYxxjzccy2d0XkHWzNE9i/l1ubEZ8P+6F7kDGmIia+l7EJznHYD/rzsMlKvdy/8W9EpGvM5kXGmG9ijxORvjHXehsYX6dp8hN3MqUngROBM0Rkcp2/j6uxTSwANxpjronZ96mIPA/8F1vjoDKU1iy0f2fVSRSiLgOis7mdKyJJG3rm1ipEv13NjJMoAODOJneVW+yPbUdOlRLg1JhEIRpDGbaaNeJuqjuD3hXu7WzsN8hNGGNWYaulwTYreDUDbq1NLrZG4T9xEoXoNaYB0Q/64oaejIicQs234GQkClFT6nwQROP7FfuGD7bmYXDM7qKY+z/Xc927gX9jm0vmJSHO5vi8TqIAgPu6veAWfdjak/rcWydRiF5jFrYaHmDLZjZFVAEnxyQK0WtHsN/w17ubzq97YgtcgE0AQ8Bp8fowGWNC7mNGH/+P0X3u/3q09sIQ50tHnP8xlYE0WWjf5sR+S47lfit82i0WYfs0JMu22A9/qPmAqc/MmPupbIp43hizNt4OY8zP2OYFgHHRN3q3w17029DbDXQ8/BDY4N73nocxptwYs7Mxpju2SSiR6LfZ3EQHicjR2GYNP/AlyUsUwLYt1+eXmPudYu7/FHP/cRHZx/0Q8RhjlhhjzjXG3GyM+TYZgTZDc55bXTMS7IsmSn5sn42mmmGMWRJvh9vx8lW3uJVbI5AMh7i3PxpjltZ3kDFmPTXNIXvG9F3YCdtnA2CyMSZcz/nzgPeTEK9KE22GaN8+amD/VzH3t8R2AEyGHWLu3ykidzbyvCFJevx4GvNa7In91twXWITt+BatcblYRC5u5GPFfR7uN8RoG+9gYCj2dd8OW5UfTbASJfHbAM9SM9lOENsbP1kWJNi3Mea+995hjPlWRGZi2+tHYKuz14jI29gmnjeNMQuTGGNzLUiwL+5ziyNeLV1UZSOvUZ9PG9j/NTVNYtsQ0++gOdykeBu3uG3MKKqGdMCOalpMTT8FqOnbUJ/PsP18VAbSZKF9a+jNZFXM/aJ6j2q6Hs08rzCJMdTV1NdiEUl8HiLSD9v0czj1J0URGq7ti36jrMLWQOyG7Xtyd7Mi3dTGBPtiP0zqNludgG37Ptnd1x3bxn4cgIh8j01y7jXGJHqMVGruc4u1IcG+WM1p1lvZwP7YIcabjEhphkKaX7tciE0WesdsW1PPsVErmvlYqg3QZKF9CzWwP3aIXrDeo5ou9u/qD9R0/GpIRcOHNFtzXovY53EDje+UV6sqVkTGYUcLFMRs3oAdiTEbO8fCW8ArwFaNuP432JElH2ETmhtFZKrbryAt3GryU0Xkr9gOmodiO41Gq6u3AW4C/iAiextjfol/pTavyXNgJFFsAhJ3rosmiv37fp2a4beNEe130tgkC+x8HSpDabLQvjX07aNnzP1lSXzc2L4B6+v2wE6T5rwWsc+jojnPw53g6HlsolCN/cB8Djs5j1Pn2I6NuOR8YF9jzFoRuRTbS70A23lwv6bGl2xuwnIrcKv7fMZghxgej/0W2h87/C7VQ2UzUUN/o71i7jf0Lb4xYv++A838P4193+hZ71FWMmpDVJpostC+bd/A/thOjcnsdPZDzP1dsdXPcYlIT2ztwwLgmxR2ftseO/67PtHXYhk11cF1n0e9RCQH+81sMbazWLQ25RRqOrvdYIz5Wz3n51K7Src+P0c7ahpjnnLnddgX2FdEfm+MebQR10gqt+17CNDTGOP1DXGbG2YCM91ZIr9wj9tHRPLr9vpXDdYq7eLeRmi4f0CDjDFVIjIPGAaMFJEsd+RDXCIyAVsDtwB4wz12dswhO2Nr0OrT0Pu
2019-10-17 17:50:31 +00:00
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
2019-10-17 17:50:31 +00:00
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgUAAAFGCAYAAAD6uOxSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3dd5xU1f3/8dciSNEIliiKiJT4UTSAXWMvib3EkqigsaAxml/UkED8GiOaqLEkYgWxRGMv0dhjFEvEgjUilo8VsERRiopI3f39cc6wwzAzOztz787s8n4+Hvu4O3PvPXPmzty5n3tqXUNDAyIiIiLtqp0BERERqQ0KCkRERARQUCAiIiKRggIREREBFBSIiIhIpKBAREREAAUFIiIiEikoEBEREUBBgYiIiEQKCkRERARQUCAiIiKRggIREREBFBSIiIhIpKBAREREAGhf7Qy0RWY2wN0nppR2Zq7rh9199zReo7UzsyOBv8WHh7r7rVXMjjShtZ4vZjYSOCM+3NrdnysjjSeAHYB57t4pZ11FeTezHYHH48NT3f3PebZZHujj7m81N/2k5BzH3u4+uVp5EZUUJMrMuprZxcDL1c6LSK3T+VJdZrYrMBE4pNp5kdqhkoJk/RU4utqZEGkldL5UiZn1BB6pdj6k9igoSNZyab+Au9el/RoiLWSZP1/cfccU034CKPT+Uz/20jqp+kBEREQABQUiIiISqfogATmtZzPPZVoOP5lkEWHKramvA34GuLuvb2abAP8HbAusAnwKPAZc6O6TmkirA3AkcDAwIO4/C3gVuBP4m7vPbyKNbYHB8fXXAlYCvgY+Ap4ELi+n1bSZrQo8AWwUnzrd3f/UjP0nA72Ai939ZDPbGTgB2BpYDZgBjAeucPfHi6SzGfBzQuvznkA98DnwLHArcK+7NxTav1xZ39cv3b2bmXUHfg3sA6wDLADeBG4DRrv7vCJpdQSOAA4CBgErEz7n14F/AmPd/dsCr5/9XIueLzkt8zcG3gBOBH4KGNARmAzcD/zF3T8r4bW+C/yGcBzXBeYCbwO3EI7jUt/3Yr0PcrarI5xPxxG+t/Ux7TuBy9z9mzz7ZL/Hxb0Pso5Jxhlmlvk8dorVDtnprAL8P2AvoB/QBZgGPEc4jx8qlO+4f3vgJ8BQYH2gG/AecBNwUbF9peWppEDyMrODCRenA4E1gA6EC9fPgFfMbEiRffsRLv5jgR9m7f9dYFdgDPCqmVmB/Tub2Z3AU8DxhB/BVQhB7MrA94FfAq+Z2eBmvq+VgIcpMyDIk955wDjCcVoLWB7oTrhIPmZmfyyw33DgecIP5feAToQf216E1uD/BP5tZiuUm7cS878t4QL+W8IPdhegK7AV4Qf7+XhRyLfvABo/5x8Bq9P4Oe8IjAI8Bj+1bHXgRULDxy0JF63OwAaE4/KWmW3aRBrbAg4Mj/t1JnxXtyQch6fjd68c7YHbgWsJn8uKhAB5M+DPwJtmtkGZaRdlZnsSLuAjgc0J76kj4bfgYOBBM7vHzL5TYP9VCAH8TcBOwJqEY7MRcC7hHFg1jbxLeRQUJGMM4W7jvqznNo5/Q6uSo8p0B64n/MBfQ7iwbwecQ7j7aQ/8Pf5gLCHedT5F+GGcD1xOuMPYAtgvpruIcAF63MzWzPP6YwgXWQg/1scQLjLbE+6Wnonr2gNXxjv/JsUL7INA5ge+ooCAcFc5HPgQOAX4ASHouQLI3I393sy2ysnH9oQf8zrCRfWouO+2hO/LG3HTXYGzKshfUzoD9xIuMKOB3WM+fgl8ErcZAJydu6OZfQ94lHBXDeGO+mDCRfDHhIsYhIvHE2a2UdbutXa+XE8INP9F43s4BHghru9G+L4Xa7R4QdzuBkJJwfbAScAXcX3mAl6OnQlB5keE0oytCefS/XF9T+ARM+taYnobE87JjCtpPP4vZp40s10I349uhJKBPxAu7FsBhxNK2wD2Be42syUaL5pZO+AhwncK4D+EEoOtgCGE4zuAUMomNULVBwlw90+BT81sRtZz/61iliqV+XEZ4u43ZT0/3sweIFQhdAQuNrNH3H1B1jZjCEHFV8Cu7v4CS7o3lgLcS7hruIisftLxYnN4fDgB2CGn+Pop4Hozu4PwQ7kCsAdwY7E3ZGadgHuAbeJTlQYEEN7n68B27j4z6/lxZjYFOC8+/hmhqDXjKEJAMJ3w/r7MWvd0fG8vA32BY8xsuLsvqjCv+SxPaIW+t7v/K+v5Z83sQWASoeTgMDP7Vc7nfAWhRABgmLv/NSftf8Y0riN8RjeZ2SB3b6jB86U7cLa7/z7ruefN7G7CnexAoD8hsJ1QJJ2fuPudWY+fMrP7CN+RzsAQM/ulu9c3M391hEBxB3f/Iuv5e83sAkKVRQ/g94SSjaLc/b9mNivrqU9zj7+ZdSacU8sBrwG7uPvnWZtMAG6MJWXDgV2AYwnnf8YxhGMGIVj6WVZ12AQzu40QPP64qTxLy1FJgRTy95yAAAB3f4ZwVwShfnHXzDozW49w1wDhRzY3IMikcT/h7gzgYDNbK2v1RsA7hBKJc4rUZ2fnrUexNxLbN9xJ+OGCZAKCjOE5AUHGVTSWFgzMWdc9LqflBAQAuPtXhB/4vwKnEaoW0nJXTkCQycMHNPZjXwnonVlnZoNo/NwfyBMQZNK4nhAUQLgjrNUROP9HKB5fQmwDcF3WU7mfY7Z/5AQEmTQ+INwtA3yHUD1UjiNzAoKMUwnnC8DRcYTCJAyh8Xs6NCcgyHYa8G78/6ScdT+Py5nACbntY9x9IWGciqXOAakeBQVSyOVF1l2b9f8+Wf/vSWO/6KYGRnkwLtsRqgYAcPe73d0Id6j35dkv49Os/zsW2W454GYai0uTDAgWEkpNlhIDhUywkFvfmmkcuYGZXWVmS10o3P1Wdx/m7pfna0SWoH8XWfde1v/Z72G3rP+vbCL9K7L+36PUTLWwx+IFKp9CxyDXg0XWvZP1f972GU14pUiAvZDw/c6knVT7jcz58qW7P19oo/j6me/Q+pkA38zWADaJz9/n7rML7D8LuDuZLEsSVH0g+cwFXiq00t0/MLOZhEZH62et2jjr/5cLtCPMp0+e11h8VxF/YPoQSib6E374fpC1ebHg9k+EluAZ00vNVAk+c/e5RdbPprGBZLbLCUWr3yHUoQ81s9cJgdSjwBMpBwLZJhdZl/1Dnv0estsHNDXe/yuEtiXLE+rta9HkIusKHYNcHxVZl/0dKec3t1iVBYRjnPF9GtvcVCJzLnfN01uhmD6E9ihG4w1CU1VDzxPaCkkNUEmB5DOthDrsTFFm96znVivz9VbOfcLMdoutmmcRSgWeAf4O/I5QdF1qkfq6cZmphjgvDvGahLx3P1kyP6ZLNFBz93cJd9vZ3Sk3BE4mNB6bbmb3mVl2Y7C0FHsP2ReD7PeQadhZT+P3IK94J5kJxMq5S24J5RyDXF+X+FrljLA4rYn12Z9BUse40nN5jaznmgrEm+zuKS1HJQWST6Gi1GyZlsbZfa+zv09bEPq7l2JxfWVs4T2WpVuhf0i4iE4kBAizCF0BS3EWoRj4esLd+ZWEqo6qcfdnzWxDQtXJAYT69r5xdUdgb2BvM7uR0ECruY3T0tTcC1vmu1JL7yFpiY8n0QzZn0fBMSWaKXMuTyQ0lC3VB3FZajAFpf9OSAtQUCD5lHK3kbmT+F/WczOy/v/Y3T+h+Y6nMSB4h9AN6t/unp12ZmCWUoxy9zPiPj8jdO/aw8wOd/cbyshfYuKF/rH4h5mtS2gMuQ+hTrc9ocHXvwmtt2tF5rNoR/geFGqElmnkmbl7nFFoOymqqfNx9az/k6oem0EoBexcZs+Q7N+F7xbcKqjVEqRlkoICyaebmfVy9yn5VsZBhzIDsbyatSp7lMOtgLsKvYCZbUm4S54MPO3umTrZE+NyEbBbbL2dzzrF3kCW7PrY4wl3Pp2Ai8zsYXdvqmg2cXGgl/WB6e7+fub5OI/8NcA1ZrYvoQslhFKDWgoKJgKHxf+3pLG/fD6bEsa7gCWrS6R0Gzaxfous/5OahnoSISjoZ2a
2019-10-17 17:50:31 +00:00
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
2019-10-17 17:50:31 +00:00
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd3jb5dXw8a8secaJM8jeA04SZigljABhlD1LocxC+0AZpS2UQmlLCzxvW0rpQ+lgFQqUMMqGMAKFQEIS9oaQnAyyIITsxPHWeP+4fxp2ZFm2Jcu2zue6fOm3dUvWOLrHuX2RSARjjDHGmOYU5LoAxhhjjOncLFgwxhhjTEoWLBhjjDEmJQsWjDHGGJOSBQvGGGOMScmCBWOMMcakZMGCMcYYY1KyYMEYY4wxKVmwYIwxxpiULFgwxhhjTEoWLBhjjDEmJQsWjDHGGJOSBQvGGGOMScmCBWOMMcakFMh1AfKViOymqh/nuhwmNRE5F7jHWz1dVf+Tw+KYFmTzfSUiEW/xRVU9MsPXvha4xlvdV1XfbMM1ZgEHAXWqWtJkX7vKLiJTgVe91V+q6h+THFMEjFHVha29fqY0eR5Hq+ryXJWlu7GahQ4mIhUi8lfg/VyXxZjuwt5XuSUihwEfA6fluiwmO6xmoePdBPwg14Uwppux91WOiMhw4KVcl8NklwULHc+f6wIY0w1l/X2lqr5s30d7qOrULF57FtDc47fPtDxgzRDGGGOMScmCBWOMMcakZM0QHaRJL93otmgP5dnZqEIUkTLgIuBEYDzQC9gAvAs8CDyiquFM36933/cC5wCqquNFZE/gV8AUoC+wBngF+LOqftrCtQqBc4FTgN288zcDHwGPAfeoan0L15gCnOnd/xDcc1EJfAHMBm5pSy9uEekHzAJ28Tb9RlV/14rzlwMjgb+q6qUicghwMbAvsAOwEZgL3Kqqr6a4zl7ABbje8MOBMLAOeAP4DzBdVSPNnd9WCa/rLaraW0QGAT8DjgNGAA3AAuBh4DZVrUtxrWLge8B3gD2APrj/83zgKeCfqlrTzP0nbsvK+6q5EQVNRgpMAj4DfgR8FxCgGFgOPAv8n6p+ncZ99Qd+jnseRwG1wCLgIdzzuN3rPdVoiCbH+XDvpx/iXrdh79qPAf9Q1aok5yQ+xthoiITnJOoaEYn+Pw72mi8Sr9MX+DFwDDAOKAPWAm/i3sczmiu3d34AOBU4D/eZ1htYCjwA/CXVuaZ9rGahmxKRb+I+AP6M+4LcASgCBuM+gB4C5ojI4A4oyym4L62TgYFAIe4L7RzgAxE5K8W543BBwT+BbyWc3x84DLgd+EhEpJnzS0XkMWAOcCHuw7EvLlDuA+wKXAJ8IiJntvJx9QJepI2BQpLr3QDMxD1PQ3D/r0G4L89XROT/NXPelcDbuA/QHYES3IfwSFzv9KeA/4pIj7aWLc3yT8F9sV+B+yAvAyqAfXAf5G97XxbJzt2N+P/5cGAA8f/zVOBmQL2gqDMbgAvGbwIm477MSoEJuOdloYh8o4VrTAEUuNI7rxT3Wp2Mex7mea+9tggAjwB34/4v5bjAeS/gj8ACEZnQxmunJCJH477YrwW+iXtMxbjPglOA50XkaRHp2cz5fXGB/QPAwbjPslLc++963HugXzbKbixY6Ei34351PJOwbZL3d14m70hEdsb9ChiK+/V8A+4DeDIuKn/KO3Q/sv8lMgj4N+6D/1+4L/wDgD/gfi0FgPu8D5Kmj2MQ7kt+AlAP3IL7RbI3cIJ33RDui+nVZgKf23FfvuA+xP8H9+VzIO7X1evevgBwh1dT0CLvOXseiH7wtytQwP0KvRJYBVyG+98cBtwKRH+9XS0i+zQpx4G4D3kf7sv2+965U3Cvq8+8Qw8D/rcd5WtJKTAd98VzG3CkV45LgNXeMbsBv296oojsCLyM+xUO7hf4KbjX60m4LzdwXyqzRGSXhNM77H2Vpn/jAtAXiD+G04B3vP29ca/3VJ0lb/SOm4YL7A8Efgqs9/ZHv9jb4hBc8PkFrvZjX9x76Vlv/3DgJRGpSPN6k3Dvyag7iD//70Y3isihuNdHb1xNwm9xX/j7AGfjaucAjgeeFJFGnSZFpACYgXtNAbyG+yzbBzgL9/zuhquVM1lgzRAdRFXXAGtEZGPCtg+zdHf3Az1wHwhTVXVpwr63gUdF5GLcl+8uwG+Aq7JUluiHzlmq+kDC9rki8hyuKaIY+KuIvKSqDQnH3I4LNrYCh6nqOzQ23as1mI77lfEXEsZ5e19CZ3urbwEHNakGnwP8W0QexX2A9gCOwj1/zRKREuBpYH9vU3sDBXCPcz5wgKpuStg+U0RW4AI+cLUxiQl7vo8LFDbgHt+WhH3zvMf2PjAW+B8RuVJVQ+0sazJFuF7xx6rqCwnb3xCR54FPcTUNZ4jIT5r8n2/F1SAAXK6qNzW59lPeNe7F/Y8eEJE9VDXSwe+rdAwCfq+qVydse1tEnsS993YHJuIC3rdSXOdUVX0sYX2OiDyDe42UAmeJyCVtaEb04QLIg1R1fcL26SJyI67pYyhwNa4mJCVV/VBENidsWtP0+ReRUtx7yg98AhyqqusSDnkLuN+rWbsSOBQ4H/f+j/of3HMGLog6J6FZ7S0ReRgXVJ7UUplN21jNQjcjIt/CtfeC++Bdmuw4Vb2VeDR/kdcvIFvuaxIoRMvwOu5XFLj2y8Oi+0RkJ9yvDHAfvk0Dheg1nsX9mgM4RUSGJOzeBViMq8H4Q4r28sSyDU31QLzn6THcBxpkJlCIurJJoBB1J/Hahd2b7Bvk3a5tEigAoKpbcR/8NwG/xjVRZMsTTQKFaBmWER+H3wsYHd0nInsQ/78/lyRQiF7j37hgAdwvyIxmUMygr3DV7I14fQzuTdjU9P+Y6PEmgUL0Gstwv64BeuKamdri3CaBQtQvce8XgB94GRkz4Szir9PzmgQKiX4NLPGWf9pk3wXe7Sbg4qb9b1Q1iMuzsd17wGSGBQvdT2KVYEuJUp73bnvh2hCz5ZYU++5OWD4uYflo4uO6030cBbgmBgBU9UlVFdwv2meSnBe1JmG5OMVxflzH0OhznMlAIYirZdmOF0BEg4im7bnRTpkTROROEdnuC0RV/6Oql6vqLck6r2XQf1PsSwxaEx/DEQnLd7Rw/VsTlo9Kt1Ad7BXviyuZ5p6Dpp5PsW9xwnLS/h8t+CBF4B3Evb6j185U/5Do+2WLqr7d3EHe/UdfQ+Ojgb+IDAT29LY/o6rbmjl/M/BkZopsmrJmiO5nUsLyxmb6/SUzhnj7fSbVAu81t1NVl4nIJlxnp/EJuxIfx/utfBxN7yP2K8T74BmDq8mYiPtA3C/h8FQB9O9wPdOjNqRbqDR8raq1KfZvI94xM9EtuCranrg2+vNEZD4uwHoZmJXlACHR8hT7Ej/gEx9DYv+DluZD+ADXd6UI1y+gM1qeYl9zz0FTX6TYl/gaacvnd6qmD3DPcdSuZOYzIfperkgyeiKVMbj+LkL8h0NLTUxv4/oimQyzmoXuZ4c2ntcno6WIW5tGG3m0SnRQwraMPQ4ROcLrZb0ZV4vwOnAfrp/GYaRfNT/Ku402Z9zgpbrNhKS/lhJEP2QbdYxT1SW4X+eJwz53Bi7FdVrbICLPiEhijVO2pHoMiV8SiY8h2qE0TPx1kJT3yzMaoLXlV3VHaMtz0FRlmvfVloySa1vYn/g/yNRz3N738sCEbS0F6C0OSzVtYzUL3U/0f7oeN/IgXV9moSzgqtdbEu35nDh2PPG1uTduvH46Yu2hXo/zf7J9r/hVuC/Xj3GBw2bckMV0/C+uOvnfuF/zd+CaTHJGVd/wRsBMBb6Na88f6+0uBo4FjhWR+3Edw7KSW6ONWvuFF32tdKbHkGkZz4fRCon/j2ZzYrRS9L38Ma6DbrqWebfpBlmQ/ueEaSULFrqfaK/wcuDjTvDFkM6vk+gvj68Stm1MWP5SVVfTehcSDxQW44Zr/VdVE68dTTiTjptV9RrvnHNww9COEpGzVXVaG8qXMd7/+RXvDxEZheuEeRyuzTiA62j2X1x
2019-10-17 17:50:31 +00:00
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
2019-10-17 17:50:31 +00:00
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgUAAAFGCAYAAAD6uOxSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3dd7zT1f3H8dcFkeEAxSqKiAj1ozjAXYsDR92r7l2taG31V7W0YocVrVtbcRXFUfeoinVrFUfFKtaJoH7UKlBxADIUkX1/f5wTbghJbtb3JoH38/G4j2+S7/menJyb5PvJWd+GxsZGRERERFpVuwAiIiJSGxQUiIiICKCgQERERCIFBSIiIgIoKBAREZFIQYGIiIgACgpEREQkUlAgIiIigIICERERiRQUiIiICKCgQERERCIFBSIiIgIoKBAREZFIQYGIiIgAsFy1CyDlMbNbgJ/Eu2u6+xdVLE7NMrPUNcKfcvc9qloYycvMVgRWd/ePE8j7FhL6vJjZusAn8e717n5yCXkcB/wt3j3C3e9J23cLZZbdzMYB3QF39w1ypNnU3UcXm3elZNTjue4+uFplWRappUBEaoaZHQa8D+xQ7bIsa8ysi5ndATxU7bJI9ailQERqgpltD9zTbEJJyp3AzsD4ahdEqkdBQZ1z9+OA46pcDJFKaJ30E9T658XdbwFuSTD/dfPsTrz+pfap+0BEREQABQUiIiISqfugDGkj2n8LXAIMAE4ENiDU7fvAfcBV7v5tQmW4heRGU/cHnot3twXGAmcBBwPrALOBN4GbgLvcvTFLNun5/QA4CdgRWAuYB4wD/kmoownNHL8GoX53AQxYFZgPfAX8h9Af/UBz5ciR94WE/yPAa8CP3H16gccOBs4BZrh7JzPrAvwK2JdQT/OA94B7gaHuPidHPqsAPwf2AXoDHYBphPfR44QR7QWVqRgZo71/DDxMaGI/FtgIWAmYCDwF/Nnd/9tMftsS/s/bA10J/6MJwAjgGnf/IM/zp/zNzFKj8Hu4+7jiX1nWst1Cjs9L2sj8K939dDPbGfgF4b2/GjAVGAn81d2foxlm1gAcQ6jLPoT/50TgSWCIu3+U5ZjjyDH7IEvaTQnv2Z2ATsBnwDOx/GNzHJN6jYtmH2TUCUD3tO+2W2OXS+brOhQ4EtiSUDffEN7jDwHXufvMXOWOefQBziDU7TqEz/DTwMVA1s+HtAy1FFRGK+DvwDBgK8KXaHtgM+BC4HUzW6d6xauI1YFXgd8B6wPtCF9EOwF3AMPNrF22A81sOTP7K/AycDywXjx+JWATYCDwgZmdlOvJzewnhBPHn4D+wJpAW2AFwpfKQYQA7CEzK6pv1Mx+T4kBQZa8tiMET78hBIcdgI7AD4ArgFfNbNUsx20aj7uA8EXZEWhDqPcdCF+W/40n3CR1IJxYbiIEb6sR6nk9QsAy1sz2ynagmbWNJ5h/E06EPQn/5xUJQc7/xeN/l+xLKJ+ZXUIIYg4iBLDLA10IAfGzZvanZrJYiXCSu5XwGVmVUBc9gVOAN81s9zKKeCjhvXo4sAbhf9SDEDS/ZWanlZF3Tma2OvAiIQDfj6a66QxsB1wGeL73qZmdRfgx8ROavku6Et4zb8V8pUoUFFTG/xG+LCYBpwM/BPYHnoj7DXjezDpUp3gVcT3hJPcmcBThJHcc4dcBwAGEL8BsbiScUCD80vpJPL4/Icj4gvCldn08+S8m/mK7hRBoTQXOA/aMeRwEXAvMjcn3JbTYFMTMTgfOj3fLCghi+R4GVgaGAnsQ3gunEn7FAWxKOPGnl6E1cD8h0PkWGExoDdma8D66MyZdFbjPzNqXWL5C/JlwEnuF8Ct3G8L/9um4vy1wS1xLIP01tALuoukX53jgNMKJoj+hJWU6oQXtgti6kvIZIYA+Me2xc+Jjm9FUdy3lMOBM4H+EX7M/BHYF/gqkfkH/IbZ85XIk4X84mvB+7Ec4kadaGFYEbjOzlUos4xDC9/e18Xm2JwTMswh1PMTMDiowrz8S6vn1eP9zmur+j6lEZrZCLH8/Qj3cQWhZ2prwebya0Hq4FvBPM9so84nM7FTgIqABmAL8mlC/uxG+Y5YnBM9SJeo+qIwuhGbw7d3907THHzazqwhBQw9gEOHLrh51ITRh/9jdUyfgUWZ2P+GX5Q+AQ83suvSmVTM7gKYTxUXunvkr8QUzuxF4nvBr8loze9Tdv0pLc17czgd2d/fXMvIYbmZPAo/E+4cQvmDyii0TqS+gcgMCCF9orYF93P3JtMdfNrPHgTGEX+JHmtkv3X1e3L8d8P14+2fufieLe9jMPiO0PnQF9gIeKKOc+XQBbgeOc/eFqQfN7GHg0fjc3wP2JnSHpBwBHBhvvwTs6e7fpO1/wcxuJfyf1wX+GP/Pr8X301tm1ikt/QR3f6uir6xwXQitNtu7+7S0x0eY2XhCVyGE9/UrefIZTugCSH1eiJ+Xxwgn0dUJgeN9JZSxETjA3R9Ne2ykmT0CvEAIUK8ys4fT3mdZxW67CWaWavKfm6PuLyB8RucTvgcezdj/pJndFp9/RUJr06LAycy+R1MA/hnwQ3dPn/74tJk9TWn1IRWiloLK+WlGQJAyEEitzHZS7I+rR9OAY9K/4ADiWIljgdQJJHMVt9/E7Vjg99kydvfJhL5bCN0Bi37px9aVtoQWgkeyBASpPB4l/BKFcOLMy8yOJvyah8oEBCnDMwKCVPk+oenX9sqEIDGlS9rtD3PkeyVwA6GbY4m+6AqaDZyeHhAAxHEaN6Q91CfjuNT/eQ5weEZAkMpjPE3TARvSjqlFZ2YEBCk30NRakFkH6eYAJ2X5vDTS9L4D2LjE8g3NclLG3f9DU9CyFqGlqWwxYEu15NyQ7bnj878GXBrvbmNm26TtPozQLQahfpdYD8HdHyAEpVIlCgoq471cA49ilH5HvNuFMOagHt3j7lOz7XD3DwndAgB7mNlysGjgXKpvcUQzAwBHEgYrQWgOTeU9y923cvfOhC6afFKDxtrmSxSbVW8hvP9fp3IBAYRBk7mkD9BLbzZ+P+3238xs59gcv4i7T3T3k9z9Ynd/uxIFzeH1XP9ncpQ/DqxMnSAfyxEcA+DuLwDvxru7Zb7OGjEfeDbbjhgopIKFfE3/r2a0dqVLD/yWGF9SoGF59qV34+1WYv6Z+hNauaApuM3l8bTbu6TdTo1FmUP+lq6biyqZVJS6DyrjpWb2v5F2ewPCgL16U8hr3IHwK3gtwmjzvoRfhAC/NLNfFvhc62V7MPXrNfZt9iAM2tqAcELaDugWk+Y70WxC6PtODUacC3xdYLkKMS7PvvQR2Ys+e+7+tpk9QWhS7k0Y4PaVmY0gdM38M9uvqoSMy7Mva/kJMxRS8jWnp6fpTRio2o3aW0HvS3efnWf/TMLJPN/3Z87AiNAak1LKd/AsQldUVu4+zsymE+p3kxLyz2aztNvDzazQ49I/y6lrLXgz9fsaoTWmXltV65qCgspobiDU5LTbXXKmqm3FvsYJhJHrpVgl8wEzW5vQFbMfOYIGQhdGc78814rbOYQWhW0JYz6uLKmkS8o3FSu9pSTzC+9wwqCxo+K+zoSBaYcCmNk7hGDmmuame5WplPJ3Trs9qYDn+DLt9qrUXlDQXP2m6iHfSWuJ7pMcSjnxTcns3smWhhAUlNoSkakSn+U14jZXCwoQuiTj+IZSB2FKGRQUVMb8ZvanT5GbmzNVbSvlNaa/v86n8MFxC9LvmNkehNH5K6Q9nJoXPZawRsEzwIMs/qs1l7cI/ZsvEb7sLogDsjLnyrcYd/8aOMbM/kgYKLkPYZBWm5hkE8Ko7V+Y2U7NrRXQwoo9saW/V5o7udWrotfKqLDU/6RSc/7TP8sHUHggNyPtdiHBVErewZGSHAUFldFcNP69tNufJ1mQBJXyGtP7pr8rZTR57K++hxAQzCOcGO8mNEE2ZqRdcckclvAxsIu7TzWzXwG3xbxvIEw7q6oYmFwKXBpfz/aEEeqHEX5pdSP0J++SM5OWl/5/Xr2A9Guk3c41fkFyW6IlLYvU5zHvr/IipP+
2019-10-17 17:50:31 +00:00
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
2019-10-17 17:50:31 +00:00
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgsAAAFGCAYAAADkcdzhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5gUVdbA4V/35EQeGEDBgB7ErIgRQV3XLGYMKBhAzNnVXT/Dqmtec8JEMGcxLSpGVEAUE+rFBIrkPDl1fX/c6upmmOlJVT2B8z7PPNO3u6runZ6Z7tM3nBtyHAellFJKqbqEW7oBSimllGrdNFhQSimlVEIaLCillFIqIQ0WlFJKKZWQBgtKKaWUSkiDBaWUUkolpMGCUkoppRLSYEEppZRSCWmwoJRSSqmENFhQSimlVEIaLCillFIqIQ0WlFJKKZWQBgtKKaWUSkiDBaWUUkollNrSDVCqrRGR6L7uU4wxB7ZoY1RCIpILdDfG/BbAtccDI91iT2PMYh+vvQnwu1t82BgztgnXGAU84RZPMMY8G/fYeJrZdhGZB/QFjDGmfx3HbGeM+bax1/ZLjefxOmPMtS3VlrZOexaUUu2SiAwHfgL2bum2bGhEpEBEngRea+m2KH9oz4JSqt0RkcHAs/UeqILyFLAvML+lG6L8ocGCUqo9Sgm6AmPMKGBU0PU0lTFmPDA+wOtvkuDhwJ9/lVw6DKGUUkqphDRYUEoppVRCOgwRsLiZ81cCtwBnAKOB/tjn/yfgBeAeY0xxEtpzAHAqsDvQAygFfgHeAu41xiwPqN6hwAducXdgDnAFcAzQBygDZgOPAU8bY5xaLhN/vd2AMcAQoBdQCcwD3sE+l3/Uc34P7O9hP0CALkAVsAL4Ajve/VJ97ajj2v/B/r4BZgH7G2NWN/Dca4FrgDXGmE4iUgBcDByGfZ4qgR+B54AHjTHldVynM3AWcCgwAMgGVmH/3t7CzrBvUJsao8bs8yOBydiu+lOArYE84C9gCnCHMebXeq63O/b3PBjojf0d/QFMBe4zxsxNUH/UEyISXRWwqTFmXuN/slrbNp46VhTErRS42xhzoYjsC5yN/dvvBqwEpgEPGGM+oB4iEgJOxj6X22N/n38B/wPuMsb8Uss5o6hjNUQtx26H/ZvdB+gELATec9s/p45zoj+jtxqixnMC0DfuNXCCO3RT8+c6DjgRGIh9bgqxf+OvAQ8ZY4rqard7je2Bi7DPbR/s//C7wM1Arf8fqvG0ZyF5wsDzwDhgF+yLZhawI/Af4EsR6RNU5SKSIyIvY19chmP/qTKwLwwDgauBX0XksKDaEKc7MBP4J7AlkOm2Yx/gSeBlEcms7UQRSRWRB4DPsUHPZu75ecC2wCXAXBEZU1flIjIS+4ZyPTAU6Il9LnKwz8vR2ADuNRFp1NiriPyLJgYKtVxrL2xQdRk2uMwGOgK7AXcCM0WkSy3nbeeedyP2BbQjkIZ93vfGvoj+6r4RBykb+4bzGDao64Z9njfDBjJzROTg2k4UkQz3jecz7Bvk5tjfcy42+DnPPf+fwf4IzScit2CDm6OxgW06UIANlN8XkevruUQe9s1vAvZ/pAv2udgcOAeY7X4IaKrjsH+rx2M/QGQAm2KD6a9F5IJmXLtOItId+AQbmB9O7LnpCuwF3AaYRH+nInIF9kPGSGKvJb2xfzNfu9dVPtBgIXnOw744LAUuBPYAhgFvu48L8KGIZPtdsYiEsZ/wjnTvegMbMOwK/A0brKwBOgCviMh+frehhoexb36zgZOwb36jsJ8mAI7AvjDW5lHsGw3YT2Yj3fOHYoOPxdgXu4fdoGAd7ie88dhAbSXwb+Ag9xpHA/cDFe7hh2F7ghpERC4EbnCLzQoU3PZNxv5OHgQOxP7NnIv91AewHTYgiG9DCvAiNgAqBq7F9p4Mwv69PeUe2gV4QUSymti+hrgD++Y2HfupeFfs7/Zd9/EMYLybCyH+ZwgDTxP7hDofuAD7BjIU2/OyGtszd6PbGxO1EBuAj4677xr3vh2JPXfJMhy4HPgT++l3D+z/3ANA9BP3VW5PWV1OxP4Ov8X+Pe6JfYOP9kjkAhNFJK+JbbwL+15wv1vPYGwgXYJ9ju8SkaMbeK2rsc/zl255EbHn/uroQSKS47Z/T+zz8CT29WkQ9v/xXmxvYy/gHRHZumZFInIucBMQApYDl2Kf379jX2PSsUG18oEOQyRPAbabfLAxZkHc/ZNF5B5sMLEp8A/si5ufzscuYwIYa4x5uMbjU0XkMeBTt51PiMjmxphKn9sRVYDtCj/SGBN9Y54hIi9iP4nuBhwnIg/Fd9GKyBHE3kBuMsbU/FT5kYg8CnyI/fR5v4i8YYxZEXfMv93vVcABxphZNa7xsoj8D3jdLR+LfeFJyO3JiL4wNTdQAPtClwIcaoz5X9z9n4vIW8D32E/uJ4rI+XG/q72ALdzbZxpjnmJdk0VkIba3ojdwMPBSM9qZSAEwCRhljIlE7xSRydiA9WAgHzgEO6wSdQJwlHv7U+AgY0xh3OMficgE7O95E+Bq9/c8y/17+lpEOsUd/4cx5mtff7KGK8D28gw2xqyKu3+qiMzHDk2C/buenuA6L2OHEqL/L7j/L29i31y7YwPKF5rQRgc4whjzRtx900TkdeAjbOB6j4hMru81wR3++0NEokMHFXU89zdi/0ersK8Db9R4/H8iMtGtPxfbO+UFVCKSTywwXwjsYYyJX6b5roi8S9OeD1UL7VlIrtNqBApRlwDRDHNj3HE8X7if0i52i2/XEigA4Ga4u8ItbkzsxToIq4CT41/43DYUY8e2o28sNbPWXeZ+nwP8q7YLG2OWYceGwQ4reD0Dbq9NBrZH4fVaAoXoNd7AfnIF+4aakIiMwH76B38ChaiXawQK0fb9TuzTeQdskBlVEHf75zquezfwCHa4ZL2xbh+VARfGBwoA7jyQR+Lu2r7GedHfczlwfI1AIXqN+cSWLYbizmmNLq8RKEQ9Qqx3oeZzEK8cGFPL/4tD7O8OYJsmtu/BWt6sMcZ8QSyY6YXtmWo2N5CL9vw8Ulvdbv2zgFvd4q4ismvcw8Oxw2tgn9/18jkYY17CBqvKBxosJM+PdU1kcqP1J91iAXZOg1+2w775Q+wNpi5vx90OcijiWWPMytoeMMb8jB1eADhQRFLBm7AXHbucWs/Ew2nYSVIQ93MYY0qMMbsYY7pih4QSiU5Wy0h0kNs9Ox77v/Ql/gUKYCdr1iV+YmB89/NPcbefEJF93YDRY4z5yxgzxhhzszHmGz8aWocv6/o9U0f73Qmd0TfON+sIrgEwxnwE/OAW/17z52wlqoD3a3vADSCiQUSiIYSZNXrH4sUHhOvNX2mgcQkeix8O/HsTr1/TUGyvGNT/mvRW3O3416ToXJdyEveMPd6olqk66TBE8nxaz+Nfxd3uj50A6Icd427/V0T+28DzNvOp/to05LnYG/upuRd29vsO2E+QAOeLyPkNrKvWnyP6adcdO90UO1msP/aNai9iAVaiN6BtsWPr0UmQFcDaBrarIeYleCx+hrj3f2yM+UZE3sZ2TQ/ATqxbISJTsUM879T2KSwg8xI8Vmv7sSsmohJ1y8cfMwA7QXZjWl/GwCXGmLIEjxdh3+QTvRbXGTBhe2+imvJ6XoId0qqVMWaeiKzGPr/bNuH6tYl/TXpZRBp6Xvz/cnQvClPP8zsL23vjW2/thkqDheSpb2LVsrjbBXUe1XjdmnheZx/bUFNjn4s/8PHnEJGNsEM/h1N3UBSh/p63Xu73cmwPxO7YuSd3N6ml60u0ZCy+Z6XmC+Hx2MlqJ7mPdcVOiDsOQES+wwY599W3LK2ZmtL+rnG3lzagjiVxt7vQ+oKF+p7f6POQ6M1svWGYOjTlDXF5zWGi2o7BBgtN7bmoyY//5R7u97p6XAA7tOnOn2jq5E/l0mAhearqeTx+iV5FnUc1Xvzv+GzsksOGKPWxDTU15bmI/zluoOGT8qrjCyJyIHa1QE7c3dF13XOwORbeA15h3U+5dfkaO376KfZF8EZ3IljNtf5JY4x
2019-10-17 17:50:31 +00:00
"text/plain": [
2021-03-10 12:30:50 +00:00
"<Figure size 555x330 with 1 Axes>"
2019-10-17 17:50:31 +00:00
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"\n",
"xlabel = {\n",
" 't_i_peak': 's',\n",
" 't_e_peak': 's',\n",
" 'p_i_peak': 'prob',\n",
" 'p_e_peak': 'prob',\n",
"}\n",
"\n",
"for cell_type in ['gridcell', 'bs_not_gridcell', 'ns_inhibited', 'ns_not_inhibited']:\n",
2021-03-10 12:30:50 +00:00
" for key in xlabel:\n",
" fig = plt.figure(figsize=(3.7,2.2))\n",
2019-10-17 17:50:31 +00:00
" plt.suptitle(key + ' ' + cell_type)\n",
" legend_lines = []\n",
2021-03-10 12:30:50 +00:00
" for color, label in zip(colors, labels):\n",
" legend_lines.append(matplotlib.lines.Line2D([0], [0], color=color, label=label))\n",
" sns.kdeplot(data=results[cell_type][key].loc[:,labels], cumulative=True, legend=False, palette=colors, common_norm=False)\n",
2019-10-17 17:50:31 +00:00
" plt.xlabel(xlabel[key])\n",
" plt.legend(\n",
" handles=legend_lines,\n",
" bbox_to_anchor=(1.04,1), borderaxespad=0, frameon=False)\n",
" plt.tight_layout()\n",
" plt.grid(False)\n",
" despine()\n",
2019-10-17 17:50:31 +00:00
" figname = f'histogram-{key}-{cell_type}'.replace(' ', '-')\n",
" fig.savefig(\n",
" output_path / 'figures' / f'{figname}.png', \n",
" bbox_inches='tight', transparent=True)\n",
" fig.savefig(\n",
" output_path / 'figures' / f'{figname}.svg', \n",
" bbox_inches='tight', transparent=True)"
]
},
2021-03-10 12:30:50 +00:00
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Stats"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"stats = {}\n",
"for cell_type, result in results.items():\n",
" stats[cell_type], _ = make_statistics_table(result, labels, lmm_test=True, wilcoxon_test=False, ttest_ind_test=True, ttest_rel_test=True)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"stats['gridcell']"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"for cell_type, stat in stats.items():\n",
" stat.to_latex(output_path / \"statistics\" / f\"statistics_{cell_type}.tex\")\n",
" stat.to_csv(output_path / \"statistics\" / f\"statistics_{cell_type}.csv\")"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [],
"source": [
"for cell_type, cell_results in results.items():\n",
" for key, result in cell_results.items():\n",
" result.to_latex(output_path / \"statistics\" / f\"values_{cell_type}_{key}.tex\")\n",
" result.to_csv(output_path / \"statistics\" / f\"values_{cell_type}_{key}.csv\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Plot PSTH"
]
},
2019-10-17 17:50:31 +00:00
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 14,
2019-10-17 17:50:31 +00:00
"metadata": {},
"outputs": [],
"source": [
"from septum_mec.analysis.plotting import plot_bootstrap_timeseries"
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 15,
2019-10-17 17:50:31 +00:00
"metadata": {},
"outputs": [],
"source": [
"psth = pd.read_feather(output_path / 'data' / 'psth.feather')\n",
"times = pd.read_feather(output_path / 'data' / 'times.feather')"
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 16,
2019-10-17 17:50:31 +00:00
"metadata": {},
"outputs": [],
"source": [
"times = times.T.iloc[0].values"
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 17,
2019-10-17 17:50:31 +00:00
"metadata": {},
"outputs": [
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAq0AAAFFCAYAAADGo9TZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdeXhU1f348fedfSZ7Qtg3BblufPmKqFWrAq1bxa1i1aLW1qpURXGv1n3ff3UrPK31647UFSsodbdVKyIoi3LZBMKafZLZZ+69vz/uJGSSmawTksDn9TzzTObec+49o+HkM2fO+RzFNE2EEEIIIYTozWw93QAhhBBCCCHaIkGrEEIIIYTo9SRoFUIIIYQQvZ4ErUIIIYQQoteToFUIIYQQQvR6ErQKIYQQQoheT4JWIYQQQgjR60nQKoQQQgghej0JWoUQQgghRK8nQasQQgghhOj1JGgVQgghhBC9ngStQgghhBCi15OgVQghhBBC9HoStHYzVVXfVlX17Z5uhxBC7ArS5wkhuoujpxuwBxg1evTo/QGzpxvSQDdMFm3xs6oy0NNNaWHffrkcOqQAu03p6aYIsavsbr/s0ud1gPR5Yg/U6V92CVr3cKOLfdiUnu8sDdNkbXWop5shhNjNSZ8nRN8lQesezqYovaIDF0KIXUH6PCH6LpnTKoQQQgghej0JWoUQQgghRK8nQasQQgghhOj1JGgVQgghhBC9ngStQgghhBCi15OgVQghhBBC9HoStIrdRiIaoezLfxGq3N7TTRFCCCFElkmeVrFbqNvyI6vnv0AiEqJK+5YDzroMm8OJ3enq6aYJIYQQIgtkpFX0eaZpsv79V0lEQtT+uIpowM+3zz2ENu8ZTLPX7CQphBBCiC6QkVbR50X9VcSCdVSuWoqhJ6heuxxT1zH3G0/UX42nsKSnmyiEEEKILspq0Kqq6krgGeBFTdN2ZPPardxzAzCincU/1TRtYrKeCwgAzlbK+zVNK+xK+0T3C5ZvIRasxzQM+h9wCIlImJr1K4nV1+IvWyNBqxBCCLEbyPZI637Ag8B9qqq+C/wf8I6maYks36epr4HNrZzPBcYlf17T5Ph+WAFrCFiaoW59l1snul2wYhvxUD2+0kE4PD4cHh+5g0YQrqmg7IuFuPMKKRy5b083UwghhBBdkO2g9TzgN8Bk4GRgClClquqLwLOapi3L8v3QNO3M1s6rqjoXK2hdAlzR5FRDIPuZpmknZrtdYtcwTZO6sjXE6mvxlQ5uPJ5TOpjqNcupWv0diqKwz0nnkT9kr3ZdD0BRlG5rsxBCCCE6LqsLsTRNe0nTtOOA4cBNwCqgHzATWKqq6mJVVS9TVbUom/fNRFXVi4BfYU0D+JWmaeEmpxuC1uW7oi2ie4SrthOs3E48VI+7YOc0AJvDScm+B2FzOKlau4INH7/V5qIsPVBFYMlb1Lz7CHq4rrubLoQQQogO6JbsAZqmbdU07X5N0w4ADgX+AlQD44EngK2qqs5VVfVEVVW7ZUhLVdVBwKPJl3/UNG1dsyINQeuK7ri/2DXqNq8j6q/ClV+EzW5POacoCoV77UesvpZwbSXh6vKM14mVr6Nm4WNE1i8mXl1GfPvq7m66EEIIITqg21NeaZq2WNO0y4GBwM+Ax4E4MBV4ByhTVfXOZJCZTfdizWddjBU0NycjrbuBUOV24qEArpyCtOcVmw1nbj6xQC2B7RszXiesfUasfC3RjUswwnXogaruarIQQgghOmFXprw6DDgJa65rDqAAOjAYuBm4VlXVOzRNe6CrN1JVdTzW3FqAqzRNM5udH4w1bUEH6lVV/VOyfT7gR+ANTdPe7Wo7RPcLVVlBq6eoX8YyrtwCYvV+Nn3+HrFAPYMnTEwZldUD1cR2rEP370Bx52BGg4TX/hd7TjGevQ/ZFW9DCCGEEG3o1qBVVdV9sBZnncvOtFQKsAwrs8CLwCDgYmA6cK+qqoamaQ918dY3Ju/zkaZp/0lzvmGU1cQaafU0O/97VVXnA+domiYZBHopPR4j4q8mHgni8OZmLOctHkDVqiUoNhvxYD3u/EJK9zu48Xxs2yqMQBU2bwHO0r2IbVuFHqwm8N0CCVqFEEKIXiLrQauqqsXA2VjB6qHJwwrgB+YAf9c07ZsmVaqAK1RVXQU8CVwOdDpoVVV1JHB68uXdGYo1BK0O4E3gAWAlUACcAdyPNSo8F/hFO++7MsOpUe2pLzouXLWDRDiIze5odbtWh9tD8T7/Q+WqJSQiYXZ89wX91INQbNbsmNi2VeiBKux5/ayR1niY+DYN++jDMU1TMgkIkYb0eUKIXS3bmwu8BZyAlf9UwRrJ/Bhrw4E3NE2LtFL9HaygtbSLzbgcsAP/1TTt4wxllgGzgO2apt3Z5HgE+IuqqsuAz4ATVVU9UaYK9E7h6h3Ew0Ec3pw2yzp9uQwYdyTlK76ifvsmqtcup2TMOIxokHjlJvRgNY7+o1BsdtwjxhPZ8A2mnsCI1GP35u+CdyOEEEKI1mR7pPWU5PMm4Fng/zRNy7z6JVUxsBr4vItt+FXy+dlMBTRNWwAsaOX8f1RV/QA4FmvUts2gNZkpoYXkaMT+bdUXHReuqSARDuJsR9AKYLPbye0/lMD2Mip++IaSMeOIbV+NHqxGcfmwOa1ZIjZPHjanFzMexgjWZD1oNQyDjz76iIULF7J8+XJ27NhBPB6nsLCQvfbai6OOOoqpU6dSXFzcrustXbqUt956i2+++YYtW7Y0XuuAAw7gxBNP5OSTT8beLLNCZ9s9bdo0lixZwsqVK3E42td9GIbB2WefzXfffYemaR2+76uvvsrNN9+M3W7n+++/b7P8Oeecw5IlS5g6dSr33HNPh+8n2kf6PNERfanfe++995gzZw4rV64kFosxaNAgJk6cyO9+9zsGDBiQsd4333zD3/72N5YuXUooFGLQoEEce+yxXHLJJeTnd+zvyB//+EfefPNNDj30UF544YU2y6uqCsB9993HL3/5yw7dqy/JdtD6Ctao6ofNFz+1RdO0b4EubVukquoEYBhWdoJXu3ItrF2yjgXazkgvekS4yhppbW0RVnOe4v4Evt9EYHsZtRtW4ajaiBGqwZ6T2lEqLi9mPIIerMHZr727BLdt2bJl3HjjjaxduxYAt9vNkCFDyMnJoaqqikWLFrFo0SJmz57NjTfeyJlnZt47o76+nltvvZUFC6zPX06nk0GDBpGTk0NZWRmffPIJn3zyCc888wxPPPEEI0Z07X08+uijLFmypMP1HnroIb777rsu3VsI0Xf1pX7vT3/6E6+99hoApaWlDB06lI0bN/Lss8/yxhtvMHv2bA4++OAW9RYsWMA111yDYRgMGDCAIUOGsGbNGp5++mkWLFjAyy+/zKBB2U6StOfJdtA6G/C3J2BVVfU4YLSmaenSUXXWqcnnhZqmVbdxfzugtLLFbEM6sHi2GieyxzTNnSOtg0e2u57D7cHmcBIL+Fm94GVGDi/AiAZxFAxOKac4PcmR1lZ/jTrkiy++YPr06USjUUaOHMnMmTOZNGkSHs/OdYAbNmzgqaee4u233+bmm2/G4XBw+umnt7hWVVUVZ511FmVlZZSWlnLZZZdxxhln4HJZc3tN0+SDDz7gwQcfRNM0pk2bxmuvvcbAgQM73G5d13n44Yd55plnOlQvkUjw0EMP8eyzz3b4nkKI3UNf6vdeffVVXnvtNRwOB/fddx+nnGJ9eRwIBLjllltYsGABM2bM4MMPP8Tr9TbWW79+Pddffz2GYXDLLbcwbdo0FEWhsrKSK6+8ksWLF3PNNdfw8ssvd+U/pSD7eVo/wcrD2h53A9n+7u6nyecPWiukquoGrGD0ilaKjU8+t/19pNjl4sE6YsF69HgUh8fXobqeolKq164gUlNOeNtazFgYxZ16DcXlxYiFs5avtbKykpkzZxKNRjnkkEN4/fXXOfHEE1M6boCRI0fy0EMPceGFFwLWVz11dam7c5mmyQ033EBZWRlDhgxhzpw5nHPOOY0dN1gbKxx77LG89NJLlJSUUFFRwR133NHhdm/
"text/plain": [
"<Figure size 750x300 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
2019-10-17 17:50:31 +00:00
"source": [
2019-12-13 10:43:57 +00:00
"cs = ['#993404', '#980043', '#d95f02', '#e7298a']\n",
"lb = ['GC 11 Hz', 'GC 30 Hz', 'NSI 11 Hz', 'NSI 30 Hz']\n",
"\n",
"fig, axs = plt.subplots(1, 2, sharex=True, sharey=True, figsize=(5,2))\n",
"ii = 0\n",
"for cell_type, ls in zip(['gridcell', 'ns_inhibited'], ['-', '--']):\n",
" for i, (ax, query) in enumerate(zip(axs.ravel(), queries)):\n",
" selection = [\n",
" f'{r.action}_{r.channel_group}_{r.unit_name}' \n",
" for i, r in data.query(query + ' and ' + cell_type).iterrows()]\n",
" values = psth.loc[:, selection].dropna(axis=1).to_numpy()\n",
"\n",
" plot_bootstrap_timeseries(times*1000, values, ax=ax, lw=2, label=lb[ii], color=cs[ii], ls=ls)\n",
" # ax.set_title(titles[i])\n",
" ax.set_xlabel('Time (ms)')\n",
" ax.legend(frameon=False)\n",
" ii += 1\n",
" axs[0].set_ylabel('Probability density')\n",
" despine()\n",
2019-12-13 10:43:57 +00:00
" plt.xlim(-5, 30)\n",
2021-03-10 12:30:50 +00:00
"for ax in axs.ravel():\n",
" ax.axvspan(0, 5, color='#43a2ca', alpha=.3, zorder=-5)\n",
"figname = f'response-probability-gc-ns'\n",
"fig.savefig(\n",
" output_path / 'figures' / f'{figname}.png', \n",
" bbox_inches='tight', transparent=True)\n",
"fig.savefig(\n",
" output_path / 'figures' / f'{figname}.svg', \n",
" bbox_inches='tight', transparent=True)"
2019-10-17 17:50:31 +00:00
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 18,
2019-10-17 17:50:31 +00:00
"metadata": {},
"outputs": [
{
"data": {
2021-03-10 12:30:50 +00:00
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAq0AAAFFCAYAAADGo9TZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdeZhU1bno/2/NPY+0yKyCLjESIxInTDSYKB5xOqJeRT3kGhOM4ohyJBGNQFQu4URB4fdL4lFROZqLQpROEI0aByIgEgTprTJ2Q9Nj9VBd8977/rGrmy66q8fqCd7P89RT1N5r7b1aYfVba6/1LptpmgghhBBCCNGf2fu6AUIIIYQQQrRHglYhhBBCCNHvSdAqhBBCCCH6PQlahRBCCCFEvydBqxBCCCGE6PckaBVCCCGEEP2eBK1CCCGEEKLfk6BVCCGEEEL0exK0CiGEEEKIfk+CViGEEEII0e9J0CqEEEIIIfo9CVqFEEIIIUS/J0GrEEIIIYTo9yRo7WFKqb8opf7S1+0QQojeIH2eEKKnOPu6AceA0WPGjDkNMPu6IY10w2TjgVqKKn193ZQWTh2UwdnDsnHYbX3dFCF6y9H2l136vE6QPk8cg7r8l12C1mPcmLw07La+7ywN0+Tban9fN0MIcZSTPk+IgUuC1mOc3WbrFx24EEL0BunzhBi4ZE6rEEIIIYTo9yRoFUIIIYQQ/Z4ErUIIIYQQot+ToFUIIYQQQvR7A34hllJqLzCqg8U/1DTtolg9N+ADXG2Ur9U0Lac77RNCCCGEEN034INWYBNQ0sb5DOCM2J+/aXZ8LFbA6ge+SFC3vtutE0IIIYQQ3Tbgg1ZN065r67xS6jWsoHULcHezU42B7D80Tbush5onBhA9UEfw2w2knjwRe0pGXzdHCCGEEM0M+KC1LUqp24HrsaYBXK9pWqDZ6cag9cteb5joV0zDoO6j/yZSXYKpR4jWlJL9g+l93SwhhBBCNHPUBq1KqSHA4tjH/9Q0bdcRRRqD1u291yrRHxl+L+Hy3QDodWVEHG1NcxZCCCFEXziaswf8Fms+62bguVbOy0irAED312BGAhj+GqI1pQCYZr/ZNl0IIYQQJHmkVSm1A3geeFnTtLJkXruT7RgP/Efs432applHnB8KDAJ0oF4p9SvgHCAN2AO8oWnaX3uxyaIPGQ1ezEiQaE0pRtiPGQ2j11fgzDqur5smhBBCiJhkTw8YCywEnlBK/RX4b+BtTdOiSb5Pex4GbMDfNU37uJXzjaOsJtZIa8oR53+mlFoL3KhpmmQQOMrpDV6MSBDdX2N9rivDv309WedP6+OWCSGEEKJRsoPWW7BGOCcBVwBTgCql1MvAC5qmbUvy/VpQSp0AXBP7OD9Bscag1Qm8CTwF7ACygWuBJ4HLgdeAf+vgfXckODW6I/VF3zH8NZiRIK7jRuNIzyO0bwvhUg09UIcjNauvmydEvyR9nhCityU1aNU07RXgldjj91uxgtixwL3APUqpL7BGX1/VNM2bzHs3cxfgAP6padr7CcpsA5YBhzRNe7zZ8SDwnFJqG/AP4DKl1GUyVeDopvuqMMMBHOn52D3p2FwpGCEfek2pBK1CCCFEP9Ej2QM0TTuINVr5pFJqAjAduAEYH3stUkr9BXgB+NuRc0676frY+wtttK8QKGzj/MdKqXeBn2CN2rYbtGqa9p3WjsdGI05rr77oG7q/lqj3IEawHntKJgD2lEyMQD3RmoO4h6g+bqEQ/ZP0eUKI3tbj2QM0TdusadpdwPHAxcAzQASYCrwNFCulHo+lqOqWWIA8Inb9P3fzco27ZJ3YzeuIfix8YDu6rxJ7SiY2lwcAuycDI1hP1Hugj1snhBBCiEa9mfLqHKx5olcA6VgLpQxgKPBrYJdSanY373FV7H2dpmnVbRVUSjmUUm2NNDf+t4l0s02inzL1KMFdG4nWHMSRNbjpuD09F72hmvChb9DrK/uwhUIIIYRo1KNBq1Lq5Ngo6m7gI+A+rJHLL2N/Ph74LrAUa6rCb5VSD3bjlhfE3t9tp117sYLRu9soNj72/lU32iP6sUjZt0SqSzDCwfig1ZOOPTWbqPcgwf1b+7CFQgghhGiU9DmtSqk84H9hLcI6O3bYBtQCK4E/aZr2ebMqVcDdSqkirOD1LuD/dOG+duCs2MdN7RT/ChgF3KyUeubIlFxKqbOBH8U+/k9n2yIGhmh9OUawDkd6LjZ7/Pc3R1oORqgew1fVR60TQgghRHNJHWlVSq0GDgJLsKYDALwP3AwM0TTtl0cErM29HXsv6OLtxwCZsT+3tzXrU1g5Ws8E/qiUaloirpS6CFiDFWi/rGna5i62R/Rzuq8KIxzA5k5tcc7mTsUMB9EbeirJhRBCCCE6I9nTA64E3EAx8DhwkqZpP9Y07VVN04Lt1M0DvsYaje2KYbH3qKZpdW0V1DTtQ6w0XDpWXtlDSqkvYtMY3seatlAI3N7FtogBwPBVYYb92F1pLc7ZXKnWhgMN3m5v6VpSUoJSCqUUs2e3P237s88+ayqfSFFREfPnz+fyyy9nwoQJjBs3jgsvvJCf/exnvPzyywSDif+5TZo0CaUUS5Ys6dLP09z999+PUop9+/Z1qt7dd9+NUoqSkpJO3/PTTz9t+u9z6NChdsvPmjULpRTTp0/v9L2EEF1zNPV7GzZsYMaMGZxzzjmcfvrpTJo0iblz57J79+4263399dfcd999nH/++U31HnvsMcrKOr9h6JIlS1BKMWnSpA6VT2Y/358kO2j9H+AS4ERN0x7TNK3Dv8k0Tduqadqpmqbd1sV7N47Q1nbwfs8A52EFydXAd7A2F/g71tSGKR0ItMUAZZomer2Vn7XVkVZXCuhhzHAAM+xP2n1Xr17N++8nSh/cMc888wzXXHMNK1as4NChQ4wYMYKxY8dit9v56KOPmDdvHpMnT2bHjkS535Nj5cqVrF27ttP1Xn75ZdatW9cDLRJC9EcDud9bunQp06dPb2r/ySefTF1dHa+99hpXX311wr5s8+bNTJ06lcLCQgzD4JRTTsHr9bJy5UquuOIKdu7cmdR2HiuSPad1OVDbkbyrSqlLgDGapj2XjBtrmvY68Hon62wCbkrG/cXAYoYaMEINmNFQ60Gr3Q4Od9Noq92TnrR7P/LII6xdu5bs7OxO1121ahXPPvssaWlpPPHEE/zkJz/B4XA0nd+1axdz5sxh69at3HbbbRQWFpKXl5e0tjd64YUXePLJJztd7/nnn2fhwoVJb48Qon8biP3ehg0bmkYqH3roIaZPn47D4SAcDrNo0SJefPFFHnzwQc444wyOP/74pno1NTX88pe/JBQKcfvtt3PvvffidDrx+XzMmTOHdevWMXPmTAoLC3G73d1u57Ek2SOtH2DlYe2I+cCCJN9fiA6Jeg9ghHzYXGnY7I5Wy9jdaZjhBvS6zj/KScRms1FRUcH8+Yl2GG7b8uXLAasDnTx5clzHDTB69GiWLVtGfn4+Xq+Xl156qdttbq68vJy7776bJ554olPTJsrKyrjrrrt46qmnuj3dQggxsAzUfu+Pf/wjAFOmTOG2225ruq/b7ebhhx9m9OjRhEIh3nzzzbh6K1asoLa2lu9973vMmjULp9MaH8zIyGDRokWMGDGC4uJi1qxZk5R2Hku6HLQqpTKVUic1f8VOpR55/IjXaKXUJOBUwJWUn0KITop6D2AE6rC3sU2rPTUTI1BHtLrz8y4TmTZtGgB/+ctfeO+99zpVt7a2lv379wNwxhlnJCyXl5fHj3/8YwC2bdvWxZa2tH79ei699FLWrVtHdnY2c+fO7VC9v/3tb1x66aWsX7+enJycDtcTQhwdBmq/N378eCZNmsS1117b4pzNZuOUU04B4ODBg3HnGoPYqVOntqjndrubrvf222+3OC/a1p3pAVlY+VZTmh0zsdJOfdPBa2zoxv1FP2Jr6HoSftOdDs6UVs/p9ZVdHpmze9ISPtYPlxah11VgS83GCDU0q3O4vD01m2jVPiJJDFpvueUWioqK2Lx5M48++ihnnXUWOTk5Har
2019-10-17 17:50:31 +00:00
"text/plain": [
"<Figure size 750x300 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
2019-12-13 10:43:57 +00:00
"cs = ['#d95f02', '#e7298a', '#993404', '#980043']\n",
"lb = ['NSI 11 Hz', 'NSI 30 Hz', 'NS 11 Hz', 'NS 30 Hz']\n",
"\n",
2019-10-17 17:50:31 +00:00
"fig, axs = plt.subplots(1, 2, sharex=True, sharey=True, figsize=(5,2))\n",
"ii = 0\n",
"for cell_type, ls in zip(['ns_inhibited', 'ns_not_inhibited'], ['--', '-']):\n",
2019-10-17 17:50:31 +00:00
" for i, (ax, query) in enumerate(zip(axs.ravel(), queries)):\n",
" selection = [\n",
" f'{r.action}_{r.channel_group}_{r.unit_name}' \n",
" for i, r in data.query(query + ' and ' + cell_type).iterrows()]\n",
" values = psth.loc[:, selection].dropna(axis=1).to_numpy()\n",
"\n",
" plot_bootstrap_timeseries(times*1000, values, ax=ax, lw=2, label=lb[ii], color=cs[ii], ls=ls)\n",
2019-10-17 17:50:31 +00:00
" # ax.set_title(titles[i])\n",
" ax.set_xlabel('Time (ms)')\n",
2019-10-17 17:50:31 +00:00
" ax.legend(frameon=False)\n",
" ii += 1\n",
" axs[0].set_ylabel('Probability density')\n",
" despine()\n",
2019-12-13 10:43:57 +00:00
" plt.xlim(-5, 30)\n",
2021-03-10 12:30:50 +00:00
"for ax in axs.ravel():\n",
" ax.axvspan(0, 5, color='#43a2ca', alpha=.3, zorder=-5)\n",
"figname = f'response-probability-nsi-ns'\n",
2019-10-17 17:50:31 +00:00
"fig.savefig(\n",
" output_path / 'figures' / f'{figname}.png', \n",
" bbox_inches='tight', transparent=True)\n",
"fig.savefig(\n",
" output_path / 'figures' / f'{figname}.svg', \n",
" bbox_inches='tight', transparent=True)"
]
},
2019-12-13 10:43:57 +00:00
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
2019-10-17 17:50:31 +00:00
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Store results in Expipe action"
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 19,
2019-10-17 17:50:31 +00:00
"metadata": {},
"outputs": [],
"source": [
"action = project.require_action(\"stimulus-response\")"
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
2021-03-10 12:30:50 +00:00
"['/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_not_inhibited_p_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_t_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_not_inhibited_t_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_inhibited_p_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_not_gridcell_p_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_gridcell.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_t_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_gridcell_p_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_ns_inhibited.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_gridcell_t_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_bs.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_not_gridcell_t_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_not_gridcell_p_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_gridcell_p_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_inhibited_t_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_t_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_not_gridcell_p_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_gridcell.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_p_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_inhibited_t_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_inhibited_p_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_p_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_not_inhibited_t_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_inhibited_p_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_bs_not_gridcell.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_bs.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_gridcell_t_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_not_inhibited_p_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_not_gridcell_t_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_gridcell_t_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_not_gridcell_t_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_gridcell_p_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_not_gridcell_t_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_not_inhibited_t_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_t_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_not_inhibited_p_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_not_inhibited_p_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_inhibited_p_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_gridcell_t_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_not_inhibited_t_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_inhibited_t_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_ns_not_inhibited.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_ns_inhibited.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_gridcell_p_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_ns_not_inhibited.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_p_i_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_p_e_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/statistics_bs_not_gridcell.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_ns_inhibited_t_i_peak.tex',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/statistics/values_bs_not_gridcell_p_e_peak.csv',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/data/times.feather',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/data/psth.feather',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-gridcell-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-ns_inhibited-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-ns_inhibited.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-bs_not_gridcell.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-ns_not_inhibited.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-bs_not_gridcell-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/response-probability-nsi-ns.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-gridcell-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/response-probability-nsi-ns-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-ns_not_inhibited-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-gridcell.svg',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-ns_not_inhibited-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-ns_inhibited-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-ns_inhibited.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-ns_inhibited.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-ns_not_inhibited-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-ns_not_inhibited.svg',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-ns_not_inhibited-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-bs_not_gridcell-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-bs_not_gridcell-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-bs_not_gridcell.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-gridcell.svg',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-bs_not_gridcell-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-ns_inhibited.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-ns_not_inhibited.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-gridcell-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/response-probability-nsi-ns-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-ns_inhibited.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-bs_not_gridcell.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-bs_not_gridcell.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-gridcell-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-ns_inhibited-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/response-probability-gc-ns.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-bs_not_gridcell.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-ns_not_inhibited-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-ns_not_inhibited.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-gridcell.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/response-probability-nsi-ns.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-gridcell.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/response-probability-gc-ns-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-gridcell.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/response-probability-gc-ns.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-gridcell.svg',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-ns_not_inhibited-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-ns_inhibited-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-gridcell-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-ns_inhibited.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-ns_not_inhibited.svg',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-gridcell-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-ns_inhibited-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-bs_not_gridcell.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-bs_not_gridcell.svg',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-bs_not_gridcell-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-bs_not_gridcell-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-gridcell.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/response-probability-gc-ns-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-bs_not_gridcell.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-ns_inhibited.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-gridcell-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-gridcell-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-ns_inhibited-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-ns_not_inhibited.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-ns_not_inhibited-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-bs_not_gridcell-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-ns_not_inhibited-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-ns_inhibited-stim-mec.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-gridcell.svg',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_e_peak-ns_not_inhibited.svg',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-bs_not_gridcell-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_e_peak-ns_not_inhibited.png',\n",
2019-12-13 10:43:57 +00:00
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-p_i_peak-ns_inhibited-stim-mec.png',\n",
" '/media/storage/expipe/septum-mec/actions/stimulus-response/data/figures/histogram-t_i_peak-ns_inhibited.png']"
]
},
2021-03-10 12:30:50 +00:00
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"copy_tree(output_path, str(action.data_path()))"
]
},
{
"cell_type": "code",
2021-03-10 12:30:50 +00:00
"execution_count": 21,
"metadata": {},
"outputs": [],
"source": [
"septum_mec.analysis.registration.store_notebook(action, \"20_stimulus-spike-response.ipynb\")"
]
},
2021-03-10 12:30:50 +00:00
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
2019-10-17 17:50:31 +00:00
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.6.8"
}
},
"nbformat": 4,
2021-03-10 12:30:50 +00:00
"nbformat_minor": 4
2019-10-17 17:50:31 +00:00
}