diff --git a/Analysis/Theory/Anscombe transformation.ipynb b/Analysis/Theory/Anscombe transformation.ipynb index b6d23b5..fa66648 100644 --- a/Analysis/Theory/Anscombe transformation.ipynb +++ b/Analysis/Theory/Anscombe transformation.ipynb @@ -41,7 +41,9 @@ { "cell_type": "code", "execution_count": 11, - "metadata": {}, + "metadata": { + "collapsed": true + }, "outputs": [], "source": [ "merfish = pd.read_csv('../MERFISH/data/rep6/middle_exp_mat.csv', index_col=0)" @@ -50,7 +52,9 @@ { "cell_type": "code", "execution_count": 19, - "metadata": {}, + "metadata": { + "collapsed": true + }, "outputs": [], "source": [ "seqfish = pd.read_csv('../SeqFISH/exp_mat_43.csv', index_col=0).T" @@ -59,7 +63,9 @@ { "cell_type": "code", "execution_count": 26, - "metadata": {}, + "metadata": { + "collapsed": true + }, "outputs": [], "source": [ "st_mob = pd.read_csv('../MouseOB/data/Rep11_MOB_0.csv', index_col=0)" @@ -68,7 +74,9 @@ { "cell_type": "code", "execution_count": 36, - "metadata": {}, + "metadata": { + "collapsed": true + }, "outputs": [], "source": [ "st_bc = pd.read_table('../BreastCancer/data/Layer2_BC_count_matrix-1.tsv', index_col=0)" @@ -193,6 +201,70 @@ "plt.savefig('Anscombe-mean-vs-var.pdf', bbox_inches='tight', dpi=200)" ] }, + { + "cell_type": "code", + "execution_count": 157, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [ + "def vst(expression_matrix):\n", + " from scipy import optimize\n", + " v = lambda mu, phi: mu + phi * mu ** 2\n", + " phi_hat, _ = optimize.curve_fit(v, expression_matrix.mean(1), expression_matrix.var(1))\n", + " print(phi_hat)\n", + "\n", + " return 2 * np.arcsinh(np.sqrt(phi_hat[0] * expression_matrix)) / np.sqrt(phi_hat[0])" + ] + }, + { + "cell_type": "code", + "execution_count": 158, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[ 0.47090384]\n", + "[ 0.31967197]\n", + "[ 0.21647288]\n", + "[ 1.73092133]\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAHHCAYAAACyWSKnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsvX+YHGd15/s505rBzAhj3HKCFzMzu5sfxBACRhAT+4LB\nYRNEgNxksws7NsIXrlAbFuUJ3GwS7eUJZJXdZNlcxAYDSmJiPBNIFkKWdUwIIeQHv+zI+Ac2Nokh\nGmGWxNYoYMsyWB6d+0d3j3t6qqrfqq7qqur+fp7nPJruqnrr7VbVt0+d97znNXdHCCGEEGKcmCq7\nA0IIIYQQeSMHRwghhBBjhxwcIYQQQowdcnCEEEIIMXbIwRFCCCHE2CEHRwghhBBjhxwcIYQQQowd\ncnBELGZ2xMweNrMdfe/fbGZuZotm9rudfU702K2d/RY7+3XfP2JmvxBxjof6jv9nPcdu6+x3npl9\n2MyOmdm3zOx2M3t133m29bX9u2b2nwr9koQQuWFmF5vZZzv3+HEz+4yZPXvINl9tZut9GvObnW2b\nNMLMXmNmd5nZA2b2j2Z2vZk9LmrfznuR2iOqgf5TxCD+Hngl8N8BzOwHgdm+fX7d3f9jQhtnufsj\nZrYT+Eszu8ndP9Gz/aXu/me9B5jZYl8b1wK3AgvAd4AfBJ6Y8rMIISqKmZ0JXAe0gD8AZoD/g/b9\nPiyfc/eLB5z/+cCvAj/u7jeb2dnAS3M4tygJRXDEIK4FXtXzejfw/iwNufth4A7gGRkOfzbwu+7+\noLs/4u43u/vHsvRDCFFJvg/A3T/g7uvu/pC7/6m73wZgZv+Xmd1pZv9kZh83s4XugWb2ok7k5Vtm\n9ptm9pdm9tqU5382bUfo5k4/jrv7Ne7+QF4fUIwWOThiEJ8HzjSzHzCzBvAKYDlLQ2Z2IfA04O6M\n/XiXmb3CzOaznF8IUWn+Flg3s2vM7MVm9oTuBjN7OfBLwE8B5wB/DXygs20H8IfAfwR2AF8BLspw\n/huAHzOzt5rZRWb2mKE+jSgdOTgihG4U50XAncDX+7a/2cy+2WPX9G0/ZmYPAZ8DrgL+qG/7H/Uc\n27+ty8/QFrX/F/h7M7slYmz+WG8/gH+X7mMKIcrC3e8HLgYc+C3gPjP7qJl9N7AX+M/ufqe7P0J7\nKOkZnSjOLuAOd/+Qu58C3gH8Q1/zF/Zp1IUR5/9r2g7UBcAfA2tm9hudB7sub+7TmNty/RJErsjB\nESFcS9tZeDXRw1Nvd/ezemx33/YdwHbgTcAlwHTf9p/sOfYnozrg7v/k7r/g7k8Fvhu4hbZjZL3n\n6e0H8HspP6cQokQ6Dsyr3f082tHef0bbYVkADvY4FscBA57U2edrPW147+sOn+/TqM/HnP9j7v5S\n4Gzg5bQ1r3eo6+19GvP0HD62KAg5OGIg7r5KO9l4F+1QcJY21t39N4BvA1cO2Z9jwNtpC9vZw7Ql\nhKgm7n4X8Lu0HZ2vAa/rc1Ie6+6fBb4BPLl7XOeh58lRbaY492l3/yTw553zixoiB0eE8hrghe7+\n4JDt/Bfg583sjDQHmdmvmdnTzGxbZ9pmC7jb3deG7I8QogKY2VPM7E1mdl7n9ZNpz+D8PPAe4BfN\n7KmdbY83s5/pHPrHwFPN7Kc607XfSIYZlmb28k6O3xOszXOA53fOL2qIHBwRhLt/pTMLKoqf76sx\ncSyhqT8G/gn4v1N2YRb4CPBN4Ku0Q9YvS9mGEKK6PAD8MHCDmT1I27G4HXiTu38E+DXgg2Z2f+f9\nF8NGRPdnaD88rQHfC3wmw/m7uvR3wP20J1P8V3dfGeZDifKw9nClEEIIMR6Y2V8Ay+7+22X3RZSH\nIjhCCCGEGDvk4AghhBBi7NAQlRBCCCHGDkVwhBBCCDF2yMERQgghxNhRu9XEd+zY4YuLi2V3QwgR\nwU033XTM3c8pux95IK0RopqE6kztHJzFxUUOH44rxyKEKBMzWy27D3khrRGimoTqjIaohBC1xszO\nMLMbzexWM7vDzN4asc+rzey+ziKtt5jZa6PaEkKMD7WL4AghRB/fob2MyAkzmwY+bWYfi1hQ8ffd\n/Q0l9E8IUQJycIQQtaazevSJzsvpjqn+hRATjoaohBC1x8waZnYLcC/wCXe/IWK3nzaz28zsQ52F\nHEUPKysrLC4uMjU1xeLiIisrWoJJ1Bs5OGLkSEhF3rj7urs/AzgPeI6ZPa1vl/8FLLr704FPANdE\ntWNme8zssJkdvu+++4rtdIVYWVlhz549rK6u4u6srq6yZ88e3Zui1tSukvHOnTtdMxvqS1dIT548\nufHe7Owshw4dYmlpqcSeiTwws5vcfWfJfXgLcNLd3x6zvQEcd/fHJ7UzSVqzuLjI6urWiSkLCwsc\nOXJk9B0SIoFQnVEER4yU/fv3b3JuAE6ePMn+/ftL6pGoO2Z2jpmd1fn7scCLgLv69jm35+XLgDtH\n18PqE+XcABw9enTEPREiP5RkLEZKnGBKSMUQnAtc04nMTAF/4O7XmdnbgMPu/lHgjWb2MuAR4Djw\n6tJ6WzFWVlYwM6Ki+fPz8yX0SIh8UARHjJQ4wZSQiqy4+23u/kx3f7q7P83d39Z5/y0d5wZ3/0V3\nf6q7/5C7v8Dd70putToUnbO2f//+SOfGzDhw4ECu5xJilMjBESPlwIEDzM7ObnpvdnZWQipqSdHO\nxyiSf+Oip+6uvDhRa+TgiJGytLTEoUOHWFhYwMxYWFhQgrGoJaNwPkaRsxYXPV1YWMjtHEKUgWZR\nCSFyowqzqPJikNaMYubR1NRU7PDR6dOnczmHZjaKuqFZVEIIUSCjSJjPK2ctaShNUVUxrsjBEUKI\nDIwiYT6PnLWQobSlpSWOHDnC6dOnOXLkSOHOjYp9ilEgB0cIITIwioT5PKIrVas9parJYlTIwRFC\niAyMamgnLroSGgXJaygtr6hL1RwuMb7IwRFCiIyMeminS5ooSB5DaXlGXeIcq9XV1VQOlIa5xEDc\nvVb2rGc9y4UQ1YR25eDSdSIPq7LWLCwsOLDFFhYWtuy7vLzss7Ozm/abnZ315eXlQs6XtS0zC+5j\nHp9J1JdQnVEERwghakaaYac8htLi1qqKex/iIyxRuUtRS0UkDVtpmEuEoLWohBCiZszPz0c6F3HD\nTktLS0MNnzUaDdbX1yPfj6K/tk53SKvbF2g7KUePHo39LJA+f0hr2oleFMERQoiakdf08dAclijn\npvt+1LGDIiz9uUvNZjOy/f7P2EVr2okQ5OAIIUTNGHbYKSpp+PLLL8fMIh2WpGUbVldXueKKK9ix\nY8eGs5Q2IhPHgw8+GOl4xTl4u3btUuKx2EBLNQghcmOSlmqoM0lOCDy6VAO0ozGrq6uReTJxJO27\nsLDAgQMHNjljcUtSdPePWvpiZWVl0zDXrl27uOaaa7TkxAQQqjNycIQQuSEHpx4kORRdms0mDz30\n0JahplCSnJx+xyPJ4Qpdd2sUa4OJaqC1qIQQQkTm2oTkqqytrWV2bqBdgiRuaOvkyZPs27dvo38n\nTpyIbSc0r0aJx6IfOThCCDGmxBXo27VrV2wCb140Gg0OHDiAmUVuX1tb48orr+SKK65gbW0tcp/+\nxOmkxGglHosthBTLyWLAk4FPAV8C7gD2RexjwDuBu4HbgAsGtVvl4ltCTDqo0F+lSCrQt7y8HLkt\nT5udnfVms5np2G4fuwwq7qfif5NDqM4UGcF5BHiTu58PXAi83szO79vnxcD3dmwP8O4C+yOEGEPM\n7Awzu9HMbjWzO8zsrRH7PMbMft/M7jazG8xscfQ9HT1JwzZLS0uJs6OimJ6eZtu28PJpwwxx9S99\nETf1/LLLLmNxcRFgJGuDifpQmIPj7t9w9y90/n4AuBN4Ut9uLwfe33HKPg+cZWbnFtUnIcRY8h3g\nhe7+Q8AzgB83swv79nkN8E/u/j3A/wf82oj7WAqDhm2iplvHsbCwwPve9z4e//jHp+rD8ePHY+vc\npCEpl6a3kGAZa4OJajKSHJzO09IzgRv6Nj0J+FrP63vY6gQJIUQsnQekbpbqdMf6p++8HLim8/eH\ngEstLjmk5vTmqZw4cYKZmZlN23vzWvrr6UxNDf5JOH78eKr+zM/Pc/DgwVTHdPvSm2czKJdGSzWI\nfgp3cMxsO/Bh4Gfd/f6Mbewxs8Nmdvi+++7Lt4NCiNpjZg0zuwW4F/iEu8c+TLn7I8C3gOHDCiUS\nlXC7srLCFVdcsZFUvLa2xiOPPEKz2YwdtumtKuwJU8e7UZKzzz47uI9mtuFMhThPXbo5FL2rlodE\nm/qniQ+74rhWLK85IYk6WY32k9THgZ+L2f5e4JU9r78MnJvU5jgk/gkxrlBykjFwFu3JDU/re/92\n4Lye118BdkQcvwc4DByen58v7osakriE2rm5uciE3WazGdRuXFJyf1uD9ulaq9XyVqu1ZaXwtNZd\ntXx5eXng+VutVuJ3FJp0rKTl6hKqM0UKjQHvB96RsM9LgI919r0QuHFQu3JwhKguZTs47S7wFuDN\nfe99HHhu5+9twDE6hU7jrMpaE+KI9FsIUT/q/WZm/pjHPGbg+RqNRi7OTfecoZ/dzHx5eTlxBtkw\n33Ho8aI4QnWmyCGqi4DLgRea2S0d22Vme81sb2ef64Gv0p4m/lvAlQX2RwgxhpjZOWZ2VufvxwIv\nAu7q2+2jwO7O3/8a+POOUNaOlZWVxGUWhqE3JyeO+fn5oJlU6+vrvPe97yXpa3aPLwbYf84ugwr3\nufvGEg5RhBb+U+HA+lPkLKpPu7u5+9Pd/Rkdu97d3+Pu7+ns4+7+enf/l+7+g+4+nnXRhRBFci7w\nKTO7Dfgb2jk415nZ28zsZZ19fgdomtndwM8Bv1BSX4eiW7gvjrg8l6RZTN08EzNj27ZtXHbZZQBc\neumlkfvPzc3x4IMPBvV30BILIbk1/cX+Qgr3ra6uDl34T4UDx4CQME+VrMphYyEmHSowRJWXVVFr\nkoZnZmdnvdVq+czMzKb3Z2ZmfHl5eWPYxsw2iuglDUnlMbQ0yLq5Qb19azab3mw2N/Wzl5BhtEaj\noRycMSZUZ0oXkbRWRdER9SFK5EV+yMEpliSno7eib4gjk5SUPErLci8m5dh0Le67SIP0opqE6oxW\nExcTQze831sNtX9VYzEcWk28WLKumJ20WnfRhA5pzczM8LjHPY7jx49z9tln88ADD/Dwww8D7aG3\n173udVx11VWbjtEK4pOJVhMXoo+4Uu8qDibqQlS+Sn+OShRlJcbOzs5yxhlnBO378MMPs7a2hnu7\nfk/XuYF2Ls+73/1utm/fvqkWTdbvQ0wGcnDExKBZEaLu9FceDl1vqazE2N27d6eufJzEgw8+uFH4\nD+K/D0AF+oQcHDE5aFaEGAd6Kw+HrrcUFenIY6WKQWtMXX/99bnfXydPnuRVr3oVZoaZsW/fPg4c\nOLDxfQDs2bNno5pzbzVkMVnIwRETg8LZYlKJinTs3bs3eKHNKBYWFjh27BiNRiN2n9XVVb7ne74n\n8zni6J1+vra2xhVXXLHhwGgoWnSRgyMmhqzhfSHGkYsuumhgUb8kug8G6+vrift98pOfzNR+Gk6d\nOrXhwGgoWnSRgyMmiizhfSHqRu8ikTt27GD79u1cdtllW4ZtAI4cOcL27dtTtT81NbVx78zNzeXe\n/yx0HRgNRYsug+ttCyGEqA395RDW1tYi9+sdtjlx4kSqc7zuda/b+Puhhx7K2NN86TowBw4ciCwH\noaHoyUMRHCGEGCOiclDiOHr0aKbclOuvv34j52XQcgyjouvAaChadJGDI4QQY0SaXJP5+flMBQBX\nV1e5/PLLufLKKxOTjEdFq9Xa5MD0D0WDpo1PInJwhBBijAjNNZmdnWXXrl2x08VnZmYSp5K7O+95\nz3u45JJLsnQzkdC8noWFBZaXl7dUOO6lO2SXZtp4bw6THKL6IgdHCCHGiEGrc0O7fs2hQ4e4/vrr\niVuuZ3p6OnZbF3fn7rvvptVqZe5vFO9973tptVqxDtbs7CzLy8tBEwXSThvP4hCJihKyYFWVrIoL\n4Akh2qDFNivB8vKyN5vNyIUtW63Wxn55rBhuZrm11bVms7mxmGaj0dj0b9wK43GLYib1K6qtuEU8\nFxYWivsPE6kI1ZnSRSSt1Vl0hBh35OBUi0GrYQ9akTvUGXH3SIcqL5udnY1dyTtupfTu/oM+Y3/b\ncQ5R15ET5ROqMxqiEkKIMWVQ3aeQ4awQVlZWuP/++4duJ46kIaVBQ1CDPmN/26qjMz7IwRFCiAml\nO6V6mJlQx48fZ//+/Zw6dSrHnm0lbYXi7vu908bjWF1d3Ugm1pIu44McHCGEmGCWlpa45pprmJ6e\nznS8mWWaaj4zM5Nq/7SRld73u5GsQU7OFVdcwb59+zh58uSG06c6OvVFDo4QotaY2ZPN7FNm9iUz\nu8PM9kXsc4mZfcvMbunYW8roa5VppzakJ2uhv9e85jWpZl/t2rUr8v00EZdBw1WnTp3aqPy8vr6+\n0Y6cm3pigy5qa8/TWwL+hbu/zczmgSe6+42j6GA/O3fu9MOHD5dxaiHEAMzsJnffmfHYTFpjZucC\n57r7F8zsccBNwE+6+5d69rkEeLO7/0RofyZJa3bs2BG7pENRmFkmp6rRaLC+vk6z2QTaQ2Rnn332\nxt/z8/OJTsnKygr79+8PjjotLCxsFAsU1SBUZ0IiOFcBzwVe2Xn9APCuIfomhBBRZNIad/+Gu3+h\n8/cDwJ3Ak4rq5DgyaucGskeMuquXr62tsba2hruztrbGQw89xLXXXjuwNk7IcFUvWoW8voQ4OD/s\n7q8Hvg3g7v8EpBs8FUKIwQytNWa2CDwTuCFi83PN7FYz+5iZPXXIvooYkqofF0n/bKhB1YhDZ5Bp\n9lR9CXFwTplZg3YtAMzsHKAaq6sJIcaJobTGzLYDHwZ+1t375yx/AVhw9x8C/jvwRzFt7DGzw2Z2\n+L777svyGWpJd7inHzPbsBCyDjvlRTfaElKNuH9RzmazuSXxWbOn6k2Ig/NO4CPAd5nZAeDTwK8W\n2ishxCSSWWvMbJq2c7Pi7n/Yv93d73f3E52/rwemzWxHxH6H3H2nu+8855xzhvgo1SB0TaWDBw9u\nmUU1PT3Ntddey+nTpzdyXAZRpnMDj0Zb4mrj7Nu3Of+8t07QsWPHuPrqq7UK+TgRUg0QeArweuAN\nwA+EHFOUjUN1USHGFYasZJxFawAD3g+8I2GfJ/LopIrnAEe7r+Os7lozqMJv1P5xVY9JqARMX7Xf\n0H3ztpmZmY0+J/Uj7vOL+hCqMyGzqC4E7vB28h5mdmZHeKLGuAtnkmY2CFE3hpxFlUlrzOxi4K+B\nL/LokNYvAfMA7v4eM3sD0AIeAR4Cfs7dP5vUbt21ZnFxMXKmUJZZQdu2bdtI7q0qc3NznDhxAoj/\n7JDt83dnXh09enTgLC1RPKE6E+Lg3Axc0PGaMLMp2t7TBbn0NCV1Fx0hxpkhHRxpTY5MTU1FDhmZ\nWeraNaNMHJ6bm+PBBx/MdGz3866srHDZZZdF7pP283fzeXqHvGZnZzV8VSJ5ThM377lL3P00sG2Y\nzgkhRATSmhzJc02l0CnVwzI9PZ3Lkg9LS0uxidNpP/+gta5EdQlxcL5qZm80s+mO7QO+WnTHhBAT\nh7QmR/JcUymqrbyjOgsLC5x55pk8/PDDkdtnZ2djnRbYOhPs4MGDuXz+tGtgieoQ4uDsBX4E+Dpw\nD/DDwJ5BB5nZ1WZ2r5ndHrNdpdOFEL1k0hoRTf806GFmBUW1tXfv3sj1pKanp2Mdkcc85jGR77da\nLY4cOcLx48dj+7B7924OHjwYec5Go8HBgwcH9nn37t3s379/Y1bZlVdeOXCWmVYXrzEhmchZDHge\ncAFwe8z2S4Dr0rZb95kNQowzDDmLqkomrRnM8vKyN5vNjRlKzWbTl5eXfXl52aempiJnOp1//vmb\n3rv00ks32ltYWIid/bSwsODu7q1Wa1Pbc3NzQTOjomaV9VvULLO0s9FE8YTqzMAIjpmdY2a/ZGaH\nOlGZq83s6gDH6a+AeHdcCCF6yKo1onji6uksLS1x7NixjR+UY8eObUSIohJ5H374Ye68885N733u\nc5/baC9p+Ojo0aOsrKxwzTXXbGrbB0yU6RKVS9NPVG5NnpEwMVpCZlF9lvYUzJuAjXmC7v7hgY23\ny6Zf5+5Pi9h2Ce3CXPcA/5v2Qnh3DGqz7jMbhBhnhpxFlVlrikBa0ybLLKKkadpR9E7djlv4s5vo\nnHXqe9yssn6yzDIToyXPWVSz7v4f3P0P3P3DXcuhj0Gl02Fyy6cLMWEUpTViCEJnEfVGedI4N7A5\nYTcpOXiYhN/QnBnl1owPIQ7OdWa2K+8Te2Dp9M72sSqfLoSIpBCtEcMR4lT0r/2Ulq5T0S2od/Lk\nSRqNBsCmIaEsCb9dx2t1dXXgzC+tPTVehDg4+2gLz0Nmdr+ZPWBm/QvZpcbMnmidq83MntPpy9a4\npBBiUihEa8RwDHIqVlZW2L1798D8ljjMjF27dm1ykgDW19c3HI7uUFjaqe/9bbr7hpOzsLBAq9Ua\nmFsTup6XqCAhmchZDPgA8A3gFO08m9fQnga6t7P9DcAdwK3A54EfCWlXMxuEqC5oFtXYkTSLKGRm\nUojNzs5umo3Va93ZU739iVszq5+4WVn9bcZ97qg+aQZVGGn+n9ISqjOhzsoTaC9Q97yuhRxXhEl0\nhKguwzo40ppqEvdjlTStG/BGozG082NmmfudddHNQY5biIM0yRQ9tT5UZ0JmUb2Wduj4POAW4ELg\nc+7+wsQDC0IzG4SoLkPOopLW1IykmUmzs7Ps3r2ba665JvPwFWRbHLNL0myupJlgg2aBaaZVMnku\n9BpFnrOo9gHPBlbd/QXAM4FvDtk/MUFoDFsEIq2pGXH5OY1Gg927d3P99ddvSRiOq3LcbDZzW1qi\nS1TOTpek9aQGzcrSTKtkqrK8RYiD8213/zaAmT3G3e8Cvr/YbolxoX92xerqKnv27JGTI6KQ1tSM\nuKTfPXv2cM0110QmDMdNAz948GDuBfW6RfriiPvBTXJgNNNqMJVZ3mLQGBbwEeAs4JeBvwL+J3B9\nyPhXEaZx8XoxTJKfqB8MkYMjraknUfk5g+77IhNQo4jrT6PRiOxDXA5OdykKkUxVcnDSCtDzgZcB\nM2mOy9MkOvUiLslvmMRBUV2GcXB6TVpTb6p232dZh2rUTti4UYVZVLFDVGZ2Zuffs7sGfBH4NLA9\nJDokRGVClaKySGvGj6rc9938v8svv5zHPvaxNJtNzGwjJ6iX/pycpaUljhw5wunTpzcSY5VLGE7/\n91fG2l1JOTi/1/n3JuBwxL9CDCRtYS4xkUhrxoxh7vu8JiX05/+tra3x0EMPce2118bOgIrLyYnK\nJbzsssvYsWOHHJ0qkxTeAQyYDwkFjcoUNq4fCvVODmQcohpGa4AnA58CvkS7eOi+mPbfCdwN3AZc\nMKhdac1wZLnv88zdSMoDSpsbmFTvR4X/Rk+ozoSIxxdDGhqVSXSEqC5ZHRwfQmuAc7sOC/A44G+B\n8/v22QV8rOPoXAjcMKhdac3oyXNSQlIeUFpHKqlgYNb+ZUEPi21CdSZkmvgXzOzZAfsJIcQwZNIa\nd/+Gu3+h8/cDwJ3Ak/p2eznw/o4+fh44y8zOHbrHIlfyrJ+SlAfUnT4eOiV9UO7QKOq7qORGekIc\nnB8GPmdmXzGz28zsi2Z2W9EdE0JMHENrjZkt0i4QeEPfpicBX+t5fQ9bnSBRMnkmJw/KA0qTBJtU\nMDBr/9LSXWW9l6RihXWg8CKwg0I8wEKUhYSHijCFjYWoLgw3RDWU1tCecXUT8FMR264DLu55/Ulg\nZ8R+e2gnNh+en58v+NsS/eRdPyXPIZ2iFt8M7WORU+/LGPoa5v86VGfSiM93AfNdCz0ub5ODI0R1\nGcbB6VoWrQGmgY8DPxez/b3AK3tefxk4N6lNaU05VD3PJG+nKfRHvqiiqUUX5YtjmM+Tm4NDu9jW\n3wEPAn8PnAbuCGm8CJPoCFFdhozgZNIa2onD7wfekbDPS9icZHzjoHalNflRdaelLNL8yBfliJRV\nbX6YiFSozoTk4PxKRxD+1t3/OXAp8PmA44QQIg1ZteYi4HLghWZ2S8d2mdleM9vb2ed64Ku0p4n/\nFnBl/t0XUSg5Np40SdVpE6OL6EOejKIYZIiDc8rd14ApM5ty908BA5cprxpa0VqIypNJa9z90+5u\n7v50d39Gx6539/e4+3s6+7i7v97d/6W7/6C7q4BgAUTp7Dgmx+ZF2h/5IqoDl1V1ehRFYEMcnG+a\n2Xbgr4EVMztIO4RcG/QEIUQtqL3W1IUiHvjidLa7ong/o5haXXWqUOm9rD4UFZHaRNzYFfAu4GJg\njrYjtA3YDbwRaIaMfxVhWcbFQ8cYNU4sxHCQIQdnnLSmDow6lyMu12JUxfGqThV+d6rQhzSE6kyS\n6OwDPgccAX4deGZIg0VbFtEJSWYqK5NciHEio4MzNlpTB4pKKk2q9js9PT3R2lo3B6LqhOpM7BCV\nux909+cCzwfWgKvN7C4ze4uZfV9iWKhihIwxapxYiHIYJ62pA0UllSblbJx55pnFDkVUmGFTJJQ/\nOgQhXlDXaFcIvRlYT3Ncnpb2qSq0OFORRZSEmBTIoQ6O11Rr6kKR9VSi2p10HU07Fbw30tNqtTSy\nEEGozoQIzTbgpcAK8A/AB4GXhzRehIWITvci6d5Y/RdWs9nccoGUVQtAiHFiGAenjlpTR6KG47s6\nOezwSdTD5DjpaJahptCH56T/l3H9PrMytIMDvAi4uiM0HwX+HTAX0miRNkh0oi6Sfms0GlsuUOXg\nCDE8WRycumpNnUl6CBx2aYRx1dGsny3U6Yt7yFZEbCt5ODh/DrwWeEJIQ6OyQaKT5iLpv0BDvPOs\nyWJKMhOTQEYHp5ZaMw4UEbkeV63L8l0tLy/7zMzMlmOmp6e3fC9JSdqK4GxmaAenqjZIdNJcJGkv\nlqwefNxxrVZrLIVATC5ZHJyq2iQ4OMo9DCfLdxXnFDWbzeB9o6zVauX62ermlE6sg5M2gpPmZs5a\nTycuRJmnDl3CAAAgAElEQVQUGk57wdXtAhXjiRyceqHcw3CyfFdpnKKycnDqOKw4sQ5OSA5O1osl\naz2dtH1Je8HV8QIV44kcnHoh7Qgny3eV1inqf1Ad9qE8hDo6uRPr4LhvvkiazeaWIlNZw30h4cYs\nEaT+CzftBVfHC1SMJ3Jw6oeiv+FkiawP40COQtvrOEw50Q5OLyEOR5ocnEEJY1lygPr7kvaCq+MF\nKsYTOThCbGYYB3IUEbY6PiCH6kzIYpu1JqQ6Z9w+/RUk4zh16hSXXXYZi4uLnH322Vm7urHAWdrV\nXctaDVYIIUQy3RXAr732WgAuv/zy4IrEo1iQsgoLfhZGiBeUxWjXtbgXuD1muwHvBO4GbgMuCGl3\nVBGcVqu1JTKybdu2gW3NzMwMHBKLs6w1eTSOLqoCiuAIsYUqaHR/Vf/egrd1G6YM1ZkiHZznARck\nODi7gI91HJ0LgRtC2s2yVENS0m//RRa3tEMaazabGxdLo9EIOqbRaGzpt2ZRibohB0eIrZQ9DLS8\nvBz54D0zM1PL34rSHZx2H1hMcHDeC7yy5/WXgXMHtZlFdPqTjpvNZqQjMOwMqK51c1/SOEt51zUQ\nogzk4AixlbLzJJNGMqqcaxNHqM6UmYPzJOBrPa/v6byXO90x0NOnT3Ps2DGOHTvG6dOnOXLkyKax\nzKgVxbMwPz+/sYLs2trapm1mxvnnn0+j0QCg0WjQarW46qqrhj6vEEKI6lF2nmRSLmp3W1mrlhd6\n3hAvKKuRHMG5Dri45/UngZ0x++4BDgOH5+fni3EJffgZUPQMeSV5zHFjrxpmEnWHEiI4DM73uwT4\nFnBLx94S0q4iOCIvys7BGRTBKat/Wc8bqjNlOjgjG6IKZdgaNlNTUxvDX4P2jap+3D8FPWR8NA+n\nSI6VyIuSHJxB+X6XANelbVcOjsiTMnV2UA5OWTlCWc9bBwfnJWxOMr4xpM0iRScuB6ebt5PksMzM\nzPjU1FSwM9Q/9hrXfjfTPerGyMPrLvvJQowXZTg4Plhr5OCIiSdpFlVZOUJZz1u6gwN8APgGcIp2\nfs1rgL3A3s52A94FfAX4IjHDU/1WtOjEORODojJRBQCTrNFobDpH0r5xDkgeXnfZ2f1ivKiwg7MG\n3Np5qHpqSJtycMSkEPI7UEQEqtYRnCKsLNEZdvgqjQOTh4VSdna/GC8q6uCcCWzv/L0L+LuEdkaS\n7ydElRgUyS8q0l90Ds7YVzLOiyKrOuYxc6uXbjZ6SFZ62dn9QhSNu9/v7ic6f18PTJvZjph9D7n7\nTnffec4554y0n0KUxaCKyVEzjE+ePMn+/fsLPe+wyMGJIWqZhlarNdI+TE9P02w2Ux93+vRpVldX\ncXdWV1fZs2dPrJMz1mW6hQDM7IlmZp2/n0Nb99aSjxJisugtp9JfQiVumnnIUkjDnHdY5OBE0K1h\n0+8kXHTRRSwvLzM1VfzXtrCwwPve9z4OHjy4xQFJS5KnXaQHXVZdBTFZmNkHgM8B329m95jZa8xs\nr5nt7ezyr4HbzexW2svDvKIT5hZCBFDbSH/IOFaVbBQ5OIMSn+JmPM3MzORSS6c//2VQXZ0sbRaN\nZmdNJqiSsRBjR9X0PFRnFMGJYFA47vjx45HbT506xd69eyO3paFbCbkb/di/fz8HDhxgeXk5MZrT\nbDZjh7RG7WnnPWaraJAQQpTDKFY1L4QQL6hKVoUITtL2PCItl156aWQkqNlseqvVil3As9FoeKvV\nqoSnnfT50lK1pwcRD4rgiBIZl6Kl4/I5iiJUZ0oXkbQ2CtEZZsrcsENUj3nMYxK3D5pSPjs7661W\nq9SbI+l7yFJfZ9S1eqLERYIThhwcURZ1eBAK0ZG4z1G2rlcJOThDMuhCjNueNoIzNze3scp52mKB\ncVZ2kb6478DMMt2Uo6zVEyUuMzMzW8qcV004q4IcHFEWVS9aGuqAJemnNKiNHJwR0uvspHVUuhdt\n3LBTVivTy0+KYvV/XyF9HKVwpXFQqyKcVUIOjiiLqhctDdWxNKMAk6pBoTqjJOMh6Z9Svra2hrvT\nbDYxs4FTytv/V7C+vp5rv9wH18ApgpWVldjPvLCwEDsFP6mPo6zVk6auQx41IIQQ+VD1qcyhtWTS\n9FcaNIAQL6hKVrWnqkFeedwqrqO0UC9/2NyTuMVKYfj1s0L7kbQwacjxiuAMB4rgiJKoeg5OqPZF\nfY48cxrHgVCdKV1E0lrVRCckLJpHHZthLCREG3VTTU9Pbxk6m56eTu0cNBqNDWdpmD5m+Qzd5LxQ\n4VMOznDIwRFlUuXJACEOWO9vRVd7FxYWKjM7tirIwRkRaSISSRGOLDY7OxtbdLDfms2mN5vN2MhG\nmvM2m81Nn6krKEnOy6DPnseTSJKDleacmkWVHTk4QsSTpCMhs3elQW3k4IyItGHRvKI5c3NzA6Mi\ng5yjqKeCUIv77HGORNJnnpmZyeVmTTtFvyrJh+OEHBwhslH1WWBVQg7OCEnrWYc4Oc1mMzF3Z3p6\neiMiMzU1lclJGaZmT6ijFlIfqDciNAx5RXBEduTgCJGNqs8CqxKhOqNZVDmQZjXU3llEcczOznLw\n4EHOPPPM2H1OnTq1MWPr9OnTmfrdvk6ysXv37sTP0F/OO2lmQO/SF8MsyRA322rPnj25zMLSchFC\niKKo+iywWhLiBVXJ6v5UlRRl6I8A5bFwZxkWl380aCZAHrMghp1FldSukvwGgyI4QmRCGhNOqM6U\nLiJpre6ikyYMGZpAXCVLuiFbrVZiNc4qj0FXuW9VQg6OENlRInEYoTqjIaoRM85hyEajsbFieNTw\nzVVXXcW1114buyJtXNGq1dXV0oeFQot0CSHqRZWGnkPTHarU50oT4gVVyer+VJUmDFnXIaqkz5RE\nmqTlUaMIThgogiNqRB2HherY57wJ1ZnSRSStjYPoFFFVN+9qyV3nKusMLUg/OypNnaBRh28lKmHI\nwRF1oswHl6zDUXn1uc7DYXJwxoA0dWZ6p21njfz0Jjq3Wq3YHKBGo7ExRX3QIqFZknkHFQ7sWl71\nc5L6kGei8iQgB0fUibKmZg/zwJRHn6vywJZVU+XgjAm9F8Dc3NyWizrqouw9Js0q5d0bZJBj1RuZ\nGeSIpKkW3E9IBCuvGjq9/arCjV9XynBwgKuBe4HbY7Yb8E7gbuA24IKQdidNayaRsiI4w5w3jz5X\nYch9GK2VgzOmpPV400Rzuhd3iGPRPf+gmV5RTxWhF3ZoBCtPqnDj15mSHJznARckODi7gI91HJ0L\ngRtC2p10rZkEynqgGSYKk0efq1BUcBitlYMj3D08j6f3BskzuTnqYk1zYYdUTM6TKtz4VaDo0HHe\nBiwmODjvBV7Z8/rLwLmD2pTWTAZlDD0P+yA1bJ+r8CA3jNbKwRHunrzCdtwNkia5OdRp6iXLhR0X\nKcp7iKoKN37ZjCJ0nLcNcHCuAy7uef1JYOegNqU1oijKHgov+/zuiuBIdDLS790nOTNxx2dZhLPZ\nbAadJ8uFvby8vGWm2PT0dO43ZBVu/LIZhfDkbXk5OMAe4DBweH5+PvN3KMQgyp60UIXzKwdHDk4q\n8vqBjrr4B0V2QodxsvZxVDdk2Td+2YwidJy3DXBwNEQlJpYq65lmUUl0UlHkEMugVcHTnKPKN92k\nM4YRnJewOcn4xpA2pTWi7oxrRDpUZwpdqsHMftzMvmxmd5vZL0Rsf7WZ3Wdmt3TstUX2ZxIockmB\npaUl9u7dG7ltZmYm1ercaVZgF6MlblX2tKuvjwoz+wDwOeD7zeweM3uNme01s+7Fej3wVdrTxH8L\nuLKkrgoxUvbv38/Jkyc3vdddTmcS2FZUw2bWAN4FvAi4B/gbM/uou3+pb9ffd/c3FNWPSWN+fp7V\n1dXI9/Pgqquu4qKLLmLfvn2sra0B0Gw2OXjwoJyUMaH7/7h//36OHj3K/Pw8Bw4cqOz/r7u/csB2\nB14/ou4IURkmfQ29IiM4zwHudvevuvvDwAeBlxd4PsFonr6XlpY4duzYRhjw2LFjlf3xE9lQhE2I\n+jPOizuHUKSD8yTgaz2v7+m8189Pm9ltZvYhM3tygf2ZCJaWljh06FDsit1CCCEmg7oNN+dNoTk4\nAfwvYNHdnw58Argmaicz22Nmh83s8H333TfSDtYRPX0LIYSY9AfewnJwgK8DvRGZ8zrvbeDuaz0v\nfxv49aiG3P0QcAhg586dnm83hRBCiPFkaWlpYhyafoqM4PwN8L1m9s/NbAZ4BfDR3h3M7Nyely8D\n7iywP0IIIYSYEAqL4Lj7I2b2BuDjQAO42t3vMLO30Z7D/lHgjWb2MuAR4Djw6qL6I4QQQojJocgh\nKtz9eto1KHrfe0vP378I/GKRfRBCCCHE5GHtEhH1wczuA7YWeimGHcCxEZ0rb+rcd1D/y2SYvi+4\n+zl5dqYshtSaOv//JzGunwv02epEkM7UzsEZJWZ22N13lt2PLNS576D+l0md+14VxvU7HNfPBfps\n40jZ08SFEEIIIXJHDo4QQgghxg45OMkcKrsDQ1DnvoP6XyZ17ntVGNfvcFw/F+izjR3KwRFCCCHE\n2KEIjhBCCCHGDjk4fZjZGWZ2o5ndamZ3mNlby+5TFsysYWY3m9l1ZfclLWZ2xMy+aGa3mNnhsvuT\nBjM7q7Nw7F1mdqeZPbfsPoViZt/f+c67dr+Z/WzZ/aoL46IdSdRZV5Kos+YkUWc9yoNCC/3VlO8A\nL3T3E2Y2DXzazD7m7p8vu2Mp2Ud76Yszy+5IRl7g7nWs23AQ+BN3/9edJUpmBx1QFdz9y8AzoP1D\nRnvtuI+U2ql6MS7akUTddSWJumpOErXVozxQBKcPb3Oi83K6Y7VKVDKz84CX0F7AVIwIM3s88Dzg\ndwDc/WF3/2a5vcrMpcBX3H1URTVrzzhoRxLSlXoxZnqUCTk4EXTCsLcA9wKfcPcbyu5TSt4B/Dxw\nuuyOZMSBPzWzm8xsT9mdScE/B+4D3tcJ4/+2mc2V3amMvAL4QNmdqBtjoB1J1F1Xkqir5iQxTnqU\nCTk4Ebj7urs/AzgPeI6ZPa3sPoViZj8B3OvuN5XdlyG42N0vAF4MvN7Mnld2hwLZBlwAvNvdnwk8\nCPxCuV1KTyeU/TLgf5Tdl7pRZ+1IYkx0JYm6ak4SY6FHwyAHJ4FOOO9TwI+X3ZcUXAS8zMyOAB8E\nXmhmy+V2KR3u/vXOv/fSzgF5Trk9CuYe4J6ep/YP0RaYuvFi4Avu/o9ld6Su1FQ7kqi9riRRY81J\nYlz0KDNycPows3PM7KzO348FXgTcVW6vwnH3X3T389x9kfYww5+7+2UldysYM5szs8d1/wb+FXB7\nub0Kw93/AfiamX1/561LgS+V2KWsvBINT6Wm7tqRRN11JYk6a04SY6RHmdEsqq2cC1zTmUUyBfyB\nu4/VlMiK893AR8wM2tfn77n7n5TbpVT8e2ClM8zzVeCKkvuTio7Avwh4Xdl9qSHSjnpSd81JotZ6\nNCyqZCyEEEKIsUNDVEIIIYQYO+TgCCGEEGLskIMjhBBCiLFDDo4QQgghxg45OEIIIYQYO+TgiNww\nM+8t/mVm28zsvnFbeVgIUS7SGhGCHByRJw8CT+sUOYN2PZWvl9gfIcR4Iq0RA5GDI/LmetorDkNf\nRdxOxdCrzezGzuJvL++8v2hmf21mX+jYj3Tev8TM/sLMPmRmd5nZinWqcQkhJh5pjUhEDo7Imw8C\nrzCzM4CnA72rKe+nXeL9OcALgP/aqZx7L/CizmJ3/xZ4Z88xzwR+Fjgf+Be018QRQghpjUhESzWI\nXHH328xskfYT1fV9m/8V7QX73tx5fQYwD/xv4DfN7BnAOvB9Pcfc6O73AJjZLcAi8Omi+i+EqAfS\nGjEIOTiiCD4KvB24BGj2vG/AT7v7l3t3NrNfBv4R+CHaUcVv92z+Ts/f6+iaFUI8irRGxKIhKlEE\nVwNvdfcv9r3/ceDfd8e2zeyZnfcfD3zD3U8DlwONkfVUCFFnpDUiFjk4Infc/R53f2fEpl8BpoHb\nzOyOzmuAq4DdZnYr8BTaMySEECIRaY1IQquJCyGEEGLsUARHCCGEEGOHHBwhhBBCjB1ycIQQQggx\ndsjBEUIIIcTYIQdHCCGEEGOHHBwhhBBCjB1ycIQQQggxdsjBEUIIIcTYIQdHCCGEEGOHHBwhhBBC\njB1ycIQQQggxdsjBEUIIIcTYIQdHFI6ZfcbMnlnSuf+bmbXKOLcQonzy1h8zu9HMnppXe6I45OCI\nLZjZxWb2WTP7lpkd7wjEs83sl8zsRMe+bWbrPa/viGnrpcAD7n5z5/Uvm5mb2b6+/fZ13v/lnvfO\nMrN3m9k/mNlJM/uimV3Rd9wRM3uo04d/MrM/NrMn9+zyduCXzGwmr+9HCFEcI9CfUz3H3WlmP913\nzJlm9g4zO9rZ5yud1zs6u7wdeFuhX4LIBTk4YhNmdiZwHfDfgbOBJwFvBb7j7r/q7tvdfTuwF/hc\n97W7xz3R7AWu7Xvvb4FX9b23u/N+tx8zwJ8BC8BzgccD/w/wX8zs5/qOfWmnT+cC/9jpOwDu/g3g\nLuBlIZ9fCFEeI9Kf3+9p52eBZTP77s75Z4BPAk8Ffhw4k7b+rAHP6Rz/UeAFZvbEfD61KAo5OKKf\n7wNw9w+4+7q7P+Tuf+rut6VtqCMWLwT+sm/T3wCz3TBv598zOu93uRyYB37G3f/e3U+5+58AbwTe\n1hHCTbj7t4EPAef3bfoL4CVp+y+EGDmj0J8N3P3jwAPAv+y89SrauvN/uvuX3P20u9/r7r/i7td3\njvk2cBPwY2n7JEaLHBzRz98C62Z2jZm92MyeMERb3wucdvd7IrZdy6NRnN1sfcp6EfAxd3+w7/0P\n03aGntvfoJnNAv8W+HzfpjuBH0rXdSFECYxKf7A2LwFmgC913v5R4E/c/cSAtqUpNUAOjtiEu98P\nXAw48FvAfWb20W4INyVn0X46imIZeKWZTQOv6LzuZQfwjYj+PQIc62zv8kdm9k3gW7Qdo//ad9gD\nnb4IISrMiPTn33T04gTt4aZfdfdvdrY1idCdCKQpNUAOjtiCu9/p7q929/OApwH/DHhHhqb+CXhc\nzDmOAncDvwr8nbt/rW+XY7RzajZhZttoOzfHet7+SXc/i3Zk5w3AX/aNjz8O+CZCiMozAv35A3c/\ny93naA9NvcrMXtfZtkaE7kQgTakBcnBEIu5+F/C7tIUmLXfTjgQ/KWb7+4E3df7t58+AF5vZXN/7\nPw18h63DUHTG7P8QWKf9FNjlB4BbU/ZdCFEyBesP7n4E+Bjw0s5bfwb8WITu9CNNqQFycMQmzOwp\nZvYmMzuv8/rJwCuJcCgG4e4P0xaM58fs8vvAvwL+IGLbtcA9wP8ws0UzmzazHwPeCfyyu38rou9m\nZi8HnkB7jLzL82mLmBCiwoxYf+ic58eB7jTza4GvAR/u9GXKzJqdKeq7OsecATwL+ETaPonRIgdH\n9PMA8MPADWb2IG1huZ12pCUL76U9I2oLnRkSf+buD0Vs+w7thL+vATcA9wO/Aex39/4cm/9lZic6\n+xwAdrv7HQBmdi7tWVV/lLH/QojRMQr9+bfdOji0Z25+hvZU9F7duYu2A3M/cCPtYfEbOse/FPgL\nd//fGfskRoS5e9l9EGOOmX0GeEO32NaIz/3fgK+4+1WjPrcQonzy1h8zuwF4jbvfnkd7ojjk4Agh\nhBBi7NhWdgegXW6fdmhyHXjE3XeW2yMhhBBC1JlKODgdXuDuxwbvJoQQQgiRjJKMhRBCCDF2VMXB\nceBPzewmM9tTdmeEEEIIUW+qMkR1sbt/3cy+C/iEmd3l7n/V3dhxevYAzM3NPespT3lKWf0UQiRw\n0003HXP3c8ruRx7s2LHDFxcXy+6GEKKPUJ2phIPj7l/v/HuvmX2E9rL0f9Wz/RBwCGDnzp1++PDh\nUvophEjGzFbL7kNeLC4uIq0RonqE6kzpQ1RmNmdmj+v+TbuyreoLCCGEECIzVYjgfDfwETODdn9+\nz93/pNwuCSGEEKLOlB7BcfevuvsPdeyp7n6g7D6JwaysrLC4uMjU1BSLi4usrKyU3SUhhBBig9Id\nHFE/VlZW2LNnD6urq7g7q6ur7NmzR06OEEJkRA+N+SMHR6Rm//79nDx5ctN7J0+eZP/+/SX1SAgh\n6oseGotBDo5IzdGjR1O9L4QQIh49NBaDHByRmvn5+VTvCyGEiEcPjcUgB0ek5sCBA8zOzm56b3Z2\nlgMHlB8uhBBp0UNjMcjBEalZWlri0KFDLCwsYGYsLCxw6NAhlpaWyu6aEEJUhtDEYT00FkMV6uCI\nGrK0tCSHRgghYugmDndza7qJw8AW7ey+3r9/P0ePHmV+fp4DBw5IY4dEERwRiaYsCiFEdtImDi8t\nLXHkyBFOnz7NkSNHMjs30u5HkYMjtqApi6JOmNmTzexTZvYlM7vDzPZF7GNm9k4zu9vMbjOzC8ro\nq5gcykgclnZvRg6O2EJVpizqSUQE8gjwJnc/H7gQeL2Znd+3z4uB7+3YHuDdo+2imCRWVlaYmor+\neS0ycbgq2l0V5OCILWR58sjbGdGTiAjF3b/h7l/o/P0AcCfwpL7dXg6839t8HjjLzM4dcVfFBNDV\nrvX19S3bik4c1nTzzcjBEVtIO2WxCGdETyIiC2a2CDwTuKFv05OAr/W8voetTpAQQxOlXQCNRqPw\n2aaabr4ZOThiC2mnLBbhjOhJRKTFzLYDHwZ+1t3vz9jGHjM7bGaH77vvvnw7KCaCOI06ffp04bOi\nNN18M3JwxBbS1rkpwhnRk4hIg5lN03ZuVtz9DyN2+Trw5J7X53Xe24S7H3L3ne6+85xzzimms2Ks\nKVO7VKNsM3JwRCRppiwWcUPrSUSEYmYG/A5wp7v/RsxuHwVe1ZlNdSHwLXf/xsg6KcaKpJzDorVr\nUL5jXtPNxwJ3r5U961nPclEtlpeXfXZ21oENm52d9eXl5aHbXVhYcDPzhYWFodsTxQMc9hFrAnBx\n57q7DbilY7uAvcDezj4GvAv4CvBFYOegdqU1IooQvStKu4rS2roRqjOlOyxpTaJTTUbhjMjhqT5l\nODhFmbRGRLGwsLDJwejawsLCWJ+7SoTqjIaoRC4UHRbVtHEhRBUocwLEoHOrdthm5OCIWqBp40KI\nKlBmEnHSufUQuBU5OKIWaNq4EKIKlDkBIuncegjcSiUcHDNrmNnNZnZd2X0R1UTTxoUQVaDMqdhJ\n59ZD4FYq4eAA+2iXVxciEk0bF0JUhTKnYsedWw+BWyndwTGz84CXAL9ddl9EdVEBKyGEiEcPgVsp\n3cEB3gH8PHC67I6IaqMCVkKIsqj6DCU9BG6lVAfHzH4CuNfdbxqwn9aHEUIIsYU8HY+4tuoyQ0kP\ngZuxds2ckk5u9p+By4FHgDOAM4E/dPfL4o7ZuXOnHz58eEQ9FEKkwcxucvedZfcjD6Q11afrePTO\nHpqdnc0UuUhqa//+/ayurm45ZmFhgSNHjmTuv8hGqM6U6uD0YmaXAG92959I2k+iI0R1kYMjRsni\n4mJujkdSW0ePHiXut7Iqv6GTRKjOVCEHRwghhEhNlEPS+36a4aukadZxM5HMrHLDVOJRKuPguPtf\nDIreiPyoesKcEEIMotFoxL6fNm8maZr1gQMHaC9avxl3Z/fu3dLPilIZB0eMjjwT5uQoCSHKYn19\nPfb9tJV9k6ZZLy0txQ5Fra+vVzLhWKDVxOtOlhW281qRdnl52WdnZze1MTs7q1W+Jxi0mrgYIUla\nZmaR24BYvUzS07hzZdVPkZ1QnSldRNKaROdRsjoYcTe+maU6f16Okhgf5OCIfrI8hKVpO04DBzkk\naR/Ios41jH6K7ITqjIaoKk7SEFDWxdXyKumttU+EEEmMon7MYx/72I2/m83mxhTxqCGnftIsRtkt\npBeX9zPJSyJUlhAvqEo2SU9VgyI0WSMxeQ0tKYIj+kERHNFDkRoRomO90aOofnRt2PMC3mw2NTw/\nIkJ1pnQRSWuTJDqDxGEY8cgjbKwcHNGPHBzRS17D4VEM0r9+jWs2m7F9WV5eTqWJy8vLke1J/0aD\nHJwxYJA4VMHBiBOFIsfdRXWRgyN6KTKCMygq0799eno61f6DtFQR7PII1Rnl4FSYQbkyZSyu1p8T\nBGxZ+6Qu67YIIYYnKU8wKg9mZmaGEydOBJeW6LZvZmzbtg0zY3FxkbPPPjvxuPbv4KOcOnUq1f6D\n8nOUg1gDQrygKtkkPVVVIUKTpT91fLJRxCkfUARnYggdpum9t5rN5pZISjdyEjdtO27m0vT0tM/M\nzCRGZYa1pKG0OurcuBCqM6WLSFqbNNGp0g9v6A1d5Lh7EVTNkawzcnAmg0FTpuN+5AdN3e6/7wbt\n32w2g5KIs1qSsyLdKA85OCJ3Qh2Xuj3Z1K2/VUYOzmQwyPGIe5gJcUR677uQ/UP7lNZCnJUqPYBO\nEqE6oxycGtI/5n3llVdmWi4h7TILofVzkkqeVxGNpdcbM7vazO41s9tjtl9iZt8ys1s69pZR93Hc\nGHRvDFNra3V1dUOL0tSWCal7M4juelOh+YxLS0tbchBFhQjxgqpkk/5UNSg0TIonj7Th1TTH1OnJ\nRhGc/KCECA7wPOAC4PaY7ZcA16Vtd9K1JomkaEmSjoToF7Tza7r5PYOiOP05PnHTwQfZqHUqr1Id\nddHZPAnVmdIdlrQ2jqKT5iINDcMO+nHO+qM+jjdUnOPWarXG7rMWTRkOTvu0LMrBGR3DFLsLXUYh\nxMws9qGr1WoFtzM1NTXS+z2P/J1JzgGSg1MT0l6kocl0gxJ665YIXDT9jlur1ZpY8RiGCjs4a8Ct\nwMeAp4a0OW5akzfDPuwsLy/71NRULo5OXDRm27ZtmdvI+37v/b4ajUamB8xeJjnyLAenJqS9SEOf\nfHpnF0SJzyTfHCHo+8lGRR2cM4Htnb93AX+X0M4e4DBweH5+vrDvadIJHaoaxrZv3z50G41GIxcn\nJ9rBpr4AACAASURBVPTzpnnAnOSH1FCdUZJxyayuria+358IvGvXroGJdDMzM9x///2Rhfa67UWd\nd1AicNqk5DqjxOPxwd3vd/cTnb+vB6bNbEfMvofcfae77zznnHNG2s9JImqh4Lw5ceLE0G2sr6+z\nZ8+ezBM5uoR+3jRJ1XktmjzWhHhBVbJxiuAsLy8PjML0F7KanZ31Sy+9NPGYuCS7ZrMZ+xRhZt5q\ntWL7GdVmUoGuuqMITjaoZgTniYB1/n4OcLT7OsnGSWuKJu1wVVF1a4qytMs4ZPm8ysEJJ1RnSndY\n0tq4iE5RIdphil5F/XiH9nPcbqxJFo9hKMPBAT4AfAM4BdwDvAbYC+ztbH8DcAftHJzPAz8S0u64\naE3RDLpXepOKu7kncTkodbI88mUajcZQOUy9D56TtJq5HJyKk3dRqq6Zmc/NzWU+dph+NpvNEr7J\n4hjHGWNFU4aDU5SNi9YUTVK0cxS5NmVZmlyXvB+YJv0BrDYODnAGcCPtp6s7gLcm7T8uolNUiDap\nBsSgGQtRTyRp+zlJTxFiK3JwJo9BmlKmEzI3N7dldmTv66w1c+L0Mok8H5gmfQi9Tg6O8egMh2ng\nBuDCuP3HRXSKuPEbjcZAJybpaSoqBydLPyfpSSKUSYkGycGZPKo83DQoypJmdlMWjSvqvp/kGVTu\nJTg4HUflMuAtndfzwHNStjELfAH44bh96i46vePR/RfpzMxMYXUhug5Q0vaoIaasIea8pleOA5MU\nTh7GwclDQ/K0umvNqCjbiUmykIhGrxMyqK1eR2WQ81Lkfa8IzugdnHcD7wLu7Lx+AvA3gcc2gFuA\nE8CvJe1bZ9FptVpbbqLemUitVmvLrKlR26BlF9I4YFlu5nGMdEySGA3p4GTWkCKszlozSsoehorT\npJmZmURHJOr9uCGr7vbu/lEzXPtnlRZ530/SQ1MUZTg4X+j8e3PPe7embOMs4FPA0/rer23xrZCy\n5N0LvmyhgHbkpTtO3X3de9OmKX+e9mYe15t2ksLJocITZXloSJ4mByeZ5eXlzBMaRmHNZjNSr6an\npyMrlU9PT0dGubdt2xa5f5Il7ZvXfT+OD4OhhOpMng7ODbQjMV2ROqdXqFK08xbgzXHb6yQ6aYZ3\nquDcDLLeBfBCLepmjrsxxzXSMa6fK4pQ4YmyvDQkL6uT1hRJ1P26vLxc6dybEF0q69zjeN+PmlCd\nydPBWQI+SrsOxQHgy8DPBBx3DnBW5+/HAn8N/ETc/nUSnTo4LaO6mZNyj7pRmqR26sy4RqaiCBWe\nKMuqIUVZnbSmKOKu3VEOpWfNS6yiAzau9/2oCdWZXAUBeArwetqFtX4g8JinAzcDtwG300kwjLM6\niU7dqnXmbWa28bQ3KJK1sLCQKEh1D8NOSjg5VHjiLIuGFGV10pqiKPshrash7u7nn39+6ZqWxYYp\n5ieiCdWZPJ2bC4HH9bw+k4TZUFmtTqJTtjgUZWnGokNykNKeu8oiMSmOTByhwhNlo9KQUKuT1uRJ\n6KyiUZm7p879G4VNTU1tDNknfVfjmGtXNqE6k6eDczM967sAU3TG0vO0OonOOFbxnJqa2jLbICmE\nXNRyFFVkkoai4ggVnigblYaEWp20Ji/GUbOKsu4sre73FheBrqpe1ZlQnclzNXHrnBgAdz8NbMux\n/dqxtLTEoUOHaDabZXclN06fPs1nPvMZjh8/vum9OE6ePEmj0ci1D1Vd0TtqxeCTJ0+yf//+knpU\nO6QhJbKyssLu3bsLX+V7XHj44YfZt28f0Nb6a665htnZ2U37zM7OcuDAgVTtrqysDLVyuXiUPB2c\nr5rZG81sumP7gK/m2H7tWFlZYf/+/Rw/fpxms0mz2cTMaDabzM3Nld29TDSbTd7znvfQ8zs0kPX1\ndaanp3Prw/z8fOL2sgQizvGqqkNWQaQhJbGyssKePXtYX18vuyuVoqvbcaytrWFmLC4uAnDo0CEW\nFhYwMxYWFjh06BBLS0vB51tZWeGKK65gdXUVd2d1dZUrrrhCTk5WQsI8IQZ8F/BB4F7gH4HfA74r\nr/a7VpewcVyot1uboeyCflkty4yGhSHXfOm13qTD0O99VMNEkzQdPA6GG6IaiYaEWl20ZliqmN9S\nBeu9b0PyCPPQmTidHLeFjIclVGdKEY5hrMqi05ucV8UpimWYmeUqoFHrZfWSxsnIOyE4L+eqzonK\nocJTB6uy1gxD7/VV5UJ9eVqr1Ur1kNWbX9P9zkKOG2bG1DiXysibUJ3JTQxo17P5JeAQcHXX8mq/\na1UVHSXnRdv555+f6/cy6P8g7rj+mQxFRXqGdU7qnqgcKjxRNioNCbWqas0wTKJOzc3Npf7sZrbl\nHt6+fXuq86a5b0P6Jh4lVGfydHA+C/wa8G+An+5aXu13raqiM65Twoe1tJGspOmWSUM9gwQiquBg\nmnOMKqoS17e61NIIFZ4oG5WGhFpVtWYYJlGn+iMxaaNWs7Oz3mq1fHp6OvW5u/frIO0Y9P+iIarN\nhOpMng7OLXm1lWRVEp2q1YsYd0v6YU8SiO4wWchTUtzSEqOKqoRcS1WO6IQKT5SNSkNCrUpakxeT\nqFW9LC8vZ3JUhkk56NeOmZkZbzabmxyepP+X/uEyUY6D85+AXXm1F2dVEZ1JDPWWaUlPMCHj47Oz\ns0Fj8FERnFEmD4c+YVc1cTlUeKJsVBoSalXRmjyZxAhOb9Qkr8kOoRbiGCVpU6PRkHMTQajO5Ong\nPACcBh4C7u+8vj+v9rtWFdGZRKEoy5IiFnk7mr2FuwaJYREVSkM/T1WrozKcgzMSDQm1qmhNnoQm\ny46bhax3V6Y1m83YKHGdJx0URajOlP6UlNaqIjqTGOodtYXc0Hk6mtu3b3f38DB2UVGUkNl44xjB\nqZpVRWvypuz7uiwLeWBpNpullPDolr/od2TqPumgKEJ1JldBAJ4APAd4XtfybN8rIjpJZbll+Vh3\n5sMg8nQ0u8Ngo6p5EcLy8vIWwa3ymHyo8MTZKDQk1KqgNXnRarWkWQH3vns5+h73wKLaWtGE6kxu\nlYzN7LXAXwEfB97a+feX82q/Kqji52h48MEH2b59e2Q14t5KxVNT0ZewmaU+Z3f5iUGVh7NUKB2G\n9v0c/3pcmBQNGTVXXnkl7373u6VZMczMzHDw4MGN10tLS4nLz+RN0nIOqo4+JCFeUIgBXwTOoDMT\nAngK8Id5td+1sp+qlHtTjvWORw/KUelO60xbdLH7VJRlGvmor7eqPsExXA7OSDQk1MrWmmEZVBJB\n1rYs07bzsmazmWkIvqr3/6gI1Zk8HZy/6fx7C/CYzt935NV+18oWHeXelGcLAUs+bN++fWMKZrPZ\nDJ410TvkFJeDU8bQUNz1NqZJxpk0hHZBwHuB22O2G/BO4G7gNuCCkP6UrTXDUKVZnr334zBthNzD\nac+RVPcqy/c3PT2dKodnkKOiHJxoQnUmTwfnI8BZtEPKfwX8T+D6vNrvWpmio9yb8bRunZz+/+te\nseyuwTXqWQx1e4JjOAcnk4bQztW5gHgHZxfwMdqOzoXADSH9qauDUyWd6s1tyRIV6f0xTzp+bm5u\nI8Ib6mAMmp0Zdf8nWW9icGj0OORBZdxnUWX5fKE6k6vzsdEoPB94GTCTd9tliU6Vnohk+VuaJ7lR\nPkGVff60kNMsKlJqCLBIvIPzXuCVPa+/DJw7qM06OjhV0qn+iGfa6Hf/j13SZ+uPwA6K5CTVl4k6\nT9KsykE1ukIqrIf+346bo5NV30J1Jg8hOrPz79lRNmz7/VaW6Ggse7wt7klqmAhKXoIUN320imJH\nBgeHHDSEZAfnOuDintefBHYOarOODk5VdCrqmkzbtyiSolP992SaoeZhFkruakfcPRnncPXnCibd\nx1GOQFTkOeq4UepE2vNl1ddQncnDwbmu8+/fA1/t/3fY9vutLNEpWzBkxVrcDZV0TBJFRl6qHNUh\nm4MztIaQk4MD7AEOA4fn5+eL/rpyp+z7COIfFtJEl5J+4NLkpfU7F1FJvcNGvZrNprdarS396n89\nNTW1sQ5W14nq3yfuPo5zBLr1c0K/7yJ1Isv5suYYhupMLk4H7bHt+TzaGmRlOThVGdOW5W9JIhH3\n/95oNBKvlyJzZ6qcl0PGISqG1BAmdIiq/4m57Htp0HXY399Wq5X6RzGP63+YiE2UFuS5f9TnyPJ9\nj1onspyv8hGcjYbgixmPezLwKeBLwB3AvqT9FcGRFWFZ/t+TGHb2U1Kot8ozqxhymvgQxy4S7+C8\nhM1JxjeGtFl1B6dK+TZdyxIhSDusMWxkoorfW5K2DFqMM+6+H7VOZDlf5XNwNhqCa4BnZzjuXDrT\nNoHHAX8LnB+3f1miM+pF2mSjs6QE46xLJcRdL6G5O0k3/ThGcHw4DfkA8A3gFHAP8BpgL7C3s92A\ndwFfoV1rZ2D+jdfAwalKxKb3+htl8n3W3JKqfW/91h8dHtTfOkdw3Gsyiwq4C3ikIyK3dYTktgzt\n/E/gRXHb5eDI8rS4p4XQGRtRDFtHZ5BQjFsOjnu+GpKXVd3BGXVNrksvvbSy110a6lDLLLS/aRci\nrloOTlZCdSZPB2chylK2sQgcpTOrIsrKEp063BSyMOv+XyY9LcQ5GXHTS0PG9JOmk4Zca70L8nX7\nMuhzjBqGc3AWoixre8NalR2cMmrddM87ylk5RZB3BCdtcb+0OThptaiXqs+iykqozuQuCsB3AfNd\nS3HcduAm4KcitpU+s6Gb/S6rv83NzQ2cnpl2pkbImH7o2HecoDWbzco/QZNDHRwyakjeVlUHp4wc\nkksvvbTsj50bab+/focjpGxDV19622k0Gt5qtTY9pITMogqJjIyD45mGUJ3J07F5GfB3wIO0p3ee\nJnCpBmCa9sJ6Pzdo3zJEZ3l5ufQfZVmx1l8oLE0UJvSJMHTsO07QhsnrGRUMF8HJrCFFWFUdnFHn\nkIyTc9MlyiGJq1VTpLMQ6pgk7ZfH0FBIP6rkRIXqTJ4Ozq1AE7i58/oFwO8EHGfA+4F3hJxn1KLT\narVK//GVjc5CQ829w0Mh++chOFWePdWF4RycTBpSlFXVwRnFcHnocOq4UaUf8VCGTSYOjRBVKXoc\nqjN5OjiH/VGRmur+HXDcxZ0v7Dbai+zdAuyK23+UoiPnRjaMdZ2gpFyZNII66lkRWWA4ByeThhRl\nVXVwio7gVG3Ycxwo0nEa9sEnRFeqpj2hOpOng/NntPNofpP2tM2DwGfzar9roxSdsn8gZfW1bhn2\npKeetE9FVXuKioLhHJyRaEioVdHBiSv7n5dFVfoVyQxyXoq+b4d1PkIcpKpFj0N1Jg/H5l20ozBz\nwBSwDdgNvBFoDtt+v41KdJR3I8tq27dv35RI2G9d4ckiTCGl58uEbEs1jFRDQq1qDk6RycVJ1bxF\nPCHOS9HRj2EdKEVwksVpH/A54Ajw68Azh20zyUYlOqp7I8tq3R+LQU89o6z8OSrI5uCMVENCrWoO\nTlFDU9u2bavM9TNK8hg2CvnhH0X0Y5jPohycMJFaAP4DcDPtgl1vAb4vr/a7NirRKftHUlZvW1hY\nyBzBaTabsWJVtSepfhi+Dk7hGhJqVXJwiooo1yWRNm/y+sEOcV6qfs+6axZVWqF6Zkek1vNuWw6O\nrA7WjeKkzcGJKhrWe0zVxsL7IYc6OF6whoRaVRyc5eXlVIXkQm2SycvpCGmnatGPcSBUZ/J0arYB\nLwVWgH8APgi8PK/2uzYK0VH+jWxY6wpcSAJi7/ZBtW6q/jTIcBGckWhIqFXFwSliaKoq10tZ5PWg\nEOq8pIl+lBkpqVKUJolQncnDsXkRcHVHkD4K/Dtgbth242wUoqP8G1moRQlldwZVFqFIqnHS/aEL\nqX5aFmTLwRmphoRaVRycvK/ZKl0vZZHng0KeTkGZ0Z46RZpCdSYPB+fPgdcCTxi2rRArWnQUvZGl\ntenpaW82m5uqomYVijjh7XdqQtbTiru+i3xCI5uDM1INCbWqODh5rjlV5afyUVLVH/OkvLwo0kaI\ny6izVYTmhOpM6SKS1ooUHTk3sqzWKwJxQtFoNIJKsvcLb1xUJ63wjELUySkHpwpWFQcnj+tzUisT\nJ1HF4ZikCG7vQ9Ty8nKmHL+k+72I/L6iNCdUZ0oXkbRWpOhoaEqWZEnXR0hRrF4bVNCvV3hDzhlC\n2ie0LD8AyMHJjaRSA2mtKknoIpnQfKuQtenS3u9FDNvF9X/YqFCozpQuImmtCNEZ9J8hm2xLO+Uz\n9Frqn2kR50yECs8gh2TYFdJDnryQg5MLy8vxC75msUlPKq4LeYwiZK2zlVe0JaQg5bAOd6jOlC4i\naS1v0SmyOqhsPKzfEclSFCvKujd5HqHmvCuqZn2aQw5OLmzfvj2367cKuSUinGFHErJGcNyLLX6Y\nRkcGEaozpYtIWstbdBS5kSVZVDJvaFGsQW2nEaJB5wxtI/QJLet4PHJwhiavRX6rlFsiwgl9QGo2\nm7nm4OTFoGFV5eCMUHTyGuOWjb+lvTGTnOe8i/eFthH6hKYITjkOTl4THeTU1Jve+7TZbMYW/wyZ\nRTXqtesG5d5oFtUIRWcUEZz+p38lM9fX0oRW457E+kUmj+S+PBME4/quHJziyUOPWq3WyPstiiXr\n0FEZUZwqzdYsXUTSWt6ic+mllxbyQ9hdUToJDY9V05KSO7NUOg0ZzhpWEIoQFc2iGr2DM+y1q+ng\nope8H3xCqUq9rdJFJK3lLTp5/jCmDf8lDY9NTU353Nycm5mfccYZpf/oT4p1HZhRC0MeglCFuh6h\nwlMHG7WDM+zDlpKJRT9VX7suK6E6M8UEs7Kykks7y8vLuDvHjh1jaWkp+Lj5+fnI9xcWFlhfX+fE\niROcPn2ahx56iFarRaPRAKDRaNBqtWi1Wrn0XzzK2WefDcCBAweYnZ3dtG12dpYDBw5EHreyssLi\n4iJTU1MsLi6mvraWlpY4cuQIp0+f5siRI0HXUf85gdRtiGrwoz/6o3zyk5/MfPzCwgKHDh3S/7nY\nRNxvzP/f3r0HR1VneQD/nn6EBHQghAZJApMIBMgDUJFyZpwpNbCFj8D62HJ8oJQCVnQGZ0vL0kK3\nHLXYsdRRqQGqlHVHIYrKTsnTZRXZWqzaHYwuzwBClEeQR1giEBI66fTZP7rDdEI6uZ2+nXv73u+n\n6hbdt2/fPt10Tp/7u7/7+8Vb7zhGqiA7LWYeVSXbFyYjI8PyeUc4hk/PSyKfj8/nuzhiaE5OziWj\nh8b7P/D7/R324/f70/48d2/AohYcANMB7ANwAMDTXTw+G0A9gG3RZU5P++yrFpxkOhYXFxf3SYyU\nnuyaJ5JlNM9YXrAkupiZdJL50TSrI59ZpxWsLiLsurR/pr0d68hIMohXKKeyP4RV59Z7YjTxmLkA\n8AKoBXAlgAwA2wEUd9pmNoA/JbLfvipwBgwY0OvvN1FPkumgbPUp73iM5hnLC5ZEFzsUOOXl5abF\nYJbuWimSSaDpvrQXosmMDGvg0ug+/wGy67l1o4nHzAXAzwBsjLn/DIBnOm0zGzYtcJL9bhOZze4t\nP0bzjOV9cETkHRE5KSK7rI7FiOLiYnz++edWh3GJrvqMiAgqKyvR2Nh48T+8vR+PW2zYsAFApI9L\nOBzu1T4OHz588XayfW3M4vpz6x3lATgSc78uuq6zO0Vkh4isEpERfRNa96ZOndqr51VWVmLJkiUm\nR0MUsWDBAjQ1NXVY19TUhMcff9wW+c8wI1VQKhcAvwJwNYBdRrY366jK4/EkfMRk9/PdRpoUE33P\n6b4YmUeqp6W9BSfeUU28YfVTeYrKrkdYsKYF5y4Ay2Luz0Kn1hoAOQD6RW8/AuCLOPuaB6AaQPXI\nkSNT+En1vu+NHVuQyVmMDoBrVc4xmmcsL3AisaIAfVjg9HYodCdwW4fknqYq8Pv9HToVd+4sHPsH\nHO+z62qk0WQ7oBthx3PkRhOPmQsMnKLqtL0XwJme9pvKU1S9zUGDBg1KWUzkbrH5JJHT+Vb0+2OB\n0/2Hk/DilAG0jB41ZmRkXPJj39Pi9/sv+aFPdulqvhWjS1dHF0aGNo/3eHf9XuxYbFjBogLHB+A7\nAIX4Wyfjkk7bDI+5fTuA/+lpv6kqcJK5aoooFZK5EMOKfn+OKnBgcrNxov+BXq/XUT9Y8a768Xq9\nHX6gjVyC3nkySiNTUSQyU3Js8ZDI/1kqigy7XrlkJ1YUOJGXxS0AvkXkaqoF0XUvAJgRvf3PAHZH\ni5/NAMb1tM9UFTi9HbiTKFXi5bbY34R4eZ0tOEkWOLFLskmnN0dPTipuVBPvv9Hdlz/ec7r7Y0hk\ngtPOp5iMPDeVow3bsd+LnVhV4KRiMbvASWYeOn7HKJWMXJVpp/zHAif+B5PQ4tRz3omcUunNF7u7\n58QrmDr/kXX1GpWVld0WOX0xkRxPRcXHAqdrK1as6PXp29zcXNPiIOqK0dZpu+S/tClwAHwA4BiA\nVkQu73y4u+2TSTq9ab2hiN58seM9J17xU1lZaeg1YvdrdLRh6hsscLqWTOd+olSzU+uMEWlT4CS6\nJJN0Eh3wjgNppY5djgTIXCxw4n4uvVqYg6ivpFNONppnJLJt+pg8ebJWV1f36rkiYnjb3NxcHD16\ntFevQ+RWIvK1qk62Og4zJJNrOksk97QrLy+35aCiRFYzmmcsH8m4ryQ64iKLGyIyw6OPPprwcwYN\nGsTihihJrilw7r//fsPbrlixIoWREJGbLF26NOHnNDQ0pCASIndxTYGTiPvuu8/qEIjIAXozV0+6\ndRsgsitXFDiJNBFXVlamMBIicpNEWo4Bth4TmckVBU4iTcScoZeIzFBSUpLQ9rm5uWw9JjKRKwoc\no9h6Q0RmqKqqQk1NTULP4YUNROZyfIGTnZ1teFu23hCRGWbNmpXQ9ux3Q2Q+xxc4P/74o6HtiouL\nUxwJEblBVVVVQgVLVlZWCqMhci/HFzhG7d692+oQiCjNVVVVJdyxuKmpKUXRELmbowucvLw8Q9uV\nl5enOBIicoPZs2cntD1PTRGljqMLnB9++MHQdhwxlIjMEAqFDG+bm5ubwkiIyNEFjhEej+s/AiIy\nwdSpUxPanldNEaWWY3/djY5B0dbWluJIiMgNNm3aZHhbnpoiSj3HFjiJjkFBRNQXWNwQ9Q3HFjhG\ncFh0IjKDiFgdAhF14uoCh8OiE1GyEpnrjq03RH3HkQWOkaMpTstARGYwOtcdixuivuXIAscITstA\nRMnq37+/1SEQURyOK3CqqqqsDoGIXKK5udnQdmy9Iep7tihwRGS6iOwTkQMi8nQy+zIyTDqTDZGz\n9JRDRKSfiHwYffyvIlJgwmsmuwsiSiHLCxwR8QJYDOBmAMUA7hERznxJRIYYzCEPA2hQ1dEAXgfw\ncl/FxwMqImtYXuAAmALggKp+p6otAFYCmGlxTESUPozkkJkA3o3eXgWgXPqgCYbFDZF17FDg5AE4\nEnO/LrouJZhwiBzHSA65uI2qhgCcAZDTJ9ERkSXsUOD0SETmiUi1iFTX19dbHQ4ROZSZuYYHU0TW\nskOBcxTAiJj7+dF1F6nqW6o6WVUnBwKBXr8QEw6RI/WYQ2K3EREfgIEA/q/zjszKNURkPTsUOF8B\nGCMihSKSAeDXANZYHBMRpQ8jOWQNgAejt+8C8IUmecTT3dN5MEVkPZ/VAahqSER+A2AjAC+Ad1R1\ndxL76/LyTSYcImeKl0NE5AUA1aq6BsC/AFguIgcAnEakCDLjtc3YDRGlgOUFDgCo6gYAG0zcn1m7\nIqI00FUOUdV/irl9AcA/9HVcRGQdO5yiIiIiIjIVCxwiIiJyHBY4RERE5DiSbv1VRKQewCGDmw8B\ncCqF4STLzvHZOTaA8SUjlbH9VFUdcX21w3JNKrjxPQPufN92e8+G8kzaFTiJEJFqVZ1sdRzx2Dk+\nO8cGML5k2Dm2dOXGz9SN7xlw5/tO1/fMU1RERETkOCxwiIiIyHGcXuC8ZXUAPbBzfHaODWB8ybBz\nbOnKjZ+pG98z4M73nZbv2dF9cIiIiMidnN6CQ0RERC7kyAJHRKaLyD4ROSAiT1sdTywRGSEim0Wk\nRkR2i8jjVsfUFRHxisj/isg6q2PpTEQGicgqEdkrIntE5GdWx9RORP4x+v+6S0Q+EJFMi+N5R0RO\nisiumHWDReQzEdkf/TfbyhjTmZ1zTap09Z1yunTJ22YTkUwR2Soi26Pv+/dWx5QIxxU4IuIFsBjA\nzQCKAdwjIsXWRtVBCMATqloM4DoAj9ksvnaPA9hjdRBxvAng31V1HICJsEmcIpIHYD6AyapaisjE\nj6ZM6piEPwOY3mnd0wA2qeoYAJui9ylBaZBrUuXPuPQ75XTpkrfNFgRwk6pOBDAJwHQRuc7imAxz\nXIEDYAqAA6r6naq2AFgJYKbFMV2kqsdU9Zvo7XOI/DjnWRtVRyKSD+BWAMusjqUzERkI4FeIzA4N\nVW1R1R+tjaoDH4AsEfEB6A/gByuDUdX/QmT27FgzAbwbvf0ugL/v06Ccw9a5JlXifKccLR3ydipo\nRGP0rj+6pE3HXScWOHkAjsTcr4NNv4giUgDgKgB/tTaSS7wB4CkAYasD6UIhgHoA/xo9hbZMRAZY\nHRQAqOpRAK8COAzgGIAzqvof1kbVpWGqeix6+ziAYVYGk8bSJteQeWyct1Mi2l1hG4CTAD5T1bR5\n304scNKCiFwG4N8A/E5Vz1odTzsRuQ3ASVX92upY4vABuBrAUlW9CsB52OQUS7Qvy0xEirBcAANE\n5H5ro+qeRi6jTJsjMiIr2TVvp5KqtqnqJAD5AKaISKnVMRnlxALnKIARMffzo+tsQ0T8iPyRVKnq\nX6yOp5NfAJghIgcRaXK/SURWWBtSB3UA6mKOIlYhUvDYwVQA36tqvaq2AvgLgJ9bHFNXTojIcACI\n/nvS4njSle1zDZnH5nk75aJdATYjjfpfObHA+QrAGBEpFJEMRDp5rrE4potERBDpP7JHVf9oCUUH\nmAAACTtJREFUdTydqeozqpqvqgWIfHZfqKptWiFU9TiAIyIyNrqqHECNhSHFOgzgOhHpH/1/LodN\nOkB3sgbAg9HbDwJYbWEs6czWuYbMY/e8nSoiEhCRQdHbWQCmAdhrbVTGOa7AUdUQgN8A2IjIj8tH\nqrrb2qg6+AWAWYi0jGyLLrdYHVSa+S2AKhHZgUjP/oUWxwMAiLYqrQLwDYCdiPx9WToCqIh8AOC/\nAYwVkToReRjAHwBME5H9iLQ6/cHKGNNVGuSalIjznXI6t+bt4QA2R3PtV4j0wbHd0CHxcCRjIiIi\nchzHteAQERERscAhIiIix2GBQ0RERI7DAoeIiIgchwUOEREROQ4LHDKNiGjsoIAi4hORejvOSE5E\n6Yu5hoxggUNmOg+gNDogFBAZFIojuxKR2ZhrqEcscMhsGxCZiRwA7gHwQfsDIjJARN4Rka3RiTJn\nRtcXiMgWEfkmuvw8uv4GEflPEVklIntFpCo6oigREXMNdYsFDpltJYBfi0gmgAnoOOPuAkSmfpgC\n4EYAr0RnAj8JYJqqXg3gbgCLYp5zFYDfASgGcCUiI4oSETHXULd8VgdAzqKqO0SkAJEjqg2dHv47\nRCbyfDJ6PxPASAA/APiTiEwC0AagKOY5W1W1DgBEZBuAAgBfpip+IkoPzDXUExY4lAprALwK4AYA\nOTHrBcCdqrovdmMReR7ACQATEWlVvBDzcDDmdhv4nSWiv2Guobh4iopS4R0Av1fVnZ3WbwTw2/Zz\n2yJyVXT9QADHVDWMyIR23j6LlIjSGXMNxcUCh0ynqnWquqiLh14E4AewQ0R2R+8DwBIAD4rIdgDj\nELlCgoioW8w11B3OJk5ERESOwxYcIiIichwWOEREROQ4LHCIiIjIcVjgEBERkeOwwCEiIiLHYYFD\nREREjsMCh4iIiByHBQ4RERE5DgscIiIichwWOEREROQ4LHCIiIjIcVjgEBERkeOwwCEiIiLH8Vkd\nABH1va+//nqoz+dbBqAU7j7QCQPYFQqF5lxzzTUnrQ6GiMzDAofIhXw+37IrrrhifCAQaPB4PGp1\nPFYJh8NSX19ffPz48WUAZlgdDxGZx81HbkRuVhoIBM66ubgBAI/Ho4FA4AwiLVlE5CAscIjcyeP2\n4qZd9HNgLiRyGP5RE5Eh3377bcbWrVuzWltbTdvnkSNHfBUVFYX5+fllJSUl4ydNmjTuvffeG2Ta\nCxCRa7HAISIAQCgUwvr16y97//33B9bX13vb19fW1vpLS0vHT5gwofTGG28cFwgEJi1fvjzpIiQc\nDqOiomL0L3/5y8a6urqdu3fv3vPRRx99d+TIkYxk901ExAKHiPDVV19lDh8+fMLdd989Zt68eYX5\n+fkTX3zxxaHhcBjTpk0r2rt3b1YwGJSmpibPmTNnvPPmzSusrq7ODIVCeOGFF4aOGDGiNCcnZ+Jd\nd9310++//95v5DXXrl17ud/v16eeeqq+fV1RUVHLggULToZCITzyyCP5paWl44uKiopfeeWVIQCw\nbt26y6dMmTJ2+vTpVxYWFpbMmDGjMBwOAwC2bNnS/9prrx1bUlIy/vrrrx9z6NAhPwC89NJLQ0eN\nGlVSVFRUfNttt12Zgo+PiGyIV1ERuVwoFMItt9xSdOrUqQ6FycKFC/MGDBjQduLEiYy2tjaJfayl\npcXz+uuvD21ubvZ8+umn2RcuXPAAwCeffJKzefPmQTU1NbuGDRvW1t3r7ty5M2vChAlNXT32xhtv\nDBk4cGDbrl279jQ3N8u11147rqKi4iwA7NmzJ2vbtm3fFRQUtF5zzTXjPvvss8tuuOGG8/Pnzx+5\nfv36A7m5uaG33347+8knn8z7+OOPDy5atOiKQ4cO7czKytJTp055u3o9InIeFjhELrdx48bLmpub\nL/nhDwaDng8//DCnq87I4XAY+/fv77d9+/bLW1paLhY/bW1tcv78ec+bb745ZOHChScSiWPWrFkj\nt27depnf79f8/Pzg3r17+69ZsyYbAM6dO+etqanJzMjI0LKysvOjRo1qBYCSkpKm2trajMGDB4f2\n79+fddNNNxW1xxcIBFoBYOzYsc2333574YwZM3687777fkzs0yGidMUCh8jlzpw54wVwSRGjqvB6\nvdrS0nLJqezMzMxwQUFBS01NTbilpaVDcRQMBj1ffvnlTwB0W+CUlZU1r169Orv9/vLlyw8fO3bM\nN3ny5PF5eXktr7322uE777zzbOxz1q1bd3m/fv0uxur1ehEKhURVZfTo0c3btm3b2/l1Nm/evP/T\nTz+9fPXq1QNfffXV4fv27dvt9xs6i0ZEaYx9cIhcbtq0aY2tra2X5IKsrKzwPffcc/qBBx44mZmZ\nGW5f7/f7NScnp3Xu3LmnwuGwdH6ez+fT0aNHN/f0uhUVFeeCwaC8/PLLgfZ1jY2NnmhMZ5YuXRoI\nBoMCADt27Oh39uzZuPlqwoQJF06fPu37/PPPBwBAMBiU6urqzLa2NtTW1mZUVFScW7x48dHGxkZv\ntKAjIodjCw6RywUCgbZnn322buHChXnBYNCjqsjKygqPGTOmee7cuaczMjJ0ypQp55csWTLs7Nmz\n3ltvvbXhueeeOxEIBNpGjRp1Yd++fVmtra0XCx2fz6dPPPFEfXevCQAejwdr166tfeyxx0YsWrTo\nisGDB4f69+/f9vzzz9c99NBDDQcPHuxXVlY2XlVl8ODBrRs2bKiNt6/MzExduXJl7fz580eeO3fO\n29bWJpWVlSfKysqC9957b+G5c+e8qipz5sw5OWTIkG77BhGRM4gqx/oicpvt27cfnDhx4qnYdZs2\nbRqwePHiQENDg++OO+5omDt37unMzMxuE8SJEye89957b+GWLVt+IiIYOnRoy1tvvXXw5ptvbkzt\nOzDX9u3bh0ycOLHA6jiIyDwscIhcqKsCJxkNDQ2exsZGT15eXsjjSb8z3yxwiJyHp6iIKGnZ2dnh\n7OzscM9bEhH1jfQ71CIiIiLqAQscIncKd3UFlBtFPwe2PhE5DAscInfaVV9fP9DtRU44HJb6+vqB\nAHZZHQsRmYt9cIhcKBQKzTl+/Piy48ePl8LdBzphALtCodAcqwMhInPxKioiIiJyHDcfuREREZFD\nscAhIiIix2GBQ0RERI7DAoeIiIgchwUOEREROc7/A+XmL7xhmSSgAAAAAElFTkSuQmCC\n", + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "i = 0\n", + "for dname in data:\n", + " i += 1\n", + " \n", + " plt.subplot(2, 2, i)\n", + "\n", + " ds = data[dname]\n", + " sds = vst(ds.T).T\n", + " \n", + " plt.scatter(sds.mean(), sds.var(), c='k', label='Genes', rasterized=True)\n", + " \n", + " plt.title(dname);\n", + " plt.xlabel('Mean')\n", + " plt.ylabel('Variance');\n", + " \n", + "plt.legend(scatterpoints=3, bbox_to_anchor=(0.0, -0.25))\n", + "\n", + "plt.tight_layout()" + ] + }, { "cell_type": "code", "execution_count": null,