Skip to content

Instantly share code, notes, and snippets.

@vals
Created March 11, 2022 03:35
Show Gist options
  • Save vals/d7e72332d3fe98fabb09d602ac5ee2e7 to your computer and use it in GitHub Desktop.
Save vals/d7e72332d3fe98fabb09d602ac5ee2e7 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import pandas as pd\n",
"import plotnine as p\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"data = pd.read_csv('Single cell studies database - Data.tsv', sep = '\\t')"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [],
"source": [
"fdata = (\n",
" data\n",
" .query('not Organism.isnull()', engine = 'python')\n",
" .query('`Organism`.str.contains(\"Human\")', engine = 'python')\n",
" .query('`Number individuals` > 0')\n",
")\n",
"\n",
"fdata['fdate'] = pd.to_datetime(fdata['Date'], format = '%Y%m%d')"
]
},
{
"cell_type": "code",
"execution_count": 43,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"C:\\Users\\vale\\Anaconda3\\lib\\site-packages\\plotnine\\utils.py:1246: FutureWarning: is_categorical is deprecated and will be removed in a future version. Use is_categorical_dtype instead\n",
" if pdtypes.is_categorical(arr):\n",
"C:\\Users\\vale\\Anaconda3\\lib\\site-packages\\plotnine\\stats\\smoothers.py:168: PlotnineWarning: Confidence intervals are not yet implementedfor lowess smoothings.\n",
" \"for lowess smoothings.\", PlotnineWarning)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA1sAAAGYCAYAAAC03/g2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdeXhcZfn/8fedZLpQQje2lh2KlB1EERCqgKiAVBCKyI4gCgjKjohSBAvI8hXlB8q+yKIgS0C2AkIr+44KZSkUpClQSmlDaZtJcv/+eM7Q6WSSTCYzPWdOP6/rmivJmTPPuWfuM5PzzLOZuyMiIiIiIiKVVRd3ACIiIiIiImmkypaIiIiIiEgVqLIlIiIiIiJSBapsiYiIiIiIVIEqWyIiIiIiIlWgypaIiIiIiEgVqLIlIiIiIiJSBapsiYiIiIiIVIEqWyIiIiIiIlWgypaIVJSZjTczN7M1446lXGa2h5m9aGafRc/l693se42Z+RIMr6al4fwox9L6vCUZzOzr0fl3cHfbRKTyVNkSqQF5/xTdzI7tYp9pZvbmko4tbcxsXeCvwFzgaOAA4NVYgxJZQszs4LzPmtytxcyeMbNjzKy+yGMeifZrNrNB3ZS5fzfHfSXa58pKP6ekM7PdzWx83HGISHWosiVSe04zsyFxB5Fi2wMZ4GfufqW7/8XdP4g7KKl5ZwEDgXfiDqRElxK+aDgQOBNYBrgIuLibx4wAju/tgczsq8D6wJvA982ssdfR1rbdgdNjOO4kwjl5fQzHFllqqLIlUlueAYYBp8YdSFKYWYOZ9a9gkStHP2dXsEypUZU6v9y9zd0XuHutdDl9PPqi4Xp3/x3wFeA94HAzW77I/u3AC8CJZrZiL491GPAJoXI3CNinD3FLidy9Izon2+OORSTNVNkSqS1NwL+Ao81s9Z52jroWPlJk+5pRl53xeds+779vZj+OuvUsMLPXzezAaJ9VzOxmM5tlZvPM7E4zG9HF4QeZ2e/NbHpUzgtmNq6LODc3s1vN7EMzazWzt8zsHDNbpmC/a6IYh5vZZWY2A1gIbN3D61AXdYF6yczmm9lcM3vYzHYqfE2AM6JNb0fHmtZd2XmPH2xml5rZB2a20MyeM7NvFuzT5RiJYmN68p7vMDO7Mnp9PjWzB6LujpjZ2KiL12dRN65fFin7m2Z2k5lNzXv+k8xstyL75o7Z4/Pp4fVoNLOLzGxGdMznu8p/tP/o6NzKHe8tMzvfzJYr2C/XJW0HMzvBzN7I27+rLrbfNrN/Rs97voXxeEeZmXXx3IueX1E5D0d5WBg9t/sstMz09HoUy29u23pmNsHM3o3KfdXM9u2pzLxy9jOzJ8zs4+i99p6Z3W5mowv2W8HMLsx7zT4ys8lm1mPlxt0/BZ4iXDeM6mK3E4FlgfG9iH05YBxwk7s/CbxIqHyVzMyGmtl50fOab2ZzLHx+XVhk3+3N7F4zm533Wp9sxbtH7mRmT0ZlzjSzq8xs+Shn15QQl5nZ0RY+++ZY+MycFr0XV4r2mQYcFP2e33Xz4Ghbl2NCrevP92PM7LXo+b1tZr8CGorsV/TzKIr7R2b2dBTzPDN73Mx2L1JG2e8JkaVFpzefiCTeicAThG5JB1ah/KOA5YHLgRbChc+1ZpYFzgYmA6cBo6N9rwN2KlLOddHPC4H+wCHA38zsEHe/JreTmX0buAP4H/BH4ANgU+A44Ktmtr27txWU/SDwMXAO4eLv/R6e0zWEb83/BfwCaIye1/1mdqC7/wWYGe3zPWAP4FjgI+DTHsrOuZ/QGvZbQpernwN3mdm67v5uiWV0V/b7hIrgCMJr80B0EXUB8CfgakKLwFlmNs3db8h7/MGEnF4HTAdWJFzgNZnZPu7+10o+HzNrAO4FvgrcBjwMrAFcBbxeZP/NCF2a6gnd194CtiV0SfuGmW3j7p8VPGwCIY9XEXJ0IHChmc1w95vzyj6UcC5PA86L9t2L0B1uU+DwIk+h0/llZmOAu4FXonJmEVpBt43Keay716QH1wIdhG56HcCRwA1m9lZUAemSme0H/CU6/hnR81sF2BFYF5gS7bd6tM9I4MboWP2AzYHvADcXll3EOtHPWcXudPeHzOx+4Edm9nt375TrIvYltGZdHf19NXCRmW3s7v8u4fEAfwO+TsjzC4TPm1HADvk7mdkPgSuifc4htKZ9lfC5tjl5LWpmtitwJ+Fz4VzCe2EP4L4SY4LQA+Es4B/RcbPAasDOhHPnA8L76jhgO8LnT87jvTjO58zsHOBk4FnCZ91A4IfAd3tRzNWE99OdQO5z5HvA7WZ2hLv/KTpWNd8TIunh7rrpplvCb4QLCQdOi/6+hdBtZ7O8faYBbxY8bhrwSJHy1ozKG1/kGDOAoXnbVwQWEC4CTyoo5/fRY9bL2zY+2vYs0D9v+1BCN6RPgGWjbQOi4z2Vv290355ROQflbbsm2nYTYCW+djtEj7kLqM/bvgLhYmd2Lp6C+NcssfxcTH8u2L5VtH1Ckdf44CLldDpuXtl/Ktj3mGh7S8H+/QmVsicK9h9U5HjLECo+/y33+XTzmvww2vf3Bdu3js6jwuc5Kdq+VcH+vybvvI+2HRxte6ng/BpEqBw/nrdtcPQavQcMy9veADwQlbNtKecX4UsDB1Ys8z1cLL+5bfcAdXnbVwNagRtLKPc2wmQuDT3sd3d0rO8WuS//2LnX9yeECvoKwCbAn6PtTxR5/CNAW/T7JoTPptuKlLl/kcc+C/wn7+/hhNbEi0p8XQdH584lPey3MjAfuL1Ibo+P4vta7vUA3o7OndXy9qsnVJwcuKaE2J4HXilhv2sAL+O+aeR9vhMqmO0UfJ4Sup43U/DZQ5HPI8L4MQeOLXK8JmAO0FiJ94Ruui0tN3UjFKlNvyD8U/1dFcq+yt0/H6/k7h8CrxH+qf6hYN9Ho5/rFinnAndfmFfObOASwsXRjtHmbxAugq4BGqMuOstbGBMyCZgHfKtI2ee6e6ljX/aMfp7leWMT3H0m8P+AIXnx9MX5+X94aJH4FPhCBcq+oODvydHPJneflnfMhYQLrcXy4e7zcr+b2SAzG06obD0MbGDFJyToy/PJveYTCsp4Angof5uZrUD4Vv9+79yKcz7hHNiTzi4uOL/mEVp88+P7JqFb2x/d/eO8fdsILQ75seYrdn7l3hN7m1mmyGP64v/cvSMvvv8R3nOlvNazCRXNsWZW9H+6mQ0jtKY84u53Ft6ff+w8lxJadT4kVGwPB26lhxYSd3+Z0IK6h5lt092+UYvmFixq1cLdZxG+GNnfShsr9xmhcraVma3dzX57Eb7cuQIYXvBZc3e0T+6zZgvCF1LXRbnIxdZOwTndg9nAqmb2tV48pi/2IFQUzy94b3xM+Kwrxf6ESulf81+j6HW6A1iORd22q/meEEkNVbZEapC7v0n4pnkn68U4mhK9VWTbbKDZ3RcU2Q7h2+hCr3SzLTfmY/3o5yWEC7v824eEi8iVipRTSveknNwF2H+L3JfrprROkft6q9jrNovir01fy57dxfbcfYsd08J4tOvNbBahwvQR4TX+cbTL0BKOCaU/n3WAj6KKeqHCPHSZHw9dB6dSPD+lxFdu7oudXxcDTxO6us42s4lm9kszW6vIvr3Vl9f6LMIsfn8HPjKzu8zs2NyYoMgowv/7F3oR0zmE7sG7Er7cmUN4LxZ+BhTzK8IF+3k97Pcjwpc4T5rZqNyN8CXAMELXNeDz8WYr598A3D0L/JTwWTI1Gqt0hZl9r2AcVu6z5m46f9ZMie7LvWa5c6LYZ1ixc6krpxBaxx4xs/fN7BYz+4mZDe5FGb1RibjXJ3Q9nE7n1yk3LX/udarme0IkNTRmS6R2/YbQr/53ZvZgF/t01frT3Xu/q5mpupuxyrq5r5hcXLkvfH5J+KddTKdZAb3z+J3u9Da2snjXM3rlH7+71rguc9JN2T3OImZmyxJawpYjjNN5mdDtrIPQ3e8HFPnircTnU47Cx+f+7u0sfaXMoNZd2V0er9j55e6zzWxrYBtCJWRbwnTdp5vZAV583Fupyn6t3f1tM9uI0CVsR0Ir4fnAmWa2s7tPLqWcIv7r7rnPlXvM7Engn4QWkgO6fhi4+3tmdhFwipkVaznEzAYSxmsZYSxlMYcRunRCmIl1jcJiouNdaWZ3AbsAYwj5ORR42sy+Fn1JlDvHD6Pr6febu3teveXuz1iYyGYnwpISXyO0sP3GzLZz99dKKaab+3pzDVfqOVBHqFjv1c0+/4WqvydEUkOVLZEa5e4zzex3hG+2u7r4+ZjwDXGh7rrbVMoGhO5HhdsgtFbAohaEBXkXdpX2JqF70IaELnb5NiqIp9pyXdmWZE52AFYFDnX3q/LvMLMfVemYU4H1zGzFIq1bGxTZFxbl4nPRBfnahByWI/e4jQhjbfJtXHD8HkXd7f4V3TCzNQjjciYQFsKORdS6MzG65brnPUMYE7Yj8Aahcr15H47xiJn9ldC975KoS2h3ziZUbM6meHfnvQhdeM+i8+cEhFatfcxsbXd/C9iP0OLSVXwfErojX2NmFh3zBMJMh9ez6LNmdgmfNblzovBchfA5UrKo4n5ndMPMvkPoJnkSoUII3VeoPo4eNyy/K2z03hjB4u+N/LgLW7KKPZdiXidMfvRC1KWzW0l9T4gkiboRitS2CwndPc4kTI5Q6DVgtJmtktsQjevo9cKjZTg+f8yFhYWYjyS0quTG7dxPmKTixFy3oHwW1jgqVjHpjduin6fmj2mJxiAcRZiw46FiD6yCtwkzkn0jf2P07fceVTpmrtWkcJrzTQmD4avh89e84JhbUzA+Lho7Nxn4lpltWVDO8YQxV38vM46JhG6TR5nZ510lo+5luSnySyo7GltW6F1C96q+nqNl6yKuVwjd+IbB52N27gW+bmadxlx1NdariDMIlbbf9rSju88lVKTWpfiMjz8ijLc6291vLbwRPtuMqELi7o+5+4P5tyj2ZaxgiYhovN3z0Z+53PyN0AVyfNTauxgzG5g3dvF5QuvXgWa2Wt4+dfRijcMucvNcQVwQzXjaxWddrvXrGwXbj6fzNdwdhIrbCWbWLy+OYYTPulJcG/38XVRpXUx+99SkvidEkkYtWyI1zN3nm9nphEHf0Plb+j8Quok9bGZ/Ily8dLnWUYUZ8LiZ3UiYYvqHhCmpD/OwZg/u/pmZHUD41vdVM7uKcHHRSBhn8j3CNMbXlBuEuz9sZtcTWv/+aWa3Ey7gDyPMtHhgLp5qc/dPo+f446iV4GFgdcLMb/8GvlyFwz5GmPHxgmgCgWmEb7l/FB1ziyoc81rCRfLPoovVhwhdwI4kjBv6YsH+xxAmRHnYzPKnft+X0OrRab2kUrj7HDP7OWFK8Gej134eoVXlq8Dl7t5VF7ZCl1mYPv1+woV4AzAWWA/4v3Liq5D7zayFMFnNu4Rzex/Ce+jqvP2OIrRs3WZmNxBaeRuibQ2EiRG65e6vmtnfCC1OO7j7wz085BLgaArOazP7AqG74y1ddQl292fN7G3gYDP7dTfdWr8ATDKzO4D/EMYjrkN4T80hzD6Iu083sx8Tlgp4zcyuJZxnwwjjlPYgfPnwiLu3m9kx0WOfMrM/E7ozf4/w+pbqVTN7ivBaTyeMwTuYUCG6Nm+/Jwnjzi4xs38QvpB5yt3fJnSj/C1wuZltSJht9GvAl6Lnmv+avWFmFxBa9B4zs5sIk4IcSugi2dWaiPll3GZmlxM+HzaJXtcPCEsGfIkw0UpuMoykvidEEkWVLZHadzVhrZZO3bDc/UkL6/CcRlgr5kPCTGHXsmhQeLUcSPiHfTzhImMK0GlNJ3efaGZfJAwm35sw+HouoVJwJZVpdTqY8I3yoYSB/1lCN6vD3f2BCpTfG7lppvciXJj8h7Dm1ZeoQmXL3T+xMInKucARhBbQlwkVmS9ShcqWu7dZWD/tt4TK/a7Aq4QK94YUVLbc/UUz+wqh29vBhBkrmwmVrN/0coxeYSxXmlkzodJ+CuFC8TVCJaDUGdogdEU7kFBpX4HQKvMGodXmym4eV22XEM6lwwjvszmElq293f2W3E7u/o6ZbUFo0duNUCGbS+hudnEvjncm4X16JuHLgi65e9bMTqVzd7LcosW39nCsWwnrCu5C6HpXzP8IXzZ9nbBe2DKELxfuBM7xvDXh3P06M3uNUBn5IaGi9TGh0nU+4X2R27fJzHYhjI39BWGii7ujeGb2EHfO+cC3CRXdoYTK0YvA0e6e/7l2E4vW+RpHaLE6BHjb3VvMbGfCjKQnEVrn7iOMTSu2FtdJhIrdUYTPuumECuYTRN1Me+Luh5vZw4QJdE4kdN/8gPBZdXTerkl9T4gkipU+e7KIiIjI0s3MHLjW3Q+OOxYRST6N2RIREREREakCVbZERERERESqQJUtERERERGRKtCYLRERERERkSpQy5aIiIiIiEgVqLIlIiIiIiJSBapsiYiIiIiIVIEqWyIiIiIiIlWgypaIiIiIiEgVqLIlIiIiIiJSBapsiYiIiIiIVIEqWymTzWZHZLPZ8dlsdkTcsVRKc3Nz/+bm5vHNzc39446lEpSj5FOOki9tOUpbfkA5Srq05QeUo1qQthyVQpWt9BkBnB79TIv+hOeUljemcpR8ylHypS1HacsPKEdJl7b8gHJUC9KWox6psiUiIiIiIlIFqmyJiIiIiIhUgSpbIiIiIiIiVaDKloiIiIiISBU0xB1ArYhmgqmFAYqjcz+z2WysgVTK0KFDB7W0tNDY2LhpNpudF3c8FaAcJZ9ylHypylEK8wPKUdKlKj+gHNWCNOUok8k8X8p+5u7VjiUVstnseMLsKSIiIiIishTLZDJWyn6qbJWoxlq2bgD2A6bEHEtFtLW1DWppaZnU2Ng4pqGhoaa/BYkoR8mnHCVfqnKUwvyAcpR0qcoPKEe1IE05KrVlS90IS5TJZGYAM+KOoyd5zcxTSj0Jkm7mzJnLAcyePfulkSNHzo07nr5SjpJPOUq+tOUobfkB5Sjp0pYfUI5qQdpyVApNkCEiIiIiIlIFqmyJiIiIiEgiWHCMmU01s2Yzu8zMBsUdV7nUjVBERERERJLi58DvWFRPOQhY1cx29RqcbEItWyIiIiIikhQnsHiDUD9gZ2D1eMLpG1W2REREREQkKZbpYntNdiVUZUtERERERJLiQaA17+8O4H1gajzh9I3GbJWoxtbZAq02nmTKUfIpR8mXqhylMD+gHCVdqvIDylEtKCVH995778Xjxo1b79NPP90YoKGh4ZNjjjnm2HPOOWfDJL0OpU7Hr0WNS5TNZscDp8cdh4iIiIhImnV0dDBlyhQWLlzI6NGjGThwYNwhdZLJZKyU/VTZKlGNtWxptfFkU46STzlKvlTlKIX5AeUo6VKVH1COakGaclRqy5a6EZYok8nMAGbEHUdPtNp48ilHyaccJV/acpS2/IBylHRpyw8oR7UgbTkqhSbIEBERERERqQJVtkRERERERKpAlS0REREREZEqUGVLRERERESkClTZEhERERERqQJVtkRERERERKpAlS0REREREZEqUGVLRERERESkClTZEhERERERqQJVtkRERERERKqgIe4AakU2mx0BjIg7jhKMzv3MZrOxBlIpQ4cOHdTS0kJjY+Om2Wx2XtzxVIBylHzKUfKlKkcpzA8oR0mXqvyAclQL0pSjTCbzfCn7mbtXO5ZUyGaz44HT445DRERERETilclkrJT9VNkqUY21bN0A7AdMiTmWimhraxvU0tIyqbGxcUxDQ0NNfwsSUY6STzlKvlTlKIX5AeUo6VKVH1COakGaclRqy5a6EZYok8nMAGbEHUdP8pqZp5R6EiTdzJkzlwOYPXv2SyNHjpwbdzx9pRwln3KUfGnLUdryA8pR0qUtP6Ac1YK05agUqmyJiIiIJJyZ7Q3sAWSB6919YswhiUgJVNkSERERSTAzOwk4mzCLtAP7m9l+7n5TvJGJSE809buIiIhIQpnZAGACi67ZLLpdFFtQIlIyVbZEREREkms4UF9k+/Jmpus4kYTTm1REREQkuT4A5hC6D+a0A6+5e0c8IYlIqVTZEhEREUkod28Dvg+0AguAhcA8YN844xKR0miCDBEREZEEc/f7zWxD4BtAG3CPuyd+ORoRUWVLREREJPHcfSowNe44RKR31I1QRERERASwYIiZZeKORdJBlS0RERERWeqZ2abAW8BsYL6Z/cbMLOawpMapsiUiIiIiSzUzGwo8CKwWbaoHfgEcGVtQkgqqbImIiIjI0m5rYCiLr2nWABwQTziSFqpsiYiIiMjSzrvYrrXMpE9U2RIRERGR1IsmvzjGzJ4xs6fN7Mi8MVmPAx8RptbPaQeuWdJxSrpo6vcSZbPZEcCIuOMowejcz2w2G2sglTJ06NBBLS0tNDY2bprNZufFHU8FKEfJpxwlX6pylML8gHKUdKnKD/Sco3XXXffIN95442AWdRX84gYbbLBhNpu9srW1lbPOOuvI884777z58+evbWatG2ywwWVPP/30s9ls9otL9IksstTlqJZkMpnnS9nP3LtqNZV82Wx2PHB63HGIiIiISO+0tbXR2NhIe3v7Ytv79+/P3LlzyZ90cMGCBfTv3x9NRCjdyWQyJZ0gatkq3Z+BpriDKMFo4AZgP2BKzLFURFtb26CWlpZJjY2NYxoaGmr6W5CIcpR8ylHypSpHKcwPKEdJl6r8QPc5euKJJ5Ztb29/tPAxCxcuZOrUqVuPGjWqNbdtwIABSyLcUixVOUorVbZKlMlkZgAz4o6jJ3nNzFNKbd5MupkzZy4HMHv27JdGjhw5N+54+ko5Sj7lKPnSlqO05QeUo6RLW36g+xztsMMOAFOBNVnUjbANeH399dd/cgmGWbKlLUdppQkyRERERGRpsDvwMWHiiw5gJrBnrBFJ6qllS0RERERSz93/Y2brAlsRpnp/0t2XitYViY8qWyIiIiKyVHD3OcD9ccchSw91IxQREREREakCVbZERERERESqQJUtERERERGRKlBlS0REREREpApU2RIREREREakCVbZERERERESqQJUtERERESnKzPqZ2QVmNs3M3jCz48zM4o5LpFZonS0RERER6cpVwDigX/T3ucAAYEJsEYnUELVsiYiIiEgnZjYM2I9FFS0IX9SfFE9EIrVHlS0RERERKWZQF9sHLtEoRGqYKlsiIiIiUkwz8C7QnrctC0yOJxyR2qPKloiIiIh04u7twK7Ah3mbpwD7xxNRMpnZd83sSTObYmYXm1lXLYKyFNIEGSIiIiIVZmb1UWWlprn7f8xsFLAR0Ab8292zMYeVGGa2O/B3FjVgrAVsZGY7uHtHfJFJUqhlS0RERKRCzGwvM5sJtEXTpW8bd0x95e6fufvT7v68KlqdnMbi19P9gK8Bm8UTTt9YsKmZbW9mK8YdTxqoZUtERESkAsxsO+CvLLr4Xg14wMw2cve34ousdpjZV4CdgIXALe4+Ld6IejSkyDYHBi/pQPrKzPoDtwG7EJ5Dq5n9wN1vjzey2qaWLREREZHK+AHhIjWnLrqNjSec2mJmBwGPA78EzgT+Y2ZbxBtVj/4JtBZsWwD8O4ZY+mo88I3odwP6Azeb2WqxRZQCqmyJiIiIVEY94SI1n0fbpRtm1ghcRrg2HUC40B8AXB1nXCU4Hngu+r2DUNHa090/ii+ksu3E4muqQTifk17hTTRVtkREREQq4w46V7YagPtiiKXWrEbnC/16YJ0YYimZu88FxgBbAd8G1nD3e+ONqmyzWbxlFkIOWmKIJTVU2RIRERGpgOgi+0gWdSubC+zh7v+NL6qa0UxoGcrnwHsxxNIr7t7m7k+5+0R3/7DnRyTW7wivea7C1UroDql11fpAlS0RERGRCnH3PwHLAiOAYe5+d8whlc3MtjCzq8zsNjP7qZlV7brR3T8BTiRc6LcRFk9uA35SrWPK4tx9IvAd4FngLcJkL9u7e+GYNOkFzUYoIiIiUkHR9Ojvxx1HX0RT1v+T0C2yHtgN+CLww2od090vNLPXCN3xFgJ/cfcXq3U86Sxqna3VbpCJFHtla+zYsT8FDgY2Bm5vamraJ+++jYArgE0INewjmpqaJufdvxehyXMl4DHgkKampulLLnoRERGRVDqPRbMpQrhmPMTMfufuU6p1UHf/B/CPapUvfWNm3yaMT5tNqAzPijmkxEtCN8Jm4Czg8vyNY8eOzQBNwO3AUOAc4M6xY8cOje4fTZih5nBgeeAN4MYlF7aIiIhIao2k83WiAyvHEIskgJn9llARPhk4lzA1/yrxRpV8sVe2mpqabmtqaroDKJwi8+vAMsB5TU1NC5uamv4CvA18L7r/AODepqamB5uamuYDvwK2GTt2bKJnrRERERGpAS8Qxk3l6wBeiyEWiZmZbQScyuJT8w8ntIBKNypS2TKzNc3sG2Y2rBLlRTYC/t3U1JQ/M82L0fbc/S/l7mhqavoYeDfvfhEREREpz1GEcWethLWjOoDD3X1GrFFJXNalc+U7A2wYQyw1pddjtszsAqDe3X8e/b0HcDPhBZ9tZt909+e6K6NEywJzCrZ9Agzu4f7GChy7k+bm5v6EWnyiDRkyZFB9fT3t7e2DZs6cuVzc8VRILqeNzc3NsQZSCcpR8ilHyZfCHKUqP6AcJV1P+Zk+fXrLhAkTtvnrX//67ba2tsZVV131qfvvv/8/zc3NSc5lzeToqquuGnjZZZdt3dbW1n+rrbZ69uKLL55ZuE+S3kMbbLDBR6+88kqmYHNbJpN5p5fnRM3kqCcjR46cW8p+5UyQsQfw67y/JwD3ELrxnUcYf7VzGeUW+hQoTN5gFi2s1tP9lfYL4PQqlV0xn3zySe7XSXHGUSWJX2ujFMpR8ilHyZfiHKUiP6AcJV0p+Tn11FM59dRTl0xAlZXoHM2YMYMrrriC6dOnU1dXx/3338/jjz/ONttss9h+SXoPTZw4keOPP55bbrkFd6e+vp4BAwY03HPPPbvRufGjFInOUYkKFzAvqpzK1ghCdz3MbB1gPWB/d/+Pmf0BuK6MMov5D3DS2LFj6/K6Em4GXJp3/6a5naOJM1aPtlfD2cCFVSq7Yk8fmBQAACAASURBVIYMGbJpfX39pPb29jGffPLJSz0/oiY0Et6Uq5KCVcyVo+RTjpIvhTlKVX5AOUq6FOYHaiRHW2211R1tbW3bApmOjg7a2tp83LhxLeeee+66+++//4LcfknL0a9+9SsmTpy4z7x5875UX1//yZe//OVr1l577d5WmmoiR5VUTmVrDrBi9PtOwMd53QZbgYG9KWzs2LENURwNQN3YsWMHAO3AI4Q+wsePHTv2D8CewNqE2QkB/gI8NXbs2B2AJ4DfAE80NTVNLeM59WjkyJELCWs+JFo2m50HUFdXN6/U5s2ky2tmbknDc1KOkk85Sr605Sht+QHlKOnSlh+onRy1tbVtQRh+k2PAcieffPLyJ5100iu5jUnL0ciRI/noo48uAy4rt4xayVEllTNBxiTgN2Z2FGHqxzvy7luPqNWrF04D5gO/BMZFv1/e1NSUBcYCexHGYv0S2D2aCIOmpqZXgUMJ63DNio69bxnPR0RERERkSfmki+2zl2gUskSU07J1LHA9Yd2r5wmVoJwDgMnFHtSVpqam8cD4Lu77N/CVbh57C3BLb45Xrmw2O4LQhTLpRud+ZrOFk8bUpqFDhw5qaWmhsbFx09y3PDVOOUo+5Sj5UpWjFOYHlKOkS1V+oHZyNGbMmCsmTZp0BtGYHzNrGzFixD3Tpk0bEV1v5ihHCZbJZJ4vZT9z94od1MyWAxa4e2vFCk2IbDY7nhqYIENEREREku3vf/87l1xyCfPnz2e33XbjxBNPpKGhnDYQiUsmkylpgoyKVrbSrMZatm4A9gOmxBxLRbS1tQ1qaWmZ1NjYOKahoaGmvwWJKEfJpxwlX6pylML8gHKUdKnKDyhHtSBNOSq1ZaukKrSZNfXi2O7u3+3F/jUhk8nMABK/kF9eM/OUUk+CpMutLTF79uyX0jCYUjlKPuUo+dKWo7TlB5SjpEtbfkA5qgVpy1EpSm2vXA5QE5iIiIiIiEiJSqpsufvXqxyHiIiIiIhIqmgkXolqbMwWaOaaJFOOkk85Sr5U5SiF+QHlKOlSlR9QjmpBmnJU9dkIzWwU8AVgQOF97n5bWYUmmGYjFBERERERqOJshNH07rcB2+c2RT8/L8jd63tVaA2osZYtzVyTbMpR8ilHyZeqHKUwP6AcJV2q8gPKUS1IU44qOhthgXMJlY7tgH8BexBWvN4f2AH4QRllJp5mI4xP2mauUY6STzlKvrTlKG35AeUo6dKWH1COakHaclSKujIe823gt8BT0d/N7j7J3Q8H7gCOr1RwIiIiIiIitaqcytaKwP/cvR2YBwzPu+9eQmVMRERERERkqVZON8L/ActHv78BjAXui/7eBlhQgbgSp8bGbIFmrkky5Sj5lKPkS1WOUpgfUI6SLlX5AeWoFqQpR1WbjdDM/gjUuftRZnYAcC3wAtAKbAlc4O4n9TLexNNshCIiIiKypM2dO5d+/foxYECnCcAlRtWcjXAZYBl3/yj6ew9gL2AgMBH4s7t39C7c5Kuxli3NXJNsylHyKUfJl6ocpTA/oBwlXaryA+nL0UMPPbTNqaee+scXXngBgBEjRvzj0UcfPXPNNdes2WauNOWo6utsSTJls9kvAs8BW6Rl5prm5ublgDnA4DTMXKMcJZ9ylHxpy1Ha8gPKUdKlLT+QrhyZWf3AgQNfz2aza7e1teU2twKXu/tPYwytT9KUo1KVM0GGiIiIiIjkMbMNzOxhM/ufmT1iZhv1obi158+fn1/RAuhHSpdYSrNeT5BhZh3kLWBcTBoXNRYRERERKcbMVgOeAJYhXF+PAJ4ws43dfVoZRXY1JCd1Q3XSrpzZCE+ic2VrGLATsBLwx74GJSIiIiJSQw4E+rPo2rqe0BJ1EHBGGeW9vcwyy7yWzWbXy5uJMEuYmE5qSK8rW+5+fhd3/dLM/gIs17eQkqnGJsgATROaZMpR8ilHyZeqHKUwP6AcJV2q8gPx5mjUqFGjpk6dWpc/F4KZ2VprrbVOND6uV1pbW3niiScuOu200y6ZPHmyAx2rrbbaHY899tjfyimvXKeeeuq6119//UELFy4cvtJKK714++23Xz1q1KjWcstL0/solgkyzOybwHXuvnLFCk0ITf0uIiIiIsXcd9997L777nR0LOrlV1dXx6233sp3vvOdPpXd2tpKQ0MDdXVLdqqFF198kW233Za2tjY6OjrIZDJss8023HfffdTXa8RQ1aZ+77Yws58Cp7v7ChUrNCFqrGVLU7kmm3KUfMpR8qUqRynMDyhHSZeq/ED8Odp0000Pe/XVV4/I/T169Og/v/zyy5f1ochYczRy5MjzP/roo6+x+IR6vt9++/3k6quvfracMuPOUSVVc1Hj44ps7gesD4wDbnT3w3pVqFSMpnJNPuUo+ZSj5EtbjtKWH1COki5t+YFk5MjM1gDWBN4pc2KMz8WdIzN7CtiyYHMrcIC7/62cMpOQoyWtnAkyio3ZWgi8B1wEnNmniEREREREapC7vwO8E3ccFfIksBmhUSWnAfh3POHUpnImyNDaXCIiIilmZtsBuxOmmf6buz8Tc0gisuSdBmwHbAy0EWZb/Lm7vxprVDWmnJYtERERSSkzO4AwvXRuNdXjzGx3d78rxrBEZAlz9xYz2wr4FjAceN7dX445rJpTUmXLzA7sTaHufl154YiIiEhczKwe+DNgQCbvriuBFWMJSkRi4+6tgL5o6YNSW7auKfg7N6uGFdkGkLrKVo3NRghaNyPJlKPkU46SL1U5SkJ+Jk+evOyIESNOmDFjxsAid6/w/vvvf3n48OHtvShSOUq2VOUHes5RNptl++23H/fqq6/uArDhhhve9dBDD92WyWQ6F5YMS12OaklFZyM0s+F5f44CbgZuBG4FPgBWIsxE+ANgH3d/qrcBJ53W2RIRkbRauHAh22yzDVOmTKHYRd2IESN45520jPmXpdUZZ5zBueeeS1tb6CFbX1/P+PHjOfnkk2OOLFnmzJnDKaecwuOPP84KK6zAr3/9a8aMGVPy492dq666iiuuuIJsNsu4ceM44YQTUrc2V9XW2TKzicDD7n52kftOBXZ09x17VWgNqLGWLa2bkWzKUfIpR8mXqhzFnZ8999xzzF133XU+UHg11GFmHXvvvffR119//dO9LFY5SrZU5Qe6z9GsWbPqR4wY8QQF57iZtX766adbJ7R1a4nn6L333mvYeOONr/3ss8/WcfcMoeeaH3zwwYdfdtllL5RSxpZbbrnfiy+++HMWrc/Vtvrqq9/25ptvnpum91E119maB+zu7hOL3PdN4HZ3H9SrQqVi4l6ToRrStiaDcpR8ylHypS1HcefHzA4GLgUGFNw1Ffiuu/+3t2UqR8mWtvxA9zkys8HAJ108dKC7L6h2fL0VR47MbCfgPgoWMgbudfddSyzjE2Bwkbsap0+fXkeK3kelKGca9w+B73dx3z7AzPLDERERkRg8y+Jr6UBYvPQv5VS0RJLG3ecAbwP54w7bgSlJrGjFaAiLZiLNMcJshD0yMwOW7eLu5foQV80qZ+r3CcCfzWwd4A5C5WtFYA9gDPDjyoUnIiIi1ebu/zGznwIXEy5AG4BHgU5DBkRq2O7AQ8DQ6O9ZwJ7xhZNIz9G5MaYVeKSUB7u7m9kzwBYsmtG0A3g/unVVEUutchY1vtzMZgC/BM6LymgDnid0NdD0kCIiIjXG3S81sweBjQi9VB53946YwxKpGHd/2cy+AGxD6Br3WNTiJRF3fyvqVnwNoUWrHpgMjO9FMT8A/gmsQXidZwO7uXtHc3NzJcOtCWUtauzudwN3m1kdsAIwUx/IIiIitc3d3wDeiDsOkWpx99nAP+KOI8nc/QYzmwRsTGj9e6Y31/nuPs3MNgS+QqhrPOPuXY2XS72yKls50Qv/QYViERERERGpeWa2FrAl0AL8093nxxxSr7j7/4D/9eHxnxFat5Z6JVW2zOwPwPnu/m70e3fc3X/W99BERERERGqLme0J3EToQlcPvGFmY9xdk8gthUpt2doNuBJ4FxhLOHm64oAqWyIiIiKyVDGzFQhrY+Uv3LU2YfKZrmbzlhQrqbLl7mvl/b5m1aJJsBpb1BhgdDabjTWQShk6dOiglpYWGhsbN81mszW9AF5EOUo+5Sj5UpWjFOYHlKOkS1V+IBk5OvDAAze/7rrr+hds7tevX78x0bpZvaEcJVg1FzUevDTO3JLNZscDp8cdh4iIiIgk0yuvvMJmm23Wafvmm2/OU089VVIZs2bN4te//jUvv/wya6yxBuPHj2fUqFGVDlX6KJPJWCn7lVPZWgDcQ+iLetfSshBcjbVs3QDsB0yJOZaKaGtrG9TS0jKpsbFxTENDQ01/CxJRjpJPOUq+VOUohfkB5ShRDjjggC3/+9//bjRo0KA5EyZMuH+77bZblRTlB5KRo2w2y5prrnneRx99NMbdG4iG3uy5555H33TTTU/09Pinn356mR133PGm1tbWldw909DQgJnNP++888YdeeSRM6r+BKosCTmqlGq2bB1N6HO6NTAPuBO4EXjA3du7e6xUX9RE/RywRaknQdI1NzcvB8wBBo8cOXJu3PH0lXKUfMpR8qUtR0nNj5kNAtYBPnL3Xi2Qk5+jfv36vRqVMyea5azmJDVHpTCzM4FTCYvT1gHv33jjjYfutddeE0nJewiSkyMzywCnAN8E5gIXuPvDJT72UOASoF/etjZ3v8jdT6hGvEtSUnK0JBWuEN0jd/+ju28LrAmcBWxIWK/gfTO71MzGVDZEERERWdLMbGfC8i4vAdPN7E/R+pq9cvjhh28ONAP/Bt41s5vNrF8PD5MKMbNNgF8SrvkGEC7iVz7++OOPiTWwFHP3rLuf6e7bufuupVa0IsOBxda0cvd6wrq2UoN6/aGZ4+7/c/ffufsXCV0FLibMVNibE0pEREQSxszWAG4HBuVt/iG9nG141qxZXHfddRcBg/M27wH8us9BLgXMbFsze8LM3jazB8zsUTObZmYPRYvGlmI9oHB2hX5z585dt8LhSmU8R16rFoCZtQPPxBOO9FXZla0cM1uR0Ez6LcKYpqVu8gwREZGU2QYoHPydISwFU7Lnn3+ejo6OgQVl9SN8OSvdMLMvExaF3ZLQm2gnYAywRvTzCTNbs4Si/sfi05ADtC2zzDLTKxWrVI67PwScDbiZtdbV1TF8+PDJwKUxhyZlKquyZWZDzOxQM5sITAfOAd4BvgusXMH4REREZMlbQPFrhF4NaB8wYAB0rrQBzC8jpqXNTwk5KJaHBqA/cHAJ5TxFmNQsS5isIQt8duyxx/6xMmFKpbn7acAW22+//RkTJ07knXfeOUHzItSuUhc1/pyZNRFasgx4ADgQuNPdP6twbCIiIhKPBwnjrEawqFXECQP3S7blllsycODAt+fPn78qi7pGdQAXVSrQFBtK91+KG9DYUyHu7mZ2APAooZXsY+CSE044YVhFopSqcPcXstmsAb+NOxbpm3JathqBo4GV3X03d79JFS0REZH0cPcWQle1Z4CFhCECFwH39aac/v37c+WVV/4YmAh8ArwL/Mjdb6xsxKn0MJ3HWuWrByaXUpC7d7j7Ze5+mLuf5O7TKhFg0plZnZkdZGZ/MLNfRUNfRJaoXrdsufv21QhEREREEiUDbEBoQRlE+KI1Q+jeVrK99trr4x/84AffqXx4qfdHYAtgf0KrohO+JHdCTs529zviCy/ZzMyA64G9o00dwE/N7IvurvFqssSUVNkysy8Cr7r7/Oj3brl7KtZrEBERWYqdByzL4tcKR5nZ5e7+UkwxLTXcvd3MDiRMlrAyYeHh/sBawDvuPjXO+GrAdsA+LN6LayhwBnBYLBHJUqnUlq1nga2Ap6Pfu1oJ2aL76vsemoiIiMRoFJ2vE7KEmfFU2VoC3N2BV6JbztsxhVNr1iKcr/3ztmWAL8QTjiytSq1sbc+iN/pS2Y0wm82OIAwUTrrRuZ/ZbHddvWvH0KFDB7W0tNDY2LhpNpvt1UxYCaUcJZ9ylHypylES87PCCiu8P3PmzNEsfq2Q+cUvftEvm8322MsF5SjpUpUfWDxHhx12mF1xxRWF61W1rbLKKjNLPH+TINU5qvX3USaTKaknn4UvTaQn2Wx2PHB63HGIiIgsCe+++y5bb701c+fOBSCbzXLccccxYcKEXpXT0dHB+eefz5VXXkl7ezt77703Z5xxBplM4dJPIpV1xBFHcM0119DQ0IC7s8IKK/Dkk0+y0korxR2apEAmkym2rEUnqmyVqMZatm4A9iP07655bW1tg1paWiY1NjaOaWhoqOlvQSLKUfIpR8mXqhwlNT+33Xbb0LPPPvs78+fPX2706NEv33rrrSXNfhcZDdyw44473jF58uTvELWQmVnbiBEj/jFt2rTfVCXoKklqjvogVe8h6JyjbDbLnnvuucPUqVPXX3bZZWdPmDDhrh133LEl7jh7IfU5ijuevqhoy5aZvU3X47Q6cfe1S91XKitqGn8O2KLUkyDpmpublyNMOzx45MiRc+OOp6+Uo+RTjpIvbTlKW34g5MjdnxswYMBCd+9fZJdl3b1mLrbSlqO0vYdAOaoFactRKUods/UPFq9s7Q4MIawB8QGwErADMBvQNKQiIiJCR0cH7t5Vf8FlgJqpbImIlKOkypa7f76mhpmdALwHbOzun+RtHwrcE90nIiIiS7n6+noaGxtfbmlp2ZAwExxAG2FGvY/ii0xEZMmo63mXTn4GTMivaAG4+2zCWhDHVCIwERGRNLPgIDO7ZtSoUb979dVX4w6pKi688MJTgDfzNs0AdnMNGhepOWa2g5n9ycwuM7NvxR1PLSinsjUMGNzFfYMJC8aJiIhI9/4PuBI4cP78+YfuvPPObL/99lvEHVSlHXTQQTOBTYEtgC8DX3D31+KNSkR6y8z2BR4EfgQcCtxrZlogugflVLYeAs41s6/lbzSzrwPnRPeLiIhIF8xsbUJPkXrAgIb29namTp16bryRdWZmu5nZvWb2qJkdZ2a9vnZw96y7P+/uz7r7gmrEKUsfM1vOzH5vZo+Z2d922WWX9eOOKa3MzIBLCZ9XddHNgD+aWX2csSVdqRNk5Psx0AQ8bGZzgJnACoRWrReAn1QuPBERkVRapXBDR0cHZtZpe5zMbB/C1NO5CtY2wLrAEbEFJQKYWX9gErA+0A9of+mll3Z7/fXX+cIXvhBvcOnUH1iuyPYBhDrAx0s2nNpRzrdTM9z9y8CuwO8JMxL+HtjV3b/k7s0VjlFERCRt3gTa8zc0NDRQX1+ftIFb57D4tUID8BMz06qwErddgQ0JFS0IrcSZSy65JL6IUixqkZ4OdORvJsxEPjuWoGpEOS1bALj7fcB9FYxFRERkqeDuM8zsaOD/AQuA+sGDB/dbf/31j4s5tELDutg+nLD0i0hclifMbJl/LVs/c+bMmMJZKuwL3J/3dx2wjya76V7ZlS0zGwmsSmg+XIy7T+pLUCIiImnn7pea2fPAdkOGDGl75JFH/m/YsGHT4o6rwHPAV1k0bTtAC2Hq9qXSuHHj1nr88cePA9YgDJ8Y7+5zYg6rZNFSPeMHDx781bFjx7L55puPPPbYY2txwdznCV3b8mW32GKLrtZ1kz5y90lmthGhVbEOuFeT3fSs15WtaFDv9cBWuU0FuzihKVdERES64e5PAU81NzcvR5idMGkOJoyLGUnoPtQO7Onu8+MMKi7vvvsujz/++GRgIOEaamvgm2b2pVp4TcxsEPAEsNacOXP63XTTTdx44403HXfccRu6+7txx9cb7v6smZ0MnAu0ApmGhobHjzrqqK/18FDpA3efCvwh7jhqSTktW5cDqwOHA68QTnARERFJGXd/J/ome3tCK8IT7v5ezGHF5vLLL4dFFS0I44XWBb5HmEikS2a2LGFikdWB14E/u/uSvobaG1iLaJxTW1sbZjaAMDPm8Us4lj5z9/PM7D5gI+DD++6777n+/ftr/FAfmFk/wjX+aOBd4FJ3bynY5wvAgcCywMPu3rTEA60h5VS2tgQOcvfbKh2MiIiIJEt0oaWLKeDDDz+EztdObYRZmbtkZo3A08DahN4/7cD3zWx7d89WIdSurEDBxCzu3gCsuARjqCh3/zfwb4CohVjKZGYZYCKh91ruPD3UzL7s7nOjfbYAJkf31wFHm9l4dz8zprATr5x1tqZT8EYVERERSbvNNtsMoLBy1J8wdquoaH2ic4BRhBal+ujnlsB+VQm0a53GOZlZW7Sdgu2rm9l3zezr0UW4JJiZrWBm3zGzb0bdRcuxD6FrbP55uhbw87x9Lom29yN88VAHnGFme5vZ2KQtX5EE5VS2TgNOMbOuZigSERERSZ3DDjuMhoaGSYQvnRcQxqmf7e6PFts/WgD6auBIOreIdQBrVi/aztz9QeACwM2sta6ujqFDhz4F/DF/PzPbi7A8wS3AQ8DjZjZ4ScYqpTOzrQn5uh24B3jFzNYso6g16dygkgGOMLNl8vYpnJvBgJuAW4E3zWzXMo6dWuVUtg4izEI4zcwmmVlTwe3OCscoIiIiErtMJsO//vWvPYGdCeOvvuLup3XzkEPouvWqjnCBvES5+0nA1ttuu+2Ef/zjH7z11ls/c/e23P3RbNM3Ei6yM1GcmwDnL+lYpWdm1gDcCTQSKvT1wAjgL2UU9ybFJ7lbARgf/f4GxXu41RHOlwHALWqUWaScMVvLsviHQ2OFYkm0bDY7gnDyJt3o3M9sdkl2A6+eoUOHDmppaaGxsXHTbDY7L+54KkA5Sj7lKPlSlaMU5gdSmqNNNtlkk9bW1lnALIBsNvvFrh6z6qqr7vLee+8Vu3jtGDJkyHNvvPHGlO4eXy2tra1ZwkUzwHrZbPbzdZL22muvrW+99dbCmPsNGDDgG32J9c033+z3rW996/jp06fvDLDKKqvcc//99184atSoik0SksL3UY/voQsvvHDkcccdVzhmMGNmW/Y2XzNnznxjjTXWeOuzzz5br+Cu+mWXXXbnbDZ788knn/yH8847b/NorB/u3q9IUQMPO+ywsdls9uXCOwpz9L3vfe9rDzzwwInZbHbYwIEDpx199NG/PvPMM5f4lxDlyGQyJS2ZYFqHrDTZbHY8cHrccYiIiEhtOPXUU7nooovIv1Cuq6tjn3324YorrqChoezlTqvm2WefZZtttum0feutt+bRR4v2lizJIYccwi233EJra6hb9evX7/PXQcr3ySefsOKKnec3WX755Wlubu51eQ8++CC77ror+fWDuro6vv3tb3PHHXcAYQmEm2++mXnz5nHzzTfz9tudl9176aWXWH/99bs91uTJk9lpp53o6OgAoL6+nsbGRl5++WVWXnnlXse+pGUymcLlr4pK3rs8uf5MbczGNJow/ex+wJSYY6mItra2QS0tLZMaGxvHNDQ0pOWbKuUo2ZSj5EtVjlKYH1COWH311Vdub2//G6FrVb2ZtdXX1398wAEH7N3Q0NDS0+OrrGh+NthgAxsyZMglc+bM2dzdM4Rxaay22mpHAU+Vc6D33nuv4YYbbniSvLVZW1tbue666zrOPffcrYYPH96pW9rpp5++zgMPPPDlhoaG7BFHHPHovvvu+1FPxyknR88999zAk08+ece5c+cO3mSTTV654oorupzsJAY9voeGDBnC6quvfvK77777PRZd13esvfba5xLGUPXKRhtt1DBo0KAr582bt14u/x0dHR0jR448DHgZYPXVV+ekk07i7rvvHnLDDTf8BBiXe7yZtQ0bNuyxUaNGHVes/PwcHXDAASd0dHTsRnRetLe3M2fOnOz3v//9cx999NHbext7UpXUsmVmxwE3uPsH0e/dcXdP4sKMS4Woyfg5YItSmzeTLprKdQ4weOTIkXPjjqevlKPkU46SL205Slt+QDnKidYkOocwq9t/gJPcfUZ1oixdd/kxs4HAb4AdgNmESUAeKvdYUXmfdXH3QHdfULD/XsDNhJkfjTAZyXbRNO9d6m2OzGw48CRh7bN2QqV4vLv/pqfHLgmlvofMrJ6wTtqehPVv/+Tu3a771p1oTbizga8CHwK/cffHC/ZZm7BA9pBoUz/gHULl8Iyu1pDLz9Eqq6zyJ8IMiPktRAuAE9z9/5Ubf9KU2rJ1PvAv4AN6HiDpgCpbIiIistRz99cJix7XDHefD5xYyfLM7HHgy4RJFCBUpJ4sUtFaBriOMFFDbuxYQ7Rt80rFFDmbUNHKH3c03sxu76lilyTu3g78LrpVorxPgaN72O3PwDAWr0usClzUi8W67yYstJ0/RrABeLjEx9eEkmYjdPc6d3867/fubsUGgoqIiIhIzMxsFTP723LLLff3XXbZhVNPPXXdJXToPYkWH468TLjQLrQaMLBgWz2wqZk9ZmY7VzCmzVm8ogWhEtj9YCMB2JjOjTb1hPXkSnUTcCZRV1VC6+fe7v5q38NLjnKmfhcRERGRGmNmQwjd5nZfsGDBmv/85z+54IILrjWz0T09tq/c/X1Cy9Ya0W3LaFuhD1h08Z3PCAvu3m1mu1QorHcovq5U72eWWPq8R1grrlDJr50HZxBayNYDlnf31IzVylFlS0RERCpu1qxZ6ulSAdE6SpUyDliRqCtfe3s70RTeR/al0MIYu4rZ3Tvc/d3oVuxCHXf/BPg14UK+sNJl0e2XfYk3z2mEMUK56SKzwB3AYxUqvyoqfE6U6+eEimo7IU/twB/cfVpvC3L3T9z99ai7aRKeW0WpsiUiIiIVY2bfNbMPR4wY8fRaa63FIYcc8qW4Y6pFZvZtM2sGsmbWbGbfqkCxg+ncGlEPHG1mk6MFjXsT4zfM7L0oxvfN7Bdm9r/o7w/MbGw5Qbr7WcC+wEMUr3BVZMFcd58CbAZcTmjlqgd2A26OJolIFDP7iplNJby+s8xs37hicfd/AV8hvHY3AD8iVMDKYmb7mtkswnOb+qcJswAAIABJREFUamZfqUyk8VNlS0RERCrCzLYCbgNWAGhububGG2+82MwKF0mVbpjZpsBdQG6xoZUJ3ec27WPRTwD9u7jvK8BEM8t0cX9hjBsC9wC5CtpKwARglejvFYHbzWzLcgJ1978CBxBm18vXCvyznDK7OM6bwHBgBOG6uAHYnTAhR2KY2WqEyuea0aZhwF/M7BtxxeTuL7j7Ee5+gLtf7WUu3mtmOwJ/YVElek3gITNbvUKhxip1TXVpZ2YGZHox04uIiMiS8n1Cd6I6AHfHzBzYgzD9eXWF/5H9CBMsDASWyfv9823tMHAONA4L3en6RTG3EbqRtRW5ZVdcZZV+c84/n8GnnLId77wzt4f9c/ctAOavCfYOtPbiYnQvQgtU7jrNor/3Al4q56UBcPfHzOxE4DwWn24bwmuxQXQr5Rh7RjEVlpP/dzvhnHi6zHjfN7PvA38lvBb1wDPASeWUV0xUudybxePuB+xuZgMKZ0us0DH7AW1ddaXswq6E1yC/ocQJa3A9WMHw4rA/4bnkcpCr9O4KXBpXUJXSq8qWmQ0gvEGvz81OKEuOmf2M8K3RMmb2OjDO3V+OOSwREZGcejpffEOxnjRh7aWVotvKhG+1u6wg9fB3/rZix//c/wNOBuYB6wC3UNp84g3TpzP8Bz+AMF11SSYRroTfI/Tfu9bs04PgE+DT6NYS3Rb7fWfY+n6oy78SN7DVYHlCV785wGeU0ZLg7heY2R0NDQ3PtrW1DSmyS6m9nkoZk2e9KK8od7/TzNYCNiU876ejqc4rpavzpc+xdyowtE7dAmxJ6C73B+CUEp9PsfdWxWOMSR2dn5uTjufWu8qWuy8ws0MoY0Vq6Rsz2w+4gEUfbusAD5vZaHfvcVV1ERGRJeAO4KeLbXHPNMG6mN1E6HK2cnRbbkkH9zfgGBYNWnob2BGYQujzVklTgW8BC6O/5wCHwLKrwbI79PDYXwH3F2wzyPwVfkK4AbRjNjcqeg7Q3e8LolAWAgsdFo4ZNuzppz7++JutbW25Q7QPgI+mwWqYrUHobtiv4Ofnv98Kq4+Dft3V9gwa7oRdCN3EFhDquPMIXQE7olt73u8dBdtbgayHVsJs9Pd2mHWqoAKzB++//8I5EyZAfelzs/x/9u48Pqryevz452RmshKSEGQHUVwAEVwrorgX0OrYaq3Wbtq611qtu+1Pbd1r3Wu11VprtXZ1mbpWq3WprSgifBVRUUEgLBK2AFnuzJzfH8+NmeyTZCazcN6v17yS3Mzc+8yc3Mmc+zzPeVS1SUT+CRxMSxn4JuAVVe1sIeYe83vQngO2xyUWhcAP/bb/NIldvEP7MvWC+zPOdY/ghow2a060/pmZ5qRWb4YRvoYb1/tSittiuvZtWl9FCuD+UR0I/D0jLTLGGLP1ESnAzW/ZjpYy3mOAMQpjfgON50BxIzAAeADkKDiph0fxgPoObluS3Nbh9jPh+jhMbz5IHFgHTTvCKRvgUdxQumBHt6bJkys2/Pzn/60855xDQh980NjJ/T5//KVwhAdf09aftfRGWHAIzPFfnvKEr59/vy8UPgic4j+JUuBeYGrr1ygAVPm3Hnt+9WpOB+73f94ZAv+AoUPh8WQefyxubN9JfhuLcRUS7sNlU0XA3cBRsFNv2tcbZQ8+SNmDD6Ku2MXaNrfaDratBdZeB+dcBncpNOfBr+GGP6bSRFx580Qh3EuYTLI1AXdeJM6pU9z5l9P83ssf4joVQriLBMer6oeZbVlq9CbZuhw3IS8KPA2spk2lGFVdm4K2mdY6u0yTF12sxhhjsohIEe5D8o64pGr7hK9j6bzIAqcB38V9sh3M5/+8NuE+AC8HVuDWUlrZ5usamhMk1ShpsFZkU0fbN0IM1Q5/12xNTc1AgM9efHHOiBEjNnZ3rL+KbAMcQ8JnLQV9Bl5D9bQuHyxS+HUYsB9UvgijjoLYIDdEsjkpG4gbmVjR2fcfQ+UPoPxdXMBuxZXda1YI/A43rHILrkJEl+MvWzTien4aj4Omr0DjKogNgS0haLwZmpaDDofNha43qwmXJBQDZf6tEPf5pbNbwL+FaJlX1/x9cXcNFKjE3bZP5gld4t/Wwfo4rKl2x3kQkeZEfRPwMfAB8D7wEaqNne3Pn19/Bq6kfiGu5+aRTu7e6nOciHwbOB8oHTt27Bvz58+nuLi4+X4xWidb7R6fq1T1DhG5Gzec+LMezmfLar3t2QL4BW7+VkdsbY3Uexg4iJbXNo57f7QeRmOMMb0jUoq72j6xzW0Hkv8Qp7iFTJcAnwKfBuHTbcLhYOzyy2/VtWsPDM6Y8Upv5helwZ+AGbT+X9pICqvbJXgauIHWE/+Fzj90t3BFsNaOgbXfcR/ye0RcovcOLjEJLYHYHuD9AI66zcWqKHbaaZP1a1+7v/h3vzu99KGHPiQhiWrzfeJXr20cg7SUH2z+Oa1dLSIBXMKW2Cs4EKiMDRo0YvOZZ94x4JZbbi3YsmUA7oN721tpZ7uuaknSduimFXFEFuMSrw/afF0OnAf8nJa/s/Nxvb+fAKNp+fzt4cqm+09NTgF+jX/uLVmyZNsTTjiBv/71r+BGlt5E678nyKOpParq4S685JXeJFvfpeOVvU163Ycb434lLm4rgK+o6upMNsoYY0wOcGsGjad1QrULrrcqmQ6NVbgP/Z+0+boEWIb7kNRKzPP2AFDYlCWJFsAfcEMgr8L1EKwCjlHVFak+kKq+LyJH4RK8Klwvzw9U9ZlUH6sDx+N6uJp7QQIK8dvh8NtUzweIe54AxA444M2CBx98qx/alBqumMRG/9bKKtf7eMems866otPeR1fsraMkbLB/G+J/HUBL8ZVKYBStq+Vt798Ob3OELZVQtL51x0MIt27Y3rgOxUm4z9L3AVck3O/HJFzkUNXQU089xV133TVSVd8SkS/hph5W45Lfc1X1Hx2/UCZb9DjZUtX709AO0w2/XOw1InID7grOut6uZ2CMMSZPiQzEze1oTqaaE6tkOhsacbUiFiTcPgA+QXVzWtrbz/z/mzeIyE30w/9SVf2niFTjPsyvT3EVva4MpP2F8QCuF2jr5kq51/i35LnqmTvghtfu7N+av0+cN1faWa34t+FnU+BPc+H9LTB7P9VP29ylw/isXLmy1DVdX/B7Lfv778n0Qa/X2RKRKlxmPhp4WlXX+aXhm/JpnGW2UTeO3ebEGWPM1sx98JsM7ErrpGpUEo+uB94D3qV1YvUJW8mHt/78X+onc7X9cawEr9J+Xp0Cr/RzO/KHaj3wf/6thZufVU1C8jURvjsftmmeeFiAW99gF9cLdvjuLY9diSuWMht4vQj+0wizaKk6qFVVVXLiiScubWlGRv6eTB/0ONkSN1b2Klz11FLcybs3sA43Dvl1kquqYowxxpjuuHlVU4A9/dseuASru/nRm2idTDXflmAXRfOaqr4sIhfg5teDG/52J/Bg5lqVp1zys8a//QfgLZHbgH/hepkpgs2/h5eCbtjuzrQMFRyGW7j3S+Amex0Eje/4vwwVFNQ/+pe/lE6aNCnlCyub/tObnq2f4tbQuBA3oXRBwu8iuEqllmwZY4wxPVR2992lXHXVrrROrCbSdbGKDXScVC3NorlSpp+p6s0i8nfc0LelqvpBptu0tVDVFSKyG+78DdXDvC+qujlkImW4Huk9/Nte/s8F1cBcKJqLq4A2JR4vrTjySCgs/B11dS/gOjReBxbbuZ07epNsnQRcpqp3+b1ciT7CLbZrjDHG9CtxRSD+Hy5JWQpcm9XrtIiU46px7zm8sHBqdOxYgh98sJyuE6sa3LCj5tvbwHL74NVzIjIeuAy30PLrwDWpXMQ2G6jqElwRE9PP1FWU/F8Hv9jsb2/5nUvA9gT2CcI+e7v1bN2Q4MZGaGycDEzeCPwMeAOaBoqsugj+Mh2eHArbrHYFOEK46oT393QuoojsjutI2QY33PR6/zmk1H333VeyatUq7r777sej0Wjz+/TCVB8nm/Qm2arGjfXuSAHt6/8bY4wxaSVuXaiXccPrCoEo8FUR2V1VF2W0cdBcuGJ3Wnqs9sTN7xAAaWoi9MEH0DrRWkbrxOotVFf2Y6vzlp9ozcHNawrgFjo+REQO0A4qKxqTVi4Be9m/OSIj4kcc8XWdOPEXBffc82b9hg0T9oOyD4AmKAzC6Jfh/J/A+Z8l7EpgVrWrPHpxsocXkb1wSzs1r3F2ILCviHwplXUYRCQUCAT+UVBQQDQaPQi3btixIrJnPidcvUm2PgC+iBuL2tbBuHUdjDHGmP50FG4oTvP/tSBujaELcIuL9h+RCtzwoMTEaseuHqIiSxtmzRpd9MorVxVs2vQaLrGypT3S5zJcUt48QqcQN/98FmCltE3mqdbEPO9FgPi1155eWVg40XNl44PgribV47ryE7uwFArWwkWbRY4uc9N9XgJeouvlDa7CnQvNF3tCwExgX/x5aCkyKxaL7RGLfV6HJ4A79y7BjZzLS71Jtm4B7hERj5aF1EaJyL64ohknpahtxhhjTLKG4BYITfy/FsQNEUsfV5m3bWLV3XD6xbTpsVqxbFkTbu7VLzpdH8ik0gjafwbycH9HxmQdzw3va/Ue11kXbBxYDzuXuWIc7mKTyIc0J14u+Vqa8JDhtB++HCX150Pz+3TiNKS262LnnV6tsyUig3CL617mb34MN5fvJ6r6l9Q1zxhjjEnK27ierERNwBspO4JbL6ltYrVdN4/6mPZDAduXHHeLsZr+8zpu6GBhwrYi3N+R6SV/WaCf4Iq6LALuAE7FLRX0EW5eXMoXkN5KzKV9Of8mXDXwwbQkMFoM9UPgBdzfeIW/fUf/dgoAIh/jJ17VsKDWVU5MPB8CtC1z33fzOngOMaBaRJ7CnX/XquqmxDuIK69/Bq7neTNwt6q+TI7o1TpbfoWb3wD74eZwrQVe0+ZKK8YYY0w/UtXXROQq3KiaRtwwmP8BN/Zqh27h0LaJVXcLAy+idWI1F9V1vTq+SbdrcFMfvoC70l4E/FhV52S0VTlMXMGXN3DrrxbiXtczcKPcQv7PXxWRyWpDZHtMVf8tIjcCF9HyHvcyrqjF87QsZB1tgCNCqi/hCtlNwc3BOhA4gJYFmLf3bycvAqaC9yGoQDQGQYGz4yme76qqbxYXF9/X2Nj4vYTNAVyhIAEOBQ4Xkamq2phwn9txf0tB/zkeLyJHqepTqWxfuvR6UWM/63w2hW3pUDgcPhs3NHFX4NFIJHJCwu8mAffiFnb8GDgzEonYgn3GGLMVUtUrROQJ3P+LlcCzmswivSJDaZ9Yje7mUR/QPrHa0PvWm/6kqltE5EDcvJQhwDxLtPrs27gKes29I20LpoVwH/TPBi7vx3blDVW9REQew/Uc1gDPqWrML/hyGO5z/b9VdZn/gBjwln+7BZECXC/jgcBBuORrcCXwNoSeBdZBaE9gV/gxIvvRMvTww1RUPY3FYjuJSNtdif+1EFfk6CvAnwBEZATub6btfW8C8jfZEpHBwLnAVNw4zxW4K4i3qepnXT22F2qAq3F/RIObN4bD4RBuXa9f4/5ojgMeD4fD4yKRiF1JNMaYrZCqvkFXQwdFhtOyflVzYtXVfAEF3qdtuXUbyZHz/KqDT2S6HXlkWBL3CSZ5P9MJVW1dNt5tW4OfnHTz4Dgw37/dgRueNxE4sBgOPNp9nh7q33sk8A3/BrACkZdxRTeeBz7uTfIVj8eHd/OwKK3/RoZ2cB/pZHtW6nGyJSL7AM/guv1ewFUpGYIrjvEDEZmhqq+nqoGRSOQRgHA4vBsJyRYuIy8FboxEInHgwXA4fB5wDPDbVB3fGGNMDnIfIkbQPrEa3sWj4rilTd6idWK1qYvHJNscwQ2VGQa8p6qL+7i/QtwQuFLgLf/DluknIjIa10PwGTCnp2sa5al5tC580JGYf7+c4NcoOAbXI/eUqr7bzf0rcIsUR4HZqlqf/lb2gfu7fde//cp/39yZlmGHB9JSZGg4cLx/A1iMyPO4xOsFkuxsCQaDb6nqDp7X6QoLRbT+G/kYaKD1nFyPHJpf2ZuerTtxQfmSJgyZ8P/AngZ+iSufmm6TgP/zE61mb/vbjTHG5DtVJBKp5Ljj9gLGJtzG4da06urKZwxYQOvEap6/3k1KiUgQeBj4qn9cEZHvq+rdvdzfYNzyK7viet62+PMX/p2iJpsuTJ48+eu4z0LgkosnROTYdCwAm2P+DjwEfAs3p6gQt1bcaFwhhxDwHG5EUtbzq2z/m5ZhkT8XkV+o6oWd3H934J/AIFzPyycicoi/sHRucMnXQv/2az/5Gofr4GhOvpqHWI/FFdtoLrjxNi7xeh54tbP30nA4fOmbb775tcWLF3u498Ni3IWuJlyidYuqvpjQpA0i8k1cz10UVzWxFvheu51nqd4kW7sAx2mbsen+i3E98OeUtKx7A3BlchOtp6XqijHGmFzm/tFvQ+tEaltgbLCkZDwiyJYtHa352FYMtwZkYmI1H9UtaWh1R84Fjva/b77yf6eI/E9Ve3N19je4RUvFv5UBj4vIyLZVvExqLVq0iNra2l/Rukz2DOBS4KeZaVV2UFUVkZOBP+Cq3i3Bze0/ENdb8inwTCoXyU0Xf5H0f9K6Oh/ABSLynKr+s839A7ghqYNo+dsYjbvIMi3NzU0fl3wt8m/3JiRfX8RN7zkEqPTvvZt/uwBoQuS/uMTrOeBNf/4Yt91225otW7aw3377nbp69epC4E3cazUCeEdVX+2gGX8XkUm46opbgKc1h4oP9SbZWkTLC9tWBa67rz9swlVeaXv8unQcrKampoj25SqzTmVlZVkgECAWi5V99tln+VJKuLz5a01NTUYbkgoWo+xnMeofsn49lZddNiT43nujC9at21YaGsaI540hGh1DPD4GGCNQ0uFj6zsenaMiqyko+JRgcIGWlLwdGzLk7U1nnfVO/XHHNba5a7C/yq2HQqFDPM9rWyygqaqq6sCamprm/5lJx0dEpqtq4odAAQZOnDhx95qamqwZopWH51H5vHnzwF1dT3z9CwOBwGE1NTW3ZKZZvZOO+CxfvhzcnMnmeZMDaLnAATAgze8/KXmf23vvvXd44403BnT0u7KysoNrampazZmaNWvW8Geeeabtmn4hYO8FCxYMrKzs7GNz17LuHHLxXY3rwXwo+N57BVXf//5ugZqag6S+/iCi0X3FnRuFtPSEXaWwjqKif8XLy/9V/u1vv84FFzB37tznafnM/lHzIWo6eV9evnz5CuAv3d2vPyW7JmJvkq0LcVfklqrqS80bReQg3NpbZ3fyuFR7B7goHA4XJAwl3A24K03HuxS4Ik37Tpn169c3f5sz6w/0wLJMNyAVLEbZz2KUIvE4BZ99RmDpUoLLlhFYupTAsmXutnQpweXLkYaGHu0yNnQosVGjiI4aRWz0aGL+1+ioUcRGjoSSkiG4ecR74aqjZdyXv/xl/v73vxOPt1zQDwaDxddcc82twK1t7t5tfMaMGcOSJe1HJt17771Z9feaj+dReXk5gUCgMBZrKXJZUFDAzJkz96f9aJuslo/xSdCn97nf/OY37L777h3+7sILL7wEuCRx22233cbOO+/c7r4DBw4MVlZW9vrvIttjFJ0wgc9eeOHzn6W+nsI33qDwlVcoeuUVQu+8g6giUEVT01cDtbVfLb/lFoqffZaGgw9e1njQQTTttRcUtu1AzBnS/V2STLZE5P9w48KbVQAviMgG3OTQbfxt64AbcHO3UiIcDgf9dgaBgnA4XIwbEvJv3IS588Ph8O3Asbj1Ah5N1bHbuA64OU37TpnKysopgUDg5VgsdsD69euz5gpnH5Xj3jhHkaaey/5kMcp+FqPkBJYulYFXXjksuGjR6IL1613PVDTa0jMVj4+WHowIUFBEVlJQ8CmBwBINhT7VwsJPdcCApbHhw5dsOfHEZc09VJWVlVMKcyRGCxYs2DUej7+IG0JYADTF4/FlK1as2A83JAY6iM/69euZNWvWwRs3btyxtLR06R/+8IdnJ0yYEN+8efPxwN20DFdqCoVCz1VUVJzYr0+sG5WVlVPmzZv38sUXX3zze++9t+DMM8988owzzuivoZvpUH7AAQcsE5H3cYtZFwLxeDzOf/7znwtxS9HkDHuf69yQIUMIBoMPRaPRI9v8atM777wzBWhVkGbAgAEUFxff3NDQ8G1aSt7HcCXuf9nbdvR3jC6++OKqp59+elY8Hi8eN27cq48//viHPXm8lpTQeMABNB5wAHVAxQUXDCp+/vmDCjZuPIympsNEdShAaMECQgsWUH7nnShsIhR6ScvKnm+aMuX5tX/846dpeXIZJMkU0BGR+2mdbHVJVU/uQ5taCYfDV9K+R+n3kUjkpHA4vCvt19nKyuy/v3ietweuu37PUCj0Vqbbkwp+V/EGoCLZLttsZjHKfhYj8BfDHIC7kDaGhPlSCbcxtJ/T0BUFlgOL/duShO8XA0tpvZBlp3ItRiLyBdxFuxG4dp+XuFRK2/j41Qvvw/XONS9g+gJwpKp6InI8cB4uRk8Cl2uSr11/2XvvvW+YM2fORUCjqhbg4j0tDUvE9IvmGP34xz8ee//99/8Mt+RMMS3FH+5S1f4a3dNnuXYOJSOV/4tEJIRbTPcbuPi+CxzbWcELf97Wpbi/iyiuEMg9falU2Z8xEpEdcRXGK3AFK4LACar691QdAJgcLy8/2ps06aeF//1vVDru9HkfV/X8GeAlsr2iYxKSSrZM7rA3z+xnMcp+ORkj94+sDHdlt90tXlW1Td3ZZ99afvPN1xds3lzkbx/Q2f1xJcV7Ko67qtw2iVrsb1tKiiq25WSMutBBsvUV4K+0LqXdBFyoqrdnoo09ISLjgA9pPcymCfizqmbF8M6eSozRyJEjz8IVxEi82BADjmhbPCFb5ds5BPa/qC9E5FVgH1onQA3ANqksvNMco/Kf/3xU+W237QXMAg7HXcxrqwG3oPKTwJOo9lddiJTq1aLGWyPP84bT9fos2WJ889cu1jDIKVVVVWV1dXWUl5dP8Twv5WWZM8BilP3SHyPPg48+KpIFC0rl00/LWL26RGpry9iwoZS6ujK2bCmlvr5MGhpKaWwso6mpFM8rJRotw/NKicXKiEbd11isFCiVLsaPF6xbR8VVV0GbuQY9oRCjsHA1JSU1lJXVMHDgCh08eAXbbluje+5ZE//qV1czbFi0k4dXABWk7vXMq/Oo7Tk0YcKEGQsXLoypamKyFRo1atRBnue1q9aVbY4++ujpkUhEVTXxb7KwpKRkP/8DZM5JjNHQoUMPWrVqVateXRGJTpkyZZbnebmy5llenUNg/4v6oqCgYLd4PN42Lyi++OKLj/A8b1GqjtMco+B1123v3XjjUuAePO+egptu2lYee2yaLF48jXXr9hTVIlzP8Uz/druWlCyO/fGP39Mjj1zf5UH6SbIJcK96tkRkKq6bdDStFxkDV/3z6PaPym2e511JDhTIMMb0E8+D9eth7Vpk3TpYuxbWrWv9/YYNUFcHmzZBXR2S8D2bNiHRzvKS9NHiYigvd7cBA9CE7ykvdz+XlX1+Hy0vh1Gj0G23hZEjIWjX6PrDvffey7nnnktTU0tHYGFhIeeffz4//Wn2VxifM2cO++67b6ttBQUFHHLIITz11FMZalXqXHTRRdx5550kfgAOhULcc889nHhiVk2dMyYpEydOZNGi9jnVsmXLGDJkSP82pr4eeeUV5NlnKXj6acRvl44dS/T990GSqkuRdqFQKHUFMhKJyA+BW3ClHz/CDQvYGvwaiGS6EUkYjyvJ+Q3conQ5LxqNltXV1b1cXl5+QDAYzJcrVRajbOB5sGBBsbz5ZoV89NFAWblyILW1A6mo2FknTz5VnnwyIp9+GqWxcSBNTRV4XjmeV0E0OlDi8bL+aKJClEBgi3/bTDC4hWDQfQ2FNhMKbaGwcDNFRVu0uHgzpaVbKCvbTHn5Fq2s3MLgwZsZMWJLdMIE3bDzzk+VDxqUWzHqXF6dR23PoalTpxYHAoEHRWSUqoZEJBqPxzfsscceX8OtKZnVJk+ezA477HDXkiVLvuAnJDFVje+yyy4nkaPxSozRLrvsUhqPx/8sImWqGhQRr7Cw8KNp06adBORKN1FenUOQw/+LOtdvMZo0adL+ixYtaq6OKkBs3LhxDw0ZMuS2VB4nqRiVlKAzZqAzZhC/6SYKfvGL0fK3v00nGPQQ+Wsq29MfetyzJSJLcRX/zs2Fhem2NjYGO/tZjNLATUyuBKpwi0o235L5OR01Z5twK9xvwFXE6s1tk/+1kRRMrs14jFIs386jjuIjIgNxE+53wRWBul5VV2aulT2zcuXKvR944IHZN95447/XrFnzEXC7qs5PvI9fCOQC4Me4eYL/Bb6hqlm3jEQH8+pGARfj5prMw8UnZz7g59s5BNn/Pici+wG/xxUYqgFOU9VnOrt/P8/Z2he3jtVI3Fyp3+AK+aS0uENNTc1Az/M2TJgw4Zf19fUn4zp+ngZOVtWsv5DUG70ZD1IKPG6JljEm5URK6DpB6ux3vVsxshsaCGySWGwNsDbhti6Jn+tTkSCZrZuqbsQlWzmpuro6dt5553Heeeed38UHxbNwVRqb56btC/xLRCZnW3XFtvyE8AeZbofJDX61v+dxF/gKcOXpnxCRfVR1TpcPTn/btgf+hVuqQ3CLyZ8N/BGYnerjXXfdddTX159OS5n8I4BHReSQVCd32aA3ydZfcFVD/pXitmQ1K5CROTbhNfu1itGGDVtk9uwymT+/gsWLB8rq1QNZu7aCjRsHyubNA6mvH0hjY0WroXnRqBua14M1mZKlECUY3EgwuJFQaAOFhRspLNxIUdEGSko2alnZRsrLN1JVtUGHDNnIqFEb41/84hB23fXXBIOnk/zQjRAw1L+RwkIQKWHnUXbLw/hAEjEqKSl6TdaQAAAgAElEQVQ5r76+vlUREGCnU0455eue583v8EEZkocxyqtzCLI7RrvuuutJ77zzToG/DAK4jl3dbrvtzvU875ZOHtYvMZoyZcq35s+fH0hoGwBjx479oed5N6XyWFVVVWUPP/wwtCRa4BLQg+67774v5lCBmfQVyBCRIuC3uKz8eToYO66qj/RopznACmSYrZIqbNgAK1ciq1e3fF2zxhWAWLu2fZGI9euReOo7vnXAABg0CCor0UGDoKoKBg1CKyvd1+ZtVVWo/zsGDXLFHrJkMq0xprVddtmFDz9svW6qiPD8888zffr0DLXKmNS75ppruP7662lsbOmwDQaDnHHGGdx8880ZbBnccMMNXH311a3aFggEOPXUU7n99tSvNFFVVcXmze1z4Q8//JBtt+2oAnx2SluBDGAisD9uMcsTOvi90npdkHxhBTIyxCa8pp7MmVMiL71ULe+/P5gVKwZRWztYNm4cxObN1dTXD6axcRBNTYPxvEGimrI5TQrxz3uZgsHEXqaNlJRs0NLSjQwYsJGqqo26zTYbGDVqo+6000adOnUjo0b1Z+m+jMco1ew8ym55GB9IIkYicqaInKSqzZ9H4oFAYENRUVEY2NJP7UxKHsYor84hyO4Y1dbWTmhsbHwA11kBQDQa1dWrV38feL2Th/VLjNauXbtzY2Pjg4lti8Vi+tlnn/0AN48yZaLRaNmMGTNefuyxxzxVbe7dipaUlCzdZpttjsPlEXmlNz1bb+GStEuBD+igGmFnq2ub9LMJr9kvbTFy852GAsMSvg7rZFtvFqxNVE/CPCUNheq2HHvskSWRyO0FW7aspPU8psTv68iB+Z52HmW/fItRvsUHkouRiISAu4Hv+puWA0eq6tv908rk5VuM8u0cguyPkYicAPwOt2ySB5yjqnd3dv9+LpBxHK54RwkQxRXCuzPVx6mpqRm4adOmDbvssss/o9HoDH/zu7jzfnGqj5cNetOzNR44pqvqKcaYFBEppGUeUHdJ1MA+Hq0BWAms8r8mft96W5uKWyv8f3Abbrrp/2XjPzhjTOr4FQTPBb7lb/o9rtJgj69Iq6oHfE9EzgfKgRpVjaWssRkkIrsAv8BVnnsHV9ktpVUWRWR34AbcuqdzgR/lUsXKfCcis4Cf4Io5vYirXlmJ+5+9WlXrM9i8VlT1ryLyOK4+QVrbNmDAAJYsWXLcyJEjC3HztVbkY2GMZr1JtuYD/by6mTF5pLY2QFMTBbfeujO33DKErpOoQX08mkdnCVP7bXVWQc8Yk4SrcB8amz9D7ApUA5f3dod+yee8KfssIjvghoYV4V6nccB+IrKrqtam6Bi74IZ4BXHTN7YHpvmVHO2iV4aJyOHAE7jqfgLsAEwGDs7WEWCq2gT0W9tUNWeKYfRFb5Kts4DfisgK4EVV7c95FBlj1QgzJ5urC31uyxaR556rkDffHCyffFLN6tXVsm5dNZs2DWLLlsHU11fT1DSIpqbBwWi00p9R+cfeHEohRii0jsLCWoqKaikurqW0tFYHDqylqqqWYcNqdbvtanXq1Fo9+OCNhELd7XIA7p9Aryvo5USMesbOo+yXVzHKlfhs2LChALiE1nOzgyJy6Zo1ayIVFRWJw4S32hhtv/32537yySdFCXPRQiIyeO+9977A87yULMo6evToS5cuXdqcaAEUisiI6dOnn+t53hNJ7CKv4gPZdR6Vl5dfV1dX15xogevBmX7GGWd8w/O8d5LcjcUoi6WzGmEdrlxjCIjj5m4kUlWt6NFOc4BVI9wKqcK6dbBqFbJqVfuKfKtWIStXwurVsHo1Euv9yBcVgcGDYehQdOjQLr9SXQ2BfKxBY4zJdps2bWLQoI473NevX09paV+ng+aHM844g9/97netBgsUFxdz2WWXcckll6TkGCeeeCJ/+9vfWm0rLi7mmmuu4Qc/sOW/Mm3y5MksXNi6pkVBQQGRSIQZM2Z08iiTS9JZjfAm8rBSSBKsGmGGpLS6kOchb7xRJq++Wi2LFlVTU1PN2rWDZOPGwWzeXE1DQzWNjdU0NVXjedXSclWyVzQQ2EhhYS2FhWsoLl5LaWmt7rRTQfy4406QDz64Ac+bp1Om1OrMmeuors7ZeQrZXAGql+w8yn55FaNcic+AAQMoLi5+pKGhYRQtPSqxkpKSJaWlpce1uftWG6P3338/rKo/IaEHsKGhQd9///0zgDdT0Z7FixcfD/yIhM9yDQ0N+uGHH56EmyPWnbyKD2TXebR58+ZLReTohIp7qKq3Zs2aI4Fkh89ZjPJAj3u2THbbaqsLiZTRvmBEZ/OgSvrYpI0kU0QCVqPa2PbBW22McojFKPvlW4xyKT4iMgF4gZb526uAQ1S11YfBrTlGfhGRXwOn4iq7BYCfqOq1qWqPiASAB4ATgRiubPePVPXWZB6fb/GB7DqPRGQA8CwwDRefKHCCqj6W7D4sRvmhT1ftjUkrt4D20G122mncxksvpfJHP/oO69ZV0nESNaCPR9tC8pX4sqZ6kDHG9DdVfU9EdgK+4G+arap1mWxTtvErq50mIr/EVQr8UFU/SPExYiLyTVzFwxHAQlX9KJXHML2nqptE5ADceVIJzFPVmgw3y2RAj5MtEXmhu/uo6iG9a47Je25NlSF0vw7UUNybE6EPPqD65JMBerqMeSPJVuJT3dSXp2WMMVsTP7n6V6bbke1UdT6uinOfiEgxcBGwF1AD3KiqH/lJ3Vz/lhVE5CDgFNwokieA+9uW9RaRctx6rbsCi4Hr8iEREZEC4GzgYFx1zTtVNaWLAndy3HHAhbik+03g56rakIL9DsfFaTvc0NRrEy+siMg+uMJ55cDzwN2aA2tp9rfe9Gxtof2crUG4cpbryaIT3vQTN5RhMMktplvdx6NFcUlSMknUBitlbowxJpeJSBB4DtdDUoj7P3iiiOyebT1ZIvJl4O+0lDs/GtgRuCzhPiXAa8BOuOfjAcf7Jetzdo0wf+jofbhhnc1F5L4pIger6qtpPO4OwFu45DYIzAQOE5FD+lIxXESGAW/j1ggLATOAo0Rkb1WtF5GDcX+XghvCehQuFzijL88nH/U42VLVIzvaLiJDcAUk/tDXRmWjra70+5YtIi+8UCFvvjlIPv7YlTJfu7aaurrBbNkyqFUhCagSd6L1ikKcYLBtKfO1lJev0UGDauNjx27Z8J3v3DJg06ajCvbbbwWlpd0lUKW49UZ6Xco8zayUa/azGGW/vIpRHsYHLEYpcfTRR09//PHH96Wl2EYQKB05cuRNnuf9rA+7Tnl8CgsL72xqakr8PBAALo1EIs8dfvjhGwCmT59+1Kuvvjq+TVn8ynHjxl3ned4dfTl+Js+jK6+8cuyVV175nYRNBYCUl5f/yvO8k3q5225jNHLkyMuXL19eSsISAMC0o48++vue573Sy+Oy4447fn/RokVVCQU+CkVkp2nTpl3ied7jZWVld23evLmAltL2QeD0m2666clzzjlneWf7zaf3umTn0aVszpaqrhaRG4DrgZSsIZFlTie3Sr8/1G6LKmzY0L6E+cqVrrT56tUtpcxXrUKifVtCTaurXcnyYcNgyJCWr0OHwrBh6JAhMGwYDB5cQCBQTRe9XgPdl3/0qUHZp32MclQwGKSqqgrg5Uy3JcUsRtkvL2KUx/EBi1GfzJo1i2effZaGhlajwgITJkw4Gtdz1Fcpi08wGKSpqand9rFjx34+BWXmzJm88cYbrZ6PqoamTp16EnBSX4+fqfNo2rRpiEjbATVSVVW1K67IRV90GqPx48ezfHnr3Ka4uDhw+OGHJ1UopTPTpk1j0aJFrbYVFRWFZsyYcTlweXl5OZs3t8+V9t577y4rd+fZe13aSr93pYDc6P3pjawt/S5z5pTISy9VywcfVGtR0WTGjz9XHnnkUfnoI6WhwZUzb2oajOcNEtXCvhxLA4E6QqGWHqiSkloGDKjViopahgyp1TFjanXSJFfKfNiwlCx4nYdlQq2Ua/azGGW/vIpRHsYHLEYpMXfu3F0bGhp+R8IHOxHxlixZ8gfgzj7sOuXxKSgoeMDf7+cl70WkIRaLHQo0ALz55pv7NzQ03ELCiBgRib7zzju308fEL5PnUU1NzWBVfYrWzz26ZcuWV4ALernbbmO0ZMmSs0Xkm4kl5hsaGvTtt98+ieSWAOjQ/PnzTwR+SOulBeJz5sz5EfBKY2PjHbihrYm5RHT16tUzcdOKOpSn73Vd6s2ixnt0sLkQmIDr+VmkqoeloG2mMyLXA9NpmQfV11UkN5FcJb7VpGDCZU/lW5lQK+Wa/SxG2S/fYpRv8QGLUSqJyFXAj3EJSwiYBxyoqr3+sJqO+IjIeOAV3ICUOC7x+KqqRhLuI8BduLL4jbjPkC8Bs1S1T+MZM30eicjJwL1AE+65rwD27W3xj2RiJG7pm5eAKbj5b8XA1ap6eW+OmbDfEPAUcJC/3yLgt8DpqqoiMgY3924bXKwLge+o6oNd7TfTMcqE3vRsvUn7AhnNV1tex508Jr12wa3b0CkVaRTVFXSfRK2iD2/WxhhjjEkvVf1/IvI0sBvu//c/VLX9eL0MU9WF/jpsR+A+nL/UtuS9/0H9TOAvuJ6bZcCTqhrr9wanmKr+TkTmAPsCdcATqprWhEJVN4vINFyBiqHA26r6Wgr264nILOBLwCjgfeCF5sqSqvqpiEzyfz8AeFVV3+3rcfNRb5KtgzvY1gAsU9VOJ8SZlPqf/7VdEhW7+urq+BlnPEpJybRQaWnOX0k0xhhjDPgfoPv8ITrdVHUNbrHlru6juIWxu11OKNekqtx/D4/ZhKsCmer9xuhiCo2qridP5mSmU2+qEb6UjoaYHlC9prNfxV2XszHG9JmIjMKVZ16uqu9nuj3GmMwTkUG49bE2AvNzqUfKHxq3G26o3by2vU4iUoErX96A6yHKyjKa/rpru+GGKr7dl+GkSRxrR1oW5l6apmNMAPbDzfX6t5+w541el+s2xhiTv0TkLGAJbqHKhSJylz/XwhizlRKRw4BPgRdxazu9KCIDMtuq5IjINsAbwGzcHKdPROQLCb/fB/jE/91sYLaIDM5EW7siItviCl/8F3gV+MBPVlJ9HBGR24EPcP8HlojIuSk+RoGI3AcsAO7BVTOvEZEvpfI4mZZUsiUidSKyMcnbhnQ32hhjTPr4H0DuwF8nxt98in8zxmyF/B6tx4AyWt4X9gFuzlijeuZeXDE3cO2vBJ4QkSK/p+gJf1vzc5voPybb/BkYk/DzUOAfabgY9i3gLP/75kWqbxaR/VN4jNOA77TZFgIeEZG8qW6e7DDCm2hfFKOt/YFDk7hfTtrqFjXOIvm0AJ7PYpT9tuoY7bbbbsfPmzfPU9WihM2BYcOGfdnzvL6uF5MqeRWjPDyHwGKU7XoUn29961t7/OEPf2hb/biwqKjocC9LpjB0FaNAIHBILBZLXP6mANjm8ssv/1I0Gi249tpr2/ZiFQYCgUMy/Nxaxai2tjaAK7eemFgFgHGPPPLIwZ7ndVpyvadGjhz5leXLl7fqlBGRpkmTJn3N87wtvd1vYoyGDBnyldWrV3fU8RM65phjjvU8L6vnKCZbxbPHpd/b7UBkP+BKXKI1H/iZqj7Sp51mIc/zriS3FjU2xpheue+++zjnnHNaLU4aCAQ49dRTuf322zPYMmNMpsydO5d99tmn3fYpU6bwxhtvZKBFPbPddtu1W/wX4OOPP0ZVGTduXLvfjRo1io8//rg/mpcUVaWyspL6+vpW20WE9evXU1JSkrJjXXjhhfzqV78iMREvLCzkxhtv5Mwzz0zJMc4880x++9vfdvi7F198kf322y8lx0mXUCiU3kWNRWQ6Lvk4GLfewzGq+lhv95cDsnZR4zbyahFJyMsF8CxG2W+rjtGoUaMGqupfcUNqgkA8Ho/HKioqvgks6o/2JiGvYpSH5xBYjLJdj+Kz/fbbFwwcOPA3dXV1kxIW0NWhQ4degpvTk3FdxWjkyJHHLV++/CL8KTQi4lVXV78yatSoCz3PY/Dgwb+ora3dP+G5xUeMGHED8Ld+fhqJWsVIRNh2221PX7hw4ffwF08Wkejo0aMfLSkpuT6VB66qqhobjUYf9o8TEJFoPB6vGzFixFfpYtHi7iTGqKqqagTu+QUS71NaWrpw3Lhx3wGifXkO2aI3ixofgOvJOgg3OfJniYvVmczKt0UkIf8WwLMYZT+LEfgLVv4SV3VsGXC+qs5ObyuTl28xyrdzCCxG2a438fGLYdyMG820EbjWvzCTFbqKkT+n6XvAOUAJ8A/gUlVt9H9fDFyLW6+qHrgNuE/7OgSsDzqKkYgUAOfi5tAGgIdxixinPDERkT2AW3BzxBYAZ6vqJ33ZZ9sYichU3GLJO+CSqyeBU1U1b2pAJN2zJSIH4XqyDsRVczlKVZ9MU7uMMcZkkKp+CoQz3Q5jTPZQ1U24ogY5x0+a7qWToheq2gD8yL9lLVWN4xLetBcmUdW3cJ/703mM/wG7pPMYmZZUsiUi/wamA68DR6jqM+lslDHGGGOMMcbkumR7tg7wv+4K/KWb6pKqqhV9apUxxhhjjDHG5Lhkk62fprUVxhhjjDHGGJNnkkq2VNWSLWOMMcYYY7KMX/zjh8CluEWnXwZOVtVVKdp/EfAL3ELHgltYuQlXJbEA+CvwA1Wt73QnPTteAfATXDGTIuAZXNGMlK0j1p96XfrdGGOMMcYYk3Fn4JKh5hLqhwHPicieqpqKFcXvxCVazYtCf8//2rwg8bdwSd7XU3AsgB8D/4+WPCUMDBGRgzJZHbK3LNlKkud5w4HhmW5HEnq0Inwu6GpF+BxlMcp+FqPsl1cxysP4gMUo2+VVfGDrjVFJSckF9fX1iWtVhYBdzzrrrOM9z1vQlwbU1tYGgJNpSaxo8z24JOyEuXPn3jFp0qSGrvaXTIyCweC50Wg0MUcpBA64+eabv+R5Xk0vnkZaJL1kQg4miBnhed6VuNL3xhhjjDHGZIXx48fz8ccft9v+4osvst9++/Vp342NjQwcOJBk8oXPPvuMioq+18gbMmQI69e3HzE4f/58xo8f38EjMiMUCnVZMbCZ9Wwl79dALize3KMV4XNBVyvC5yiLUfazGGW/vIpRHsYHLEbZLq/iA1t1jM4VkRNUNeT/HA8GgxtDodBRwJa+NKCoqIjy8vJ7Nm7cOJmWvKE582pONqJlZWUfVlRUfLO7/SUTo5KSkp9t2LBhRsLziRUWFq6urKz8Mm7h45xiPVt5pjcrwme7rlaEz0UWo+xnMcp++RajfIsPWIyyXb7FB7beGPkFLH4PHO9vWgl8yV+UuM9EZCjwD2Bvf9M7gAfs7v/8f7h1eJd1t69kYiQiA3BFN2b5mxYDh6tqTl4UsJ4tY4wxxhhjcpSqNgIniMgPgXJgSYoKYzTvf5WITAVG43qzPsX1bo3Gzd/6VFXjKTzeJhE5Alcrodjff871aDWzZMsYY4wxxpgc55d6T0m59w72HQeWtNn8aTqO5R9PgawphtEXbauJGGOMMcYYY4xJAUu2jDHGGGOMMSYNLNkyxhhjjDHGmDSwZMsYY4wxxhhj0sCSLWOMMcYYY4xJA0u2jDFJE5GjReT3InKfiBya6fZkKxEJi8j9/uv0xUy3x5hsJSKjReQWEfmTiFzsrxeUc0RkJxG5Q0QeFpFzRCSQ6TZ1RUQGT5w48ZwTTzyRyZMnn+6va5SK/ZaLyM9E5M8icr2IVKdiv/6+jxORB0TkHhGZnqr9mhYiEhSR8/y/49tFZFy6jnXqqacOE5Gb/HP/MhEpTtexMs1KvxtjkiIiPwBuxa2xocBJIvINVX04sy3LLiLyfeAO/8fm1+nbqvpgBptlTNYRkTHA28AAIAR8BZglIoepaiyjjesBEdkFeAP3HILAscB0EfmaX746q/gJ0NyPPvpo6KJFixCRk4F9RWQfVa3vw35Lgf8BOwCFQBPwdRHZTVXX9bHNlwDX0PL/57sicoyqPt6X/ZoWIiLA34HDcX/LUeBkEdlLVd9P5bFWrFjB008//RowkJZz/wgROTiV64NlC0u2kuR53nDc4mrZbnzzV8/Lj7/Xqqqqsrq6OsrLy6d4nrc50+1JgZyL0bJly4IicrOqNveGC0AgELjL87z3LUZObW1tQERuVVXxNzW/Tr/yPG9BitvYIxaj7JaH8YFuYjR69OhLly5dWk7LZ5FCYHo4HD7b87xX+qmNSessRtXV1bfU1tYWAs29WSHg2NNOO+3bnuf9X0Ya24UJEyactnDhwqGqGgJQ1ZCIjJ86deplnuc92tv97rPPPsfOnj17x+b9AoUiMnzixIlXe573297ud/bs2aW4RKvV/59QKPRrz/OWJt43D8+jfnuf+/a3v737Aw88cBT+64s7L0uGDBlyh+d5l6TqOFVVVWVXXHEFQAWtz/2pM2fOPNfzvH+l6ljpFgqF3krmfpZsJe904IpMN6IHHsp0A1IlGAxSVVUF8HKm25JiOROjYDBIRxdoY7FYRSwWm2MxcuLxeGevU7nneXNCoVAHj+ofFqPslsfxgU5iNGHCBJYubfVZmeLi4sCMGTNu7ZdW9VBnMRoxYgS1tbWt7hsKheTQQw+9v/9al7xp06axcOHCVtsKCwtDhx122E+An/R2v4cccghvv/02jY2Nn28rKCgI7b///mcBZ/V2v9XVHY9EFJGhwJzEbXl8HqX9fe6www7jT3/6E01NTYmbA2PHjv0ikLLh8MFgkNWrV6OqrXKQ4uLiwMyZM3+equP0E+n+LpZs9cSvgUimG5GE8biT8hvAwm7umxOi0WhZXV3dy+Xl5QcEg8F8uVKVUzEaOHCgBAKBf8discRx/fGioqKaQCBwtMXIKSsrk0Ag8EIsFiun5U1YCwsLV4ZCoSPT0M6kWYyyWx7GB7qJ0SeffHK2iHwzoSeEhoYGffPNN0897bTT5vZjO5PSWYzWrFlzpYjMSnwenufx7rvvHvuVr3xlcUYa24W5c+d+HTiXhM+AjY2NsdmzZ18K9LpXYfbs2TMaGxuvpqWHj3g8Hn3rrbduAv7S2/0GAoFQQUHBv+PxeOKcnlgwGPwEOD7xvnl4HvXb+9z8+fN3bGpqepiEBEJEvBUrVjwOXJeq40Sj0bIdd9zxZRHx2pz78dmzZ5+FG5KbVyQLhxObPvA8bw/clZ49k+3ezHY1NTUDgQ1AxYgRIzZmuj19lasxEpGZuAsOcdybsQccpqqvW4xaiMhhwBO4eQWCG/f+RVX9b8ob2gMWo+yWb/GB7mPkF2V4FZiIO0+KgNtU9Uf92tAkdRYjEdkGN1dpFBADioHLVPX6jDS0GyISAp4FDiguLg40NDTEccnQiX2ZYyYiBcCfgGNw87VCwEvA4X2dhyMiRwN/w/2dCNAIHKCq8xLvl2/nUX+/z4nIFbhRXA24ZPwTYF9VXZuqY9TU1Aysr6/fsPPOO8+NxWKTaDn371bV76fqONnEeraMMUlR1WdFZBIwA5dwPamqn2a4WVlHVZ/3X6eZuNfpKVVdkuFmGZN1VHWTiEzFfTgfBsxT1ZyZr9FMVT8Tkd1wz6MKeENV/5PhZnVKVT0RmTFz5szzDj300J//5z//+dHjjz9+e1+LeahqXESOB44EdgQ+BR5NRbETVX1cRCYDh+Iu9P1DVWv6ul/Tmqr+VEReAPYC1gCPqGrKewhLSkq49957v3jyySfPwNVDeAd4LtXHyRaWbBljkqaqHwIfZrod2U5VFwGLMt0OY7KdqjYAf8x0O/pKVeuA32e6HclS1WhzIYJzzz33lVAolJJhTn7C9o9U7KuDfb8HvJeOfZsWqvoKkPYCNTNmzPC2lmrGts6WMcYYY4wxxqSBJVvGGGOMMcYYkwaWbBljjDHGGGNMGliyZYwxxhhjjDFpYMmWMcYYY4wxxqSBJVvGGGOMMcYYkwaWbBljjDHGGGNMGliyZYwxxhhjjDFpYMmWMcYYY4wxxqSBJVvGGGOMMcYYkwaWbBljjDHGGGNMGliyZYwxxhhjtmrinCciy0SkVkQeFpHKTLcrW/ivz7kJr8+fUv36iEiliPzZ3/9SEfmhiEgqj5EJwUw3wBhjjDHGmAw7B/g5LZ+NjwFGiciBqhrPXLOyxg+AXwAB/+ev4F6fA1Lx+ohIAfAksBdQCAwCbgJiwC/7uv9MsmQrSZ7nDQeGZ7odSRjf/NXzvIw2JFWqqqrK6urqKC8vn+J53uZMtycFLEbZz2KU/fIqRnkYH7AYZbu8ig/0LUaFhYUXNzU1JX4uLgT2v/LKK7/sed7ilDY0eVkTI//1CSRuAva7/PLLe/T6dBajyy67bPtrr712Wpu7B4qKii72PO+1vrU+PUKh0FvJ3E9UNd1tyQue510JXJHpdhhjjDHGmNQaPXo0q1atarf99ddfZ/fdd89Ai7LLmDFjWLlyZbvtqXp93nrrLaZOndpu+/Dhw1myZEmf958OoVAoqSGO1rOVvF8DkUw3IgnjgYeAbwALM9yWlIhGo2V1dXUvl5eXHxAMBvPlaqLFKLtZjLJfXsUoD+MDFqNsl1fxgb7FKBAIXCEih6tqyN8UDwaD60tKSo4EGlPf2qRkTYwKCgouF5Ej+vr6dBajwsLC4mAw+EQ0Gq3ArykhIl4gEPgHcE0qn0t/s56tPON53h7AHGDPZLs3s11NTc1AYANQMWLEiI2Zbk9fWYyyn8Uo++VbjPItPmAxynb5Fh/oW4xEZADwCPBFf9NK4AhVnZvaViYvm2Lkvz5/B2b4m1YBh/f09ekqRiKyB/A0MMTf9CxwrKrm9MUN69kyxhhjjDFbNVXdJCIzgbFAKbBIVTPVo5V1/NdnFu71KQM+TPXro6pvicgYYEdgM7BY86BXyJItY4wxxhiz1fM/2H+S6XZkq/54ffwE7p10HqO/2TpbxhhjjDHGGJMGlmwZY4wxxhhjTBpYsmWMMXiNuMwAABNDSURBVMYYY4wxaWDJljHGGGOMMcakgSVbxhhjjDHGGJMGlmwZY4wxxhhjTBpYsmWMMcYYY4wxaWDJljHGGGOMMcakgSVbxhhjjDHGGJMGlmwZY4wxxhhjTBpYsmWMMcYYY4wxaWDJljHGGGOMMcakgSVbxhhjjDHGGJMGlmwZY4wxxhhjTBpYsmWMMcYYY4wxaWDJljHGGGOMMcakgSVbxhhjjDHGGJMGlmwZY4wxxhhjTBpYsmWMMcYYY4wxaWDJljHGGGOMMcakgSVbxhhjjDHGGJMGlmwZY4wxxhhjTBoEM92AvgqHw2cDJwG7Ao9GIpETMtsiY4wxxhhjjMmDZAuoAa4GDgMGZ7gtxhhjjDHGGAPkwTDCSCTySCQSeQxYk+m2GGOMMcYYY3pGRA4UkbdEZJWI/EtExmW6TamSDz1bxhhjjDHGmBwkInsCz+M6gQqA6cBrIrKLquZ8Z4olW0mqqakpAooy3Y7uVFZWlgUCAWKxWNlnn302MNPtSZHy5q81NTUZbUgqWIyyn8Uo++VhjPIqPmAxynZ5GB+wGOWCdjEqKio6q7GxUWgZcRcCKgcNGnRCTU3NAxloY1JGjBixMZn7WbKVvEuBKzLdiO6sX7+++duXM9mONFmW6QakgsUo+1mMsl8exygv4gMWo2yXx/EBi1Eu+DxGs2bN4vHHH2/1y6KiosIf/vCHdwB39HfDekCSuZMlW8m7Drg5043oTmVl5ZRAIPByLBY7YP369fMy3Z4UKcedlKOAugy3pc8sRtnPYpT98jBGeRUfsBhluzyMD1iMckG7GL3yyivfwCVVgeY7NTY2xv/85z8feMopp8zPSCtTKOeTrXA4HMQ9jyBQEA6Hi4FYJBLxUnmcESNGNAKNqdxnOnietxmgoKBgc7Ldm9kuYShAXT48J4tR9rMYZb98i1G+xQcsRtku3+IDFqNc0FGM1q5dezewI3Ce/7socOq77777av+3MPVyvhoh8BOgHvgxcJz//T0ZbZExxhhjjDGmW+r8CNgWmAoMV9X7M9uq1Mn5nq1IJHIlcGW6j+N53nBgeLqPkwLjm796Xko79zKmqqqqrK6ujvLy8inNV3lynMUo+1mMsl9exSgP4wMWo2yXV/EBi1Eu6CpGTU1NAB4wxvO8MRlpYA+EQqG3krmfqGq625IXPM+7khwokGGMMcYYY4xJr1AolFSBDEu2kpRjPVsPAd8AFma4LSkRjUbL6urqXi4vLz8gGAzmy5Uqi1F2sxhlv7yKUR7GByxG2S6v4gMWo1yQTzFKtmcr54cR9pdQKLQCWJHpdnQnoZt5YbJ/BNmueW2JdevWzcuHCaIWo+xnMcp++RajfIsPWIyyXb7FByxGuSDfYpSMfCiQYYwxxhhjjDFZx3q2kpRjwwhhK5lMmaMsRtnPYpT98ipGeRgfsBhlu7yKD1iMckE+xcgKZKSYFcgwxhhjjDHGgBXISLkc69myyZTZzWKU/SxG2S+vYpSH8QGLUbbLq/iAxSgX5FOMrEBGilmBjMzJt8mUFqPsZzHKfvkWo3yLD1iMsl2+xQcsRrkg32KUDCuQYYwxxhhjjDFpYMMIjTHGGGOMMSYNrGfLGGOMMcYYY9LAki1jjDHGGGOMSQNLtowxxhhjjDEmDSzZMsYYY4wxxpg0sGTLGGOMMcYYY9LAki1jjDHGGGOMSQNLtowxxhhjjDEmDSzZMsYYY4wxxpg0sGTLGGOMMcYYY9LAki1jjDHGGGOMSQNLtowxxhhjjDEmDSzZMsYYY4wxxpg0sGTLGGOMMcYYY9LAki1jjDHGGGOMSQNLtowxxhhjjDEmDSzZMsYYY4wxxpg0sGTLGGOMMcYYY9LAki1jjDHGGGOMSQNLtowxxhhjjDEmDSzZMsYYY4wxxpg0CGa6ASZ7hcPhIuBO4FBgMPApcG0kEnnI//0k4F5gMvAxcGYkEnnF/91w4NfAXsBwYEIkElnYZv+TgVuBLwD1wD2RSOSyLtrzVeDnwFDgP8DJkUhkuf+7g4HLgT2A+kgkMiyJ55fS/WVCOmMUDoefBqYnHC4ENEYikYFdtMdi1EaaY1SEO4eOAQqBN4DvRyKRD7toj8WojTTHqBy4GQjjzqE/AOdHIpFoF+2xGCXoY3y+BFwC7Ao0Af8Czo1EIqsS9n81cAbuM9HDwDmRSMTroj0WnzbSGaNseE0tRt3G6ELgW8BYYAPwAHB5JBKJddGerSZG1rNluhIEanAnZgVwOvCrcDi8bzgcDgER4FGgCrgeeDwcDlf5j40DzwBf7mjH4XB4EPAc8HvcSb8t8KfOGhIOh8cDvwNO8+//IfDHhLtsBu4DfpTME0v1/jIobTGKRCKHRyKRAc034DHgL501xGLUqbTFCPfc9wd2A4YA7+E+zHfIYtSpdMboZlxsdgAm4uLV1UUli1F7fYlPhb9tJLA9EAPub95xOBw+BTgBlyzviPsw9pPOGmLx6VTaYkSGX1OLEdB9jAqA7wLVuPe4I4ELOmvI1hYjUdVMHNfkqHA4/BTuisZ83Ie6EZFIJO7/bg7wq0gk8ts2j1HaX+29GtguEol8I8njXgPsGIlEvub/PAhYBYyPRCIfJdzvIOBP3V21SPX+skmqYtTm94Nwb9KHRCKR1zq5j8UoSSk8j+4B1kUikYv8n/cEXvKT446OazFKUgpj9BlwdPN5Ew6HvwP8LBKJbNvJcS1GSehNfPzffQH4VyQSKfd//g/wUCQS+ZX/81H+Y0d3clyLT5JSFaOE7QeRgdfUYpR8jBJ+fxEwPRKJHNXJ77eqGFnPlklaOBwuw139eweYBPxf80npe9vfnox9gdpwOPxaOBxeEw6Hnw+Hw7t0cf9JwLzmHyKRyFpcF3iyx0v3/rJCimOU6ERgcWeJls9ilIQUx+heYHo4HB4VDoeLgZOBp7q4v8UoCSmOkfi3xJ/HhMPhik7ubzHqRh/jc6D/uGatXh//saMsPn2T4hj1lMUoCWmOUXe/36piZMmWSUo4HBZcF+1s4J/AANy43ETrgQ6vcnRgNHASrkt3BPAq8Jjfld2Rvh4v3fvLuDTEKNF3/X13xWLUjTTEaCGwGFgKbAJm0MXQjRQcL937y7g0xOj/t3fnMVZWZxzHvz9EXBvUWFdajVaJVesStUoXqSJprD0tqG1qFde0tom11aICMU7VGonaRXHXahe7pK6vdUMliBXFXaRVLFiNGETLJiBCldM/zrnycrkz996B15mG3ye5Yd7tvM99nwwzz5zl3gecE0LYLISwPfCjvH/jTs53jrqwJvnJf40fw6rfI/XXL8j/dvZ8nJ8mKshRu5yjJqrMUQjhNNLcrku7CGGdypEXyLCm8jflNaSxukOLooghhMVA/UIJ/YFFLTb7HnBXURRP5HtcQJp8uWsIYW/ShHOA14ui2J30i2S37hdC+O7abK83qihHtbb3Ik2Y/Vpp31p9ps7RKtp5j9eQfphsRfrBchLwSO4lHoZz1JaKcnQ68GvSfLolwPWkv7a+7e+j9qxJfkJakKkATimK4rHSofrraz1ai5yf9lWUo67u5xy1qcochRCOBUaRphzMzfvW+Ry52LIu5W/KK4F9gCFFUSzJh6YBZ4UQ+pS6nfcGrm6x6amdHSjSyji31O2eBuxVimtz4NO0MNRgbbfX21SYo5qTgPuLophd2+EctafCHO0JdBRF8U7evjaEcBnwWeeoPVXlKA9nOa50n1OBp4q0Spdz1KI1yU9Iq6yNJ60CeWtd07XnM7l07ayiKBbi/LSlwhx1yv/PtafKHIUQjiH1Zg0pz1t1jlxsWXPjgAOBQ4uieLe0fyLwPnBmCOFy4EjSCjV31E7Ic0hq+uXtZUVRROBGoAgh7EcaFzyKNBxqeidx/AGYEkI4BHgcOB94vDbxMYTQh7T0db/SvWNRFMs+pvZ6UlU5IoTQjzRf6/stxOEcda6qHE0BRoQQJpCGUJyYz5vRSRzOUecqyVEIYSfSX1fnAoOAc4ERXcThHDXWrfzkXt6HgFH5l7R6NwMjQ1ooYAlpJcLfdBGH89O5SnLUC56pc9Q8R98h9eAPLYqilQJnncqR52xZp0IIOwA/JC1X/EYIYXF+jS7SZ5AE4CjS8KUxwDfzX3FrluYXpImLS0lLvFMUxSPAmcDtpF9CBufrG372TFEULwEnkxYEmAsMJBUBNV/O7T9A+oyFpXReuK319npKlTnKaisJ3d0sFueosYpzNDJf9xIwj/RZQcOLolhAA85RYxXnaB/SH5QWk4a+/KQoioc7i8U5Wt0a5uenpGG2V5SuW1xq/gbSR1o8A8wk5erCzmJxfhqrOEc9+kydI6B5ji4CNgMeLR2/r7NY1rUceel3MzMzMzOzCrhny8zMzMzMrAIutszMzMzMzCrgYsvMzMzMzKwCLrbMzMzMzMwq4GLLzMzMzMysAi62zMzMzMzMKuBiy8zMzMzMrAIutszMzMzMzCrgYsvMzBqS1CEp5tcKSQslTZU0TtJu3WxzsKTRazvWtUXSREl/a+G81ySNK213SFrcxn0G5+e6X3djbeEe4yS9VlX7ZmbWnIstMzPrylLgIGAQcBRwMzAEeF7Ssd1obzDQa4utNXAD8JWeDsLMzHqXvj0dgJmZ9WorYoxPlLYflHQVcA9wo6TJMcZXeyi2XiPGOAuY1dNxmJlZ7+KeLTMza0uM8X3gNKAfcEptv6QRkv4uaZ6k+XlI3gGl4x3AecAmpeGJE0vHd5N0Vx6uuETSPZJ27ioWSTvmdo6XdGO+dp6kX0jqWzqv4TA/SYtzXPX7R0iaKWlpfh8Dm8SxSvuS1pd0iaTXJS2TNFvS3ZL61126haQ/SlqUzz2rQdsHSZqQn8nCfP5WdedsJ6mQ9J6kNyWN7CpeMzP7eLhny8zM2hZj/KekN0lDDGt2BH4HzCQVYscAkyR9Lsb4Cmmo3YC8/5B8zbsAknYCJgPTgBOAFcAY4GFJA2OMy5qEdBEwHvgWsC9wPrAcOKcbb29fYOfStRcCD7QYR80o4FTgbOAfwJbAUGCDuvOuBn4PDAOGA2MlTY0x3g+p0AImAvcC3wY2yfEUwIGldu4iPdsfAAvy/QcAH7QYr5mZVcDFlpmZddcbwDa1jRjj+bWvJfUBHgT2JxVPo2OMsyTNYvWhiZB6vOYDh+WeMyRNBv4NnAxc1SSWmTHGE/PXD0jaGDhD0tgY4/w239fWwMExxn/lOJ4HXgaOB65rsY0DgPExxnLctzU477YYY0e+zwTgCNLcuPvz8YuBp4HhMcaYz5sGvCjp8BjjvZK+CuwHHBpjnJDPmUTKz39ajNfMzCrgYYRmZtZdAuJHG2kY4B2S5gAfAv8FBgK7ttDWUFLvzAeS+uYhgPOBF0gFWzN31G3fDmwM7NnCtfWm1QotgPz1NFbtSWrmWeDwPLxw/1x8NjK+dJ8VpKJuAEAuGL8A/BVYr/RcpgOzWflcPg8srBVaua35wEfbZmbWM1xsmZlZdw0A3gKQ9AlS4bADcAbwJVIx8AKwYQttbQn8mFSglV+DgE+1cP3bnWxv28K1zdqq7WunrZ8DY0m9YU8Cb0k6T5LqzltQt72clc9rc2A94Jes/ly2Y+Vz2RZ4p0EMc9qI18zMKuBhhGZm1jZJuwPbk5aChzR3awBwRIzxhdJ5/Wltlb55pBUOGw0XXNTC9Vt1sj07//s+sH75BEkbkHq/mrVV2/dMC3EAkOd2dQAdkj4DnJS3XyXN0WrFAlLP4UXAnQ2O14YIzgY+2eD41q3Ga2Zm1XDPlpmZtUXShsAVwDLSohcAG+V/l5fOG0RaNKNsOasvEgHwELAH8FyM8em61/QWwhpWtz0ceA94MW/PAvrVrW44hDQUst4eknYpvY9dcmxTWohjNTHGGTHG0aSCsuUPg44xLgEeB3Zr8EyejjG+lk99EugvqbboCJI2Z+UiJGZm1kPcs2VmZl3pI6k2V2lT0hyo7wE7ASeUfuF/AlgMXCnpYlKvVwfwZl17LwF9JZ1OWn3w3VxMnQc8RVrc4jrSELhtgIOBR2OMf2oS586SbgL+TFpN8GzgV6XFMe4DlgDXSxpL6oU7nVJxWDIHKCSdSyrGLsjv47dNYviIpDtJPWHP5ft+HdiC9udRjQQmSPpLfm/zc+yHATfFGCeSFtN4FrhF0tmkHrHRrD5E0czMPmYutszMrCsbkXpXIqmYeh14GBgWY3y5dlKMcY6ko4FLSQtdvMLKpc/L7iYNFRxFGpo3CRgcY5yRP5Prwnx8U9LwuEnA1BbiHAMMJi0m8WFuY0wpvrmSjgQuIw3Jex44DnisQVvPklYOvIQ0H2oKcGptlcQWPUZahv5M0s/a6cAxMcaH2miDGONkSV8EfgbcRFpSfxYpBzPyOVHSN4BrgGtJBdnl5GGd7dzPzMzWLuWVZM3MzP7vSNqRtDz80THGW3s2GjMzs1V5zpaZmZmZmVkFXGyZmZmZmZlVwMMIzczMzMzMKuCeLTMzMzMzswq42DIzMzMzM6uAiy0zMzMzM7MKuNgyMzMzMzOrgIstMzMzMzOzCrjYMjMzMzMzq4CLLTMzMzMzswq42DIzMzMzM6uAiy0zMzMzM7MK/A/DICyKvYEXbAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 1000x400 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"<ggplot: (-9223371872006454920)>"
]
},
"execution_count": 43,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"p.options.figure_size = 10, 4\n",
"\n",
"(\n",
" p.ggplot(fdata, p.aes(x = 'fdate', y = 'Number individuals'))\n",
" + p.geom_smooth(color = 'red', method = 'lowess')\n",
" + p.geom_point()\n",
" + p.scale_y_log10()\n",
" + p.theme_minimal()\n",
" + p.labs(x = 'Date published', title = 'Number of human donors in scRNA-seq studies')\n",
")"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"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.7.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment