{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# One drug analysis" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In the introduction we have analysed a IC50 input file. We performed a full analysis that is we analyse all associations between all drugs and \n", "all features. This may take a while and sometimes one is just interested in a single drug or a sub set of drugs. In this notebook, we show how to \n", "restrict the analysis." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "source": [ "%pylab inline\n", "matplotlib.rcParams['figure.figsize'] = (10,6)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "As before, we create an ANOVA instance. Pur input IC50 will be the test file" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from gdsctools import ANOVA, ic50_test" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "an = ANOVA(ic50_test)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Drugs are stored in the **drugIds** in case you forgot the drug you are interested in.. eg. the first drug identifer is " ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "drug_name = an.drugIds[0]" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " [-------- 21% ] 10 of 47 complete in 0.0 sec[---------------- 42% ] 20 of 47 complete in 0.1 sec[-----------------63%---- ] 30 of 47 complete in 0.1 sec[-----------------85%------------ ] 40 of 47 complete in 0.2 sec[-----------------100%-----------------] 47 of 47 complete in 0.2 sec" ] } ], "source": [ "results = an.anova_one_drug(drug_name)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now, you can create an HTML report or simply call one of the volcano plot" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnkAAAGXCAYAAADPr0bOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl8VPW5x/HPM1kggSQQFknYFQVcwFQEvVYEF8TixbrU\nvVWrVVvvra120bZqa63X9nbTWu3V1qq1itq6tS6tioqoKKiAgoCIKDEokhAISwhJnvvHOYlD1slk\nkkkm3/frldeQc37nnGd+M+LDbzV3R0RERERSSyTZAYiIiIhI4inJExEREUlBSvJEREREUpCSPBER\nEZEUpCRPREREJAUpyRMRERFJQUryRHoAM5tmZm5mdzY4fm54/McJek5C75dKmqsbM/txePzcNt7v\nm2a2zMx2htc/n8BwRSQFpCc7ABERaRszOwm4EdgEPAZsA1Z04vNHAe8DL7j7tM56roi0jZI8EZHu\n54vh6ynuPjepkYhIl6XuWhGR7mdY+LomqVGISJemJE+kGzKzWWZ2h5m9Y2ZbzGybmS0xsx+YWa9O\neP5hZvaMmVWYWbmZ/cvMprRQ/vlw3NgoMzvTzBbUXRueb3LMYNT1d4bnpzVx7ggzmxveb5OZPWFm\nk1oYA9fXzK4M62uzmW01s/fM7EEzOzbG929mdoaZzTGzVWH9V5jZa2b2DTPrkL9b68bvAdPDQ++H\n73G3ujGzdDP7upm9En4/tpvZYjP7lpk16sExswPN7Bdm9rqZfRqO81tjZreYWWHDGAi6agGOiHp+\n/ecXfs7NjhNsbhyima0Nj5uZ/Xf4GW03s8UNyo0PvxPrwlg/CT+L/Zp53hfM7N9m9lFYvsTM5pvZ\nNc1WtkgKUHetSPf0JyALeBtYCuQBk4GfAUeZ2Qx3r+mIB5vZ8cDDBH9/vEbQmjQRmAfc2crlVwIX\nAC8B/wSGtzOWk4AHgDRgAbAWOACYD/y5ifJpwDPAFGAj8DxQGcbxBYKxbf+K4dG9gHuBUmA58AYw\nEDgU+D3BZ3FunG+rJYuBu4CZwB7A34Gt4bmPAcwsC3icIBEsI6iXSoL3/Btgupmd6O61Ufe9AjiZ\n4Ls0Pzx2IPB14ItmNsndS6Ji+HtY/hPgqaj7zCcx/gCcB7wAvANk1p0wsy8Ccwg+g8Xh+xsOnAr8\np5kd5+7zospfAtwM1BB8714g+KzGAz8GfpKgmEW6HCV5It3TRcC/3X1H3QEzyyFIPI4HzgLuTvRD\nw2fcQfB3x1fd/c/hcQP+B/h+K7f4CnCku7+QgFhygdsJEryz3P3eqHPXAlc1cdlUgmRnITDV3Ssb\n3G/vGB9fDZwIPO7uu6LuMQh4AjjHzO6ITjYSwd0fAR4JW8j2AL7j7msbFPslQYJ3P3CRu28OY8sh\nSI5mAxcSJFJ1/g+41N0/iXovEeBHBEnQdcBX62IIW9ZOBla4+7mJfI+hk4Aid18WfdCCCR/3ALuA\n4939mahzMwkmodxjZmPcvSo89T3AgUPcfVFUeQOmdUDsIl2GumtFuiF3fzQ6wQuPVQDfDn89oYMe\nfQowCJhXl+CFz3aCpGpdK9f/KREJXuhUIB94NjrBC/2EoFWvoUHh60vRCR6Au29x99djebC7V7v7\nI9EJXnj8U4LWSui4z6BZZjYY+BrB53BeXYIXxlYBnA9UEbTQEXXuuegELzxW6+7XAh/R+e/l5w0T\nvNC3gD7AldEJHoC7PwXcStCqNyvq1CCgPDrBC8u7uz+X2LBFuha15Il0U2a2N0EX4xiC//FFAAtP\nx9oi1VaHh69zGp5w911m9neC/xE357EExnJY+PpgE7HUmNlDwGUNTi0GaoHzzGw58JC7l8YbgJkd\nCMwARgLZBPWfE57uqM+gJdOADOCphv8IAHD3j83sXeAAM8tq0BI8gKCVb3+gH0ELKeH98s0s393L\nOvoNhJr7nswIXx9q5vyLwDcJussfDo+9DnzezP4E/Mbd305YlCJdnJI8kW4m7Gb6JUGrnTVTLKeZ\n4+1VNwj/g2bOr23l+g8TFwoF4WtzrYeNnuXuq8zsu8ANwG3AH8zsbeBZ4E53XxrLg80sk2D84Rkt\nFOuoz6Alo8LXr5nZ11opm0/QSoeZnUFQH31bKJ9DMMavMzT3PRkVvn4U/GfQrIFRf74EeISgu/mr\nZvYJwbi8h4C/ddTYVZGuQEmeSPdzGkELVTFBq9krwKdhS1omsJPmk7/2au99K1sv0qSmhpbUxeJt\nuZG7/9rMHiRYa+4YgtbJbwPfMrNvu/uNMdzmMoIE7y2CMV9vAJvCz2AfYCUd9xm0pK6eFgNLWim7\nE8DMRvLZhJlvEUza+Kiulc/MXiaYUJLI99PiUKGGXelNXHdXK/d/NepeS81sX4LJKl8gaO08Nfx5\nxcymu/vOWIIW6W6U5Il0PyeGrxe7++MNzu3Zwc+um2E5spnzzR1vTd0g+eZakpqahVsXy4g2XAOA\nu68Dfgf8LlxS5HSC2bi/MLO73X1TK/HWfQZnNDF2rKM/g5YUh6/z3f2/Y7zmCwSzV3/ZTIIbz/uJ\n5/OMRTGwF3B5W7rZw6TxkfCHcKmVewmS1/OBW+KMR6RL08QLke6nf/ha3MS5Uzv42S8295wwWTo5\nzvuuD1/3aeK++cDnmrjmpfC10TPDpVJOiuXB4SSKewhm3GYS21i6ZH4GLXmOYKmQ480sI8Zrmn0v\nZjaVYBZvQ3VJXHMNBRsJZsCObrguXxjXETHG1tDT4euJLZZqRZiY/z78df/23EukK1OSJ9L9rApf\nL7SogUlmdjjw3Q5+9oMEa8NNM7Nzop5tBDNam2tVa5G7v08wDusAM6ufyWlmfQjGiuU2E0sZcIyZ\nnd7g3I+A0Q0vMLPpZna0NVis2MxGE6yb5jSduDVU9xlc3OA+pxAsE5MU7v4RwRI3o4D7zKxRgmbB\nwsenRR2qey9nh/VdV24ouy+zEq0uidsrTKgbxlFFsH5dPsGYuLp7pgO/oonPJka/AnYAvwzXSNyN\nBQtdn2Nmw8Lfs83sm2bWr0G5CEH3LbQ+I1yk21KSJ9L93ESwaO83gLfN7D4zm0cwmLy5/yknRNQy\nHDXAnRbsXHEvwaLM3yVYty5ePw5f/27BDhaPAe8BE2hitmW4PMjXwljuM7OXzexeM3sL+AFBcgif\ntTpBsGjz08DHZvakmd1jZv8iWHC3H3Bz1KK/LflF+NwbzGxR+NyFBInnb9r2thPuUoL3eDLwXriz\nwxwze9bM3gfeZPcJI48By4BJwGoz+5uZ/ZMg+dsEvNzwAWES9xQwBFhiZneb2R/N7LyoYtcSzGT+\nbfjZPETweZ5B62PqmuTuq8PrMwi+J++a2aNm9pCZLSJYnPlOPpt4kQncCHxiwe4f94UzwNcStAau\n5bPviUjKUZIn0s24+yrgYOAfBP8zm00w9ukid+/oljzc/VGCxXafI+jqmkXQ3XoETSQEbbjvnwl2\nOXiHYHmUyQTv8VCCZKOpax4CjibYuWJCGEsJwWSKuhma0WO3/kmwsO8qgoTvS8B+BDs1nEKQIMUS\n6zzg88BcgjFrxxMkkyfzWTdgUoQTJo4DziGYgDCeoOt6HMHndA3BZJG68lUE9XUrwcSY48Nrfkcw\nMWW3tQCjXAD8BRgAnEmQ/Nd3w4br2M0m6Ab/XHhuAcF3d2073t+jBJ/1LQQtrzMIvgN9CRLWUwl2\nIYFgN5BLCD73QWE8RxJ8n64BDmrPEjoiXZ0Fa5iKiKQWM3sKOJZgp4NXWysvIpJqNLtWRLqtcNxY\ndRPbcV1KkOC9S7C/rohIj6MkT0S6s8MJ9ip9k2CB5l4EXcijgO3A+a7uChHpodRdKyLdVri125UE\nyd4eQG/gY4Ixeje4+/LmrxYRSW1K8kRERERSkGbXioiIiKQgjckDBgwY4CNGxLWGKwCRSITa2toE\nRtQ9qR4CqgfVAUBZWRkA+fn5SY4kufRdCKgeVAd1ElEPixcv3ujug1orpyQPGDFiBHPnzo37+pyc\nHCoqKhIYUfekegioHlQHAHPmzCESiXDqqcnc5Sz59F0IqB5UB3USUQ/5+fkfxFJO3bUiIiIiKUhJ\nnoiIiEgKUpInIiIikoI0Jk9ERKSNampqKCsrY9eu5rb23Z2Z0dOXLFMdBNpSDxkZGeTn55OWlhbX\ns5TkiYiItFFZWRn9+vWjf//+mFmr5dPS0qipqemEyLou1UEg1npwd8rKyigrK2PQoFYn0jZJ3bUi\nIiJttGvXrpgTPJF4mBn5+fkxtxY3RUmeiIhIHJTgSUdr73dMSZ6IiEg31KtXLw466KD6n7Vr1/LC\nCy8wYMAAJk2axH777cf06dN5/PHH66+59tprGTlyJAcddBATJkxgzpw5Td77gw8+YMaMGRQVFXHU\nUUdRXFxcf27WrFkMHDiQE044odnYLr/88vq49t13XwYOHFh/7oorrmDixIkccMABfOtb36ofn3bU\nUUex33771V+3YcOG3e7597//nYyMDBYtWgTA2rVrycnJqS//jW98o77sAw88QFFRERMnTuSKK66o\nP37XXXdRUFBQf82f/vSn3Z6xZcsWRo4cyTe/+U0Atm/fzuzZs9l///2ZOHEiP/jBD+rLfvjhhxx9\n9NFMmjSJoqIinnzyyWbrA+DTTz/liCOO4IADDuDRRx+tP37SSSdRUlLS4rXx0pg8ERGRbigrK4vX\nX399t2MffPABn//85+uTiMWLF3PKKaeQlZXFkUceCcCll17KZZddxrvvvsuUKVM4+eSTycjI2O0+\n3//+9zn77LP5yle+wnPPPccPf/hD7rrrLiBI4LZv387tt9/ebGy/+tWv6v988803s3jxYgBefvll\nXn75Zd544w0AjjjiCObNm8cRRxwBBEnYpEmTGt2voqKCm2++mcmTJ+92fK+99mpUB6WlpVxxxRW8\n+uqrDBo0iPPOO4+5c+fWv/8vfelL3HTTTU3Gfc011zB16tTdjl122WVMmzaNqqoqZsyYwVNPPcXM\nmTO5/vrrOeWUU7j44otZvnw5s2fPZvXq1c3WyZw5c/jyl7/MmWeeycyZMznhhBP45z//SVFREYWF\nhc1e1x5qyRMREUlRBx54ID/84Q+55ZZbGp3be++9yc7OZtOmTY3OvfPOO/VJ0bRp0/jHP/5Rf+7I\nI48kJycn5hjuv/9+Tj/9dCDofqysrKSqqoqdO3eya9cuBg8e3Oo9rrnmGr7zne/Qu3fvVsuuWbOG\nvffeu36ywlFHHcVDDz3U6nWvv/46GzZs4Oijj64/lp2dzbRp0wDIzMykqKiovlXTzOp3rti8eTMF\nBQUt3j8jI4MdO3awc+dOIpEI1dXV3HTTTVx++eWtxhYvJXkiIiLtlJ6R0eyPRbV42e23t1i2LXbs\n2FHf7XjKKac0W66oqIiVK1c2Ov7GG28wZsyYJpOsCRMm1CdGjzzyCBUVFZSWlrYpPghaFteuXcv0\n6dMBOPTQQ5k2bRrDhw9n+PDhzJgxg/Hjx9eXv+CCCzjooIP42c9+Vt+N++abb1JcXMysWbMa3f/9\n999n0qRJHHnkkcyfPx+AMWPGsHLlStauXUt1dTWPPfbYbt3NDz/8MEVFRZx22mmsW7cOgNraWr73\nve9xww03NPteysvLefzxx+uT36uvvpq//vWvjBo1itmzZ/Pb3/62xbo444wzePrppznuuOO4+uqr\nufXWWzn77LPJzs6OpSrj0uWSPDO7w8w2mNnbzZz/rpktDn/eNrMaM8sPz601s7fCc4s6N3IREZHO\nU9dd+/rrr/O3v/2t2XIN12S78cYb2W+//TjssMO4+uqrm7zm5z//OS+++CKTJk1i3rx5DB06lPT0\nto/weuCBBzjppJPq13lbvXo1K1asYO3atXzwwQc899xzvPjiiwDcfffdLF68mOeff5758+dzzz33\nUFtby3e+8x1+8YtfNLp3QUEBa9asYdGiRfzv//4vX/7yl9myZQv9+/fn5ptv5swzz2TatGmMHDmy\nPvbjjz+e1atX8+abb3LkkUfy1a9+FYBbb72V4447juHDhzf5Pqqrqzn77LO55JJL2HPPPYGg+/Wc\nc85h7dq1PPbYY5x33nnU1tY2Wxd5eXk89thjLFy4kKKiIp544glOOukkLrroIk477TReeeWVNtdv\na7pckgfcCcxs7qS7/6+7H+juBwJXAi+4e1lUkenh+cad+iIiIh2geteuZn/8a1+rL+df+1qLZTvC\n4sWLGTduXP3vl156KcuWLePee+/lvPPOo7KystE1hYWFPPjggyxatIif/vSnQJCktNX999/Paaed\nVv/7o48+ypQpU+jbty99+/Zl5syZvPrqqwAMHToUgJycHE4//XQWLlxIRUUFy5Yt4+ijj2bMmDG8\n+uqrnHTSSSxatIhevXoxYMAAAA466CD23HNPVq1aBQTJ3Msvv8z8+fPZZ599GDNmDAADBgygV69e\nQNBqWDc2cMGCBdxyyy2MGTOG73//+9xzzz27TbK4+OKLGTNmDJdeemn9sTvvvLO+BfXQQw+lsrKS\njRs3xlQv1113HVdeeSVz5szhc5/7HLfffjtXXXVVm+u3NV0uyXP3eUBZqwUDZwD3dWA4IiIi3dbS\npUu5/vrr+frXv97o3IknnshBBx3E3Xff3ejcxo0b61ulfv7zn3Puuee2+dkrV66kvLycQw89tP7Y\n8OHDmTdvHtXV1ezatYt58+Yxbtw4qqur6xOkXbt28cQTT7DffvuRl5fHxx9/zOrVq1m9ejVTpkzh\noYceYtKkSXz66af1iwqvWbOG1atX17ey1c3M3bRpE3/4wx/qW+zWr19fH8s//vGP+uT3L3/5S/09\nfv7zn3P22Wdz/fXXA0G37JYtW/j1r3+92/sbPnw4c+fOBYIxjJWVlTEtWvzuu++yfv16pk6dyvbt\n24lEIvVjFROt286uNbNsgha//4o67MC/zcyB/3P321q4/kLgQoARI0a0aRBpQ2lpae26PlWoHgKq\nB9UBUP8Xd0+vh1T9LphZm7eaindrqrbcMxKJMH/+fA4++GC2b9/O4MGDufHGGznmmGOAIO7o2K++\n+mrOOussLrroIiKRz9p9XnzxRX7wgx9gZhx++OH8/ve/r79m6tSprFixgq1btzJq1Cj++Mc/cuyx\nx3L11VczadIkZs+eDQRdtaeddtpu3bynnnoqzz//PEVFRZgZxx57LCeccALbtm1j1qxZ7Nq1i5qa\nGo466iguuuiiRu+vLva0tDReeuklrrnmGtLT00lLS+PWW2+tT7Iuv/xylixZAsBVV11VP+7v97//\nPf/4xz9IT08nPz+fO++8s8k6rHtOcXEx//M//8O4cePqZ/ZecsklXHDBBfzqV7/iwgsv5KabbsLM\n+POf/xxTl/aPfvQjrrvuOtLS0jjrrLM48cQTufnmm/nJT37S5HekPX+PWFfcR87MRgH/dPf9Wyhz\nGnC2u/9n1LFCdy8xs8HA08B/hy2DLSoqKvK6bDweOTk59TNsejLVQ0D1oDqAYLxOJBLh1FNPTXYo\nSZWq34WSkhLGjh0bc3lt6aU6qNPWeli5cmWjJVby8/Nfj2VYWpfrrm2D02nQVevuJeHrBuBhYHIT\n14mIiIikvG6Z5JlZHnAE8GjUsT5mllP3Z2AG0OQMXREREZFU1+XG5JnZfcA0YKCZFQPXABkA7v6H\nsNiJwL/dfVvUpXsAD4f7vKUD97r7U50Vt4iIiEhX0uWSPHc/I4YydxIstRJ9bA0wsWOiEhER2Z27\nt3sDeZGWtHfeRLfsrhUREUmmjIwMysrK2v0/YZHmuDtlZWWN9hVuiy7XkiciItLV5efnU1ZWFvPi\nt2bW4xNC1UGgLfWQkZFBfn5+3M9SkiciItJGaWlpMS18WydVl5JpC9VBoDPrQd21IiIiIilISZ6I\niIhIClKSJyIiIpKClOSJiIiIpCAleSIiIiIpSEmeiIiISApSkiciIiKSgpTkiYiIiKQgJXkiIiIi\nKUhJnoiIiEgKUpInIiIikoKU5ImIiIikICV5IiIiIilISZ6IiIhIClKSJyIiIpKClOSJiIiIpCAl\neSIiIiIpSEmeiIiISApSkiciIiKSgpTkiYiIiKQgJXkiIiIiKUhJnoiIiEgKUpInIiIikoKU5ImI\niIikICV5IiIiIilISZ6IiIhIClKSJyIiIpKClOSJiIiIpCAleSIiIiIpSEmeiIiISApSkiciIiKS\ngpTkiYiIiKSgLpfkmdkdZrbBzN5u5vw0M9tsZovDn6ujzs00s5VmttrMrui8qEVERES6li6X5AF3\nAjNbKfOiux8Y/lwLYGZpwO+B44B9gTPMbN8OjVRERESki+pySZ67zwPK4rh0MrDa3de4exUwBzgh\nocGJiIiIdBNdLsmL0aFmtsTMnjSz/cJjQ4F1UWWKw2MiIiIiPU56sgOIwxvASHffamZfAB4B9gas\nibLe3E3M7ELgQoARI0aQk5MTd0BpaWntuj5VqB4CqgfVAUAkEsHMenw96LsQUD2oDup0Zj10uyTP\n3bdE/fkJM7vFzAYStNwNjyo6DChp4T63AbcBFBUVeUVFRdwx5eTk0J7rU4XqIaB6UB0A1NbWEolE\nenw96LsQUD2oDup0Zj10u+5aMxtiZhb+eTLBeygFFgJ7m9loM8sETgceS16kIiIiIsnT5VryzOw+\nYBow0MyKgWuADAB3/wNwCvB1M6sGdgCnu7sD1Wb2X8C/gDTgDndfloS3ICIiPVRFRQXr169n586d\n5ObmUlBQQGZmZrLDkh6qyyV57n5GK+dvBm5u5twTwBMdEZeIiEhzampqWLhwIcuXL8fMiEQiVFdX\nk5mZyWGHHcbEiROTHaL0QF0uyRMREeluXn31VZYvX05hYSGRyGcjoaqqqnj22Wfp168f+fn5SYxQ\neqJuNyZPRESkK9myZQsrVqxolOABZGZmkp+fz4IFCwhGFol0HiV5IiIi7VBcXEwkEmmU4NXp06cP\nmzZtory8vJMjk55OSZ6IiEg77Nixg7S0tBbLRCIRqqqqOikikUBcY/LMrA9wGHAoUAgMAnoTLGXy\nKfAO8IK7r0pQnCIiIl1Sbm4uu3btava8u1NbW0tWVlYnRiXShiTPzDIIli+5kCDBi/5nS91uE97g\nmk+Ae4Hb3X1l+0IVERHpeoYNG4aZUV1dTXp64/+tlpeXM3ToUHJzc5MQnfRkrXbXmlmWmf0I+Ai4\nBziCIDmsAZYD84CHCZK5fxEsSryRIPEbAnwbWG5mT4eLF4uIiKSMrKwsDjnkENavX09lZeVu5zZv\n3kxVVRX/8R//kaTopCdrsSXPzM4FfkrQJWvA28B9BInd6+5e2cK1o4ApwAnAbOAo4EgzewD4nruv\na3/4IiIiyTd+/HgyMjJYtGgRZWVlmBnuzh577MHRRx/NgAEDtKWXdLrWumvvAKoI9nj9g7svifXG\n7r4WWAvcb2bZBF293wNOIxizd20c8YqIiHRJY8aMYc8996S0tJRdu3aRnZ1Nv379kh2W9GCtJXn/\nB1zf3lY3d98O3G1mfwG+RLhNmYiISCqJRCIMGjQo2WGIAK0kee7+9UQ+LNxj9oFE3lNEREREGtM6\neSIiIiIpSEmeiIiISApqc5JnZplmNsLMhjRxrq+Z/dLMlpjZG2b2UzPT6o8iIiIinSyeHS8uAH4H\n3AV8tcG5x4HP89niyBOBw81sumtnZhEREZFOE0937bHh673RB81sNnA4wa4XfwX+COwKj325HTGK\niIiISBvFk+SND19fb3D8TIIE7+fu/mV3vxD4FkGr3pnxhygiIiIibRVPkjcI2O7umxocnx6+/jHq\n2F/C14lxPEdERERE4hRPktcHqI0+EG5hNghY5+7v1x13921AOZAff4giIiIi0lbxJHllQF8zi96r\n5cjw9eUmyqcDW+N4joiIiIjEKZ4k743w9XwAM4uEf3bgueiCZjYI6At83I4YRURERKSN4kny7iSY\nTHGDmT0JvAYcStBa92CDsoeHr+/EG6CIiIiItF2bkzx3fwD4M5BGsJzK54BK4GJ3L29Q/DSaaOET\nERERkY4Vz2LIuPv5ZnYH8B8EEyueiZ5wAcHOGMBm4G7gifYGKiIiIiKxiyvJA3D3l4CXWjhfBVwY\n7/1FREREJH7xjMkTERERkS4u7pY8ADObQDAubySQ5e7nR53LIFg7z919fbuiFBEREZE2iSvJM7M8\n4A7gi3WHCCZYnB9VLANYAvQ3s33cfU17AhURERGR2LW5u9bM0gkmUnwR2A48TjC7djfuvp0gEYzw\nWTIoIiIiIp0gnjF55xOsi7cGGOvuswlm0Tbl7+HrcXE8R0RERETiFE+SdyZB1+y33b2klbJvEuxz\nu28czxERERGROMWT5B1AkOT9u7WC7r6LoJVvQBzPEREREZE4xZPk9QEqwnXwYpEJVMfxHBERERGJ\nUzxJ3kYg18z6tFbQzPYmSAq1hIqIiIhIJ4onyVsYvs6Koey3w9dmd8YQERERkcSLJ8m7k2BdvGvN\nbEhzhczsu8DFBOP37oj15mZ2h5ltMLO3mzl/lpktDX9eNrOJUefWmtlbZrbYzBbF+kwRERGRVNPm\nxZDd/REze5ygJW+Rmd0D9AYwsy8RTMz4ErBPeMkcd5/XhkfcCdwM3N3M+feBI9x9k5kdB9wGTIk6\nP93dN7bheSIiIiIpJ95tzU4H7gFOAL4bdXxO+Grh68PsvgtGq9x9npmNauH8y1G/LgCGteX+IiIi\nIj1BPN21uPs2dz8R+E/gMaCUILEzYAvhjhjufrK7N9oNI4HOB56MDg34t5m9bmYXduBzRURERLq0\neFvyAHD3xwm2Navb7izN3XcmIrDWmNl0giTv81GHD3P3EjMbDDxtZiua6yoOk8ALAUaMGEFOTk7c\nsaSlpbXr+lShegioHlQHAJFIBDPr8fWg70JA9aA6qNOZ9dCuJC+au1fTSevhmdkE4I/Ace5eGhVD\nSfi6wcweBiYDTSZ57n4bwXg+ioqKvKKiIu54cnJyaM/1qUL1EFA9qA4AamtriUQiPb4e9F0IqB5U\nB3U6sx7i6q5NJjMbATwEfNndV0Ud72NmOXV/BmYATc7QFREREUl1bW7JM7Op8Two1hm2ZnYfMA0Y\naGbFwDUUSySaAAAgAElEQVRARniPPwBXE2yTdouZAVS7+yRgD+Dh8Fg6cK+7PxVPrCIiIiLdXTzd\ntc8TTHBoC4/1We5+RivnLwAuaOL4GmBi4ytEREREep54x+RZ60XaVV5ERERE2qHNY/LcPdLSD9AP\nOBZ4ASgjWLi42439ExEREenOEp58ufsWd38aOJJgz9rHzGx0op8jIiIiIs3rsBY2d3fge0AecFVH\nPUdEREREGuvQblR3X0mwA8YxHfkcEREREdldwhZDboqZZQBZQO+OfI6IiIiI7K6jJ0R8kWCNuw0d\n/BwRERERiZLwljwzywSGAycDPyBYI+/JRD9HRERERJoXz44XNW0pDnwE/KStzxERERGR+MXTXWsx\n/lQC9wCHuHtJQqIVERERkZjE0107vZXz1cAmYJW7V8dxfxERERFppzYnee7+QkcEIiIiIiKJo+3G\nRERERFKQkjwRERGRFNRid62ZTU3Ug9x9XqLuJSIiIiIta21M3vME69y1l8fwLBERERFJkFgSL0vA\ncxJxDxERERGJUYtJnrtrzJ6IiIhIN6QkTkRERCQFKckTERERSUFK8kRERERSULtmvJpZJnAgMAzo\nQwsTLNz97vY8S0RERERiF1eSZ2a9gJ8BFxIkd61xQEmeiIiISCdpc5JnZunAv4DDCVruNgCDgVqg\nBBgI9A6LbwVKExKpiIiIiMQsnjF55wNTCRK6Se4+JDy+wd1HAH2B6cDLBEnkj9x9dCKCFREREZHY\nxNNdewZB9+sP3f2NhifdvRZ4wcyOAJ4E7jCzd5oqKyIiItJd1NbW8sknn7By5Uq2bNlCdnY248aN\no6CggLS0tGSH10g8LXn7h69/a3B8t3fn7jXAZUAG8J04niMiIiLSJVRXV/PCCy/wxBNPsH79eqqr\nqyktLeWpp57iqaeeorKyMtkhNhJPkpcDbHH37VHHdhJ00+7G3d8GKgjG74mIiIh0S2+++Sbvvfce\nhYWF9O/fn+zsbPLy8hg2bBgbN27kpZdeSnaIjcST5G0AshscKwN6m9ng6INmZkAmMCi+8ERERESS\nq7KykmXLljFkyBCC1GZ3gwcPZu3atZSXlychuubFk+QVA+lmNiTq2Nvh68wGZacBvYDNcTxHRERE\nJOlKS0upra0lPb3pqQxmhpmxYcOGTo6sZfEkeS+Er4dFHXuYYDmVX5rZl8xsbzM7BbiLYJLG3PaF\nKSIiIpIc7t5kC140M8PdOymi2MST5D1GkNCdFXXsTwSteQOBOcAK4H6CnTC2AT9pX5giIiIiyZGb\nm4u7t5jE1dbWkpeX14lRtS6eJO8VYDTw33UH3H0XcBRwH8EkjLp0dz4wzd1XtDNOERERkaTIzc1l\n5MiRlJY2vb/Dli1b6NevH4MHD27yfLK0OcnzwAfu/lGD45+6+1lALjAUyHX3qVofT0RERLq7Qw45\nhN69e7N+/XqqqqqAYFmVDRs2UFVVxfTp04lE4mk76zhx7V3bEnevBtYn+r4iIiIiydKnTx9mzZrF\nypUrWbZsGVVVVaSlpTFu3DjGjRtHbm5uskNsJJ69a38A3OPuH3ZAPCIiIiJdUlZWFgceeCATJkyg\nurqa9PT0Ltd6Fy2eyK4D1pjZXDM718xyEh2Umd1hZhvM7O1mzpuZ3WRmq81sqZl9LurcOWb2bvhz\nTqJjExERkZ4tEomQmZnZpRM8iC/JWxdeN41gVu3HZvZXM5tpZol6t3fSeM29aMcBe4c/FwK3AphZ\nPnANMAWYDFxjZv0TFJOIiIhItxHPxIuRBAneHcAWIAs4HXgcKDazX5rZxPYE5e7zCHbRaM4JwN3h\nJJAFQD8zKwCOBZ529zJ33wQ8TcvJooiIiEhKimviRZiEzTOzS4DZwFcIEqwhwLeBb4ddrXcD97p7\noidiDCVoUaxTHB5r7ngjZnYhQSsgI0aMICcn/l7ntLS0dl2fKlQPAdWD6gCC7hwz6/H1oO9CQPWg\nOqjTmfXQrtm17r4TeBB40MwGAmcAZwMHAwcAvwBuMLNn3P249gYbpallp72F440Put8G3AZQVFTk\nFRUVcQeTk5NDe65PFaqHgOpBdQDBwqiRSKTH14O+CwHVg+qgTmfWQ8JGDLr7Rnf/nbtPAcYBPwM+\nBNKAGYl6TqgYGB71+zCgpIXjIiIiIj1KR00LySFYFDm7g+7/GPCVcJbtIcDmsEv4X8AMM+sfTriY\nER4TERER6VESthiymQ0n6Kr9MjC27jBQRTApoy33uo9gcsdAMysmmDGbAeDufwCeAL4ArAa2A+eF\n58rM7KfAwvBW17p7SxM4RERERFJSu5K8cI28UwgSu6kESV3duLhXgb8Ac9qaaLn7Ga2cd+CSZs7d\nQTDzV0RERKTHimfHiwjBTNovEyxl0pvPErsPgHsIljd5N1FBioiIiEjbxNOS9xEwOPyzEayV9zeC\nxG5eogITERERkfjFk+TtAdQA/yZYB+9Rd69MaFQiIiIi0i7xJHmXA3919w2JDkZEREREEqPNSZ67\n/6YjAhERERGRxEnIEipmNpLPxultcPcPEnFfEREREYlP3Ishm1mhmf3OzDYAa4AF4c8aM/s0PDcs\nUYGKiIiISOziSvLMbAawDPgGMJDP1ser+xkQnnvbzGYmJlQRERERiVU86+SNBR4hWB+vDPgDMJdg\naRWAocB04CKCBPAhMyty95UJiVhEREREWhXPmLyrCBK8pcAx7v5pg/MrgblmdiPwDHAA8COCxZNF\nREREpBPE0117FODABU0kePXcfSPwNYLu26PjC09ERERE4hFPktcP2Orui1or6O4Lga3hNSIiIiLS\nSeJJ8tYDaW18xvo4niMiIiIicYonyXsCyDKzI1sraGZHAdnAP+N4joiIiIjEKZ4k76fABuBPZrZP\nc4XMbG/gdoJWvOviC09ERERE4hHP7NqxwJXAb4AlZvYgny2h4sAwgiVUvgRUApcB48xsXMMbufu8\nOOMWERERkRbEk+Q9T5DM1Tkr/GlKL+BPzZzzOJ8vIiIiIq2IN8myBDw7EfcQERERkSa0Oclz97j3\nuxURERGRzqHuUhHptnbt2sW6detYtmwZFRUV9O3bl/HjxzNy5EgyMzOTHZ6ISFIpyRORbqmyspJn\nn32W9evX069fP3Jzc9m5cyfz5s1jwIABzJgxg+zs7GSHKSKSNOp6FZFu6dVXX2Xjxo0MGzaMvn37\nkp6eTp8+fRg6dChbt25l/vz5yQ5RRCSpWkzyzOw7ZpaVyAea2cFmdlwi7ykiPcvWrVt57733GDRo\nUJPnBwwYwLp16ygvL+/kyEREuo7WWvJ+Aawxs2+bWbv2nzWzz5vZP4EFwMHtuZeI9GybNm3CzIhE\nmv4rzMwwM0pLSzs5MhGRrqO1JO96IBf4JbDezP5mZieb2eDWbmxmGWGr3U/N7D3gBeALwELgkfYG\nLiI9m7u3eN5MqzSJSM/W4sQLd/+Rmd1KkOydCZwEnAhgZuuAJcCnQBmwE+gP5AN7AhOBuultBrwH\nXOXucxL/NkSkJ8nLywOgtra22da82tpa+vfv35lhiYh0Ka3OrnX3j4BzzOxK4ELgqwRbl40If5r6\n53TdP6GrgceB/wP+5a3901tEJAa5ubmMHj2adevWMXhw446F0tJSCgsLyc/PT0J0IiJdQ8xLqLh7\nCfBj4Mdmtj8wFZgCFAKDgN5AKUHL3nJgHvCSu1ckOGYREaZMmUJ5eTnFxcXk5+fTu3dvqqqqKCsr\nIycnh8MPPzzZIYqIJFVc6+S5+9vA28AtiQ1HRCQ22dnZzJo1i/fff5+lS5eyYcMGsrOzOeSQQ9hz\nzz3p3bt3skMUEUkqLYYsIt1WZmYmY8eOZezYsckORUSky9FiyCIiIiIpKK6WPDPrBYwnGI+XEx6u\nAEqAd9x9Z2LCExEREZF4xJzkWbDo1LnAV4DP03wrYK2ZzQfuAu7SjFoRERGRzhdTd62ZDQfeBP4I\nHAGkESyT0tRPWljmT8Ab4bUiIiIi0olabckzs74Ey6GMDA89DzxJMLv2I2B7eDwbGArsD8wEphMs\niPy8mU10960JjVxEREREmhVLd+13CBK8EuAUd1/QQtmlBAng/5rZFODvwCjgMuDaWIMys5nAjQSt\ngn909xsanP8NQRIJQXI52N37hedqgLfCcx+6++xYnysiIiKSKmJJ8k4i2NXinFYSvN24+6tmdg7w\nNHAyMSZ5ZpYG/B44BigGFprZY+6+POre344q/99AUdQtdrj7gbHGKSIiIpKKYhmTNxrY7u7PtvXm\n4TXbCfayjdVkYLW7r3H3KmAOcEIL5c8A7mtrbCIiIiKpLJaWvBogM56bhzNyI0BVGy4bCqyL+r2Y\nYPu0pu4/kiAJnRt1uLeZLSLYN/cGd3+kmWsvJNiLlxEjRpCTk9NUsZikpaW16/pUoXoIqB5UBwCR\nSAQz6/H1oO9CQPWgOqjTmfUQS5K3CjjIzE5297+38f4nE+xp+3YbrrEmjjW3DMvpwN/cvSbq2Ah3\nLzGzPYG5ZvaWu7/X6IbutwG3ARQVFXlFRfxb7Obk5NCe61OF6iGgelAdANTW1hKJRHp8Pei7EFA9\nqA7qdGY9xNJdO4cg8fqjmX0p1hub2SnA7QQJWlu6U4uB6GVXhhFM+mjK6Q3v7e4l4esagpnARY0v\nExEREUltsbTk/Q44GzgQmGNm1wNP0fISKscCexEkh2+E94jVQmBvMxsd3v904MyGhcxsLNAfeCXq\nWH+C8YM7zWwgcBjwizY8W0RERCQltJrkufsuMzsSuJdg/bu9gG+0clldl+uTwFnuXh1rQO5ebWb/\nBfyLYAmVO9x9mZldCyxy98fComcAcxrsqDEe+D8zqyVopbwhelauiIiISE8R07Zm7l4OfCFM9s4h\nWN5kSDPFPyZYNuUud5/bTJnWnvcE8ESDY1c3+P3HTVz3MnBAPM8UERERSSUx710LECZtc6G+a3Qo\nUDdFpAL4yN03JTRCEREREWmzNiV50cJkTgmdiIiISBcUy+xaEREREelmlOSJiIiIpKAOTfLMrJ+Z\nvW9mjRYjFhEREZGOE/eYvBilASNpfscKERERkW6ptLSUVatW8fHHHxOJRNhrr70YPXo0ffr0SXZo\nQMcneSIiIiIp56233uK1114jMzOTnJwcampqWLhwIW+++SbHHnssgwcPTnaIGpMnIiIi0hbFxcUs\nWLCAPfbYg0GDBtG7d2+ys7MpKCggKyuLp59+mh07diQ7zNZb8swsrgWNQxntuFZERESky1m6dCl5\neXmkpzdOo/r27cuWLVv44IMPGDduXBKi+0ws3bXTCMbUWSvlRERERFLazp07+fjjjykoKGi2TG5u\nLmvXru0WSV4NQbfuvwi2LGuLXsDpbQ1KRERS065duygpKaGiooLMzEwKCwvp27dvssMSiZm7Y2aY\nNd/2ZWbU1tZ2YlRNiyXJWwHsCzzo7ne05eZmNgAleSIiAqxZs4aXXnqJqqoq0tPTqa2txd0ZP348\nkydPJi0tLdkhirSqV69e5Obmsn37drKzs5ssU1FRwZgxYzo5ssZimXjxevh6UEcGIiIiqau4uJhn\nn32W3Nxchg4dyh577EFBQQFDhgxh2bJlvPrqq8kOUSQmZsaECRMoLS3FvfEKcVVVVQCMHj26s0Nr\nJJYk7w2C8XhFHRyLiIikIHfntddeIz8/n169eu12LhKJUFhYyIoVK9iyZUuSIhRpm7322otx48ZR\nXFzM1q1bAaitraWsrIxPP/2UqVOnkpubm+QoY+uufSN8nWBm5k2lrc2rBO5GiyGLiPRY5eXllJeX\nU1hY2OT5SCRCJBKhuLiYfffdt5OjE2m7SCTCYYcdxrBhw1i6dCklJSWYGaNHj2a//fZj0KBByQ4R\niCHJc/f5xLmenrtvA86N51oREUkNVVVVLQ5SB0hPT+8S64qJxCoSiTB69GhGjx5NTU0NZkYk0rWW\nH9aOFyIi0qGysrJw9/pZiU2pqqoiJyenkyMTSYyuOmmoa6WcIiKScnJzcykoKKC8vLzJ89XV1ZgZ\nw4cP7+TIRFKbkjwREelwU6ZMoaqqis2bN+92vLKykvXr13PIIYeQlZWVpOhEUpO6a0VEpMPl5+dz\n/PHHs2DBgvpB6u5OdnY206ZNY++99052iCIpp81JnpnVtPGSnUA5sAx4Evizu29q63NFRKR7GzBg\nALNmzaK8vJzt27eTnp7OwIEDu9xgdZFUEU9LXlv3sO0NDAl/jgS+a2anuvuLcTxbRFLIli1bWLly\nJe+99x41NTUMGTKEfffdlyFDhrQ6G1O6r379+tGvX79khyGS8uJJ8qYDo4BfA9nAA8DzwEcECWAh\nMA04FdgGXAZsAQ4Gzgf2AB41s/3cfX27oheRbqu4uJhnnnmGSCRCv379iEQibNiwgTVr1jBx4kQO\nPvhgJXoiIu0QT5K3HLgX2Awc6u6rmijzZzO7DngK+CnwOXd/xMx+C8wD9gG+CVwZX9gi0p1t27aN\nuXPnkpeXt9tg+/79+5OXl8eSJUsYPHgwo0aNSl6QIiLdXDwDIa4i6Hq9oJkEDwB3fxf4GjAc+EF4\n7FPgcoIWv5lxPFtEUsDatWupqalpcjZlJBKhf//+LFmyJAmRiYikjniSvFlApbvPba1gWGY7cELU\n4WeAaiD5O/eKSFJ8+OGHLS5827dvX8rKyti5c2cnRiUiklriSfIKCZK0WNUAQ+t+cfcqgjF6vZq9\nQkR6vLodEkREJD7xJHnlQF8zO7C1gmGZHILxe3XH0oA8oDSOZ4tIChg+fDgVFRXNnt+6dSsDBgyg\nVy/9W1BEJF7xJHkvEIypu83M8porFJ67DXCC2bd1RgFpBLNxRaQHGj16NJFIhMrKykbnamtr2bRp\nExMnTtTsWhGRdognybsO2AUcBKwws6vM7Agz29vMxpjZVDO7GngHmETQtfuzqOtPCV/ntSdwEem+\n+vTpw1FHHUV5eTkbNmygqqqKmpoaysvLKSkpYeLEiYwcOTLZYYqIdGttXkLF3d82s9OBvxCseffj\nZooasAP4iru/FXV8O3AjcE9bny0iqWPYsGGceOKJ9YshV1dXU1BQwNSpUykoKFArnohIO8W1d627\nP2xmBwA/BE4CGi5dXg48BFzv7msaXPu7eJ4pIqknLy+PyZMnM3ny5GSHIiKScuJK8gDc/X3gAuAC\nM9sTGETQerehYWInIiIiIp0r7iQvWpjUKbETERER6SLimXghIiIiIl1c3EmeBU4yswfN7H0z22Zm\nW8M/P2BmX7Q4R06b2UwzW2lmq83siibOn2tmn5rZ4vDngqhz55jZu+HPOfG+PxEREZHuLK7uWjPb\nA/gb8B91h6JOjwRGACcDL5nZqe7+cRvunQb8HjgGKAYWmtlj7r68QdH73f2/GlybD1xDsHSLA6+H\n126K/d2JiIiIdH9tTvLMLBP4F3AAQXL3GvA0QUIGMAw4GpgCHAY8aWZTwu3MYjEZWF03ecPM5hDs\nfdswyWvKscDT7l4WXvs0MBO4L8Zni4iIiKSEeFryvg5MINh/9mx3/2cTZa4ysy8A94ZlLwZuivH+\nQ4F1Ub8XEySMDZ1sZlOBVcC33X1dM9cObeJazOxC4EKAESNGtLhZemvS0tLadX2qUD0EVA9dpw7q\nFlh2d3Jzc8nMzOy0Z0ciEcysS9RDMnWV70KyqR5UB3U6sx7iSfJOJegKvaSZBA8Ad3/CzC4hWDT5\nNGJP8poax9dwl/J/APe5+04zuxi4Czgyxmvr4ruNYNs1ioqKvKV9NFuTk5PT4j6cPYXqIaB6SH4d\n1NbWsnz5cpYsWcLOnTsxM9LS0hg/fjwTJ07slGSvtraWSCSi74L+ewBUD6A6qNOZ9RDPxIvxBNua\n3R9D2fuBqvCaWBUDw6N+HwaURBdw91J33xn+ejvBFmsxXSsiqc3deeWVV3jllVfIycmhsLCQgoIC\n8vPzeeutt3jmmWfYtWtXssMU6TG2bNnCmjVrWLVqFRs2bMC9ybYX6QDxtORlAdvdvbq1gu5ebWbb\nw2titRDY28xGAx8BpwNnRhcwswJ3Xx/+Optgn1wIxgpeb2b9w99nAFe24dki0s198sknvPPOOwwd\nOpRI5LN/x6anp1NYWEhxcTFr1qxh7NixSYxSJPVVVVXxyiuvsHr1asyM3r17s2PHDvr378/06dPp\n379/6zeRdoknyfsEGG5mI9z9w5YKmtkogi3PWiwXLUwM/4sgYUsD7nD3ZWZ2LbDI3R8Dvmlms4Fq\noAw4N7y2zMx+SpAoAlxbNwlDRHqGFStWkJ2dvVuCF23AgAEsXbpUSV4XUFlZyQcffEBJSQmRSIQR\nI0YwbNgwMjIykh2atFNtbS3PPfccJSUlFBYWYmZkZWWxY8cONm/ezOOPP87s2bPJzc1NdqgpLZ4k\nbx5wNvAbMzvFm2l3DdfI+zXBmLgX2vIAd38CeKLBsauj/nwlzbTQufsdwB1teZ6IpI7S0lKys7Ob\nPZ+VlUVJSQk1NTWkpaV1YmQSrbi4mLlz51JdXU2fPn1wd1avXk2fPn045phjGDBgQLJDlHZYv349\nxcXFDB3aeO5jXl4eGzZsYPny5RxyyCFJiK7niGdMXl3i9kXgOTM7yszq/9llZhlmdjTwXFjGgd8k\nIlgRkdb06tWL6urmR5NUV1eTnp7ebEufdLxNmzbx9NNPk5ubS2FhIXl5efTr14+hQ4eSlpbGU089\nxY4dO5IdprTDu+++2+I/tvLz81m5cmWL/61K+7X5bzl3XwxcTjCT9XDg38BWM/vIzIqBrQRdrVPD\nS74TXiMi0uHGjh3L5s2bmz2/adMmxo4dS5wb8kgCvPPOO6Snp9O7d+9G53JycqiqquL9999PQmSS\nKNu2bWtxFnt6ejq1tbWaBNXB4vqnrLv/lmDCwwqCZC8DKAAKwz8bweLF/xmWFRHpFHXrXm7a1Hij\nm23btlFTU8O4ceOSEJlAMPv53XffJT8/v9kyeXl5rFy5shOjkkTLycmhsrKy2fPV1dWkpaV16tqV\nPVFc25oBhGvk/dPMDiDYRmxweGoDwQSJtxIQn4hIm/Tq1YuZM2fyzDPP8NFHH5GZmUkkEmHnzp30\n7t2b4447jry8vGSH2WO5e6vjIdPT06mqinWTJOmK9tlnH959913cvclW840bN7L//vtrXGwHizvJ\nqxMmc0roRKTLyM3N5Ytf/CIff/wxJSUl1NbWMnjwYIYOHaqZm0kWiUTIz89n69at9O3bt8kyFRUV\nDBs2rJMjk0TaY4892HPPPVmzZg1DhgypT+bcnbKyMrKyshg/vi1L6Eo82p3kiYh0RZFIhMLCQgoL\nC5MdijQwYcIEnn32Wfr06dOolae2tpbt27erS72bMzMOP/xw+vTpw7Jly3D3+iVUhg4dyqGHHkqf\nPn2SHWbKU5InIiKdauTIkey1116sWbOGgQMH1k/A2L59O6WlpUycOJHBgwe3chfp6tLT05k8eTIT\nJkxg48aNZGVlkZaWRr9+/ZIdWo/RYpJnZnMT9Bx396MSdC8REenG0tLSOOKIIxgyZAhLly6lrKwM\nMyM3N5fp06ez1157afZzCunduzfDhg3T3rVJ0FpL3rQEPUcb1YmI9DCVlZWkp6c3Ofg+LS2Nfffd\nl3HjxrFt2zbMrMnuWxGJX2tJ3k86JQoREUkZxcXFLF68mE8//ZRevXqRkZHBxIkTGTNmTKNFqCOR\nCDk5OUmKVCS1tZjkubuSPBERidlbb73FggULyMvLo6CggKysLEpLS3nhhRdYv349hx9+uHYbEekk\n+i9NREQSoqysjNdee42CgoLdNp7Pzs5m2LBhrFq1irVr1yYvQJEeJmFJnpm9ZmbvJep+IiLSvaxa\ntYpevXqRnt64k8jM6N+/P0uWLElCZCI9UyJb8oYDoxJ4PxER6UbWr1/f7ALHAH379mXTpk3alF6k\nk6i7VkREEqJu0/nm1NbWYmYakyfSSfRfmoiIJMRee+3F5s2bmz1fXl7O6NGjleSJdBL9lyYiIgkx\natQoMjIy2Lp1a6NzVVVV7Nixg/322y8JkYn0TEryREQkIbKzs5k5cyY7d+6kpKSELVu2sHXrVj75\n5BPKysqYPn06gwYNSnaYIj1GIveufQXon8D7iYhINzNw4EBOPvlkPvzwQ9auXfv/7d15dBzXndj7\n76+70VgbxEYC4AqQBLiIAkmJpETJsqh9i6nEsWc0jrdnTTyO40wyfhM/z7Fnotgz73nsZ8/kHdsT\ne7zInkRR7HhjtMsWSYuSKIqkuEMUKRAksRE70ECj0eiu+/6o6mYD6AaaWBpo8Pc5p04DVbeqfn1R\nQP9Qt+69eL1eVqxYQXV1tQ56rFSazViSZ4z54EwdSymlVObKycmhtraW2tpana9UqTmkzbVKKaWU\nUguQJnlKKaWUUguQJnlKKaWUUgvQTHa8UEoppaZteHiYCxcucPLkSQKBAAUFBWzatInq6mq8Xu9c\nh6dUxtA7eUoppeaNwcFBnnnmGV5//XU8Hg+LFy9GRDhw4ADPPfccgUBgrkNUKmNokqeUUmreOHDg\nAENDQyxdupTc3FxcLhd5eXksW7YMv9/PwYMH5zpEpTKGNtcqpdQMsiyL9vZ2Ojs7ycrKIhwO4/Ho\nn9pU9PT00NTUxNKlSxNuLysro7Gxkf7+fgoLC6/p2IFAgM7OTizLYtGiRRQXL6xhXbu7u+nv78ft\ndlNaWkpeXt5ch6TmAf3Lo5RSM6SpqYkDBw4QCAQ4c+YMIsLTTz/Njh07qK2tnevw5r2enh5EBBFJ\nuD26vqenJ+Ukb2RkhMOHD/POO+8AYIzBGENlZSXve9/7rjlZnG/6+vo4cOAAV65cia0TETZs2MC2\nbdv0H4zrnP70lVJqBjQ1NfHCCy9QUlJCUVERRUVFuFwuCgsL2b9/P5FIhA0bNsx1mPNasuRuqizL\n4ve//z0XL16koqICl+vqE0rd3d0899xzfOADHyA/P39Gz5suAwMDPPvss4TDYQD6+/txuVyUlJTE\nOiOMhFAAACAASURBVK3cddddM16vKnNokqeUUtNkWRZvvPEGJSUl45rJvF4vFRUVHDp0iNWrV5Od\nnT1HUc5/xcXFsTttiRITy7Ji5VLR1tZGY2Mjy5YtG7etpKSEtrY26uvr2bZt2/QCnyNnzpyhvb2d\n1tZWwuEwWVlZGGO4ePEiPp+PkZERbrjhBsrLy+c6VDVHNMlTSi14w8PDXLx4kfr6eoLBIMXFxWzc\nuJGlS5eOurszVV1dXfT39ydMJoDYs3mtra1UVVVN+3wLVVFREatWraK1tZXFixeP297R0UF1dXXK\nTazvvvvuhM+mlZWVcebMGW666aYZuQ7SKRwOc/jwYRobGykqKiIrK2vUdr/fz/nz56mvr9ck7zqW\nWVe1Ukpdo/7+fvbs2cOBAwcIh8Pk5ubS09PDCy+8wL59+2JNXdMxPDw8aZLgcrkIBoPTPtdCd9tt\nt+Hz+WhubmZwcJBwOMzAwABNTU0UFxezc+fOlI/l9/snvHPq8XiIRCKEQqGZCD2hcDjM4ODgjJ9j\nZGSES5cukZOTMy7BA/D5fIRCIRobG2f0vHMlWo8jIyNzHUpG0Tt5SqkFyxjDvn37GBkZGdVj0+v1\nUlhYSENDA0VFRdx0003TOo/X6401JSZjWRY5OTnTOs/1IC8vj4cffpiLFy9y+vRp+vr68Pl87Nq1\ni1WrViVMaJLx+Xy0trYmvZsXDodxu93XdMxUBQIB6uvrOXPmDOFwmJycHCorK9m8efOM9Ow1xtDd\n3c3y5cuTlnG73fT29k77XHPJ7/dz5swZ3nnnHSzLQkSoqanhxhtvzPhOM+mgSZ5SasFqb2+no6Mj\nYTOqiFBRUcHp06fZtGnTtGZSKCsrw+fzEQgEEiYU0WSisrJyyue4nni9XmpqaqipqZnWcWpqanjv\nvfeSJlVdXV1s3LgRt9s9rfOMNTg4yPPPP8/AwABlZWV4PB6ys7Npbm6msbGRBx98kIqKimmdQ0Qo\nLy9nYGCARYsWJSwTiUQSNntniv7+fp555hnC4TClpaWxO68NDQ00NDTwyCOPUFJSMtdhzmvaXKuU\nWrCuXLky4RASHo+HcDg87bsdLpeLnTt30tXVNa5JNvos3rZt27TTRZpVVlayYsUKWltbx91p7e3t\nJSsri/Xr18/4eQ8dOkQgEKCioiJ2/blcLkpLSyksLGTv3r3TfkwgOzubqqoqLMsaNwuIMYaenh58\nPh9r1qyZ1nnm0oEDBxARlixZEqtHt9vN4sWLyc7OZt++fRhj5jjK+W1e3skTkQeB/wK4gR8YY742\nZvvngT8GwkAH8CljzEVnWwQ46RS9ZIzZnbbAlVLzSrJemonKTdeKFSu47777eO211+jq6qK3txcR\nobu7m9tuu42NGzdO+xwLXVdXF62trYyMjFBSUsLSpUun1ZTqcrnYtWsXb775JufOnYuttyyLJUuW\n8P73vx+fzzcToccMDg5y4cKFpHfq8vLy6Onpoa2tbcKm1sm4XC62b99OIBCgq6uL7u7uUdsrKysp\nLCzM2PEZe3p6aG1tTdqZqbCwkObmZjo6OliyZEmao8sc8y7JExE38B3gPqAJeEtE9hhjzsQVexvY\nZowJiMi/Ab4O/KGzbcgYsyWtQSul5qWysrIJH9SORCKISNLmrmtVVVXFihUraGtrY3BwkKysLB57\n7LFpNQVfD4LBIK+++iqXLl3C7XbjdrsZGRkhOzubXbt2Jf2gT4XX6+WOO+5gy5YtdHV1YVkWhYWF\nlJaWzsr4cX6/HxGZsCNOVlYWXV1d00ryANatW8fFixcpKCggOzubUCiEiJCVlcXAwADbtm3L2ObM\naD1ORETw+/2a5E1g3iV5wA7gvDGmAUBEngYeBWJJnjFmb1z5g8BH0xqhUiojRO9mJJsGq729nQ0b\nNsxohwi3282yZctYsmQJLpdLE7xJWJbFK6+8Qnt7+7hkbmhoiBdffJHdu3dTVlY2rfP4fL4Zv2uX\nSCpDsViWdc0zUYyMjHD58mXOnTvH8PAwZWVl1NbWcv/993P8+PHYjB6RSIS8vDzuvvvujG6qTaUe\njTEZN/RNus3HJG8ZcDnu+ybglgnKPw48H/d9jogcxm7K/Zox5tczH6JSKhO4XC7uuecenn/+ea5c\nuUJJSQlZWVkEAgG6u7uprKxk69atcx3mda2trS1ps1xubi7BYJDjx49zzz33zEF0166kpASv10so\nFEqa4EcikWvqeNHf389LL70U62ns8XhoaGjgzJkz1NXVsX37drZs2cLg4CAulwufz5fxyU9ZWRki\nQiQSSdgxJtrTdrrJ/0I3H5O8RPdnEz4wIyIfBbYBd8atXmmMaRGR1cArInLSGPNegn0/DXwaYOXK\nldP6D8/tdqflP8T5TuvBpvUwv+rA5/Px0Y9+lHfeeYdTp07h9/spKirioYceYs2aNbMyfAbYCaaI\nzJt6mCuTXQttbW0UFxeTm5ubcHt2djYtLS1kZWVlzBA0O3fu5LXXXmP58uWxZMvlcpGbm0tbWxu1\ntbUpD4ptWRYvvPACLpeLtWvXxtaXlpYSiUQ4e/Ysy5cvZ926dZSWls7G25kx1/J3wefzsX37do4d\nO8by5cvHNd02NzezefPmUUMjZYp0/n2cj0leE7Ai7vvlQMvYQiJyL/Al4E5jzHB0vTGmxXltEJF9\nwFZgXJJnjPk+8H2ArVu3Gr/fP+WAfT4f09l/odB6sGk9zM86WLduHevWrRvVGSMYDM7IAMWWZdHb\n2xtrKsvPz8eyLFwu17yrh3Sb7Fro7u4mHA4zNDSUtEwoFKKnp4eCgoLZCHHGVVVV0d7ezunTp/F6\nveTm5uJ2u+nq6qKyspLt27enfF00NzfT3NzMsmXLEtaRz+fjwIED4+bmnY+u9e/C+vXr6ezs5Pz5\n8+Tm5pKdnc3w8DDBYJBVq1ZRV1eXkb9f6fz7OB+TvLeAGhGpBpqBx4CPxBcQka3A94AHjTHtceuL\ngYAxZlhEyoDbsTtlKKUUwIw/bN/Q0MCRI0diD4obY6iurmZoaChjJ75Pp6KiIlpaWpIObBsdYzCT\nhp9xuVzceuut1NTUcP78eXp6eigtLWXHjh1UVlZeUzLW1NQ04R3MnJwcuru78fv9M9aBaL7weDzc\neeedbNy4kXPnzuH3+ykvL6empib2zKua2LxL8owxYRH5HPAi9hAqPzLGnBaRrwCHjTF7gG8ABcDP\nnT/Y0aFSNgDfExELewzAr43plauUUjOmvr6eAwcOUFpaGms2siyL5uZmjh07Nu2ZNK4Ha9as4fjx\n47E7n2N1dnZyww03zFqz+mwqLS2NNaFO9e5NOByeNJkRkUlnXMlU0XHytAft1My7JA/AGPMc8NyY\ndX8V9/W9SfZ7HbhxdqNTSil7PLSDBw9SWVk5qqeky+WirKyMSCSyYOYNnU3FxcXU1dVx7NgxysvL\nY50VLMuiq6uL/Pz863qMwfLycs6ePZt01o5wOIzH48mYpmyVXvMyyVNKqfnu8mV7EIBkQ2H4fD46\nOjqSTnWmrtq2bRsFBQUcO3aMzs7OWJN6dXU127dvv66bvVeuXElWVhZDQ0MJO6e0t7dTV1eXkXc6\n1ezTJE8ppaagu7t7wufEoonK0NDQdZnk9fT00N3dTV5eHi6Xi8WLFydtdhQRNmzYQG1tLd3d3ViW\nRUFBwXWd3EV5vV7uuusuXn75ZQKBAMXFxbhcLoLBIJ2dnVRWVlJXVzfXYap5SpM8pZSagpycnAln\n04i63u6w9PX18fzzz3Py5Ek6OzvxeDzk5eVRW1vLI488MuGQF9F5SRcaYwzNzc28+eabtLS04HK5\nWL16NevXr09p2JPly5fz6KOPcurUKRoaGjDGkJuby+23387atWuvu2tMpU6TPKWUmoKVK1fy9ttv\nJ90eDAbJz8+/rsbJ6+/v58knn+Sdd96hsLCQ5cuXk52dTW9vL8eOHaOpqYlPf/rT05qmLNMYYzh6\n9GhsOJXFixdjWRYXLlzg7Nmz7Nq1i9WrVzM4OBibyi1Rs2xJSQnvf//7ue2224hEIni93lmZlk0t\nLJrkKaXUFJSWllJdXU1jYyOVlZWjPnBDoRCBQIDNmzdfVx/E+/fvp76+flRnFJfLRX5+Pnl5ebS2\ntvKb3/yGz3zmM9fN8BfNzc0cOXKEtWvXEgqFgKudc0KhEM8++ywrV66kq6tr1BA8W7dupaioaNzx\nPB7PNU+Jpq5feqUopdQUiAjve9/7ALhw4QIejwe3200oFMLtdrNhw4Z5PwPBTPL7/Rw5ciQ27dZY\nIkJxcTHnz59POo3ZQnTixAmKiooSTs0VCAQ4deoUAwMDbNmyJTYUSktLC01NTTzyyCOUlJTMQdRq\nodAkTymlpsjr9XL33XfT09PD5cuXGR4epqioiJUrV/KrX/1qrsNLq/7+fnp7eycckDc3N5f29na6\nurquiyQvEonQ1tZGZWXluG2WZXHq1CmKi4sZHh6O3fF1uVyUlpbS19fHgQMH2L17d7rDVguIJnlK\nKTVNxcXFSccxu15EkxRjEk41HjPZ9oVqZGSElpYW+vr6EBE8Hg+Dg4P4fL6ETfqLFi2iubmZrq6u\n6+qOsJpZmuQppZSatkWLFlFSUkJ/f3/SoWUCgQCLFi1akD1oE3G73SxbtozGxkaampoIBAJ4vV6M\nMXR0dNDV1cXKlStZv359wv1FhMHBQU3y1JRdH0++KqWUmlX5+fls376dgYGBhEPLWJZFT08PtbW1\nlJeXz0GEc2PlypUcOnQIr9dLSUkJBQUF+Hw+SkpKGBkZ4fLly5SVlSXdP9GzfEqlSpM8pVRGMMYQ\nCAQIBALzoskv1Xgsy2JwcJChoaE0RgfDw8MMDg4SDofTds477riDuro6rly5Qm9vL+FwmHA4jN/v\np7W1lcrKSh599NHrpmct2INmr127loGBAfx+P5FIJFYvAIsXL2ZgYGDcftHpyq6Xu55qdmhzrVJq\nXrMsi4aGBo4fP05/fz/GGIqKiti8eTOrV69O+xAlxhgaGho4ceIEPT09iAg+n4/NmzezZs2aWAIT\niURobW3lZz/7GUNDQxhjWLJkCVu2bGH58uWzFl9nZycnTpygsbEx9uzXDTfcwIYNGxKOvzaTCgoK\n+NjHPsbevXs5cuQI7e3tscGQb731Vh566KHrKmmxLItz585RV1dHOBzm3LlzdHd3xwZDLi8v5/Ll\nyzQ1NY0aJNqyLFpbW9mxY0dsLl+lpkKTPKXUvGWM4eDBg5w+fZrS0tJYL8XBwUFeeeUVOjs72bFj\nR9oSPWMMb731FsePH6ekpCT2wRwIBNi3bx/t7e2xwWrr6+vp7u6moqKCrKwsXC4X/f39PP/889x+\n++1s3LhxxuNramripZdeIjc3l4qKClwuFyMjIxw/fpwLFy7w0EMPzfoUaz6fj927d7Nr1y56enpi\n05pN97myvr4+/H4/brebsrKyjJjlwbIsIpEIbrebgoICNm3aNGp7JBKJJXTt7e14vV6Gh4cJh8Ns\n3ryZG2+8cY4iVwuFJnlKqXmrubmZ06dPs2zZslFNfPn5+eTm5nLy5ElWrFgx4VRZM6mtrY0TJ06M\niycvL4+cnBzOnDnDqlWrGBgYoK2tjVAoxOHDh2OD3LpcLlasWMFrr73G0qVLEw52O1UjIyPs37+f\n4uLiUXfssrKyqKyspL29naNHj8bG9ptthYWFFBYW4vP58Pv9Uz5OX18fb7zxBi0tLbFk3uPxsGXL\nFm644YZ53fTr8XgoLCwkEAgkvIvqdrtZtWoVmzdvpqKigsHBQQoLC6murp5wKBqlUqVJnlJq3jp9\n+jSFhYUJP8ijMynU19enLcmrr68nPz8/aTyFhYWxOVujMxisW7cuVsayLBobG8nLy+PcuXNs3779\nms4/ODjI5cuXGRwcJD8/n+XLl1NQUADYd/GCwWDSO2ZlZWWcO3eOm2++edabbWdKf38/zzzzDCIy\nalaRkZERDh48yNDQEDt27JjjKCdWV1fHq6++mvDnYlkWgUCA++67jyVLlsxBdGqh0yRPKTVvXbly\nZcKehz6fjytXrqQtnra2NgoLCyeMJzpbgWVZ4+atjTZbtrS08O6776ac5BljePvttzl27BgiQlZW\nFiMjIxhj2LJlC1u2bKG7uzvp0CXRc4sIAwMDGZPkHT9+HMuyxj3Hl5WVxbJlyzh58iS1tbUzekd0\npq1ZsyY2hEphYWGsmTkYDNLR0UFdXd119ZyiSi9N8pRS85bb7SYSiSRtkotEImmdx9Pj8Ux4zmis\nLS0t5OTkjNoWDofp6emht7eX3t5ejh49yiOPPJJSgnLq1CmOHj1KZWXlqCE1LMviyJEjuN3uWGwT\niTYZZ4Lh4WHOnz+fdLgVl8uFx+OhoaGBm266Kc3Rpc7j8XD33Xdz4cIFDh06FEvOCwoKuPPOO6mp\nqbmu5jdW6aVJnlJq3qqtraW+vj5pU1Zvby91dXVpjef48eOUl5cTCoXo6+ujv7+fnp4egsEgfX19\nVFRU4Pf76e3tpbCwEGMMQ0NDvPfee1iWhdfrxbIs+vr6+OUvf8nOnTvZsGFD0nOGQiHefvttKioq\nxo2Z5nK5qKys5NixY9x7770TDpcSDAbJzc0dl1RGH/73eDzzKtkYHh6eNCnNycmhv78/jVFNTVZW\nFjfffDPV1dUMDg7GemRnSsKtMpcmeUqpeWvdunWcPn2aoaGhcU2MgUAAsJvD0mXNmjW8/vrrHDt2\njK6uLizLwu12k52djTEGj8fDypUraWhooKuri4GBATweD36/H5/PR0FBAcPDw3i9XtasWcOSJUs4\ncOAARUVFCec3BWhvb4+NmZZI9A5eJBJh6dKltLe3j0uKLcuio6ODO+64I5Yotre3c+rUKRobGwEo\nKiqirq6O1atXz4vkI/p+jTFJk89QKDTrvYVnUlZW1rxuWlYLz9z/JiulVBKFhYXcd999+P1+Wlpa\n6O/vp7+/n5aWFgYHB3nggQcmfEZuJo2MjPDuu+8yPDzMpUuXCIfD5Obm4vV6CQQChEIhtm3bxuLF\ni9m0aROlpaW4XC66u7tpb2+npaWFjo4OACoqKlixYgUejwefz8epU6cmPO9kRIRgMMimTZtwuVw0\nNjbS09PDwMAA7e3ttLa2snXr1lgnkPfee489e/Zw5coVKioqWLp0KcYY9u7dy/79+ydt9k2HvLw8\nVqxYQW9vb8LtxhhCoRDV1dVpjkypzKF38pRS89qyZcv48Ic/zMWLF2lqagJg8+bNrFq1atxzb7Ol\nr6+P3/72t/T19VFbW8uaNWtob2+ns7MTgKqqKhYvXhzr+LBq1Spyc3PJz8/HGBO7GxVNnsrLyykp\nKQHsOV8vX77MyMhIwrHfJrtTZVkWly5divW4jUQiBINBRITly5ezatUq1q5dS3FxMQB+v59XX32V\nJUuWjBpoNz8/n7y8PM6fP8+yZcuora2dfsVN09atW9mzZw/Z2dmj6sEYQ2trK9XV1RN2zFHqeqdJ\nnlJq3svNzWX9+vVJJ3KfTb29vTz77LO43e7YUC1ut5vly5cnnLkiHA4zMDBAUVERnZ2dZGVlUVBQ\ngNvtxhhDOBwelcxFmyKTTY22ePFiCgoKYklcPMuyOHz4MG1tbWzZsiWWZIbDYTo6OhARtmzZMup8\nDQ0NAAlnUhARSktLOX78+LzoEFBWVsZDDz3E/v376enpweVyxeqppqaGW2+9dc5jVGo+0yRPKaWS\nGBwc5IUXXsDj8aQ0OG1zczPvvvsu4XCY3t5e3G43Xq+XcDjMihUrKCgowOv10t7ezvnz56mtrWVw\ncJDi4uKkMzi4XC7uuOMOnn/+eSKRyKjm6cuXL9PY2Mg999wzavgUj8dDZWUlTU1NXLhwYdRdudbW\n1tjYeonk5eXR2toa66gx1yorK/nwhz9Ma2srAwMDuN1uysvL09ZMr1Qm0yRPKaUSMMbEhrxIZRyz\ntrY2Tp06RVFRER6Ph9bWVlwuF1VVVRw5coRQKBRLxEpKSrhw4QKlpaUEAgHuueeeCe9IVVZW8sgj\nj/Dmm2/GZn4wxtDR0cEtt9ySdADkkpKS2FhyUS6XC8uyJnzfE3V2mAvRO6dTFQqFYglvJBKhsrKS\n6urqCZNdpRYCTfKUUiqBixcvcv78+ZSSC8uyePfdd1m0aNG4XrDRpuazZ8/GhjCJDmj86quv8gd/\n8AdUVVVNeo7y8nJ2795NX18fwWCQ7OxsXn755QkHQM7Ly6OlpWVU0lZdXc2rr76a9E6Y3++noqIi\nbc87zraenh5eeuklBgYGKCgoiI1jeOTIEe6444609s5WKt00yVNKqTEsy+LQoUOUlpamdEerv7+f\n4eHhcc/MRZWVlREMBmO9bQEKCgooKyujoqLimoYsWbRoUazpODs7m5GRkYTP1wGxbfHvYeXKleTk\n5MSSnnjhcJi+vj5uu+22lOOZz4aHh3nxxRdxuVwsW7Ystt7n8xEKhdi7dy8FBQVJB1xWKtNpkqeU\nUmNcuXIFv98/KjGYSDgcnjQZzM/Pp6amhsrKytiYesFgkFOnTrF+/fopjU23fv16Dhw4kDS57O7u\nZtOmTQwNDREMBmOdQO6//35efPFF+vv7KSwsxO124/f7GRkZ4dZbb51W0+h8cunSJQKBQMK5jb1e\nL/n5+Zw8eVKTPLVgaZKnlFJjnD179poG2fV6vUl7x0YZY8jOzh41a0VOTg7d3d10dHRMKdFYtWoV\nx48fp6enJzZESpTf7ycUCtHe3s6JEydiMVRUVHDzzTfzwQ9+MNYkHYlEWL9+PTU1NeOOk8nOnz8/\nbv7geIsWLeLSpUtJh69RKtNpkqeUUnGMMTQ3N6fUmzbK5/ORn59PMBhM+CxbdMaKRAmUy+Wiq6tr\nSklednY2DzzwAHv37qW5uZmsrCxEhFAoFOtg0dnZSXl5eWz4Eb/fzzPPPMP9998/Z8PSpEs4HB43\nFVw8l8uFiEzYEUWpjBbtSXU9LzeDMUmWgW99y3R3d5vu7m4z8K1vJS1nIFauu7vbjGzenLRc8OMf\nj5Xre+WVCY/Z98orsbLBj388abmRzZtHnX+iY87WexoZGVlw72kqP6fI448vuPe0EH9Oyd7TAJhv\nL1tmnnrqqdiyB5IuP3v8cfPUU0+Z73znO+bxnTvN/wnmy2Ae3rrVPLx1q/kCmE+D+QcYdcxfVVWZ\nPWB+CuaFab6nzs5Oc+bMGfP6ww+bt8BccmL7cVycv6qqip37Rz/6kfkhmFAG/5xSufbeAPOPY35O\nTz31lPnZ44+bPWCeBvPfwFhpeE8jBw9el79P8e8p+hmxkN7TXP2cgMOp5Dc6rZlSSsUZAKYyeEhx\ncTE7qqoQoBMYCAYZCAYZAW4Gkj3dl+uUnw6Xy0VFRQU3lZWxDVgENAElScrn5OQwDDRP87ypMEA7\n8GZHB7///e85ceIE/YODaTgz1AIhIJxkeydwE1P7eSuVCWSy50iuB1u3bjWvvPLKlPf3+Xz4/f4Z\njCgzaT3YtB4yuw5aW1t54YUXqKysnNL+0SbRAwcO4PF4uPfeeyfsVBEOh+nv7+cjH/nIVEMeJ5X3\ncOXKFbZu3UpdXd2MnXesUCjE4cOHqa+vx+v1kpWVRTAYxLIsbrrpJjZv3jzr4/GdOnWKN954g0WL\nFsWGjQkGg3R2dlJVVcWuXbvGDXszGzL5d2KmaB3YZqIeSkpKjhhjtk1WTp/JU0qpONNNOkSEwsLC\n2Jhsk/WaNcZMqWftRCZ6Di3KsqwJx9ibCa+//jotLS3jeilHIhEOHTpEbm4u69atm9UYNm3aRFFR\nEceOHYsNJJ2bm8vtt99OTU1NWhI8peaKXt1KKRUnOzubdLZwjIyMJB0CZapKS0vxer0MDw8nTOSM\nMViWNeW7lano7e3lvffeY82aNQSDwVHbolOTHT16lLVr16aUlE5HdJ7hoaEhLMsiNzd3xhNrpeYj\nvcqVUiqOz+ebdOqvmRQIBKioqJjRY7rdbm6++Wba29uJRCKjthljaG1tpaamZlbnf41O65bszmh2\ndjZDQ0OxwaHTITc3l/z8fE3w1HVjXl7pIvKgiJwVkfMi8sUE27NF5H86298Ukaq4bX/hrD8rIg+k\nM26lVOaLDnUy9u7TbBkeHmbJkiUzftx169axfft2rly5QltbGz09PbS3t9PS0kJVVRU7d+6c8XPG\nC4VCk96hExHC4WTdIpRS0zXvmmtFxA18B7gPu4PYWyKyxxhzJq7Y40CPMWatiDwG/C3whyKyEXgM\nuAFYCvxWRGqNMaP/lVVKqQnU1tZy8ODBaxoQeSosy8LlcrF48eIZP7aIsGXLFtasWUNjYyO9vb3k\n5uZSXV1NaWnpjJ9vrOLiYkZGRpJujw7xMNNN1Uqpq+ZdkgfsAM4bYxoARORp4FEgPsl7FHjC+fp/\nAd8Wu03gUeBpY8wwcEFEzjvHeyNNsSulFoBVq1bxxhtvxJKw2dLT08OaNWtmNZn0+XzceOONs3b8\nZCorK2NNsol0d3ezYsWKWW0yVup6Nx+ba5cBl+O+b2L8EFOxMsaYMNAHlKa4r1JKTSgvL49169bR\n0dExa+ewLIuhoSE2btw4a+eYS1lZWdx55510dnbS19cX68xiWRYdHR2ICLfccsscR6nUwjYf7+Ql\nekp3bFe3ZGVS2dc+gMingU873w6UlJScTTnC8cqY/nimC4HWg03rQesgquy73/3uhPXwhS98IV2x\nzJWk18LnPve5NIcyp/R3QusgaibqYVUqheZjktcErIj7fjnQkqRMk4h4sAd4705xXwCMMd8Hvj8T\nAYvI4VQGJVzotB5sWg9aB1FaD1oHUVoPWgdR6ayH+dhc+xZQIyLVIuLF7kixZ0yZPcAnnK8/BLxi\n7LaAPcBjTu/baqAGOJSmuJVSSiml5o15dyfPGBMWkc8BLwJu4EfGmNMi8hXsCXn3AD8E/snpWNGN\nnQjilPsZdieNMPBvtWetUkoppa5H8y7JAzDGPAc8N2bdX8V9HQQ+nGTfvwH+ZlYDHG9Gmn0XQF6z\nLQAAEcpJREFUAK0Hm9aD1kGU1oPWQZTWg9ZBVNrqQdI5fY9SSimllEqP+fhMnlJKKaWUmiZN8qZI\nREpE5GUROee8Ficp93UROS0i9SLy/0myiRwz0DXUwUoRecmpgzPx09AtBKnWg1O2UESaReTb6Yxx\ntqVSByKyRUTecH4fTojIH85FrDNtOtMwLiQp1MPnnd//EyLyOxFJaQiITDJZHcSV+5CIGBFZkD1N\nU6kHEfkD53o4LSJPpTvG2ZbC78NKEdkrIm87vxMPz0YcmuRN3ReB3xljaoDfOd+PIiK3AbcDdcAm\nYDtwZzqDnGWT1oHjp8A3jDEbsGcgaU9TfOmSaj0AfBXYn5ao0iuVOggAHzfG3AA8CPy9iBSlMcYZ\nFzcN40PARuCPnOkV48WmYQT+DnsaxgUlxXp4G9hmjKnDnqno6+mNcnalWAeIiA/4U+DN9EaYHqnU\ng4jUAH8B3O78PfgPaQ90FqV4LXwZ+JkxZit259HvzkYsmuRN3aPAT5yvfwL88wRlDJADeIFsIAu4\nkpbo0mPSOnAubI8x5mUAY8yAMSaQvhDTIpVrARG5GSgHXkpTXOk0aR0YY941xpxzvm7BTvZnftLW\n9IpNw2iMCQHRaRjjxdfN/wLuWUh39B2T1oMxZm/c7/5B7HFMF5JUrgWw/9H7OhBMZ3BplEo9/Gvg\nO8aYHgBjzEL7xz+VOjBAdE6/RSQZ03e6NMmbunJjTCuA87pkbAFjzBvAXqDVWV40xtSnNcrZNWkd\nALVAr4j80rkt/Q3nv5yFZNJ6EBEX8E3gP6Y5tnRJ5VqIEZEd2P/8vJeG2GbTdKZhXEiudUrJx4Hn\nZzWi9Ju0DkRkK7DCGPNMOgNLs1SuhVqgVkReE5GDIvJg2qJLj1Tq4AngoyLShD2ayL+bjUDm5RAq\n84WI/BaoSLDpSynuvxbYwNX/WF8WkfcbY34/QyHOuunWAfY1dgewFbgE/E/gk9hjHWaMGaiHzwLP\nGWMuZ+pNnBmog+hxKoF/Aj5hjLFmIrY5NJ1pGBeSa5lS8qPANhbWoyswSR04/+j9Hfbfv4UslWvB\ngz1ZwS7sz8dXRWSTMaZ3lmNLl1Tq4I+AJ40x3xSRndhj/26a6b+JmuRNwBhzb7JtInJFRCqNMa3O\nh1ai283/AjhojBlw9nkeuBXImCRvBuqgCXjbGNPg7PNr7DrIqCRvBuphJ3CHiHwWKAC8IjJgjJno\n+b15ZQbqABEpBJ4FvmyMOThLoabTdKZhXEhSmlJSRO7F/qfgTmPMcJpiS5fJ6sCH/Wz2PucfvQpg\nj4jsNsYcTluUsy/V34mDxpgR4IKInMVO+t5KT4izLpU6eBz72WSMMW+ISA72nLYz2nStzbVTFz+1\n2ieA3yQocwm4U0Q8IpKF/Z/rQmquTaUO3gKKRST67NXd2DOSLCST1oMx5l8ZY1YaY6qAPwd+mkkJ\nXgomrQOxpyn8FfZ7/3kaY5tN05mGcSGZtB6cpsrvAbsX4DNYMEkdGGP6jDFlxpgq5+/AQey6WEgJ\nHqT2O/Fr4C4AESnDbr5tSGuUsyuVOrgE3AMgIhuwn9/vmPFIjDG6TGHBfqbmd8A557XEWb8N+IHz\ntRv7j1o9dmLzrbmOO9114Hx/H3ACOAk8CXjnOva5qIe48p8Evj3Xcae7DoCPAiPAsbhly1zHPgPv\n/WHgXeznC7/krPsK9gc42H+8fw6cx55Le/VcxzxH9fBb7I5n0Z/9nrmOOd11MKbsPuzexnMe9xxc\nCwJ8y/lcPAk8Ntcxz0EdbAReA447vw/3z0YcOuOFUkoppdQCpM21SimllFILkCZ5SimllFILkCZ5\nSimllFILkCZ5SimllFILkCZ5SimllFILkCZ5Sl0HRGSfiBgReWKuY0k3EXnCee/75uj8N4rIz0Sk\nVUTCTizHxpS5Q0SeFZEOEYk4ZX49F/HONRHZ5bx/HfpBqWnSGS+UyhDOpPYfAj4C3IQ9P2wEe+yx\nVuwx2F4FfmeM6Z+rONVVIlKNPRaWz1nVjT1WYGdcmVuBV7D/HhugC/vn2pPWYOOISBXO9FvGmCfm\nKg6l1PRokqdUBhCRIuxR4uPn+wwDAWAlsBq4Hfgz4P/AHnQ63iXgLHHJhUqLP8FO8M4DdxljmhKU\n+Q/Yf4tfwx4odT5MeVYF/Cfn6yfSfO4A9rWqlJomTfKUygw/xU7wIsDfY8+k8p4xxnLmQ92IPQ/i\nRxLtbIz5eLoCVaPc6Lz+JkmCF1/m6XmS4M0pY8whYP1cx6HUQqBJnlLznIjUAB9wvv2yMeZr8duN\nMWHsaeNOAF8Xkdw0h6iSy3NeB6ZZRimlrpl2vFBq/tsS9/VvJitsjBkau26yjhcikiUifyYix0Rk\nUES6nX0+NNn+0YfknQfmfSLy1yLyjogMiUiXiDwjIrdMFrez//8QkUsiEhSRPhE5JCJfEJH8SfZ9\nSEReFpFeERkQkePOflmTnTcVIlIhIl9zjtvnxNcgIj8QkY0Jyjc6HQd2Oav+U1w9Resq2rmgyinz\n4zFlqsYcM0dE/lRE9otIp4iERKRNRH4tIg+m8B5uEZEfi8h552fcLyJnRORHInJ/fOzA3rjvzZjl\nyWusu1tE5L+LyAWn3gZF5KLzPv5SRJaPKZ+w44WIVCWIJdnyRJJYpnyNKZWJ9E6eUpllOVA/kwd0\nPtyeA97vrIoAw873d4rI15LtO0YlcBRYCwQBCygBHgHuF5EPGGNeTHB+D/APwB/HrR4A8oHtzvIp\nEXnAGHMxwf5PcPX5MYBe7Obrv3XO/VqK8SckIv8M+B9AgbNqBAgB1cDjwMdE5F8bY34at1sHkIP9\n/rOAQUbfqQthd5gBWIz9D3c/EJ+gR+JiqAGeBWqcVcYpXw48CjwqIv9gjPlsgvjd2JPB/2nc6kHA\nDWxwlg8CRXGxFwLFzvdXGK1v7DmSEZFPAD/GnpAe7OsqjP0c6Ursa+wy458hTSTaySjp6bA7IyWK\nY1rXmFIZyxijiy66zOMF+06Phf3BfgKoncIx9jn7P5Fg2391tkWALwAFzvoy4L8423om2N84Szdw\nGrgLO2kR7A/Pd5ztjYArwf5/72xvA/4NUOKsz8K+E3bU2X5k7P7A7rjz/wxY4azPBT6LnVREY983\nhXrb4RzDOPW0HnA721YC33G2jQDbrqXe48o0OmU+mWR7EXDBKfM74A4g29m2CLuzjd/Z/u8T7P+3\ncXX0w/jrBzspehT7ecD4fXZF95nGdZuHnYga4J+ANXHb8oGbga8DD8/EuYH/29lvGHjfTF1juuiS\nycucB6CLLrpMvgDfj/ugtpwPpe8AnwI2ATLJ/gmTDSdRiTjbvpxk3yfjzv1Egu3Rbe3AkgTbb4wr\nc/uYbZuc9zMI3Jjk/D7suz0G+Odjtp121u9L9OGM3bs1eu59U6j3Q86+X5mgTDQR/nWq9T6mTCMT\nJ3nf4GqC50lS5l84ZTriywC1cT/fv72G9z2lRGvMMXY4xxhIFvdMnRu7R3n05/zJMdumdY3poksm\nL/pMnlKZ4bPAV7E/qATY6qz7IXASaBORb4lI+TUe919i33ULAH+XpMxXUzzW940x7WNXGmNOYt+J\nAqgbs/lx7PfzrFNuHGOMH3v4GIAHoutFpA67WRbgr40xVoLd/xFoTjH+UURkM/adyBHgmxMUjTbT\n3us0jc4YERHsRB7gm8buZJPIr7HvmpVh3yGL+gT2z7eL0U3a6dDrvHqB0tk6iYjchd3bHOD/McY8\nOabIlK8xpTKdPpOnVAZwPtz/SkS+id3T9k7sBGQD9ofoEuxmu4+JyCPGHoYiFTc5r4eNMYNJzv2e\niFwGVkxyrDcn2NaC/QxbyZj173NeHxKRtgn2jz4Ptypu3TbnNYw9CPQ4xh5iZh/wryY4djLR2FzA\nWTvfSiia2OVjJzPjEt1p2MjVOntSRBIlslHxdRT9WdzmvL5sjAnOYFypeA+7qX498KaI/APwInDS\nGBOZcM8Uich64BfYza6/AL6UoNh0rjGlMpomeUplEGNMH/DfnAURycH+EPtT7OSvDPiFiNSk+KG+\n2HltmaRcM5Mnef4JtkXvQI3t7brUeS3g6ofsRPLivo4+ZN9pjBmeYJ9k49NNJhqbG7uDQyryJi8y\npRjg6s/qWmKocF7T3pnAGBMRkceAX2En+F9zloCIvA78EviJMSYwleOLSBl2Z5Ri4DDwMWNMoqnQ\npnONKZXRtLlWqQxmjAkaY35rjNkN/MRZvRx7YORURG9PTTZPaNLbWNMUvQv2RWOMpLDsSnCM2Zrj\nNBrbOynGJsaYxlmKAaAixRieTHCcOZkH1hhzHPtO3r/Efq70FHanmHuB7wLviMiNyY+QmIhkYzev\nrsZ+lm63STB0kGMmrjGlMpImeUotHN+P+3pdivtEmxaXTlhq8u1TFW0+u+YPeq7Gvtj50E9m2RSO\nDVdjWz2HY6jFNy9OpY5andeq6YcyNcaYkDHml8aYPzHG3Ih9R/Iz2L2xV3D1n5Nr8WPsafwGgA8Y\nY1onKDuda0ypjKZJnlILR/w4bBM1X8Y76rxuS5bIiMhqJm+qnaroGHaPiEgqTWnxDjuvHq4+dzWK\niLi4OiDxtYrG5sXuvToXTmF3qAB4bAr7v+683uc07acq9uyfTPAw4lQYY7qMMd8D/i9n1VYRSblj\nhoh8BfgjJ8aPOHcLJzKda0ypjKZJnlLznIhUi0htCkU/Eff10aSlRvsl9odlPvDvk5RJ9DD7TPlH\n7KbEIuyhQpISe1aO2Ie0MeYEVweG/pKT0I31Kezm66k4DLztfP03IjLhM3EiMrZTybQ5HW5+5Hz7\nCRFJmMxOEMOT2EOolAL/+RpO3R/3dVHSUhPHMtHdVUgy8PMkx/wY8JfOt39ujPnfKew25WtMqUyn\nSZ5S898NQL2IPCsiH4+f7sr5UNoqIj8GPu+sPgQcSOXAxh7d/4fOt18RkT+PfsiJSKmIfAs7UepN\ndozpMMYcwx6oFuAzIvJzEdkSvXskIm4R2Swif4ndW3PLmENEE9C7gKeiU2SJPQXYZ4BvTzV25yH+\nz2DfFV2J3UP0QyISezBfRJaJyEdF5GXsQYdnw1ex37sHeEFEPh+fcIrIIhF5UER+wphexsaY81xN\nbL4g9jRsNXH7LhaRPxSRX40557vYs3IA/PEU7+Y9JiKvicifOHeDo+d0i8gD2J0wAN4wxkz6MxKR\n24EfON9+zxiTbMifUWbgGlMqc831QH266KLLxAv2uF1mzDKMPfaZNWb9EWBpgmPsI/lgxgXYyUH0\nGGHs56Wix/4qsN/5+osJ9o/ut2uC9zDR+d3YY/TFv48hoBN7jLr49bcn2P+vx5Tpjtvv91ydCWHf\nFOv/PieW+PrpxB6zMP68/3gt7zuuTCMTDIbslKkGjo05Xw/2FGPx684lqd9vjynnHxN/b4L9fhC3\nfRC7h24j8P+mWG+fHHPOoFNvkbh1zcD6Mfvtim4fs/6JuP06sJ+1S7b8+UxeY7rokqmL3slTap4z\n9nyvNdjNqT/HbqIcxm5+CgDnsKf0egzYboyZbDiUsccfAO4B/iP2tGkh7N60+4EPGmP+kqtNdjN+\nR88YEzHG/Bn2mH3fB85iJwKLsBOZ17A/4LcYY8bNQ2uM+TLwz4BXsJsZs7Hr6IvO+wqN3eca43sZ\nez7ev8C+Q9qHXR8WcAb7Tuhu4N9N5zyTxHABe1zAjwPPYHeoyMd+XvAC9jAlnwJ2Jtg3Yoz5HPZz\ni/8duIQ9lE0Ie8aQH2L3fh3r32LX+ynn+5XYY8iVpRj2HifeHwPHsettEXaCeQi72fUGY8w7KR4v\nXhn2sDbJllFNrtO9xpTKVGLMnPSsV0plCKf5tgs7oXi/MSbhwMNKKaXmF72Tp5SazOexE7xu4K05\njkUppVSKNMlT6jonIj4Redp5eL8obv0qEfkGdjMWwN+b9E+NpZRSaoq0uVap65yT2PXErYpOT+aL\nW/cL4DFjD+mhlFIqA2iSp9R1TkQ8wJ9g9yLdhD0jQS52z8PDwE+BXxj9Y6GUUhlFkzyllFJKqQVI\nn8lTSimllFqANMlTSimllFqANMlTSimllFqANMlTSimllFqANMlTSimllFqANMlTSimllFqA/n/N\n8i/ErOHeBAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "results.volcano()" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "### DRAFT version to use new JS volcano plot (v 0.12)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from gdsctools.report import ReportMain\n", "from gdsctools.volcano import VolcanoANOVAJS\n", "v = VolcanoANOVAJS(results)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": true }, "outputs": [], "source": [ "html = v.render_drug(999)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Created directory report\n" ] } ], "source": [ "report = ReportMain(template_filename=\"volcano_standalone.html\")" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": true }, "outputs": [], "source": [ "report.jinja['volcano_jsdata'] = html" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": true }, "outputs": [], "source": [ "report.write()" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from IPython.display import HTML" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "display(HTML(''))" ] } ], "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.5.3" } }, "nbformat": 4, "nbformat_minor": 1 }