{ "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": { "collapsed": false }, "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": false }, "outputs": [], "source": [ "from gdsctools import ANOVA, ic50_test" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "WARNING: column named 'MEDIA_FACTOR' not found\n", "TISSUE FACTOR : included\n", "MEDIA FACTOR : NOT included\n", "MSI FACTOR : included\n", "FEATURE FACTOR : included\n" ] } ], "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": false }, "outputs": [], "source": [ "drug_name = an.drugIds[0]" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " [-----------------100%-----------------] 47 of 47 complete in 0.3 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": 14, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmkAAAGbCAYAAACfwwddAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3el30+eZP/73JW+ysGVZ3hfAG8sAIdiYNiGBJhhD2s5k\npiEkf8A3QOd5E5Izz78JNM87QOecOWcefNNAfqftmWkKNiTjtE0DtiFhDWDZBG9gW5LlTbYl3b8H\nn4+MF0m2bFmL9X6doyP7s+i+LBH7yr1ctyilQERERETxxRDrAIiIiIhoISZpRERERHGISRoRERFR\nHGKSRkRERBSHmKQRERERxSEmaURERERxiEkaERERURxikkaURETkjIj4RKQi1LEwX/M9/f4DkYky\ncaz0/RSRN0WkQ7/+jdWLlIgSEZM0IqLY+RRABYAzANpXuzERqdMTwndXuy0iWrnUWAdARJSMRKRK\n//KsUupfo9w8t5ohSgDsSSMiiq3hGLQpMWiTiMLEJI0owYmIRZ8H5Z/bZBeRVhE5sgrtnNfbsYvI\npWBt+Odl6V8fFxGHiHyqf39ej9Mc4L4m/32L/IytItIwaz5cxbzr3xSRNr1dn4g8XOoQXzTeTxE5\nD+Ch/q3/Z3hj1vnjevz+9kO910uKV0SaALTq356aPY8w3M8k1Oe7jPiX/VkRrXVM0ogSmD5kZgdw\nDNof/TMAmgHUATgvIsci1E4dgE4ARwD4ADQBqAJwHsDxEPcdB/Dv+j0ds06FGm6bc05ELADaoP2M\ndmg/o0WP4WCANt/Ds7lel/TrrdASk49CtBu19xPA/wVwUv+6Cdp72K7H0IRn75m//YN6+3PiDzPe\njwCcntdm66zzS/5MZrW/4PMNM/5lf1ZESUEpxQcffCToA8ApAF4A78w7Xgntj+TFecfP6McrQh0L\n0E6bfs2v5h3/SD/uA3AgwGv6ALwx757z+nFzgHaaAHiDXD//Z/QnAt55P48DwNDs1weQAy2Z8Qb7\nGaP5furXVenXfTjr2HH92P+Zd20OtITKB6B2BfHWBfkcw/1MAn6+y4h/2Z8VH3wkw4M9aUSJ7RMA\nR5VSv519UCnVqX9pXWkDem9NLYA2pdTH89p5H4AtxO1nlFL/3wratkDrvWsL8DP+MshtOfp516xr\nh6ElKLsXaXLV389FnITWI3VBH8a06O+BADirX/N2HMU7//MNN/6VfFZEax5XdxIlMKXUdQDX/d/r\nw5L1AI5GsJk6/fl3Qc43I/iQZ9MK2/avgGwOcv46gF3zjp0FcFxE7PrXzQBalVJdizUWpfczlEr9\n2RHkvIKe2ABxEe/8zzes+LGCz4ooGTBJI0pwInIG2pwf/x/INv3REKEm/L0xziDnQ/WkhTq3FP4k\nbSjE69fOPqCU+qWItAE4AeA9/QERaQZwUk9sgorC+xmsXf/P2oZn89UCmfOexire+bEsJ/6VflZE\nax2TNKIEpk/SboA2R+iMUurGrHNBJ/SHyf9HNTfI+eoQ9wZL7IKpmve9//68INdbAh1USp0DcE5f\nrbgHWs/ScQBtIlI9azhwjii9nwEppWwi4v/6ylLuiVK88z+T2WY+3+XEr1+7rM+KKBlwThpRgtLn\n+jRAm6/1r7P/QEeYf1Xm20HOL1hhuRz6z1M57/A1/bkxyG0NmLXyUEQqReSUv9yDUsqllLqsz1/z\n9+4ELWWB6LyfoQwDqBORnPkn9J/LJyK79O9XPd4gn0ko4cS/7M+KKFkwSSNKUEopfy/Ggt4kETkV\nwXY6oZd1mF+/Sm+nEuFVsPcPXc5P+hbErE8i97c9p/yFPsw3nwPAu9B6ZuYnCjX6c8Ah2Gi9n4v4\nv9Am2V+e134dtJ/L4U/GVhjv/J7JJX8mi1hy/FjBZ0WULDjcSZTYLgB4U0QuQUtm8gC8CS1psgHY\nLSIf6aswV+IEtLlGp0TkBLTJ6nXQErQLeptL5a+tdkZEGqHVXzsIrVaWv6ZWoLbP6G134lkpCf/1\ndkBLXETEH0+niFzWz/nnbLUtsto0Wu9nQEqpX4vI29CSUju0ZMffY6awcJ7ZcuM9KtrY5Bk9CQ/3\nM1lx/BH4rIjWvljXAOGDDz5W9oBWq+whtKTlGvQaWNCGiuyYVYcKWm2x+XXFFhwL0k4OtMKjD/S2\nHgB4A9ofXi/m1kkL+Zp6bP66WV5oK0cr/PeFaNuuX38RgBkBanjp17876z3x6rH+CgHqgMXw/VxQ\nJ21e/P73Z8j//qw0Xv34p7Pel13L+UyW8PmGE/+yPys++FjrD1GK++wSUfwSkYMAhlSAlX4i0gHA\nopQKtrCAiChhcU4aEcW7MwAuz5+3JCJvQhsWC1ZDjYgoobEnjYjimr767zy0cg9noQ3h7YE2POcA\nsFux+CkRrUFM0ogo7olIA7SyDPXQJqK3QZvzdFLN2lKIiGgtWXNJmoisrR8oyR0+fBgXL16MdRhE\n/LdIcYX/HtcWpZQEOr4mS3DY7fZYh0AR8vHHH+P//b//F+swKMl98sknaG9v5+8Wihv83bh2WK3W\noOe4cICIiIgoDjFJIyIiIopDTNIoru3fvz/WIRABAHbv3h3rEIhm8HdjcmCSRnGNv4goXjBJo3jC\n343JYU0uHCAiIiJNqInpFF3hLj5ikkZERLTGrbVyW4lIJGCVjZA43ElEREQUh5ikEREREcUhJmlE\nREREcYhJGhEREcWEwWBY8MjLywt4rr6+HtevX59zf25u7oLr3nrrLQwPDwdsr7m5Gbt374bBYMCh\nQ4fQ2dk5c85/3P94++2359zrdDrR2NgIq9WK999/f87xQ4cOReotmYNJGhEREcVMe3s7nE7nzGNo\naGjBOZvNhrfffhu7d++ek1iJCJqbm2fubWtrg81mw4cffrigHX8y9W//9m9wOp2oq6vD0aNHZ853\ndnaivb0dNpsNNpsN586dm3P/hx9+iMOHD8Nut6O9vR2XL1+eOT47aYskJmlEREQUMxaLBWazeeYR\n6FxFRQXeffddHD9+HKdOnZpzjdVqnbm3trYWb7/9Nmw224J2/L1ob7zxBsxmM95//320t7fPnHc6\nndi1axcqKipQUVGxIBYRmbNKVkTgdDpx/fp1HDhwIBJvxQJM0oiIiJKdSODHUq9fgXDKgxw9ehSt\nra1B73c6nWhqago4/Pjmm2+iubl55vvW1lZUV1cDAGw2GywWCw4dOoSampqAQ6YffPABmpqaYLVa\nUV1djQMHDuDkyZM4ffr0kuMPF5M0IiIiipnq6uo5c8GuXLkS9Nrc3Nw5vWRKqTlzyfyFe995552A\n9+fk5AAALly4gLfeegtnzpwBAAwPD6O6uhrvv/8+mpqaAGDOUKj/3kuXLsFut+M3v/kNbDYbnE4n\nKisr0djYiJqaGnz22WfLfyMCYDFbIiKiZBdusdsIFsdtbm5GVVXVzPcVFRVBr7Xb7XOuFRFcuHAB\ndXV1AACHw4GGhgZcvnwZDQ0NC+53Op04evQonE4nrly5gl27dgEAamtrce3atZnrzp07h9zc3JBx\nnz59Gh999BHee+89HD58GEeOHEFjYyOOHDmypJ97KdiTRkRERDFTVVU1Mw8sVIIGAOfPn8eePXuC\n3l9bW4u33nprzlyz2RoaGlBfX49r167NJGgAcP369Tn3+HvcgrHZbHA4HKisrERbWxuOHz+OyspK\nVFVVoaurK/QPHAb2pBEREVFccjgcsFgssNvtuHDhAn7729+io6Mj5D3V1dUBr7lw4QJEBMePH58z\nZFpVVQWlFA4ePIjLly+jsrISJ0+eXDDcOdv7778/MxetqqoKv/vd73D06FHYbLZFE81wsCeNiIiI\nYmKx/Sx3794Nq9WKmpoanD9/Hm1tbYsmQVVVVTh//vyC462trWhvb0d1dTVqampQU1ODTZs2AQDq\n6upw6tQpHD16FFarFU6nc0EJDr/r168jLy9vJo5z587hzJkzqKqqivgiAllrm66KiAp3l3mKX0aj\nEW63O9ZhUJL75JNPUF9fj5qamliHQgQgvN+NVquVG6zHARFBoPxE/3wCZqvsSSMiIiKKQ0zSiIiI\niOIQkzQiIiKiOMQkjYiIiCgOMUkjIiIiikNM0oiIiIjiEIvZEhERrXGL1SOj+MQkjYiIaA1j7dDE\nxeFOIiIiojjEJI2IiIgoDjFJIyIiIopDTNKIiIiI4hCTNCIiIqI4xCSNiIiIKA4xSSMiIiKKQ0zS\niIiIiOIQkzQiIiKiOMQkjYiIiCgOMUkjIiIiikNM0oiIiIjiEJM0IiIiojjEJI2IiIgoDjFJIyIi\nIopDTNKIiIiI4hCTNCIiIqI4xCSNiIiIKA4xSSMiIiKKQ0zSiIiIiOIQkzQiIiKiOMQkjYiIiCgO\nMUkjIiIiikNM0oiIiIjiEJM0IiIiojjEJI2IiIgoDjFJIyIiIopDTNKIiIiI4hCTNCIiIqI4xCSN\niIiIKA4xSSMiIiKKQ0zSiIiIiOIQkzQiIiKiOMQkjYiIiCgOMUkjIiIiikNM0oiIiIjiEJM0IiIi\nojjEJI2IiIgoDsVtkiYipxY7JyLHohcRERERUfTEZZImIscBHAlxyTEReQCgI0ohEREREUVVXCZp\nSqmzAGwhLjmmlNqklLoSrZiIiIiIoikuk7QlsIpIg4i8G+tAiIiIiFZDQiZpSqlzSqnLAPJEpCHW\n8RARERFFWsIlaSJyTET889WGAFTFMh4iIiKi1ZAwSZqIWPQvbQCa9a/zAFyLTUREREREqyc11gEE\nIiJvAqgXkXeUUr/VDzcDqFdKXRaRIyICAINKqRuz7z18+DA+/vjjme/379+P/fv3Ryt0irDU1FQY\njcZYh0FJrr6+HqWlpfy3SHGDvxsTV0tLC1paWpZ0rSilVjmc6BIRZbfbYx0GRYjRaITb7Y51GJTk\nPvnkE9TX16OmpibWoRAB4O/GtcRqtUIpJYHOJcxwJxEREVEyicvhTiIioljzer3o6+tDZ2cnPB4P\nCgoKUFFRgaysrFiHRkmCSRoREdE8o6OjaG5uhsPhgNFoREpKCh49eoTW1lbs3bsXmzdvjnWIlASY\npBEREc3i9XrR1NSEiYkJlJSUzBw3m82Ynp5GS0sL1q1bh7KyshhGScmAc9KIiIhm6evrg8PhgNVq\nXXAuLS0NZrMZ3377bQwio2TDJI2IiGgWm82GzMzMoOfNZjOePHmC8fHxKEZFyYhJGhER0SzT09NI\nSUkJeY2IwOv1RikiSlZM0oiIiGYpLCzExMRE0PPT09NIS0sL2dtGFAkrWjggImZoe2daAVgAOAHY\nAdiUUq6Vh0dERBRdlZWVaG1thcfjQWrqwj+TAwMDqKurC3iOKJLC+hcmIjkA3gLQCKABWmIWqEqu\nEhEntK2cmgB8yqSNiIgSQVZWFvbu3YuWlhZYLBZkZ2cD0HrQBgYGUFhYiG3btsU4SkoGS0rSRKQB\nwAkAb8471YlnvWdOaEmbv1etEsBR/XFGRC4AOKOUuhKZ0ImIiFbHli1bkJWVhevXr6Ovrw+AtrKz\nrq4O27ZtQ3p6eowjpGQQMkkTkVoA5wDU6Yf8PWPNSqnri724iNQBOAit5+0ogKMi0gbg2PyN0YmI\niOJJWVkZysrKMDY2Bp/Ph8zMTA5xUlQFXTggIv8OoA1ar9gJALlKqUNKqV8vJUEDAKVUu1LqtFKq\nEVoP2y/153YR+c3KwyciIlpd69atQ3Z2NhM0irpQqzv3ADiqlKpRSp1TSg2vpCGllFMpdVYpVQ3g\nkP76RERERBRA0P8tUErtXq1GlVLNAOpX6/WJiIiIEl1E+2711Z/1AHIAtCuluiL5+kRERETJIuxi\ntiKSIyIficg1EamYdbwW2mrPSwAuAOgQkd9FKlAiIiKiZLKcOmmd0BYTAHNrpJ2fdfw6gFpoqznt\nSql/XWmgRERERMkk3J60D/BsZ4GTSqlOYKaOWpV+TbU+n80/5+y4vjMBERERES1RuEnaUf+zUurX\nAY5f8CduSql2AJeh9bZxkQARERFRGMJN0ioBKKXU5XnHG/Xn+XPQ2vXnOhARERHRkoWbpA0DMxur\nQ/+6EnryBm1Hgtly9GfncgMkIiIiSkbhJmk2aMOXB2cdO6E/Xw9Q8PYgtOStc3nhERERESWncJO0\nM/rzeRH5lYi8C+C9eecgIpUi8im0xQQCoHXFkRIRERElkbBKcCilzorIUQANAE7POmVTSp0DABGp\nAvBw1rnTK91SioiIiCjZhL3jgFKqUUTeg7ZYwArgGoCTAS61ATjlT96IiIiIaOmWtS2UUuo05vak\nzT5nwzJ2MiAiIiKiZ5hMEREREcWhZW+wru/VeRBaoVoLtHlp/6pvHVWplLoRoRiJiIiIks6ykjR9\n5eab8w636c97AFwSEQeABiZrREREROELe7hTRFrxLEFrBvDreZfYoBW9zQXQJiIHVhQhERERURIK\nK0kTkWN4tsVTo1LqkFJqzspOpZRNKZUL4DNoNdLOgIiIiIjCEm5Pmn93gZMB9u+cQyl1FNpOA9Xs\nTSMiIiIKT7hJWjW0bZ4uLPF6/16eVWG2Q0RERJTUwk3S/BumD4V5nyXM64mIiIiSWrhJ2nUs3GA9\nlAb9mRusExEREYUh3CTtd/rzORGpCHWhiHwEfZhTKfVZ2JERERERJbFwN1g/LSJvA6gF0CEiZ/Fs\n3hlEZBe0Omkn8GwV6IkFL0REREREIS2nmG0DgPP68wk8S8J2QytoK7OuPcEN1omIiIjCF3aSppRy\nAmgUkYMAjkLbFqoK2qKCLmjFbNsAfKiUGo5cqERERETJY9l7dyqlmjFrqJOIiIiIIifsbaGIiIiI\naPUxSSMiIiKKQ2ENd4qID9qOA0u+BYBSSqWEFRURERFRklvOnDRZ/JIZzmW8PhEREVHSCzdJsy7h\nfBW0shxHAFxTSh1eTmBEREREySzcYraL9Yw5oZXgaNZLdFwSkT8rpV5bboBEREREyWjVFg7oJTo+\nA3BI34mAiIiIiJZotVd3XtOf61e5HSIiIqI1ZbWTNP8cNssqt0NERES0pqxakiYitQDeg1ay4/pq\ntUNERES0FoVbJ+0SllYnzQptw3UAGFZKXQ43MCIiIqJkFm4JjoNhXu8E0BDmPURERERJL9wk7ZdL\nvM4OoEMpxWFOIiIiomUIt07a2dUKhIiIiIie4QbrRERERHGISRoRERFRHAo63CkirVjaSs5FKaX2\nROJ1iIiIiJJFqDlpdVGLgoiIiIjmCJWkvRWhNiLSG0dERESUTIImaUqpC9EMhIiIiIieWdWFAyKS\nIyJviEjOarZDREREtNas9urO4wAuADi6yu0QERERrSnh7jgAvVfsAwC10PboDMW/f2duuO0QERER\nJbNwN1jPAdAJwBLGbU5ovWlEREREtETh9qR9AC1BcwL4EFrCdgLaJuqnAbTq1+0B8C6ANtZIIyIi\nIgpfuEmaf27ZMaXUZwAgIjZoyVnOrBWhF0SkCcAlEXlHKfXbyIRLRERElBzCXThghVb3rNl/QCnV\nDmAYQP3sC5VSzQAuAzgjIuYVxklERESUVMJN0nIAQCk1PO+4DYF3KGgCIAAOhh8aERERUfIKN0m7\nDkBEZNe846368Yp5x236c1X4oRERERElr3CTNP/CgFPzjnfozyfmHff3oM3veSMiIiKiEMJdOHAG\nwDEAjSIyBOCoUuoKtBIbHwF4Tz/eDq1G2nH9vuZAL0ZEREQUKxMTE7DZbLh9+zYmJiZgNpuxfft2\nVFRUID09PdbhhdeTpi8SeAtaz1gu9GFMpZQNwDlo889OQ5uL5u9tO6uU6oxUwEREREQrNTIygv/+\n7//G1atXkZ6ejoKCAiil8Je//AV//vOf4Xa7Yx1i+NtCKaUuKKVyARyCtnrTf/wEgPeh1U5zQetN\nO6GU+mWEYiUiIiKKiJaWFkxPT6OkpARGoxEGgwEmkwmlpaVwOp345ptvYh1ieEmaiNhF5EMR2aWU\nap7fQ6aUOq2UqlZK5Sql6pVS5yIbLhEREdHKDA4O4smTJ7BaA+9umZ+fD5vNhtHR0ShHNle4PWkW\nACcBtInIAxH5VYAVnURERERxy+FwQESCnjcYDBAROJ3OKEYVII4wr/8ltCFOAVANbf5Zh4hcE5F3\nWLSWiIiI4p2IQCkV8prFzkdDuAsHziqlGqHtPDA7YdsN4CwAh4hcFJE3Ih4pERERUQQEG+b08/l8\nEJFFr1ttYS8cAACllNOfsCmlDNBWfH4GLWFrhLZ3p1dEficiB5bThojMr8U2+9wREWkQkWPLeW0i\nIiJKXlarFaWlpRgcHAx4/unTp6ipqYHJZIpyZHMtK0mbT1/xeTRAwnYUWjmOsIjIcQBHgpyr09u8\nrH9fu9y4iYiIKDnt27cPJpMJvb29GB8fh8fjwejoKLq7u1FYWIgf/ehHsQ4xMkman54w1QOYnTgF\nn5kXhFLqLJ5tKTXfWwAc+tc2cF9QIiIiCtO6devwj//4j9i3bx9SUlLgcrmQmZmJgwcP4tChQ3FR\nzDbcHQcWEJEGaD1mR6Gt/pydlF0A8LuVtjGPBYB91vd5EX59IiIiSgLp6enYtGkTNm3aFOtQAgo7\nSRORHGi9V2/j2ZDkgsRMKfXZysMLHsYqvjYRERFRzIWVpInIJQANmJskDUObd3bGP09slTmhrS4F\ntK2phqLQJhEREVFUhduT5p//NQxtGPN8lBIziIhFKeXU262HVv6jEstYmEBEREQU78JN0s4B+HS1\nEzMReRNAvYi8o5T6rX64GUC9Uuq6iNTrc+GcSqkbs+89fPgwPv7445nv9+/fj/37969muLSKUlNT\nYTQaYx0GJbn6+nqUlpby3yLFDf5uTFwtLS1oaWlZ0rUSDxV1I0lElN1uX/xCSghGoxFutzvWYVCS\n++STT1BfX4+amppYh0IEgL8b1xKr1QqlVMC59ssuwSEitSLyqb4l1JBevPaBvuPAh9wiioiIiGj5\nwk7SRKRSRFoBtAJ4E9qWULl4tp9nI7RN2B0i8qsIxkpERESUNMJd3ZkDoA1arTJAmyd2HlrCNgyg\nCkAdgA8A5AA4pW9i+nGAlyMiIiKiIMLtSfsAzxK0RqXUIaXUOaXUdaWUTSnVrJQ6rZTKxbOtoU5x\n6JOIiIgoPOEmaW/qzycWW+GplDoKrXdNoG3lRERERERLFG6SVgVAQRvmXIpP9WdLyKuIiIiIaI5w\nk7Rh/XmpdTty591HREREREsQbpL2KbThyxOLXagvMmjUv11qzxsRERERIfwk7SS0XrH3QpXXEJFK\naKtAcwCcVkp1Lj9EIiIiouQT7rZQdQDegbY91GkR+QBaL5kN2kbnNdD21azTr3cCGAyW0LE0BxER\nEVFg4SZpzdDmo/m3L8gFcDTE9RYAp4OcUwCYpBEREREFEG6S9lkE215bm4YSERERRVBYSZpe+4yI\niIiIVlnQhQOrvUsAdyEgIiIiCi7U6k6niHwY6WRKRHJE5CMAjki+LhEREdFaEipJ+yW0khsOEfmN\niBxYSUMi0iAi/w4tOXtPf30iIiIiCiDonDSl1FkR+RTAKWjFa0+IiAPaCs8mAK1KqRvB7heRWmjl\nOBoBNODZ7gNnAZxUSnEXAiIiIqIgQi4cUEo5oSVnJwEcB/ABtJIbRwFARPwrNJ0A7ACseLZPp8x6\nKSe0XrmzTM6IiIiIFrek1Z16snYaWgHbOgAHofWQ1UPbVSAXz3rKAG1XgmvQetyalVLXIxk0ERER\n0VoXbp00KKXaAbRjVpFaEbFA60Wz6wkdEREREa1A2ElaIHpixuSMiIiIKELC3WCdiIiIiKIgrJ40\nEcmBNg+tFkAeni0ScALogLZ4oF0p1RXBGIkoibhcLjx+/Bjj4+Mwm81Yv349TCZTrMMiIoq6JSVp\nejmNU9BKaQBzV27Op0TECeBDpRQ3UCeiJfF6vfjmm29w7949GAwGpKamYnp6GgDw4x//GNu2bYtx\nhERE0bVokiYi70JL0GbrBGDDs3lo/oUDVXi22vO0iJwAsFsp5YpYxES0JrW2tuLOnTsoKyuDyLP/\nD/R4PPjb3/6GtLQ0bNq0KYYREhFFV8gkTUQa8CxBa4dWhPbyIvdUQSvRcQpANYDLAPasPFQiWqtG\nR0dx584dlJaWzknQACA1NRUFBQVoa2tDVVUVUlJSYhQlEVF0LbZw4KT+fEEpVb9YggYASimbUuos\ngEpo9dLqRGTXCuMkojWsr68PAGAwBP6VlJGRgYmJCQwNDUUzLCKimFosSTsIQAE4Fu4L62U5zkKb\nv1YffmhElCwmJyeDJmh+IgKPxxOliIiIYm9JJThWsJXToP5sCXkVESU1s9m8aALm8/m4ypOIkspi\nSVonABGRd5b5+m/Peh0iooBKSkqQkZGBycnJgOeHh4dRXFwMi4X/v0dEyWOxJO2M/1lEfrXUFxWR\nShFpBVAHwKmU+my5ARLR2peWloZ9+/ZhYGAAo6OjM8eVUnA4HJiamsILL7wQwwiJiKIv5OpOpdRp\nEXkbWvHa0yJyCkAztJWe/uK1gFZ+wwIgH8Cb0BYN+JdoHV2FuIlojdm4cSNee+01tLa2ore3d2aV\nZ3l5Oerr65GbmxvjCImIomvROmlKqd16cvYutMSrUX8sph3AMaXU9ZWFSETJory8HOXl5XA6nZie\nnobRaER2dnaswyIiiokl7TiglDopIh9CW+3ZCG21pr+ALaD1qPmL214F0MzkjIiWi3PPiIjC2LtT\nL6lxQX8QERER0SpaUgkOIiIiIoouJmlEREREcYhJGhEREVEcYpJGREREFIeYpBERERHFoaCrO0Wk\nAcDxSDSilHp78auIiIiIyC9UCY5cRGa3ABWB1yAiIiJKKkGTNKXUBRGpAXASwLFZp07j2ZZPS8Ek\njYiIiChMi+3daQNwQkSc0LaFciil3o9KZERERERJbKkLBz5c1SiIiIiIaI4lJWn6llDDCG+Yk4iI\niIiWacl7d0Kbl5a7WoEQERER0TPhbLDOjdWJiIiIooTFbImIiIjiEJM0IiIiojgUzpw0iEhFmK9v\nV0q5wryHiIiIKOmFlaQBsCHM4rQiMwtCnQCaAXyolLoRZrtERERESSXcJO26/ly7jLb820y9KSIn\nlVIfL+MGOicxAAAgAElEQVQ1iIiIiCLO5/Ohu7sbN2/exNDQENLT07F161bU1NQgKysrJjGFNSdN\nKbUbQJP+rQ3ACQDVSimDUsoAwAqgHsA5/ZomaMnZ7OMC4JSI7Fp5+EREREQr4/V68eWXX+LSpUsY\nHx9Hfn4+TCYTvv32W/z+97/H4OBgTOIKK0kTkSMA3gPQrpSqUUqdU0p1+s8rpZxKqXal1AkAbwFo\nBPDpvOPvQ0vUPojcj0FERES0PHfu3IHNZkN5eTmys7NhMBiQnp6OoqIiGI1GNDc3w+PxRD2ucFd3\n+hOrYyGvwkxdtcsAGmf3mimlTutfHgyzbSIiIqKI8ng8uHnzJgoLCwOez8rKwvj4OHp7e6McWfhJ\nWh20hQMdS7y+XX+en5B1ArCE2TYRERFRRI2OjmJychLp6elBrzEajeju7o5iVJpwFw4MAzBDm192\nZQnX+xcYdM47XgnAEWbbREQUZwYHB2Gz2eByuZCdnY2qqirk5+fPXtlPtCYYDNEvLRtuktYM4AiA\nUwD2hLpQn7/m70Frn3XcP1TaFmbbREQUJ7xeL77++mt8//33yMjIQHp6Ovr7+3Hz5k1UV1dj3759\nSE0N908MUfRlZ2fDaDRiamoqaG+a2+1GWVlZlCMLf7jzpP68W0QeiMg7ImKefYGI1IrIvwM4rx86\nq5TqFJFKEfkIwBn9+Knlh01ERLHU3t6Oe/fuobS0FPn5+TCbzcjPz0dZWRk6Oztx7dq1WIdItCQp\nKSnYuXMnnj59CqUWloJ1uVzIyspCSUlJ1GMLtwSHDcAhaMOe1QDOAnCIiFdEhkTEC6AVwHH9lmal\n1C/1rw9CWxkKaInb5RVHT0REUed2u3H79m2UlJQsGNYUERQXF+P777/H+Ph4jCIkCs/WrVuxefNm\n9PT0YHh4GB6PB263G/39/fB4PDh48GBMeobDHmBVSjUrpXKhldLohFZOQ6DVQ/N/3Q7ghFLq0Kxb\nHQA+A9A4K3EjIqIEMzg4CJ/Ph5SUlIDnDQYDfD4fBgYGohwZ0fKkpKTg5ZdfxmuvvYbc3Fy4XC54\nPB7U19fjX/7lX2C1WmMS17LTQr2UxmkAEJEqAFUAOmbXTZt3/QUAF5bbHhERxQev17vowgARgdfr\njVJERCtnMBhQXl6O8vLyWIcyIyJ9d/owqC0Sr0VERPEtKysr4Nyd2Xw+X8y20iFaK1aUpOkrOCuh\nzU+zQkvUOgC0chN1IqK1KS8vDwUFBRgeHkZOTs6C8y6XC1arFQUFBTGIjmjtWFaSJiLvQludqaDN\nQZtPiUg7gGNM1oiI1p6XX34Zf/rTnzAwMIC8vLyZeWh2ux1KKTQ0NLBWGtEKhZ2kicinAN70fwtt\nnpkDgBPaLgL10IrY7gbQJiK7magREa0tubm5+Kd/+ifcvHkTDx48mDleXV2NnTt3BuxhI6LwhJWk\n6cOb/gTtpFLq10Guq4NWJ61Sf960kiCJKHFMT0/D5XJBKYWcnBykpaXFOiRaJWazGS+99BL27Nkz\ns61ORkZGrMMiWjPC7Uk7oT+fDpagAYBSql1EdkMr0VEtIgeUUkvZRoqIEpTX68WtW7dw8+ZNTE9P\nQ0SQkpKC7du3Y+fOnaw+v4alp6eH3PeQiJYn3Dppe6DNQzuz2IVKKSe0baQAbWN2IlqjfD4fWlpa\n0NraipycHJSUlKC4uBi5ubm4ceMGvvjiC5ZjICIKU7hJmn+SwdASr7eH+fpElIB6enrQ0dGBsrKy\nOcObqampKC0txaNHj/D48eMYRkhElHjCTdL8OwwcXOxCXcOs+4hojbp9+zbMZnPQ8xaLBTdv3oxi\nREREiS/cJM2/Y8A5EakIdaG+mXqV/m1zqGuJKLHZ7XaYTKag500mExwORxQjIiJKfOEmaR9C21w9\nF0CHiPxGRA6ISIWImEWkVkSOiUgrnm2mflopNRzJoIkovmRkZMDj8QQ97/F4OLGciChMYSVp+mKA\nBmiJmkBb7dkEbZcBB4BWaIsK/AsFziql3o9YtEQUl7Zu3Rqyp8zhcOAf/uEfohgREVHiC7cnDUqp\ndqVULoBfArgMwAUtYfM/OqENi+5WSv0ygrESUZyqqqpCZmYmhocXdpq7XC6kpaWhpqYmBpERESWu\nZRcuUkqdBXDW/72IWPSeNiJKMpmZmXjttddw+fJl9Pb2Ii0tDSKCqakpZGdn47XXXsO6detiHSYR\nUUKJWHVJJmhEyc1iseAXv/gF+vr60N/fDwAoKipCSUkJUlJSYhwdEVHiYQlwIooYg8GAsrIylJWV\nxToUIqKEFzRJ01doqkg0opTaE871+h6hTgBVSqlzAc6fUkqdFJFjgc4TERERJbpQPWkx2cpJ35wd\nSqnLIlIlIrVKqevzLjsmIm/g2V6iRERERGtKqCTtrQi1EW5v3FsALulf26DtbrAgSVNKfbbSwIiI\niIjiVdAkTSl1Idi5VWbB3D0/8wJcYxWRBgB1SqlfRycsIiIiougJu05aMCISfOO+ZbxcqJNKqXNK\nqcsA8vRkjYiIiGhNiUiSJiIPoe04EAlOAFb961wAQ/PaOqYvLIB+rgpEREREa0ykSnCE7PkK0+8A\n1EPbzaAS2rZTs4vl2qBtPwVoQ6FNs28+fPgwPv7445nv9+/fj/3790cwPIqm1NRUGI3GWIdBccDl\ncuHhw4d49OgRRASVlZWorq5GVlbWqrddX1+P0tJS/lukuMHfjYmrpaUFLS0tS7pWlFp5lQ0R6QBQ\noZSKSMVKETkGLRmbKcEhIq1KqXr9a39PWqVS6uN59yq73Q5aG4xGI9xud6zDoBiz2WxoaWmBwWCY\nScpGRkYAAAcPHlz1umyffPIJ6uvrubUVxY1Y/26cmpoCAKSnp8cshrXCarVCKRWwsysui9nOqn12\nedax+llfc2UnUZIYGhrCl19+ifz8/Dl/EEwmE9xuN5qamvDGG2/AbI7ktFgiCuTx48e4ceMGBgcH\nAQB5eXnYuXMnNm7cCJFIDqoREMGFA0REq+Hu3bvIyMgI+H/sRqMRBoMBHR0dMYiMAnG5XOjt7cWT\nJ08wPT0d63Aogtra2nDx4kW43W6UlJSgpKQEU1NTaGpqwrVr12Id3poUqZ40GyK0OwER0Ww2mw35\n+flBz+fk5KCjowO1tbVRjIrmGxkZwdWrV/Ho0SMYDAYopZCWlobnn38e27dvh8HAPoFE1t/fj+vX\nr6O0tHTOZ5mVlQWTyYTvvvuOW8KtgogkaUqpxki8DhHRfD6fL+QwioggEnNraflGR0fxpz/9CR6P\nByUlJTOf1/T0NP7+979jbGwML7zwQoyjpJW4e/cu1q1bFzDZ9s8VvXPnDpO0COP/2hBRXCspKYHL\n5Qp6fmRkBCUlJVGMiOa7c+cOJicnkZeXNyehTktLQ1lZGW7fvg0u6EpsfX19yM7ODnrebDajr68v\nihElByZpRBTXtm/fjpGREfh8vgXnvF4v3G43tmzZEoPICNB6y+7du4e8vECbw2i9LOnp6Zw3mOBS\nU1MD/jfo5/P5kJISkQIPNAuTNCKKa+Xl5di5cyd6e3tnym4opTA8PIy+vj786Ec/QkFBQYyjTF6T\nk5OL/oE2Go0YHh6OYlQUadXV1XA6nUHPOxwObNq0KYoRJYe4LMFBRDTbnj17UFJSgm+//XZmSKW0\ntBT79u1DaWlpjKNLbmlpaVBKwefzBV0cMDU1hczMzChHRpG0adMm3Lx5E263e0ERXbfbDa/Xi82b\nN8courWLSRoRxT0Rwfr167F+/Xp4PB6ICIdW4kRGRgYqKirQ19cHq9Ua8JrJyUlUV1dHOTKKJLPZ\njMbGRly5cgUOhwPZ2dkQEYyMjEBE0NDQAIvFEusw1xwmaUSUUFJT+Wsr3jz//PN49OgRxsfHYTKZ\n5pzr7+9HeXk5ioqKYhQdRUpZWRmOHDmCzs5O/PDDDwCALVu2oKKiIirbsyUj/rYjIqIVsVqt+OlP\nf4ovvvgCTqcTKSkpM0OglZWVeOmll1iNfo0wmUzYvn07tm/fHutQkgKTNCIiWrGioiIcPXoUfX19\ncLlcSElJQXFxMXJycmIdGlHCYpJGREQhORwOPHjwAD/88ANEBBs3bsTmzZsX7JeakpKC8vLyGEVJ\ntPYwSSMioqAePHiAr776CmlpaTNJ2Z07d3Dr1i288sorqKioiG2ARGsY66QREVFAg4ODaGlpQUFB\nAQoKCpCRkYGMjAwUFhbCarXiiy++YP0zolXEJI2IiAK6d+8ejEYj0tLSFpxLT09HSkoK7t+/H4PI\niJIDkzQiIgqos7MzZO0ri8WCrq6u6AVElGSYpBERUUBKqZClM0Qk5H6ORLQyTNKIiCigsrKykHPO\nhoeHuZqTaBUxSSMiooC2bduG0dHRgL1lHo8Hk5OT2LJlSwwiI0oOTNKIiCigkpIS7NmzBz09PXA6\nnTO7CDgcDvT39+PFF19Efn5+rMMkWrNYJ42IiILatWsXCgsLcfv2bXR3dwMANm7ciO3bt3M/TqJV\nxiSNiIhCKi0tRWlp6aILCYgosjjcSURES8IEjSi6mKQRERERxSEmaURERERxiHPSiIhoVU1MTMBm\ns+Hu3buYnJyE1WrF9u3bUVZWhpSUlFiHRxS3mKQREdGqGR4exueff46JiQnk5uYiMzMTo6OjuHTp\nEiorK/GTn/wEqan8U0QUCIc7iYhoVfh8Ply+fBkigpKSEhiNRqSkpCA7Oxvl5eV49OgRvvvuu1iH\nSRS3mKQREYUwPT0Nu92O7u5u9Pb2wuv1xjqkhNHf3w+n04mcnJyA5/3116ampsJ+bZfLha6uLnR1\ndcHlcq001LiilMKTJ0/Q2dmJx48fY3JyMtYhUYywj5mIKIj79+/j6tWruHXrFkZHR9HZ2Qmj0YiX\nXnoJGzZsiHV4cW9gYCDkUGZqaiq8Xi+Gh4dRUFCwpNecmJjA119/ja6urjnHN2zYgL1798JkMq0k\n5Jh7+vQpvvrqKwwPD8+UPBERPP/883j++edhMLBvJZnw0yYiCuDevXtoaWmB2WyGxWJBQUEBSkpK\nkJmZiYsXL+Lx48exDjHuRbr47dTUFC5duoSenh6UlJTMFNktKSlBf38/Ll26tKxeuXgxNDSEP/3p\nT5iYmMDU1BQGBwcxPDyM1NRUXLt2Da2trbEOkaKMPWlERPNMTU3h6tWrKC4uXtATZDQakZ+fj6+/\n/hplZWXs2QihsLAQ09PTQc97PB4YDIagw6HzPXr0CENDQygtLZ1zXESQn5+P3t5edHV1YfPmzSuK\nO1ZaW1sxODiI/v5+iAjS09Ph8Xjw+PFjWCwWXL9+HVu3boXZbI51qBQlTNKIKOG4XC48fPgQ9+/f\nx/T0NAoKCrBjxw6UlZVFpOemv78f09PTQYfqMjMz4XQ6MTg4iMLCwhW3t1YVFxcjJycHw8PDAROx\np0+fYufOnUhPT1/S692+fRsWiyXoeYvFgtu3bydkkjY6Oopvv/0WfX19yM3NXVCaZHh4eGYe3s6d\nO2MUJUUb/xeQiBLKkydP8Pvf/x63bt2CyWRCXl4eXC4X/vznP+Ovf/0rfD7fituYmJhYUg8ZJ3SH\nZjAY0NDQAKUU+vv7MTk5CZ/Ph9HRUfT09GD9+vV4/vnnl/x6Y2NjyMjICHreaDRidHQ0EqEHNT4+\nDofDgfHx8Yi+7uTkJLq7u5GdnR2wdlxOTg5GR0fXzDD76OgoHA4HJiYmYh1KXGNPGhEljMnJSTQ3\nNyMrK2vOBHGz2Yzs7Gzcu3cPBQUF2LJly4rayczMXDTZU0qFTBhIY7FY8Prrr6OjowN3797F8PAw\ncnNz0dDQgPXr14dVzHbdunWYnJwMujjA7XYjKysrUqHPYbfb0d7ePpMkKaWwceNG1NXVITc3d8Wv\n73a7MT4+jqKioqDXKKVWPQldbU+fPkVbWxv6+vpgMBiglEJVVRVqa2s5jBsAkzQiShj+cgR5eXkL\nzvnnJX377bfYtGnTiuaKFRcXIy0tDR6PJ+CQ58TEBLKyspCfn7/sNpKJyWTCc889h+eee25Fr7N9\n+3Z89dVXQZM0p9OJl19+eUVtBDIwMID/+Z//QUZGBoqLiyEiUErh6dOn+OMf/4if//znK/63kJmZ\nidzcXIyNjWHdunULziulACCh/8319PTg4sWLWLdu3cy8Qp/Ph8ePH6Onpwc///nPlzw/MVlwuJOI\nEkZ3d3fAP2B+RqMR4+PjGBsbW1E76enp+NGPfoS+vr4FqwXdbjeGhobw4osvctFAlG3cuBF5eXkY\nGBiYSVoALYEZHByE1WpFRUVFRNtUSqGlpQVZWVnIzc2dUxYjNzcXJpMJX3311Zx4lsNkMqGiogKT\nk5Nwu91zzvl8PtjtdhQVFWH9+vUraidWPB4P/vd//xdWq3VOImYwGJCfnw8RwTfffBPDCOMTe9KI\nKGH4h0dCWekfS7+tW7fCYDDg6tWrcDqdM4mByWTCoUOHEvaPZbRMTU2hu7sbP/zwA5RSKC0txYYN\nG5CZmbns10xPT8ehQ4cC1knbuHEjXnzxxSUvQliqp0+fwul0oqysLOD57Oxs9PT0YHBwcMm13gIx\nGo2ora2dKWRrt9tnzhkMBmzYsAEWiyVh6/P19/fD7XbDarUGPJ+bm4vu7m64XC4Oe87CJI2IEkZZ\nWRkePnwYdIWff05SqN62cGzevBmVlZVwuVzYtm0bNm/ejKKiIm4KvoihoSE0NTVhYmICJpMJIoKu\nri5cvXoVr776KsrLy5f92pmZmThw4ABcLhccDgeUUrBarav2h318fHzRHlODwYCxsbEVJWkAsGvX\nLvT29iIrK2tmuN1gMCAlJQXj4+PYv38/jEbjitqIFafTueh/NwaDAaOjo0zSZmGSRkQJY/369TAa\njQHn7SilMDAwgJ/85CcRHYZMS0uD1WpFeXn5gvpctNDExAQuXryItLQ0lJSUzBw3m81wu91obm7G\nP//zP694sr3ZbI7KH/OlJORKqbATd4fDgbt376KjowNerxfFxcXYsWMHfvazn+HWrVu4e/culFLw\neDzIy8vDq6++Ouf9TDTp6elLWozD/wGai0kaESWM9PR0NDY24s9//jNGR0dhsViQkpKCsbExjIyM\nYPv27di0aVOsw0xqXV1dmJycDDisZTQakZqairt372Lv3r0xiC58BQUFEBH4fL6Ayb/X64XBYAir\nF+3x48doampCenr6TE00fxmZHTt24Mc//jF27doFt9uNlJSUiPUMx1JxcTGUUkF3oZiamkJGRkZC\nL4xYDZz1SkQJpaCgAL/4xS+wa9cuTE1NweVywWq14qc//SlefPHFiG5DROG7f/9+yBV6VqsVDx8+\njGJEK5OZmYkdO3agt7d3QU+Qz+dDX18fnn/++SUPQ46OjuLy5cvIy8tDfn4+UlNTISIwm80oLS3F\nzZs30dXVhfT0dJjN5jWRoAFaz+emTZvQ19e34JzP58PTp09RW1vLnrR52JNGRAknKytrZsPp1eBy\nuTA4OAiv1wuz2RyxxQjJYHp6OuTkfYPBAK/XG7RnKh7V1tZiamoK9+7dQ2pqKtLT0zE1NQWPx4Md\nO3aE9e+ws7MTAALW2DMYDMjNzcW3336LysrKiMUfL1544QV4PB50dnYiLS0NaWlpmJychNfrRX19\nPbZu3RrrEOMOkzQiIt3U1BS+/vprdHR0zPTI+Xw+3L17l8OoS1RYWIjHjx8HrGUHaD1JVqs1YRI0\nQJuXtnfvXmzbtg1dXV0YGRlBdnY2KioqQm5TFUh3d3fIgrtZWVkzuzOstWLJaWlpePXVV7Fz5050\ndXVhfHwcOTk5qKio4GKBIJikERFBS8a+/PJL9PT0oKSkZM6w6a1bt9DS0oLNmzfzj8kitm7divv3\n7wftKXM4HDh48GAMIls5i8WCXbt2rXo7weZtrRV5eXlBk3iaK3H+V4aIaBX19/fj8ePHCxI0QJuX\npJTC7du3YxRd4igsLMSuXbvQ09MzZ39Lt9uN7u5u1NTUYOPGjTGMMLYqKipCbu00MjKCgoKCiNd7\no8TEJI2ICMCDBw+CbjcEaBtcf//99/B4PFGMKjHt3r0bBw4cAAD09fWhr68Pk5OTePnll7F///6E\nGuqMtI0bN87UPZvP5/NheHh41eZaUuLhcCcREbSipaF6L/yrzqanpwPu57mWjYyM4Pvvv4fNZgMA\nbNmyBdXV1UHnVokIqqurUVVVNbNFl8lkWnPJ2cTEBMbGxpCSkoKcnJwl/XwmkwkHDx5EU1MTRkdH\nZ+4bGRnB+Pg4du/ezd0saEZy/aYhIgrCbDbDbrcHLXng9XqRkZGRdMNQN27cwB/+8Af09fUhLS0N\nAPDFF1+gtLQUR44cwbZt24LeKyIhJ8knqrGxMVy/fh0PHjyY2Ww9KysLdXV1qK6uXvT+kpIS/OIX\nv0BHRwcePnwIr9eLsrIybN26FUVFRVH4CShRMEkjIgJQU1ODe/fuBZ207XA48MorryRVHae7d+/i\nv/7rv+Dz+VBWVjbTU+TxeDA4OIj//M//xIkTJ5aUmKwV4+Pj+PzzzzE+Po6ioqKZ98TtduOLL77A\n2NgYqqur8eTJE0xNTSErKwvFxcULel+zs7Oxa9euqCxEoMTFJI2ICNqE95qaGjx8+BDFxcVzkrGx\nsTEUFxeH7DVaazweDy5evIjp6ekFvTupqakoLCxEf38/mpqaUFlZueaGMoO5desWxsbGFrwnRqMR\nRUVFuHDhAkpKSpCRkQGDwQCfzwej0Yh9+/ZxGJPCxiSNiAja0NzLL7+MdevW4datWzMFbP1DWfv3\n718z1d+Xor+/H52dnSF3DzCbzXj48CEGBgaSYphuenoa9+7dC1o+4uHDh+jr60NRUdGcfV7dbjcu\nXbqEn/3sZwm9/yZFH5M0IiJdSkoK6uvrsWPHDgwODkIpBbPZjM8//3xNzq0KZWJiAhMTEygsLAx6\njdFohMPhwMTERBQji52JiQn4fL45Q5c+nw8igvHxcXR3d6OwsHBmsYSf0WiE2WzGtWvX8Prrr0c7\nbEpgTNKIiOYxGo0oLy+PdRgxlZqaipSUFHg8nqCrWf3nkmW1a2pqKpRS8Pl8ePLkCbq6ujA2NgYR\ngcFgwPT09Mx2R/NlZ2ejt7cXLpeLBZFpyZJjEgEREYWloKAAJSUlGBkZCXrNyMgISktLkZ+fH8XI\nYsdkMqGoqAh///vfcfPmTSilkJubC7PZPJO0ORwOFBcXB7zfYDBgamoqylFTImOSRkREC2RlZeHH\nP/4xxsfHAw5njo+Pw+12Y+/evTAajTGIMDaysrLQ3d2NnJycmb01DQYDrFYrPB4PnE5nwKLI/h64\ntbYfJ62u5OijJqK44vP50N/fD6fTiZSUFBQWFiI3Nzcu4jEYDCgoKAg4Odzj8aCvrw/Dw8NIT09H\ncXFxVIaunE4nnjx5gunpaVgsFpSUlESlFMgrr7wCu92Ov/zlL3C5XMjMzATwbG5WY2MjXnzxxVWP\nI154vV709PTgpZdewv379+HxeJCWlgafzwePx4Ps7GwUFxdjYGBgwXC50+nE+vXrkZ2dHaPoKREx\nSSOiqBocHMSVK1cwOjqKlJQUKKXg9XqxYcMG7Nu3L+q9Mna7HVeuXIHL5ZoTT3l5Ofbt2zfTK9Lf\n349vvvkGk5OTSElJgc/ng1IKmzdvxgsvvLAq87Kmpqbw17/+FZ2dnTAYDBAReL1emEwmvPrqq6u+\nojIjIwNHjhzBc889h6+++gr9/f0zuwns378fNTU1STMfDXjWe1hWVoaioiIMDAxgZGQEqampyM/P\nh91ux40bN/D06dM5Sdrw8DCmp6dRX18fw+gpEYl/mflaISLKbrfHOgyKEKPRCLfbHeswKEKGh4fx\nhz/8ASaTacFqyYGBAeTl5eHw4cNRKxjrcrnwxz/+Eenp6Qt6xAYHB2E2m/Gzn/0M//Ef/wGHw4HS\n0lI4HA64XC6kpaWhuLgY09PT2LRpE1555ZWIxqaUQlNTE3p7exfMcRofH4fL5cLrr78Oq9Ua0XaD\n8fl8mJiYgIggMzMzYMHfpVBKYWBgAA8ePIDdbofJZMLmzZtRWloa94WCR0dHcf78+ZBlNO7evQun\n04mKigqICHw+HwoLC/Hiiy8GLd2xHPzduHZYrVYopQL+B5U8/wtERDF3+/ZtGAyGgOUsCgoK0NPT\ng76+vqitrPz+++/h8/kCDlnm5+ejp6cHvb296OzsnEkuUlNTkZGRgcnJSdy8eRMmkwlTU1PYuXPn\nshIm/1BZWlranMTnyZMnePz4McrKyhbcYzKZ4Ha78d1330U8OQzGYDCsuE6cz+fD119/jbt37yIz\nMxOZmZkYHBzEo0ePUFBQgIMHD84MqcajdevWwWw2Y3x8POC8MwDIycnB66+/jnXr1sHj8WDdunUx\nHcqnxMYkjYiiwuPx4P79+yHrbq1btw4PHz6MSpLm8/lw7969kIlVdnY2bty4gZ6eHqSnp6OysnJO\nIpWZmYmxsTE8fPgQnZ2dYSVpLpcLd+/exffffw+v14vMzEzs2LEDmzZtQkZGBjo6OkIO/VosFnR1\ndWFycjJhJqPfvHkTd+/eRVlZ2Zz3MScnBwMDA2hpacHhw4djGGFoIoJdu3bhypUrMBqNC3ZZcLlc\nyMrKwoYNGwKW4SAKF1d3ElFUeDweAAi5fVB6evqCQqCrxev1hqwB5o/H6XRiaGgIWVlZC4b4vF4v\nAG1fz66uriW3PTg4iD/84Q+4f/8+8vLyUFJSgszMTFy9ehWff/75TCHZUJu5+9/H6enpJbcbS9PT\n07h58yaKiooCDpX6e1KHhoZiEN3SVVVVoba2Fn19fRgaGoLb7cb4+Dj6+vpmFlMwQaNIYU8aEUVF\neno60tLSZgp+BjIxMRG17YVSU1ORmZmJqampmWTI4/HMlJyYmpqCy+WCyWTC+Pj4nHv9q0H9uxI4\nHA588cUX2LBhA3bt2hXyj7TX68WVK1eQmZk5Z9g3IyMDpaWlGBgYQFtbG3JyctDX1xd0NaDH44HB\nYMW8lzYAABXWSURBVEiYXjSHwxHyswe0HR/6+/sjOncr0kQE9fX12LhxI+7fv4+nT58iLS0Ne/fu\nxcaNG5OqHAmtPiZpRBQVBoMB27dvx40bNwIW+1RKwe12Y9OmTVGJR0SwY8cO/P3vf0dmZiZ++OGH\nBb04w8PDqKqqgtvtxp07dzA0NITCwkLY7fY585ImJydRUlKCmzdvwul04sCBA0Enwff392NkZCTg\nXDMAyMvLw8OHD9HY2IjvvvsOPp8vYO+j3W7H1q1bZ5IepRSePn2KgYEB+Hw+FBQUoKioKG42Pvdv\nnxSKf/VqIigoKEBBQUGsw6A1jkkaEUXNli1bZjbkzsvLm0kgPB4P+vv7sXXr1qj+4cvOzsb9+/fh\ndDqRl5cHi8UysyLP6XRiy5Yt2Lp1K/72t79BKYXJyUncv38fdrsdxcXF8Pl8GB8fh8ViQWlpKUpK\nSvDo0SN0d3dj48aNAdscGhoK2Zvkf09SUlKwY8cOfPfddyguLp6TjNntdhiNRmzfvh2ANhfqyy+/\nnFnYAGg9dhaLBQcOHIDFYonk27Ys/sUZwZJOQPt3EK3VqkSJgEkaEUVNZmYmfvrTn+Lq1avo6uqC\niEAphdTUVOzZswfPPffcsks7hGNqagrt7e24desWamtr0dPTM1PMFtASpaqqKlRWVsJgMKCwsBBe\nrxcpKSkYGBiAiKC3txcmkwk1NTWwWCwzWyNZLBbcunUraJJmMBjg8/lCxjcxMYHbt29jZGQEU1NT\nuH37NiwWC4xGI3w+HzZu3Pj/t3dvv01c+x7Avz9jyJU4IWlCwqW5VSgFmksLtFWklk2kShReWvbu\nSx9PqM77KWj/BYV2v++ye/6Bkva8VSLhoqpSJVRCgTSUIkIINIQEgu1c6lycZJ2HmXHGjmc8duyJ\n7Xw/0lJC1oxnjWdq/zprrd/C4cOHUVJSgoWFBfT29mJpaWnN07mpqSlcunQJJ0+eXPfMzPUqLi5G\nY2MjRkZG4k4emZ2dRWlpqW16C6LNhkEaEbmqpKQER48exezsLKanpyNL6tgNkk+n+fl5XLt2DRMT\nE6irq4PH40FFRQWam5sRCoUgIti+fXvUhILKyspIEtuKigoUFBRARDA3N4exsbGocWhFRUWRYC+e\nnTt32nbp/fnnn7hz5w4KCwvh8/kikwrm5ubQ1taGpqamqHFqIyMjmJ2djRvc+Hw+TExM4OHDh3jj\njTdSebvS6tChQ5icnMSzZ89QWVmJbdu2YWVlBX6/H8vLyzh+/HjW50ojchODNCLaEKWlpXHzpWVS\nOBzGtWvXMDk5uSaoKSwsXDPoOxQKYXR0NDIL0+PxIBQKQSmFgoICVFdXw+fzYWRkBFVVVSgrK0M4\nHLYdzF9VVYXa2lpMTk6uWZh8amoK169fR2tra1T7iouLMTc3h4GBATQ3N0ft8/vvv8Pn81ker6Ki\nAvfu3cuKIK2oqAjHjx/H/fv3MTg4iMXFRXg8HjQ3N+PAgQO250G0GTFII6JN49atWxgfH0ddXV3C\nbWdnZ9Hf34+VlRUUFRWhoqICZWVlWFpagt/vx759+yLBWCgUwm+//YYjR44gEAjg7bfftn3t9957\nD1euXMHY2BhKS0uxbds2hEIh3LlzB3v27FkTiAFagDM9PY3h4eGogGthYcE22N26dSuyaRWWwsJC\ntLa24uDBgwiHw/B6vSk9PZufn0cgEACgdTFncxJcolQxSCOiTWFiYgIDAwOOAjSlFAYGBiJpOkZH\nRwFog/nr6+sxNTWF4eFhtLS0ANCedAUCAQwMDKCpqQn19fW2r19SUoIPP/wQT58+xf379xEKhVBX\nV4fJyUk0NTVZDqwvKytb03VZXl6O6elpy4Xe5+bmsmLiQKxU04cY4wn/+OMPGMsaigj27duHN998\n07VucyI3MEgjorynlMLPP/+M8vJyRykppqam8Ndff8Wdaej1erF//37cvXsX4+PjkS7S5eVljI+P\n49NPP3X0VMfr9eLVV1+NTDBYWVnBo0ePbJ8qeTyeSFJgw/79+9Hb22sZpAUCAdeWjso0I8fcs2fP\notKLGKtHTE1Noaura1Mt+k75jXcyEeW958+fIxgMOnqKBmjdl3bBXGFhIRoaGtDY2BgJqnw+H8Lh\nMMbHx7F3796k2+jxeFBTU2P7VGxmZgZNTU1Rf6urq0NDQwMeP36M6urqqBQcz58/R21treVM01zz\n5MkTjI6Orlk2zOPxoLa2FqOjo3jy5AkaGxs3qIVE6cUgjYjy3r1795LKBG+kBkm0TXl5eVR2/KWl\nJdy7dw9tbW0pdbsdPHgQvb29KC0tXRMkLi0tYX5+HiUlJejr68PMzAzKysrQ0tKCzs5O+Hw+3L17\nNyppbEtLCzo6OvLmydLg4KBt162R/oRBGuWL/Pgvl4jIwsrKCh4/fpxUklyfz2cbpBmBUOySTV6v\nN5JSIt6qCons3r0bBw4ciAQjpaWlUEphenoafr8fIoIbN25g+/btKCgogN/vx6VLl7B37168//77\nOHjwYCT9h8/ny7sliqanp22DtOLi4qxf+5MoGQzSiCivzczM2Ga5j6e4uBg1NTV48eJF3KAgGAyi\nvr4+7tMyEUEgEEgpSBMRHDlyBLW1tbhz5w7GxsYgIti1axe8Xi9mZ2ej1jYtKCiAz+fD06dP8csv\nv6Czs9O1tU83QmFhIcLhsOW4vcXFxbwLTGlzY5BGRHltZmYmpf1aWloQDofh9/uxsLCAxcVFzM7O\nYnFxETU1NWvGhhmKioowMTERmfmZLBGJTCgIh8MQEYRCIXz33XeW2fhramrw4MEDtLe3u7KyQDgc\nxuTkJJaXl7F9+3bX8pu9/vrruH79uuX7EAgEcPjwYVfaQuQGBmlElNfC4XBK+23duhXt7e3w+/14\n/vw5tmzZgsrKSuzatSuyxmc8Xq8Xc3Nz62lyVBsAbb1PpZTlMY2nhC9fvsxokKaUwt27d3H79u3I\n+6qUQl1dHd55552MB2sNDQ0YHBxEIBBARUVFVF0wGERJSQnHo1FeYZBGRHkt0QQAOx6PB1VVVdiz\nZw9ee+21DVtXUinlqLt2PefqxI0bNzAwMICampqoReIDgQB++OEHnDhxwnJmajoUFhbigw8+wI8/\n/oinT59Guj2Xl5dRVVWFo0ePMqkt5RUGaUSU17xeb8aDF7Pl5eW0J1T1+XwJF2VXSmU0aW0wGMTg\n4GBkvVOziooKTE5OYmBgAJ2dnRlrA6Al9D158iRevHiByclJANraqtXV1ZZPGolyVVYGaSLyMYAg\ngEal1DfJ1hMRGUpLS1398p6bm0tqJqkTlZWVqK2tjdvNB2jdnLt3785od+PIyAi8Xq/lE70dO3Zg\naGgIhw4dSmklgWSICKqrq1FdXZ3R4xBtNOfTnVwiIh0AoJS6qv+7PZl6yi8//fTTRjeBclxZWZmj\nvGeJDA0NOdpuZWUlKndaunR2dsLj8WBiYiIyHmxxcRHj4+PYtm0b3n333bQf0ywYDNoGX0bwtrCw\nkNF2kIafjZtD1gVpAP4BIKD/PgygK8l6yiP8IKL18nq92LlzZ8qzPA1OgjSjSzLe0671Mrr52tra\nMD09jWfPnmF2dhYdHR04ceKE7SLr6VBSUoLFxUXLeqUUlFJRY9Uoc/jZuDlkY3dnOQC/6d+x/0ua\nqJ6IKMqBAwfQ19eX0UHtgPa0qaGhAcXFxRl5/eLiYrS2tqK1tRVLS0uuriRQX1+PgYEBy/pgMIg9\ne/Zw4D5RGmXjkzQASDSAhKNDicix2tpaFBUVYX5+PqPHmZubSzk/WrLcXurplVdeQWNjI8bGxtZ0\nHYdCISwsLKCtrc3VNhHlO3Fz1pMTInIOwGWl1FUROQWgQSn1VRL12XVCRERERDaUUnEfPmVjd+e3\nAN4CcBVAA4DLACAi5UqpoFW9wepEiYiIiHJJ1nV3KqVuAYCIHAMQVErd1quuJKgnIiIiyhtZ191J\nRLTRmKuRsomD+/G8UuqsiHTzfswvWfckjQjQPpRE5JiIdFvUn9d/xq0nShVzNVI2cXi/dYvIAwAP\n3WwbZR6DNMo6/FCiDcZcjZRNnNxv3Uqp15RS19xrFrmBQRplI34o0UZirkbKJk7utx16z8PnLrWJ\nXMIgjbIRP5RoozFXI2UT2/tNKfWN3vNQqU+qozzBII2yFT+UaKMEAezQf68A8DLJeqJ0sr3fRKRb\nn1gAva7RxbZRhmVjnjTaBCwG/PuVUt/DwYeSaVvjQ+lqBptLm8u6cjUSpVmi+3EYQL++bSV4P+YV\nBmm0IRJME+eHEm0YpdQtEXnLIlfjWzb1RGnn4H68qs+GB4BJ3o/5hXnSKCvpT8uGYcoLJCL9Sqm3\n9N+Nx/sNSql/bVAziYiIMoZBGhEREVEW4sQBIiIioizEII2IiIgoCzFIIyIiIspCDNKIiIiIshCD\nNKIcIyLlInJBRG6KSEBEVkTkoYhctFtwXkRO6dv2udleN+nvyYpp9m+mjnPe9N53O63LF5vhXiLK\nBgzSiHKIiJyBtmRWN4B2AGUAFIB6AKcAXBCRIRFpsHmZzTClO2PnqF+DzwH49OMoJ3VuEJEzInJO\nRHwuHXIz3EtEG4ZBGlGOEJHTAM7p/7wMbeH5Cr28BeBLva4R8RP8BqAlwPw1sy3Ne5/oPy8opbYo\npf7XYZ0b/gngDLQk0JnEe4nIBVxxgCh3nNd/9iilPompuwXgloh8C+AmgEYR6Tav7KCvdcrls9av\nXP/Zk2SdWzL+dIv3EpE7+CSNKAeISCNWu9C+sNpOKXULwH/0f55yoWmbmaRY54aNPj4RpQGDNKLc\nUG76fSrBtpehdUMNm/+YaLC3iDSKSI9p0Hu/MQDfNCC/3rR9l/63i6Z/G9sFRKRPX2/Qkj6GyrxP\nf6LB9iLSEdPOIRH5PMF7su426BMCVrDalXjZmBxgV5eG8zUmijw07XdRRNpN2/Tox/dBC9BSmkAh\nIqfjTEhZM8bN6l4SEeO8bYvFsZN+b4jynlKKhYUlywu0L98VvZxL8TVO6fv3xqnrMr3+sl4ix4PW\nhboMoD7OPhehjYMy72ve/+M4xyvXX9PqmP0W53DaZp8+AP367x85eD+SagO0yRr9pvoHAG4A+Miu\nbp3n2wFt/JfVft36duf048Ue/29J3B+XbY7jB+BLdC/p98JLixJ57XTcCywsm6FseANYWFiclZgg\noA/AsST3twvSjEDgBoC/QZs12h5zTKsgbUj/+W8AZXpdg2lff5zjGQHBSz3IKdP3+dx0vIsx+zTG\nBCFGO4/p78eKqTgJ0pJug77fQ71uTQBkVZfi+ZbHXhfTe2s+X3PwFNCvU1uS94Y5yP7IdB2PQQvQ\nov7nwO5esjmG8d78Ox3XgYVlM5QNbwALC4uzon9x+REdjBgBWzeAhgT7Wz39ML6gH1jsZwRhVkHa\nCoBvLdpr7FcWZ7+o1zPVHzO9brvp7xeMdppfz1RvDihtg7RU26DXJRWkreN8nV6Xz01/SzVIMwKl\n/4lTZwRLNxLdSzav3xP7Guu9Diwsm6FwTBpRjlBKPYIW+JyFNt7MmMXXBS2AeZji+KzP9J/nLeqt\n/h5pmt6meO017DD9/nf953+UUiNx9rsKLb0DoJ2bwRifdFYpNR2nHWvaYCPVNqQi1WMZ1+WCxete\ngPbeN66zfcDq9Vkz4UAp9RWAJqyeR1L01DEfQ2tr7Gu4eR2Icg6DNKIcopSaUkp9pZRqhvbF+hmA\n76A9QTG+sM87SGhr1qDve9Gi/orF383tGnF4LEDL6ZbodW+Zt9Vnt+qHUv9n0YZkUkIk3YZ1SPVY\nxnX5Lt4O+n2wRSn13+tsH6B1pwLavXMu9t5RSj1K8hoD0CZ5APgaeoAW5zXcvA5EOYd50ohylFJq\nCsA3eoE+k/IzaF1RjdC6mGy/2EzBDyyeTkEp9UjENqPDsF1lHMYxexK8LrA6q9XYJ5hge+NpYyba\nkKqkjyUikWOmEhyl4Cy0J1WN0LpZz4jII2jdoD1JBsAAIudg7PelRXDt5nUgyjkM0ojyhP5FelXv\n7jwPoF1E2pWWO82K0+AnnYx0DjcdbGsEgMYXtD/B9k7PI5U2pCqVY6WjC9MxPeBv1lNefAZt0kgD\ntNm0p0VkGNqTMLt7KVYPtHO/qZT6p8U2bl4HopzDII0oB4jITWhfnGf1MUKWlFJfich5aOOL3sRq\nd1E8D/Wflk8pzE910uQRtLVGu5VStx3uY7QzUfDiNLhJpQ2pWs/5QkTKrJ5yppvSVqj4Rs+L1gVt\nmSvjyexVEWnQAzpb+v13DFo3vF2uPDevA1HO4Zg0otxgjBk6lOR+tk+ezIP7zclRY6T7qc4wtADS\n8lz0xLpdpiSq5nbG7c7Ug0mnC4un0oZUJX2smECoyWKfDj157NfrbJ9xfHPX95RS6nul1D9Mx/fB\nwbgwEemCNiPUGIdmF2C6eR2Icg6DNKLcYCyYfirRhAARMZaDUnC2APav0L4oP7Oot/p7qoyuLbvX\nvQkttUgDACilglj9QreaxXk6k21Yh1SP5eS6HEN61uocAhB3sokeyN/S22IbKOmBco/epi+VUtcS\nHNfN60CUcxikEeUApdT3WB2Tc9NquSX9KYYxS/OKw0HnRoqH07HLCOkBX7qX5jHWHu0QkXOxlSJi\njGUKxHSBmdsZu9xSF7Ss+5luQypSPZaxX7zr0gHtuiisBvBRmyTZRuPeWpNuRT9WO5wF/VehncsV\nm3FoZm5eB6Lcs9GJ2lhYWJwVrE1mG4D2lKEH2he1efmgl4hJDgr7FQfMGez7oH1ZGwlOe2GflDZu\nslV9G+M1Y9vSbaobghaAXcBqMti4CVkRnbC2X9/nsmkf4zycrDiQahtSWXEg1WOZr0u/fl16TH+L\nTUxsJLh1lODYpn1fx9wDUQmL491L0JawMrb16/djvOKHaamwVN8bFpbNUDa8ASwsLM4LtKcKFxG9\ntqG5LEN7OhEvI//H8b7YTfVfx3m9XqyuGxq75qKRDT5RkGaVTb7b4hwexAuATOffF+f8X0JbJuqc\nXpcwSFtHG4b0Y8QL0uzqkj6WzXVZhmkZLtO252K2c/Q+xOwb+94uY+1STmvuJWhP24ztE5X/Ssd7\nw8KS70WUSsdwBiJyk2n2XQW0mZnDRlHrmAmoj0nqgPbU7opS6rY+zsgPLZHslnU3fu0xP9aPF4TW\n/kTjmIx2dkEL2oaVRYLbTLbBzWPFni+AX5VFV7beFd6uv36PcjAb07SvMTmgAdp99Su0Bc5dmV3q\n5nUgygUM0og2OX3cWQW0oOxRnPoOaF1tAaVUpdvtIyLarDhxgIg+gTYGKNHszoTLQxERUfowSCMi\nI8/WmTizJs9gdRbhF7E7EhFR5rC7k4igJ0Q15xkzr4GpoK108C/XG0ZEtIkxSCMiAJEB52ehLSVV\nDm3w9hUAXyjmqCIich2DNCIiIqIsxDFpRERERFmIQRoRERFRFmKQRkRERJSFGKQRERERZSEGaURE\nRERZiEEaERERURb6f9Apz9/mpSzPAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# These 4 lines are for the notebook only if mpld3 is enabled, \n", "# the matploltib figure does not show up\n", "try:\n", " mpld3.disable_notebook()\n", "except:\n", " pass\n", "\n", "results.volcano()" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": false, "scrolled": true }, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\n", "
\n", "" ], "text/plain": [ "" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import mpld3\n", "mpld3.enable_notebook()\n", "mpld3.display(results.handle_volcano.current_fig)\n", "# 1. hover the mouse on circles to get more information.\n", "# 2. in the bottom right corner you can find small icons to click on to focus in the images. Click again \n", "# to be able to hover the mouse agains" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.5" } }, "nbformat": 4, "nbformat_minor": 0 }