Skip to content

Instantly share code, notes, and snippets.

@vals
Created April 16, 2018 19:04
Show Gist options
  • Save vals/316eeaf9b8c796f8c32037488c1519a6 to your computer and use it in GitHub Desktop.
Save vals/316eeaf9b8c796f8c32037488c1519a6 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"%pylab inline\n",
"import pandas as pd"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"klein = pd.read_csv('Klein_RNA_control.csv', index_col=0)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"gene_info = pd.DataFrame({'mean_': klein.mean(1), 'var_': klein.var(1)})"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAECCAYAAAD6oXArAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvFvnyVgAAG/JJREFUeJzt3V2MXPWZ5/Hvz233DLbZkJidVdwGN1kjJ82MdjNYcDHSCMnOrNnBIcpGM3hoBwaLxt1y5L3acdRzsytZzF7sSkG8eJyNMdAnJhaKRk7iiBVIiL3gwk00ozF4PfRYadMGycZ5UXhRWMOzF13FVBdd3VV1zqn6V/XvI7XU59/n5amiqMf/d0UEZmZmVau6HYCZmaXFicHMzBZwYjAzswWcGMzMbAEnBjMzW8CJwczMFnBiMDOzBZwYzMxsgY4lBkl3SPo/kg5LuqNTzzUzs9bkSgySjkq6JOlMXflOSeckzUg6WCkO4F3gd4G5PM81M7PyKM+SGJL+mPkv+6cj4vcrZQPAPwFfYT4BnAZ2A/83Ij6W9G+A/xkR9+YN3szMiperxhARLwO/qCu+DZiJiPMR8SHwLHB3RHxc+fsvgd/J81wzMyvP6hLuOQS8WXM8B9wu6evAfwCuAx5tdLGkMWAMYN26dbd+8YtfLCFEM7P+9eqrr74TEf+63evLSAyLiogfAj9s4rwjkt4Gdn3+85+/dXp6uvzgzMz6iKTZPNeXMSrpInBDzfGmSpmZmfWAMhLDaeBmSTdJGgTuAU6W8BwzMytB3uGqx4FXgK2S5iTtjYirwH7geeAscCIiXmvlvhHxo4gY+8xnPpMnPDMza0OuPoaI2N2g/BRwqt37StoF7NqyZUu7tzAzszYluSSGawxmZt2TZGKQtEvSkV//+tfdDsXMbMVJMjG4xmBm1j1JJgbXGMzMuifJxOAag5lZ9ySZGMzMrHuSTAxuSjIz654kE4ObkszMuifJxGBmZt3jxGBmZgskmRjcx2Bm1j1JJgb3MZiZdU+SicHMzLrHicHMzBboaGKQtE7StKS7OvlcMzNrXt6Neo5KuiTpTF35TknnJM1IOljzp78CTuR5ppmZlStvjeEYsLO2QNIA8BhwJzAC7JY0IukrwOvApZzPNDOzEuXdwe1lScN1xbcBMxFxHkDSs8DdwHpgHfPJ4gNJpyLi4zzPNzOz4uVKDA0MAW/WHM8Bt0fEfgBJ9wPvNEoKksaAMYAbb7yxhPDMzGwpZSSGJUXEsWX+fkTS28CuwcHBWzsTlZmZVZUxKukicEPN8aZKmZmZ9YAyEsNp4GZJN0kaBO4BTpbwHDMzK0He4arHgVeArZLmJO2NiKvAfuB54CxwIiJea+W+XhLDzPpVlmUMDw+zatUqhoeHybKs2yF9St5RSbsblJ8CTrV7X0m7gF1btmxp9xZmZsnJsoyxsTHef/99AGZnZxkbGwPg3nvv7WZoCygiuh1DQ9u2bYvp6eluh2FmVojh4WFmZ2c/Vb5582Z+/vOfF/YcSa9GxLZ2r09yrSQvu21m/ejChQstlXdLkonBfQxm1o8azc1Kbc5WkonBNQYz60eHDh1i7dq1C8rWrl3LoUOHuhTR4pJMDK4xmFk/uvfeezly5AibN29GEps3b+bIkSNJdTyDO5/NzPqOO5/NzBLRC3MUmpFkYnBTkpn1momJCfbs2cPs7CwR8ckchV5MDkkmBjOzXpJlGYcPH6a+af79999ncnKyS1G1z4nBzKxF9U1GBw4c+FRSqEptjkIzOr7stplZL1tsWYulpDZHoRlJ1hjc+WxmqZqcnPwkKSxHUnJzFJqRZGJw57OZparZpiFJ7Nu3L7k5Cs1IMjGYmaXqc5/73LLnbN68mWeeeYbHH3+8AxEVr2N9DJK+BBwArgdejIgnOvVsM7NOkVToSqndkHejnqOSLkk6U1e+U9I5STOSDgJExNmI2Af8GfBHeZ5rZtZJWZZx/fXXI4krV64seW4vdjbXy9uUdAzYWVsgaQB4DLgTGAF2Sxqp/O2rwE/IsYmPmVknTUxMMDo6umxCgDQXxGtHrsQQES8Dv6grvg2YiYjzEfEh8Cxwd+X8kxFxJ9B7vTFmtuJUJ641Y8OGDUkuiNeOMjqfh4A3a47ngCFJd0h6RNLfskSNQdKYpGlJ05cvXy4hPDOzpVUnsI2OjjacuFa1efNmpqameOedd/oiKUAHO58j4iXgpSbOOyLpbWDX4ODgrWXHZWZWa2JigieeaG5sTNFbcqaijBrDReCGmuNNlTIzsyRVawiSmk4Ka9as6Yv+hMWUkRhOAzdLuknSIHAPcLKVG3iCm5l1SnWJi+WWtqi1bt06nnzyyb5pOqqXd7jqceAVYKukOUl7I+IqsB94HjgLnIiI11q8r5fEMLOOaGWJi2p/wrvvvtu3SQFy9jFExO4G5afwkFQz6wHNLnGxYcOGvuxPWEySS2K4KcnMyjYxMcHq1auXHXUEMDAwwHe+850ORJWGJJfdlrQL2LVly5Zuh2JmfeiWW27h9ddfb/r8p556qq+bjuq5xmBmfWuxPZhbSQqSmJqaWlFJAUDNVKM6rabG8OAbb7zR7XDMrAfVb6jTjl5NCpJejYht7V7vGoOZ9aVWRhstZnx8vCeTQhGSTAxmZnnl2Wt5+/btPbuXQhGcGMysLzWzoU69VatWMTU1xQsvvFBCRL0jycTgCW5mlkeWZU0tk13v6aefXrHNR7WSTAzuYzCzdu3YsYPR0dGWrpG0ovsU6iWZGMzM2rFjxw5efPHFlq7ZsGFDT+/PXIYkJ7iZmbVqYmKi5aQwMjLCa6+1tJTbipBkjcF9DGbWiizLml4uu+q6665zUmggycTgPgYza1a7fQq//OUvS4qo9yWZGMzMlpNlGatWrWq5+QjgmWeeKSGi/tGxxCDpa5K+K+kHkv6kU881s/6TZVlT+zHXW79+fc8uc9FJeTfqOSrpkqQzdeU7JZ2TNCPpIEBE/F1EPAjsA/48z3PNbGVrtemoOnHtN7/5jZNCE/LWGI4BO2sLJA0AjwF3AiPAbkkjNaf8deXvZmZNy7KM9evXI6ml6zZs2OCJay3Ku4Pby5KG64pvA2Yi4jyApGeBuyWdBf4G+GlE/CzPc81sZZmYmGh51BHMb8W5UnZdK1IZfQxDwJs1x3OVsm8BO4BvSNrX6GJJY5KmJU1fvny5hPDMrJe0mxTWrFnDoUOHSoio/3VsgltEPAI80sR5RyS9DewaHBy8tfzIzCxV7cxPqHryySfdfNSmMhLDReCGmuNNlTIzs6a1W1O45pprcu3DYOU0JZ0GbpZ0k6RB4B7gZCs38AQ3s5WlfgvOtWvXtjWTOSKcFAqQd7jqceAVYKukOUl7I+IqsB94HjgLnIiIluade0kMs5VjYmKCPXv2MDs7S0QwOzvLBx980NI9Nm7c6JnMBco7Kml3g/JTwKk89zaz/pdlGYcPH255olotT1grnvL8Bynbtm3bYnp6utthmFlJhoeHmZ2dbft6r466OEmvRsS2dq9Pcq0kNyWZ9b8sy5wUEpVkYnDns1l/a2dF1Krt27cTEU4KJfJGPWZWqizLOHDgwCd7MK9evZqrV6+2da+Um777SZKJQdIuYNeWLVu6HYqZ5ZBlGXv27Fnwhd5uUhgfHy8qLFuGO5/NrDTXXnst7777bu77bN++nRdeeKGAiFaGvux8NrP+UERSmJqaclLosCQTg0clmfWm2hnM69evz3Wv6kxmz1HovCQTg0clmfWeLMsYGxv7ZAbze++91/a9PJO5u5JMDGbWew4cOFDIOkXj4+NcvOh1N7spyVFJZtZbJiYmPhmO2i53MKcjyRqD+xjMekeePROqRkZGnBQSkmRicB+DWW+YmJhoewZz1fbt2z2LOTFuSjKzpmVZxuTkJBcuXMg9C9kb6qQryRqDmaWhdvjptddey+jo6CejjvLYvn27k0LCOpYYJH1B0vckPdepZ5pZ++qHnxYxWQ08Ya0X5N3B7aikS5LO1JXvlHRO0oykgwARcT4i9uZ5npl1zuTkZKH/qp+amvKEtR6Rt8ZwDNhZWyBpAHgMuBMYAXZLGsn5HDMrWf2+y3n2SqjnhNBbciWGiHgZ+EVd8W3ATKWG8CHwLHB3s/eUNCZpWtL05cuX84RnZk2qbzYqKils3LjRS2X3oDL6GIaAN2uO54AhSRskHQa+LOnbjS6OiCPAfwV+Njg4WEJ4ZlavqFnLtSLCM5h7VMc6nyPiSkTsi4h/GxEPL3Ou5zGYdUiWZblnLdeqLn5nvauMxHARuKHmeFOlrGme+WzWOQ899FBh95qamvLid32gjMRwGrhZ0k2SBoF7gJMlPMfM2pRlGevXr0dSrlVQq0ZGRtzB3EfyDlc9DrwCbJU0J2lvRFwF9gPPA2eBExHR0nx3NyWZFa866kgSo6OjhSQEmK8leEmL/pJrSYyI2N2g/BRwqt37es9ns2JlWcb999/f9n7Li/FqqP3LS2KY9Zn6+QhlJAXPXu5vSnn0wLZt22J6errbYZj1jOp8hLLWIRoZGXGzUQ+Q9GpEbGv3etcYzHrUYjWDopexqOW+hJUjyRpDTR/Dg2+88Ua3wzFLTtk1g1rj4+M8/vjjpT/HipO3xpBkYqhyU5LZ4opey2gx7lzuXW5KMluBLly4UOr93bm8siW5g5uHq5otbc2aNXz44YeF3zflFgTrnCRrDJ7gZjZvsQ7mW265pfCk4PWNrFaSicHMFl8Ke3R0lNdff73Q54yPj3t9I1sgyaYkMyt+B7V6Hm1kjSRZY/DqqtaPFmsWWqq8rFFH1QXvnBSsEQ9XNeuAxeYdrF27lvvuu4+nnnrqU+URwQcffFBoDFNTU179dIXwPAazHtBo3sHAwAAfffRRqc8uawSTpStvYnAfg1kHNJp3UHZSSPkffpaujvUxSFon6SlJ35Xk+qytKDfeeGNHn7dx40YnBWtb3o16jkq6JOlMXflOSeckzUg6WCn+OvBcRDwIfDXPc816zaFDh1i7dm3pzxkcHGRqaoqLF1vaTddsgbw1hmPAztoCSQPAY8CdwAiwW9II83s/v1k5rdz6s1lCalc9lVTKMyQxNTXFb3/7W3cwW255d3B7WdJwXfFtwExEnAeQ9CxwNzDHfHL4e5ZISJLGgDHofPXbrEhZlnHgwAGuXLnySVnRzTsbN2507cAKV0bn8xD/UjOA+YRwO/AI8KikPwV+1OjiiDgi6W1g1+Dg4K0lxGdWuk4si+0+BCtLxzqfI+K9iPjLiBiPiGyZc71WkvW0sjfMcVKwMpWRGC4CN9Qcb6qUNc0zn60XVWcwSypl1nJ1xrL7EKxsZSSG08DNkm6SNAjcA5ws4TlmyciyjG9+85ulJITt27cTEd5W0zom73DV48ArwFZJc5L2RsRVYD/wPHAWOBERLX2i3ZRkvSTLMkZHR/n4448Lve/mzZu9YY51Rd5RSbsblJ8CTrV7X2/UY6mqDj29cOFCaUtNeEtN6zavlWTWpCzLeOCBB0pbd+i6667zvghWiL7c89mdz5aiffv2lZYUvFmOpSTJxOA+BktJlmVI4t133y383tWhp94bwVKSZGJwjcFSUE0Io6Ojhd97w4YN3h/BkpVkYnCNwbolyzKuueaa0hJCdejpO++846RgyfJ+DGYVQ0NDvPXWW6Xc2yONrJckWWNwU5J10o4dO5BUSlIYHx8nIpwUrKckmRjclGRlqC5ZsWrVKoaHh8myjIGBAV588cXCn+VOZetlbkqyFaF+tdPZ2dnS+hBcO7Be58RgK0KZq52CJ6dZf0myKcl9DFakLMtKWdwO5hNCRDgpWF9JMjG4j8Hyql0Cu4wmo+oCd04I1o/clGR9Z2JigieeeKKUe7vJyFYCJwbrK5/97Gf51a9+Vcq9U15w0qxIHWtKkvQFSd+T9Fynnmkrx8TEBJIKTwpr1qwhIpwUbEVpKjFIOirpkqQzdeU7JZ2TNCPp4FL3iIjzEbE3T7BmtaprGUkqpelofHy8tNVUzVLWbFPSMeBR4OlqgaQB4DHgK8AccFrSSWAAeLju+gci4lLuaM0qduzYUcrEtJGREW+haSteU4khIl6WNFxXfBswExHnASQ9C9wdEQ8DdxUZpFmWZRw4cIArV66U9gw3F5nNy9PHMAS8WXM8VylblKQNkg4DX5b07SXOG5M0LWn68uXLOcKzfpFlGXv27CklKbgPwezTOjYqKSKuAPuaOO+IpLeBXYODg7eWH5mlqsxhp+AaglkjeRLDReCGmuNNlTKz3MrqQwAnBLPl5GlKOg3cLOkmSYPAPcDJYsKylao60qjopLBx40Y3GZk1qdnhqseBV4CtkuYk7Y2Iq8B+4HngLHAiIgoZzuElMVae6p4IRS9fUd0x7eJFV2bNmtXsqKTdDcpPAacKjYj5RfSAXVu2bCn61paALMuYnJzkwoULrFq1io8++qiU57h2YNYeL6JnHVXdF2F2dpaIKDwpVGsITgpm7UtyrSTXGPrXQw89VMq+CJ6YZlYc1xgst8W2zFzMLbfcwnvvvVfos0dGRogIJwWzAinFKndNjeHBN954o9vh2BLqt8wEGBwcZM2aNZ8kgdWrV3P16tVCn+saglljkl6NiG1tX59iYqjatm1bTE9PdzsMW8Lw8HBpu6MtJuXPq1kq8iaGJPsYrHdcuHChI89xQjDrnCT7GLznc++48cYbS72/RxiZdV6SicGdz73j0KFDrF5dfMXTCcGse9yUZG0bGBjg448/LvSeTgZm3ZdkjcHSNDExwerVqz/ZNa3IpOAaglk6kkwM7mNIR3WOQnX7zKJnKjshmKUnycTgPoY01C5fUTQnBLN0uY/BFjU0NMRbb71V+H2dDMzSl2SNwbqjtg+h6KTgGoJZ7+hojUHS14A/Bf4V8L2I+N+dfL4trqzaAbiGYNaLmq4xSDoq6ZKkM3XlOyWdkzQj6eBS94iIv4uIB5nf+/nP2wvZijA0NPTJ6KKymoycFMx6Uys1hmPAo8DT1QJJA8BjwFeAOeC0pJPAAPBw3fUPRMSlyu9/XbnOusD9B2a2lKZrDBHxMvCLuuLbgJmIOB8RHwLPAndHxD9GxF11P5c0778DP42InxX3MqwZ1T4EJwUzW0rezuch4M2a47lKWSPfAnYA35C0b7ETJI1JmpY0ffny5ZzhrWy1+ySsX7/e8xDMrCkd7XyOiEeAR5Y554ikt4Fdg4ODt3Ymsv5Tv09CkRvkOBGY9be8NYaLwA01x5sqZdYFtTWE++67r/AtNF07MFsZ8tYYTgM3S7qJ+YRwD/AXuaOyltXXEIpqMnIiMFt5Whmuehx4BdgqaU7S3oi4CuwHngfOAiciIvd+i14So3WTk5Ol1BDMbOVJcmtP7/nculWrVhX6RZ7i58LMmpN3a88kl8RwjaF1ReykVu1DcFIwW9mSTAxedrt1eVZAdTIws1pJJgbXGDrHCcHM6iWZGFxjWFp1jaPan2bVNhc5KZjZYpJMDK4xNNZKEqhyIjCzViSZGMzMrHuSTAxuSiqOawlm1qokE8NKbkoqsv/AScHM2pFkYlipGiWBdvoVzMza5cTQYxrVAlw7MLOidHTZbSuGk4CZlSnJGoM7n83MuifJxLCSO5/NzLotycSwUrn/wMxS0LE+BklfAg4A1wMvRsQTnXp2Nyw2kqiZL3gnATPrtqZqDJKOSrok6Uxd+U5J5yTNSDq41D0i4mxE7AP+DPij9kNOn4edmlkva7Yp6Riws7ZA0gDwGHAnMALsljQi6Q8k/bju5/cq13wV+AlwqrBXYGZmhWqqKSkiXpY0XFd8GzATEecBJD0L3B0RDwN3NbjPSeCkpJ8A3283aDMzK0+ePoYh4M2a4zng9kYnS7oD+DrwOyxRY5A0BoxBMbuSmZlZazrW+RwRLwEvNXHeEUlvA7sGBwdvLTsuMzNbKM9w1YvADTXHmyplK56HnZpZL8uTGE4DN0u6SdIgcA9wsoig+mGCm1c6NbNe1exw1ePAK8BWSXOS9kbEVWA/8DxwFjgREa8VEZSXxDAz655mRyXtblB+Cg89NTPrK0kuidEPTUlmZr0qycTgpiQzs+5JMjG4xmBm1j1JJgbXGMzMuifJxOAag5lZ9ySZGMzMrHuSTAxuSjIz654kE4ObkszMuifJxGBmZt3jxGBmZgs4MZiZ2QJJJgZ3PpuZdU+SicGdz2Zm3ZNkYjAzs+7paGKQtE7StKS7OvlcMzNrXrMb9RyVdEnSmbrynZLOSZqRdLCJW/0VcKKdQM3MrDOa2qgHOAY8CjxdLZA0ADwGfAWYA05LOgkMAA/XXf8A8O+A14HfzReymZmVqdkd3F6WNFxXfBswExHnASQ9C9wdEQ8Dn2oqknQHsA4YAT6QdCoiPm4/dDMzK0OzNYbFDAFv1hzPAbc3OjkiJgEk3Q+80ygpSBoDxiqHv61vvmrSZ4Dlxro2Oqe+fKnj5X6/Hnin6aiXj62Zc/LGX1tWVvzNxr5YWbPxtxv7UvEt9/d+j7/Ra+mXz37t770e/9bmQ15ERDT1AwwDZ2qOvwH8r5rjPcCjzd6vyWdOt3ndkXbPqS9f6ni533sx/rqyUuJvNvY88bcbu+NvXNbotfTLZ9/x/8tPnlFJF4Ebao43VcpS8KMc59SXL3XczO/t6Gb8eWNv5h7Nxr5YmeNfXlnxN3ot/fLZb/b5S+n1+AFQJbssf+J8H8OPI+L3K8ergX8CtjOfEE4DfxERrxURWOUZ0xGxraj7dZrj755ejh0cf7et9PibHa56HHgF2CppTtLeiLgK7AeeB84CJ4pMChVHCr5fpzn+7unl2MHxd9uKjr/pGoOZma0MXhLDzMwWcGIwM7MFnBjMzGyBnk4Mvboon6QvSTos6TlJ492Op1WSvibpu5J+IOlPuh1PqyR9QdL3JD3X7ViaVfmsP1V53+/tdjyt6sX3vFYvf+bb+r7JMwmi3R/gKHCJmglzlfKdwDlgBjjYxH3+G/BfgLt6LfbKNauAqV587yvXfBb4Xg/H/1wnY8/zWpifQLqr8vsPuhl3nv8W3X7PC4i/45/5AmNv+vumWy/qj4E/ZOFM6gHgn4EvAIPAPzC/rtIfAD+u+/k95hfvuwe4v8OJIXfslWu+CvyU+bkfPfXe11z3P4A/7OH4u50YWnkt3wb+feWc73cz7nbiT+U9LyD+jn/mi4i91e+bPGsltS16eFG+ImKv3OckcFLST4Dvlxfxp55bxHsv4G+An0bEz8qNeKGi3v8UtPJamF+LbBPw9yTSBNxi/K93NrrltRK/pLN06TO/mFbf+1a/b5L4gFUstijfUKOTI2IyIv4z8y/yu51ICktoKXZJd0h6RNLfAqfKDq4JLcUPfAvYAXxD0r4yA2tSq+//BkmHgS9L+nbZwbWo0Wv5IfCfJD1BQcselGTR+BN/z2s1ev9T+8wvptF73/L3TVdqDEWKiGPdjqFVEfES8FKXw2hbRDwCPNLtONoVEVeAVP/nXlREvAf8ZbfjaFcvvue1evkz3873TUo1hpQX5VtOL8cOjj8lvf5aHH/3FBZ7SonhNHCzpJskDTLfsXyyyzE1q5djB8efkl5/LY6/e4qLvUs96seBt4H/x3w72N5K+X9kfsXWfwYmu93z32+xO/60fnr9tTj+/o3di+iZmdkCKTUlmZlZApwYzMxsAScGMzNbwInBzMwWcGIwM7MFnBjMzGwBJwYzM1vAicHMzBZwYjAzswX+PwTqvfJWjAVXAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.scatter(gene_info.mean_, gene_info.var_, c='k')\n",
"plt.loglog()\n",
"plt.xlim(1e-4, 1e3)\n",
"plt.ylim(1e-4, 1e5);"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"from scipy import optimize"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [],
"source": [
"phi, _ = optimize.curve_fit(lambda mu, phi: mu + phi * mu ** 2, gene_info.mean_, gene_info.var_)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"mm = np.logspace(-3, 3)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x116bea630>]"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAECCAYAAAD6oXArAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvFvnyVgAAIABJREFUeJzt3Xd8VFX+//HXZyaZNCCBRJCWBEhAENYCgg1FiiJSLKiAWFY0SBYXdYu67NpWfuKuuwpfC7KKCIyioigoiBRZsC7Y6EIIEEILBAiQQtr5/ZEEkpBJJpnMzJ3weT4ePB6ZM7d8Jsa8c+659xwxxqCUUkqVsfm7AKWUUtaiwaCUUqoCDQallFIVaDAopZSqQINBKaVUBRoMSimlKtBgUEopVYEGg1JKqQp8Fgwi0kdEVovINBHp46vzKqWUqh2PgkFEZohIhohsqNQ+UER+FZEUEXmstNkAJ4BQIN2T8yqllPIe8WRKDBG5ipJf9rOMMV1L2+zAVmAAJQGwBhgJbDHGFItIC+Dfxpg7PC1eKaVU/fOox2CMWQUcrtTcE0gxxqQaY/KBucAwY0xx6ftHgBBPzquUUsp7grxwzNbA7nKv04FeInIzcB0QBbzsamcRSQKSACIiIrqfd955XihRKaUarh9++OGQMeacuu7vjWCokjHmI+AjN7abLiL7gCEtW7bsvnbtWu8Xp5RSDYiI7PJkf2/clbQHaFvudZvSNqWUUgHAG8GwBkgUkXYi4gBGAAu8cB6llFJe4Ontqu8C3wKdRCRdRMYYYwqB8cASYDPwvjFmY22Oa4xZaIxJioyM9KQ8pZRSdeDRGIMxZqSL9kXAoroeV0SGAEMSEhLqegillFJ1ZMkpMbTHoJRS/mPJYBCRISIyPSsry9+lKKVU4Cgugs//4vFhLBkM2mNQSqlaKi6Cj5Phu1c8PpQlg0F7DEopVQvFRfDxOFg3F/r+1ePDWTIYtMeglFJuKi6C+Q/Auveg79/gqj95fEifPfmslFKqnhUXwfyxsP4D6PcE9P5DvRzWkj0GvZSklFI1KCr0SiiARYNBLyUppVQ1KoTCk/UaCqCXkpRSKrAUFcL8JNjwIfR/Cq58uN5PocGglFKBoqgQProPNs6H/k/DlQ955TSWvJSkYwxKKVVJUQF8OKYkFAb83WuhABYNBh1jUEqpcooKYN69sOljuHYSXPF7r55OLyUppZSVFebDvN/Clk/huufgsmSvn1KDQSmlrKp8KAx8Hi59wCen9emlJBGJEJG1IjLYl+dVSqmAU3gS3r+rJBSu/4fPQgE8X6hnhohkiMiGSu0DReRXEUkRkcfKvfUo8L4n51RKqQavIA/euxO2LoZBL0CvsT49vac9hpnAwPINImIHXgGuB7oAI0Wki4gMADYBGR6eUymlGq6CPHjvDti2BAa/CD3v93kJnq7gtkpE4is19wRSjDGpACIyFxgGNAIiKAmLXBFZZIwp9uT8SinVoBTkwrsjIXUlDJkK3e/2SxneGHxuDewu9zod6GWMGQ8gIvcAh1yFgogkAUkAsbGxXihPKaUsKD8H3h0BO1bBsJfhotF+K8XndyUZY2bW8P50EdkHDHE4HN19U5VSSvlRfja8czvs/ApufBUuHOXXcrwRDHuAtuVetyltU0opVdnJ4+C8DXZ/Bze9Dhfc7u+KvHK76hogUUTaiYgDGAEs8MJ5lFIqsOUdgzm3wO7v4ZY3LBEK4Pntqu8C3wKdRCRdRMYYYwqB8cASYDPwvjFmY22Oq1NiKKUaKqfTSXx8PE3DbPz8h3iKd6+BW9+Crrf4u7RTPL0raaSL9kXAoroeV0SGAEMSEhLqegillLIcp9NJUlISIcU5LL0zgi7NChkxv5hhiSe4o4u/qztNjDH+rsGlHj16mLVr1/q7DKWUqhfx8fGcyEhj6Z3hdDnHxs3v57JoWyFxcXHs3Lmz3s4jIj8YY3rUdX9Lzq6q024rpRqivMzdfHl3OOfF2Bg6N4dF2woBSEtL83NlFVkyGHSMQSnV4Bzbx1djmtC+qY3B7+bwxfaiU29Z7ZktSwaD9hiUUg1KVjrMHERc02BunFfEih2nQyE8PJxJkyb5sbgzWTIYtMeglGowjuyEt66H7EyC7/2Me56cTlxcHCJCXFwc06dP54477vB3lRXo4LNSSnlL5nZ4e0jJk813fQytLvLJaXXwWSmlLKLsGQWbzUa/C9qS8+o1UJgH93zqs1CoD5YMBr2UpJQKNMnJydx5553s2rWLbs2Fd689yrGso3waMw7O7ebv8mrFksGglFKBxOl0Mm3aNIwxXNLKxsq7I8grhN5vZTP+mVf8XV6taTAopVQtlb9kFB8fz4QJEzDGcGWsnWV3RXAkz3DVW9mkHC623DMK7vD5tNtKKRXIyqa1yMnJAWDXrl0A9Gtn55MR4ew+Vky/WTnsPV5yY4/VnlFwhyV7DDr4rJSyqokTJ54KhTKDEoP4dFQ4248Uc/XM06EgIpZ7RsEdlgwGHXxWSllV5UtDw7sE8fHtYaw/UMw1b+eQkX06FB544AHLPaPgDksGg1JKWVWzZs1OfX3PhcHMvSWM79KL6D87m8O5JaEQFxfH7NmzefXVV/1Vpkd8NsYgIp2BCUAMsNwY85qvzq2UUvXtwZ4Opl4fypKUQm5+P4ecgpJ2EanXmVL9wdOFemaISIaIbKjUPlBEfhWRFBF5DMAYs9kY8wBwG3CFJ+dVSilfcjqdxMTEICJkZmbyl94lofDR5gKGzj0dChCYg82VeXopaSYwsHyDiNiBV4DrgS7ASBHpUvreUOAzPFjERymlfCk5OZnRo0eTmZkJwOT+IUzqG8qsX/K57YNc8k/Ph2eJCfGKij2f5sijYDDGrAIOV2ruCaQYY1KNMfnAXGBY6fYLjDHXA4E3GqOUOuuUPbgGIMArg0J59IoQXlmTzz0f51FU7ndwdHS0XyfEKyo2fPLzHgb8+78eH8sbYwytgd3lXqcDvUSkD3AzEEI1PQYRSQKSoGF0yZRSgcfpdDJx4sRTzygE2eDtG8MY1S2YyV+d5PHlJ09tGxcXx6RJk/wWCMXFhk/X72Pq8m2kZJygU4vGHh/TZ4PPxpiVwEo3tpsuIvuAIQ6Ho7u361JKqfKSk5N57bXT98aEBcEHt4ZxQ8dg/rw0j39+k3/qvfpekrM2iosNizfsZ8ryrWw9cIKOLRrxyqiLub7rudgf8ezY3giGPUDbcq/blLYppZQlVe4hlIkMgYUjw7ki1s79C3N548fTo8zBwcF+GU8oLjYs2bifKcu3sWX/cRKaN+L/Rl7EDd1aYrNJvZzDG8GwBkgUkXaUBMIIYFRtDmCMWQgs7NGjx/1eqE8ppU6pPMVFmeYRwpLR4XQ5x8bt83KZt6nw1HsRERG8/vrrPr18ZIzhi00HeGnZNjbvO0b7cyKYMuJCBv+mFfZ6CoQyHgWDiLwL9AFiRCQdeNIY86aIjAeWAHZghjFmYy2POwQYkpCQ4El5SilVo6qmuIiLFJbeGU6rxjYGv5PD0tSSW4/8MZ5gjGH55gxeWr6VDXuOER8dzou3X8DQC1rXeyCU8SgYjDEjXbQvQm9JVUoFgMpTXHRtbuPzO8IJCxb6z87hu/SSUIiOjvbpeIIxhpW/HuSlZVv5JT2LuOhwXrj1Am68sBVBdu9OWmHJKTF0riSllLclJycTFBRE+eWNr4y1s/q3EQBc9Vb2qVCw2+1MmTLFJ3WVBEIGN736Db+duYbM7Hz+cctvWPbI1Qzv3sbroQAWnXZbLyUppbzp/PPPZ9OmTRXahnYKYu4tYezKKua6OTmkZZ0OjLffftvrl4+MMazedogXl23lp7SjtI4KY/LN3bilexuCfRAG5WmPQSnVYFVeUMfpdFYZCmMuCuaj28JYd6CYK2ecDgURYc6cOV4NBWMMX6cc4tZp33LXjP9xICuPSTd15cs/9mFEz1ifhwKAlO9GWUW5HsP927Zt83c5SqkA5Opuo8r+0tvBpL6hLN5WyPAPKs575O1Q+HZ7Ji8u3cr/dh7m3Cah/O6aDtx2SVtCguweHVdEfjDG9Kjz/lYMhjI9evQwa9eu9XcZSqkAFB8ff8ZzCeXZBKYMDGV8Twezf8nn3gV5FBaffn/cuHFemzb7+9RMXly2le9SD9O8cQi/uyaB2y9pS2iwZ4FQxtNgsOQYg1JKeaq6tZZDg+Cdm8O4qXMwL3xzkj8vPUn5P5H79evnlVBYu/MwLy7bytcpmZzTOIQnBndhVK/YeguE+qLBoJRqkJo1a3ZqRtTyosOEBSPDuLSNnQmf5zH1+9NTXNhsNmbNmlXvl49+TDvCi0u3snrbIWIaOfjrDZ25o1ccYQ5rBUIZSwaD3pWklPKE0+msMhTaRQmfjw4nNtLGrR/k8tHmwgrv13co/Lz7KC8u3cp/tx4kOsLBXwadx+hL4wh3WPJX7ymWrE6nxFBK1VX//v1Zvnz5Ge09Wtn4dGQ4QTboNyuHb3afXkihvtdnXp+exYvLtrJiSwZNw4N5dOB53HVZHBEhlvyVe4bAqFIppdzgKhRuSAziveFhZGQbBjpz2Jp5epQ5OjqaKVOm1EsobNiTxUvLtrFs8wGiwoP503WduPvyeBoFSCCUCaxqlVLKheTk5CpD4cGeDl68LoSf9hcz+J0cDmSfHmbu0qULGzfWaiq3Km3ae4yXlm3li00HaBIaxB8GdOSeK+JpHBrs8bH9wZLBoGMMSqnacDqdFdZQALALvDgwlAd7Opi/uYDR83MrPKMQFRXlcShs2X+MKcu2sXjDfhqHBvFQ/0TuvbIdTQI0EMrocwxKqYBW1eWjxg6YOzyMQYnB/PObkzxa6XZUEaG4uJi62nrgOFOWbeOz9ftoFBLEvVfEM+bK9kSGWyMQ9DkGpdRZyel0cuedd1L5j9u2TYRPR5Wso5C0MJf/lFtcp8zs2bPrdM6UjONMWZ7Cp+v2Eh5sZ/w1CdzXux1R4Y46Hc+qfBYMInIjcAPQBHjTGPOFr86tlGpYnE4no0ePPqO9RysbC0aEEx4sXO/MYVlqUYX3GzVqxLRp02o90Lz94AmmLt/Ggl/2EhZs54GrO3B/7/Y0i2hYgVDG04V6ZgCDgQxjTNdy7QOBKZQs1POGMWayMeZj4GMRaQq8AGgwKKXqpKpQGNUtiDeHhrHvuKHfrGw2Hzp9qaiuD67tOJTN/y3fxsc/7yEkyE7SVe1J6t2e6EYhHn8GK/O0xzATeBmYVdYgInbgFWAAkA6sEZEFxpiy6Qz/Wvq+Ukq5zel0MnbsWLKzsyu02wQm9Q3hsStDWLmzkOHv55KZe/ryUl1uR03LzGHqim3M/2kPwXbhvt7tSbqqPTENPBDKeLqC2yoRia/U3BNIMcakAojIXGCYiGwGJgOLjTE/enJepdTZJTk5+Yy7jqBkkPmdW8IY3DGY19bm8/vFFSfCi4uLq9Wqa7sP5/DyihTm/ZhOkE245/J4xl7dnuaNQ+vhUwQOb4wxtAZ2l3udDvQCHgT6A5EikmCMmVbVziKSBCQBxMbGeqE8pVQgcRUKHZoKC0aG0zHaxrjPcpm2tuIgc3BwMJMmTXLrHOlHcnjlyxQ+WJuOzSbceWkcyX060LzJ2RUIZXw2+GyMmQpMdWO76SKyDxjicDi6e78ypZRVVfV8AsCA9nbmDg/HGMOA2Tms3Fl0xjZvvfVWjZeP9h7N5ZUvU3h/7W4E4Y5esYzrk8C5kWdnIJTxRjDsAdqWe92mtE0ppdxWVU9BgMeudPBs3xA2ZhQzbG4OO45WvF01LCysxsV59mXl8uqX23lvzW4MhtsvaUtynwRaRYXV98cISN4IhjVAooi0oyQQRgCjanMAnURPqbOL0+lk4sSJpKWlERsbS0ZGBrm5uRW2aeyAt28sWUPhnfUF3L/wzCeZjxw5Uu15DhzL47WV23nn+zSKjeHWHm353TUdaNM03BsfK2B5ervqu0AfIEZE0oEnjTFvish4YAklt6vOMMbU6rlznRJDqbNHcnIy06ZNO/WgWlWrrnWOsTH/9jA6NLOdsYYCQKtWrdizx/WFiYzjeUxbmYrz+10UFhuGX9yG8X0TaNtMA6Eqnt6VNNJF+yJgkSfHVko1fE6ns0IoVGV4lyDeGhbGiXxD37dzWJ1WcTyhunWZDx4/yev/3c6c73dRUGS4+aLWPNg3kdhoDYTq6FxJSim/qW5d5mAb/GNACA9dGsI3u0ueT9h3ouLvK1ezo2aeOMn0VanM+nYXJwuLuPGi1vy+byLxMRFe+RxW0yDnStJLSUo1fE6n02UoxEcJ7w0Pp2drOy99V7Imc0GlOe+qCoUj2flMX53K29/sJK+giGEXtubBvgm0P6eRtz5Gg2TJYNDBZ6UaNlcL6gAM61Ry6UgEbnovh4+3VFx+s1+/fixbtqxC29GcfP6zOpWZX+8kp6CIoRe04sG+iSQ010CoC0sGg1Kq4XA6nUyYMOHUGsxBQUEUFhaesV2wDZ4fEMLDl4awZk8Rt88781bUype+s3IKePOrVGZ8vZPs/EJu6NaSCf0SSWzR2Hsf6CxgyWDQS0lKNQxVTY1dVSi0byq8c3M4vdrYmfr9Sf609CT5lZ5ZGzdu3Kmvj+UVMOOrHbz51Q6O5xUyqNu5TOjXkU7naiDUBx18Vkp5TePGjTlx4kS124z+TTCvDgqlyMB9C3L5cPOZwVF2+eh4XgFvfb2TN1anciyvkOvOb8FD/TvSuWUTb32EgNQgB5+VUg1DdaHQJAReuyGMUd2CWbWrkDvn55KWdeYfqnPmzGHY8Nt5ecU2/rN6B1m5BfTv3IKH+ifStXWkN8s/a1kyGPRSklKBqfwTzOHhrp8VuKyNHefNYbSNFP66Io/nvsqnuFImREVFkb7/IG9/u5Pez6/gSE4B/c5rzkP9O9KtjQaCN+mlJKVUvXA6nSQlJVU7T1GQDSb2DuFvVznYedRwx0e5fL/nzAnwWsXG88SsZby+KpXD2fn06XQOD/XvyIVto7z5ERoMTy8laTAopepFTEzMqTuPqtI5xsasm8Lo0crOrF/yGb8oj+MVZ7ZAgkK49nd/51DMhRw6kU/vxBgeHtCRi2Obern6hkXHGJRSfpecnOwyFGwCD13qYFLfEI6fNNz8Xg7zKz2bIEEOLr7199DlOracOMmV5zbh4QGJdI9r5ovyVSWWDAYdY1AqcLhaMwGgXZQw88YwrooL4uMtBYz9NI+M7HJXKezBNL7gOqJ7j+JQaBMua96IV++4mJ7tNBD8yZLBoE8+KxUYXK2uJsD93YP517WhFBXDXfNzmb2u3BzZ9iAa/eZaIi+7jaDGMVzSrhmPDOjIpe2jfVe8csmSwaCUsqbydx25Gp/s0FT4z5AwrmkXxNLthdy7IJf0Y6Xb2oJo9JsBRF52K0FNmnNJfFMe7t+RyzpEIyI+/CSqOhoMSimXygdBREREtc8llI0l/P2aEAqKSh5We/On0l6CzU6jbv2JvOx2giKbE3piL2+M6cUVCRoIVuSzYBCR9sBEINIYM9xX51VK1U3l20+rC4Uu59iYMTSMXm3sfLKlgORFeew9bkoCoWvfkkCIOpeTe7dw3wXhPDn2Pg0EC/N0BbcZwGAgwxjTtVz7QGAKJSu4vWGMmWyMSQXGiMg8T86plPKNiRMn1rh2ssMOj18Zwl96O8jKM4yYl8N7GwtBbER07Uvk5SMIbtqSk/u2cu8FjXj6uUc0EAKApz2GmcDLwKyyBhGxA68AA4B0YI2ILDDGbPLwXEopL6q87rKrtRLKXB1n5/XBoXSKsfPO+gImfJ7HoVwh4vxriLx8JMHNWnFyfwoZ854me9v3GggBxNOlPVeJSHyl5p5ASmkPARGZCwwD3AoGEUkCkgBiY2M9KU8p5abKl42qC4XoMOGFa0O450IH2w8Xc+3sbJbuMISfdxWtrhhJcHQb8g9sJ+PDv9M0Zzc51azFrKzJG2MMrYHd5V6nA71EJBqYBFwkIo8bY56ramdjzHQR2QcMcTgc3b1Qn1KqkgkTJtR42Qjg7guCeeHaECJDhEmrT/Lsqnxsib1pee8IHDGx5GfsIOOjSeRu+w5jims8nrImnw0+G2MygQfc3FafY1DKR5xOZ7VTWQB0a27j5UGhXBUXxFdphYz99CQ7m11K07tG4TgnjvyDuzj48XPk/PoNUVGR5GgoBDRvBMMeoG25121K29ymTz4r5Ttjx451+V6TEHi6Twjjezo4mme4b0Eec7O702TwKM5p3o78Q2kc/GQyOVu+Bgxz5szhjjvu8F3xyiu8EQxrgEQRaUdJIIwARnnhPEqpOnI6nYwdO5bs7GyX29zRreSyUfMI4fUfCnh2ZzdMz7uIad6OgszdHFzwD3K2fAWmmC5durBx40YffgLlTTZPdhaRd4FvgU4iki4iY4wxhcB4YAmwGXjfGFOrnxhjzEJjTFJkpM65rlR9cTqdxMfHIyKMHj3aZSj8poWNlXeHM+fmMNKyDFd9Hs8zMZMJHvwkEuTg0MIX2Pvm78jZvApMMXPmzNFQaGA8vStppIv2RcCiuh5XLyUpVb+cTif33HNPlestl4kOE/7eN4Ski4M5kmf43XctmN/iARzXdMJ2ZC+HPvs32RtXQun4Qdlym6rh8ajHoJSynrKegc1mIz4+vsZQCLLB+J7BbHuwEfdfHMy0X5tyZcZjfHbhi0hYJIcWvcTeN8aRvWHFqVCYM2eOhkIDpgv1KNWAuLOKWnn92tmZMjCU85vbWb43jGcK72ZX86sozDpA1jfvcWLDcig+vcKajiUEBk8X6tEeg1IBqqqegTvTWAB0jLbxyYgwlt0VQVhoCHfvG86YZm+wPfQ8Mj//P/ZMH8uJdV9UCAUdSzh7WLLHUG6M4f5t27b5uxylLKe2PYMyzcKEJ652kNzDQW6RjalZVzO78T1kH8si69v3ObF+KRRVvOQ0btw4Xn311fosX3mZrvms1FkoPj6+xrmMygu2QfIlDp64uuSp5TnHLuT/Qsey/3hhSSCsW3JGIOjgcuDSNZ+VOgulpaW5ve1N5wUxuX8IHaPtfHkinucKxrGpMIJjy+Zy/OfPoajgjH30QbWzmyWDQW9XVap6wcHB5OfnV7vNpW3svDAghCtig9hyMpp78u9l+ck4sv77ISd+XowpPHN/K19BUL5jycFnfcBNqRJVDTCff/751YZCh6bC+8PD+HZMBAktI3ms4D4GZj/N/OU/sef1+zm+9pMzQiEqKkpDQZ1iyR6DUqrqqbBHjx7tcvuYcOGvVzkY1yOEAoJ5sWAo045fyf7/LeL4j69iCvKq3E8Hl1VlGgxKWZS7t56GB8PDlzr485VhRATDe0XX8K8T15Hy/QqO/zgek59b5X4aCMoVS15KEpEhIjI9KyvL36UoVW+quixUXXtNdx0F2SCpezApEyJ5tm8o39h70O/4UySvCOGn1x7h2HcfVBkKXbp0wRijoaBc0ttVlfKBqp47CA8P5+677+btt98+o90YQ25u1X/pA9zcOYjnBjSiY9Ni/lfcif+XfSMrv1vPsbULMPlV9zL0TqOzhz7HoFQAcPXcgd1up6ioqIo9qnZNvJ3nBzbhkhZFbC1uzeTcm/joux0cW7MAc7Lq2VLduYNJNSz6HINSAcDVcwfuhsLFLW1MHhjJgNgi9hRH8UjOUN7+7hBH10ylOO+Ey/2s/Iefsi6fBYOIRACvAvnASmOM01fnVsrfYmNja/WkcpnEZjYmDYzi1sRCDpswns67gde+z+Xg929RnHvM5X6tWrViz55aLZyo1CmeLtQzQ0QyRGRDpfaBIvKriKSIyGOlzTcD84wx9wNDPTmvUoFm0qRJhIeHu719mybCm8ObsWl8YwYl2Hnp5FAu/Ooanpn6HgdWznEZCg6Hgzlz5mgoKI942mOYCbwMzCprEBE78AowAEgH1ojIAkrWfl5fupn7F1WVCnDlZz0VkWov78SECxP7NeWBCw0iMKugH5N/CGX7VwspznF9l56IMHv2bB1cVvXC0xXcVolIfKXmnkCKMSYVQETmAsMoCYk2wM9U01MRkSQgCUq630oFKqfTyYQJE8jMzDzV5ioUGjvgj9c045FLIMxWxLzCK5j0UyQbVy2mOPuoy3PoJSPlDd4YY2gN7C73Oh3oBUwFXhaRG4CFrnY2xkwXkX3AEIfD0d0L9Snlde5Oix0WBA/2bsafL7MTHXySzwp78Pcfz2Htl19QlH2k2n11YFl5i88Gn40x2cBv3dx2IbCwR48e93u3KqW8o6anloNtMPbyZkzsHcy5jlxWFnbm6Z9bsXr5CopOZLrcD/R5BOV93giGPUDbcq/blLa5TWdXVYGobCyhuruP7AJ39WzKk31CiQvN5vui9oxZH8uSpf+l6PjX1R5fl9VUvuKNYFgDJIpIO0oCYQQwygvnUcoynE4nd911F8XFxVW+L8Dt3aN4um8EHcOPs67oHIZvvpSPF6+m6Hj1D3HqgjnK1zy9XfVd4Fugk4iki8gYY0whMB5YAmwG3jfG1OrPHJ12WwUSp9PJ6NGjXYbCsN9Esu6PbXl3cDEmtAl3b+3DxVP38eH78yk6fsjlcePi4pgzZ46GgvI5T+9KGumifRGwqK7H1UtJyqrKLhelpaXVONXEtZ0j+X8Dm9K9yWF2FIcydvsFzPz0f+Qf3VztObSHoPxN50pSyk1Op5N77723xnmHrk5swqRBMVwRdYj04mj+uasz0z/9kbzD+6vdLyoqiiNHqr8TSSl3eDpXkk67rZSbHnjggWpD4dL2TVj+YAIrR0FiZAGP7rqMTq/lMXXWohpDYdy4cRoKyjIsGQw6xqCsxOl0IiKcOFH1ZHUXxzZm8e8S+fZOuKBpLk/svpSEacX8Y+YScg7tq/bYc+bM0bURlOVYcnZVHWNQVlA2qOxK19YRPDu0DcOa7yPLHGfS3p78a+EWjuz/osZjR0dHM2XKFH0eQVmSjjEoVY7T6eS+++4jL6/q9ZEBOp4bwbPD2nJLi31kE8rrB87n+YXbOLR3t8t9yujAsvIFXY9BqXrSunVr9u7d6/L99s3D+fuwOG5vuZeTZPLywYv4fwtTOZC+osZjayCoQGLJMQYdfFa+1L9/f0TEZSjExoTx1pgubHnAwc2Ps10SAAASiklEQVQtM/jPoQvo+HY4E15byYH0qhfgKTNu3DiMMRoKKqBYMhh08Fl5g9PpJD4+HpvNRnx8PE6nE7vdzvLly6vcvnXTUKbfcz7bkkMZ2Xo/sw5347zZjRn36ir21LDojg4qq0Cml5LUWaHybKe7du1yObDcIiqEvw1NYEz8Puzs5d0jXXhm8X62p6yu8Tx6yUg1BBoM6qxQ02ynAOc0CeEvQxMZ234fwezhg6OdeWpJBlt/rX5yO9CH01TDYslg0NtVVX1yOp3VznjarJGDx4Z1JLnDfsLYzYfHOvPUkkNs2vxtjcfWQFANkY4xqAapbDxBRFxeMooMD+a5EV3Z+XAUf+iwm6XH23PRh8257aXv2LQ5pdrjl01wp6GgGiJL9hiU8kRycjKvvfaay/cbhwXxxyHn8ftOB4mypfHp8U48uSyLH9f9r8Zjaw9BnQ00GFSD0rRpU44erXqN5IiQIP4wtDMPnXeQprY0Pj+RyFPLT/D9z2vcOraVHwZVqj75LBhEpD0wEYg0xgz31XnV2aG6XkKYI4iHh3TmkS6HiLbtYll2Ak+uyOGbH3+o8bg1Ta2tVEPk1hiDiMwQkQwR2VCpfaCI/CoiKSLyWHXHMMakGmPGeFKsUuWVTW4nIlWGQqjDzp9v6srOP7dgUtddrM9rwdWL2jDghR/55sctNR5/3LhxGgrqrORuj2Em8DIwq6xBROzAK8AAIB1YIyILADvwXKX97zXGZHhcrVKl+vfv7/LBtJAgG7+74Xz+1O0I59rTWJ3bjpH/jWLF9z/XeFxdV1kpN4PBGLNKROIrNfcEUowxqQAiMhcYZox5Dhhcn0Uq5XQ6mTBhApmZmS63CbbbGDeoC49ekEUr+y6+zYvj7hWRfPHNL26dQ8cQlCrhyRhDa6D8dJLpQC9XG4tINDAJuEhEHi8NkKq2SwKSAGJjYz0oTzUUTqeTO++80+Uv7iC7jaSBXXj8wmO0CUpjTV4s930ZyeKv19d4bB1DUOpMPht8NsZkAg+4sd10EdkHDHE4HN29X5myqppuO7XbbNx3XWcevzibuKA0fjjZluTVkSxctQGo+a9/7SEoVTVPgmEP0Lbc6zalbUp5rLoxBLsI91zbmb90z6V98G5+zm/Nw/+NY/5KDQSl6oMnwbAGSBSRdpQEwghgVL1Upc5a1a2aZhPhrv6d+esleXQITmdDfkuGr47nwxXrqSkQWrVqxZ49+neLUu5w93bVd4FvgU4iki4iY4wxhcB4YAmwGXjfGFMvt3PolBhnn7I1EaoKBQFG9+vM5sfa89bl6ZwkiJFfx9Nt8jY+XLGO6kKhX79+GGM0FJSqBXfvShrpon0RsKheK0In0WvonE4nEydOJC0tDZvNRlFRUZXbCTDimvN44tJCznPs4deCFoz+ph3vLFuPMcU1nkcvGSlVNzqJnvKpsnURdu3ahTHGZSjc3uc81j+WyDtX7SXIZrjnu/Z0mbwd59Jfqg2Fsh6ChoJSdWfJuZK0x9BwjR07ttp1EW65qhNPXm7oFrKX1MJzGLOmAzOXrKPYRYCU0QfTlKo/2mNQHqtqycyqnH/++WRnZ1f53rArOvLTo+cx75p9NLYXkrS2Ax0n72TGop+qDYUuXbpgjNFQUKoeiRW73OV6DPdv27bN3+WoalReMhPA4XAQHBx8KgSCgoIoLCyscv8hlyfyVG87F4fuZVdhNJPXRfHG4vUuty+jPQSlXBORH4wxPeq8vxWDoUyPHj3M2rVr/V2GqkZ8fHy1q6O5MqhXAk9dHcwlYXtIL2rK5HXNeH3xegoLqg8EK/+8KmUVngaDJccYVOBIS0ur1fbXXdKBp/qEcGl4OnuKmvL7XxKY9tk6CgqqDxcNBKV8x5LBoIPPgSM2NtatHkP/7h14qm8oV4TvZn9RJI+sS+CVz9aTn6+BoJTV6OCz8sikSZMICnL990Wfi9qx6o9dWTr4IAmhWfxpYwLt/nmAF+f/SH5+gcv99JZTpfzHkj0GFRjsdjvFxVU/U9D7gnie6d+YPo12kVFcwGObEpm6cAO5eenVHlPDQCn/s2SPQVlTcnIyQUFBp1ZNqyoUrvhNHMsf6caqGw/TNTyTiZsTaffCQZ7/4Ady8066PLb2EJSyDkv2GHSMwTrKpq+oaRzh0q6xPHNtFAMa7+RwcT5P/NqRfy/cSHb23mr30zBQynr0dlXlUlXPKFR2SZe2PHNdMwY22cGR4gimbm/NvxZs5PiJ3GqPbeWfO6UCnd6uqryidevW7N3r+q/97p1a8/T1MdwQuYOs4gKeTenIPz/ZxLET+6o9rgaCUtanwaBOSU5OZvr06S4ntgO4qGMbnrk+hsFRqWSZQp5L7cjzn2wm61j1PTsNBKUCh0+DQURuBG4AmgBvGmO+8OX5VdVq6h0AdEtoyTODWnBj01SOm3ye39GR5z/ZwpEsDQSlGhq3g0FEZgCDgQxjTNdy7QOBKYAdeMMYM9nVMYwxHwMfi0hT4AVAg8FP3AkDgPM7nMszg1pyY9NUctjDCzs7MvmTX8k8qoGgVENVmx7DTOBlYFZZg4jYgVeAAUA6sEZEFlASEs9V2v9eY0xG6dd/Ld1P+YE7odC53bk8fUNLbm6WSh7pTElL5LkFWzl42HUgaBgo1TC4HQzGmFUiEl+puSeQYoxJBRCRucAwY8xzlPQuKhARASYDi40xP9a1aFU37owhdIprzlOD23Br9HZOsptX0hOZ9MlWMjK1h6DU2cLTMYbWwO5yr9OBXtVs/yDQH4gUkQRjzLTKG4hIEpAEJfPwqLorv4RmeHi4y7UQABLbnsOTQ9pye0wqhexi2p5Enl2wjf0HNRCUOtv4dPDZGDMVmFrDNtNFZB8wxOFwdPdNZQ1P5WcQXIVC+zYxPDUkjhHnpFLETv6zN4FnP0lhbzWBoGGgVMPmaTDsAdqWe92mtE35Qfkegs1mq/aSUXyrGJ4cGscdzUsCYca+Djy7IJX0AxoISp3tPA2GNUCiiLSjJBBGAKM8rkrVWuUegqtQiG3ZjCeHtmN0i1QMO5i1vz3PLNxB2r6qA0HDQKmzj9tTYojIu0AfIAY4ADxpjHlTRAYBL1FyJ9IMY8yk+ipOp8RwX00rqbU9tyl/G9qeu87dgWBwHujA0wt2smvfIZf7aCgoFZga5NKeuuZz7dlstip/kbduHsXfhnXg7pY7sGN4J6M9Ty/cxY49rgMBNBSUCmQNMhjKaI/BfZV7DOfGRPK3YQnc23onQRQx91AHnl6YRsrugy6PYeWfBaWU+xrkJHo67XbtlYVCi+gmTByWyH1tduEghfcOJfD0Z7vZuusHl/tqICilytMeQwPRPDqSicMSub/tLkIoYF5mB576bC9bdu6vdj8r//dXStWN9hjOQiUPkJeIiWrM4zd2Ysf4EMJI4aMjHXjys31sSq36wXINAqVUTbTHEGDKQqFZkwgeu7EzD8SnEcFJPj7agacW7Wd9yplzIFn5v7FSqv41yB6Dci2qSQSPDTuPce3SaSJb+eRIB578/CC/bNWpp5RS9cOSwaCXks507GgmGz96np0PNSJStvHZ0fY8seQQP275qdr9tLeglKotvZRkMeXHDwCaNApjydQH6bJrNk3IZnFWO55Ycpi1m3efsa+V/1sqpXxHLyU1IOVDoXFEKH8Y2oXfJ+yj6a5p/Bx+GY2u+yuDLuztxwqVUmcDDQaLaRQewiNDz2dC4j6a2VJYerwdTy7N4pt1nwMlvYLKvYqydqWUqg8aDBaRm32cv9x2MQ932k+MLYXlJ+J5cmkIX6/75YxtNQSUUt5kyWA4mwaf83JO8PP8f5Ow7U0mdT7KyhNxPLE8lNU/r/N3aUqps5QOPvtJXm42P3/8Egm/TieGo2wIuZDfz9vDlz9ur3J7K/93UkpZiw4+B5i83Gx++WQq7bZM51IOs9HRjQN9XqPr5YNY8fiZdyWBhoJSyrd8Fgwi0hmYQMl6DsuNMa/56tz+UPkXvMMRxKrZk4jfNI1eZLIpuCsZfabS9YohFbbTEFBK+ZvNnY1EZIaIZIjIhkrtA0XkVxFJEZHHqjuGMWazMeYB4DbgirqXbH3lQyE4KIgHh15I6p9b0WvTJI4Et2B931l0fnz1GaGglFJW4G6PYSbwMjCrrEFE7MArwAAgHVgjIgsoWcntuUr732uMyRCRocA4YLaHdVteUJCdsdd349ELjtDWnsravNaMXRHGwtVfIza38lgppfzCrWAwxqwSkfhKzT2BFGNMKoCIzAWGGWOeAwa7OM4CYIGIfAa8U9eira7lOZF8m9SUuKBUfsprxfiVYSz4ajOAhoJSyvI8GWNoDZSflyEd6OVqYxHpA9wMhACLqtkuCUgCiI2N9aA8/9l3MIsvMxNYsDGU+au3+LscpZSqFZ8NPhtjVgIr3dhuuojsA4Y4HI7u3q7LW347zfWKaUopZWWeXNfYA7Qt97pNadtZz9WdRXrHkVIqEHgSDGuARBFpJyIOYASwoD6KMsYsNMYkRUZG1sfh/MIYc8Y/pZQKBO7ervou8C3QSUTSRWSMMaYQGA8sATYD7xtjNtZHUSIyRESmZ2Vl1cfhlFJK1YK7dyWNdNG+iGoGkpVSSgUeS9472RAuJSmlVKCyZDDopSSllPIfSwaD9hiUUsp/LBkM2mNQSin/sWQwaI9BKaX8x5LBoJRSyn8sGQx6KUkppfzHksGgl5KUUsp/LBkMSiml/EeDQSmlVAUaDEoppSqwZDDo4LNSSvmPJYNBB5+VUsp/LBkMSiml/MenwSAiESKyVkQG+/K8Siml3OfuQj0zRCRDRDZUah8oIr+KSIqIPObGoR4F3q9LoUoppXzDrYV6gJnAy8CssgYRsQOvAAOAdGCNiCwA7MBzlfa/F7gA2ASEelayUkopb3J3BbdVIhJfqbknkGKMSQUQkbnAMGPMc8AZl4pEpA8QAXQBckVkkTGmuO6lK6WU8gZ3ewxVaQ3sLvc6HejlamNjzEQAEbkHOOQqFEQkCUgqfXmy8uUrN0UCNd3r6mqbyu3Vva7p6xjgkNtV11ybO9t4Wn/5Nm/V727tVbW5W39da6+uvpreb+j1u/osDeVnv/zXgV5/J/dLroIxxq1/QDywodzr4cAb5V7fCbzs7vHcPOfaOu43va7bVG6v7nVNXwdi/ZXavFK/u7V7Un9da9f6Xbe5+iwN5Wdf6z/9z5O7kvYAbcu9blPaZgULPdimcnt1r935ui78Wb+ntbtzDHdrr6pN66+Zt+p39Vkays++u+evTqDXD4CUpkvNG5aMMXxqjOla+joI2Ar0oyQQ1gCjjDEb66Ow0nOsNcb0qK/j+ZrW7z+BXDto/f52ttfv7u2q7wLfAp1EJF1ExhhjCoHxwBJgM/B+fYZCqen1fDxf0/r9J5BrB63f387q+t3uMSillDo76JQYSimlKtBgUEopVYEGg1JKqQoCOhgCdVI+EeksItNEZJ6IjPN3PbUlIjeKyH9E5D0Rudbf9dSWiLQXkTdFZJ6/a3FX6c/626Xf9zv8XU9tBeL3vLxA/pmv0+8bTx6CqOs/YAaQQbkH5krbBwK/AinAY24c5xngz8DgQKu9dB8bMCcQv/el+zQF3gzg+uf5snZPPgslD5AOKf36PX/W7cl/C39/z+uhfp//zNdj7W7/vvHXh7oKuJiKT1Lbge1Ae8AB/ELJvErdgE8r/WtOyeR9I4B7fBwMHtdeus9QYDElz34E1Pe+3H7/Ai4O4Pr9HQy1+SyPAxeWbvOOP+uuS/1W+Z7XQ/0+/5mvj9pr+/vGk7mS6swE8KR89VF76XEWAAtE5DPgHe9VfMZ56+N7L8BkYLEx5kfvVlxRfX3/raA2n4WSucjaAD9jkUvAtax/k2+rq1lt6heRzfjpZ74qtf3e1/b3jSV+wEpVNSlfa1cbG2MmGmMeouRD/scXoVCNWtUuIn1EZKqIvA4s8nZxbqhV/cCDQH9guIg84M3C3FTb73+0iEwDLhKRx71dXC25+iwfAbeIyGvU07QHXlJl/Rb/npfn6vtvtZ/5qrj63tf6941fegz1yRgz09811JYxZiWw0s9l1JkxZiow1d911JUxJhOw6v/cVTLGZAO/9XcddRWI3/PyAvlnvi6/b6zUY7DypHw1CeTaQeu3kkD/LFq//9Rb7VYKhjVAooi0ExEHJQPLC/xck7sCuXbQ+q0k0D+L1u8/9Ve7n0bU3wX2AQWUXAcbU9o+iJIZW7cDE/098t/Qatf6rfUv0D+L1t9wa9dJ9JRSSlVgpUtJSimlLECDQSmlVAUaDEoppSrQYFBKKVWBBoNSSqkKNBiUUkpVoMGglFKqAg0GpZRSFWgwKKWUquD/A/Rt4uxPIXvgAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.scatter(gene_info.mean_, gene_info.var_, c='k')\n",
"plt.loglog()\n",
"plt.xlim(1e-4, 1e3)\n",
"plt.ylim(1e-4, 1e5);\n",
"\n",
"plt.plot(mm, mm)\n",
"plt.plot(mm, mm + phi * mm ** 2)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
"gene_info['drop_rate'] = 1 - (klein > 0).sum(1) / klein.shape[1]"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEACAYAAABfxaZOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvFvnyVgAAFrtJREFUeJzt3X9s3PV9x/HXO3ZDcNw1KPakxQk2ndKqzjYNaqFOlTakZFtAG0xqN8HuWEBoBmetqDatAmWqKiar+6FNBQmbGIFgcVpgaJrcNRVTWFGlqVSY9YeaIKqU4SSmE8YlFbUDJOa9P3yhl8ud/f187/u9733v+3xIJ933e5/v3fsu59d98/1+vp+PubsAAMWwIesCAACtQ+gDQIEQ+gBQIIQ+ABQIoQ8ABULoA0CBEPoAUCCEPgAUCKEPAAVC6ANAgXRn9cJ9fX0+NDSU1csDQC69+OKLb7h7f9ztMwv9oaEhzc7OZvXyAJBLZjbXzPYc3gGAAiH0AaBACH0AKBBCHwAKhNAHgAIh9AGgQNYNfTN71MxeN7MfNnjczOwBMzthZj8ws2uSLxMAkIQoe/qPSdq7xuPXS9pZuY1Kmmy+rPZmZpfcACAP1g19d/+WpJ+u0eQmSf/iq56XtMXMfiWpAttNo4CPEvxRfiwatam3PuptYGCguTcNoGMkcUx/QNKpquXTlXWoEuXHopkflLW89tprwT8U+/fvb+o1AbSnlp7INbNRM5s1s9mFhYVWvjQCTU5ONvxBOHz4cNblAYgpidCfl7Sjanl7Zd0l3H3K3UfcfaS/P/Z4QchYuVzmhwDIqSRCf0bSn1V68XxC0s/c/ScJPC9ypN4PQW9vLz8GQJuJ0mXzq5K+LemjZnbazO4ws7vM7K5KkyOSXpF0QtLDkjr6YLC7B60vsqWlpYt+DABkb92hld39lnUed0l/kVhFORAn4N29bvBVP9dabTohNC+8h7GxMU1MTGRcDVBMXJHbQu5+yS1qm3rrG93aXe1JYgCtQ+h3oJAfiHb4kaCbKNA6hD7a5gfhwv8A9uzZ09LXBYqE0MeasvghePbZZ+kGCqSE0EewVv0IXOj5AyA5hD6alvaPgJlp165diT8vUESEPhKXxg/A8ePHOeQDJIDQR6qS/gEol8u64oorEnkuoIgIfbRMUuF/5swZjvUDMRH6aLmkwp++/UA4Qh+ZSSL8L/TtBxANoY/MJRH+BD8QDaGPtuHuGh4ejr09wQ+sj9BHWzl27FhTe/0EP7A2Qh9tyd31gQ98INa2ZqahoSH69AN1EPpoW++++27svf65uTmNjo4S/EANQh9tL27wLy8v68CBAwlXA+QboY9caGaPH8AvEPrIDXfX5ZdfHrwdJ3eBXyD0kSvLy8vasmVL8HZdXV0pVAPkD6GP3HnzzTeDt3nvvfcYqA0QoY+cinOM/8yZM4zVg8Ij9JFbcYJ/cnIyhUqA/CD0kWtxgp8TuygyQh+5Fyf4OcyDoiL00RFCg39ycpLgRyER+ugY09PTQe0JfhQRoY+OUSqVtHv37qBtOLGLoiH00VGOHj0avA2DsqFICH10nNDj++VyOaVKgPYTKfTNbK+ZvWxmJ8zsnjqPX2lm3zSz75rZD8zshuRLBaILDX66caIo1g19M+uS9KCk6yUNS7rFzGrntPsbSU+5+9WSbpY0kXShQKjQqRc5qYsiiLKnf62kE+7+iru/K+kJSTfVtHFJv1S5/yFJryVXIhDPsWPHgtpzUhdFECX0BySdqlo+XVlX7YuSymZ2WtIRSZ9NpDqgSaGHedjbR6dL6kTuLZIec/ftkm6QdMjMLnluMxs1s1kzm11YWEjopYG1hQzFzN4+Ol2U0J+XtKNqeXtlXbU7JD0lSe7+bUmbJPXVPpG7T7n7iLuP9Pf3x6sYCBQ6FHNPT09KlQDZixL6L0jaaWZXmdlGrZ6onalpc1LSbkkys49pNfTZlUfbCLlo6+zZsxzmQceyKMc8K10wvyypS9Kj7j5uZvdJmnX3mUpvnocl9Wr1pO7n3f0/13rOkZERn52dbfoNAFGFdsuMOy8vkCYze9HdR+Ju3x2lkbsf0eoJ2up1X6i6f1zSJ+MWAbTC9PQ0F2Kh8LgiF4VRKpWCJlbn2D46EaGPQlleXo7c9uzZs4zLg45D6KNwQo7VczgInYbQB4ACIfRRSGNjY5HbMhgbOgmhj0KamGBMQBQToY/CYm8fRUToo7BC9/a5ShedgNBHoYUMz8BgbOgEhD4KLXROXfrtI+8IfRReyAxbt912W3qFAC1A6KPwQmbYOn/+fIqVAOkj9AFJ27Zti9yWQzzIM0IfkDQ/XzsvUGMMzYA8I/SBipC++OztI68IfaDi0KFDkduyt4+8IvSBilKpFNSevX3kEaEPVOnq6orcdt++fSlWAqSD0AeqPP7445HbrqyspFgJkA5CH6hSKpU0PT0duT2HeJA3FjKLUJJGRkZ8dnY2k9cG1hPSkyervyEUk5m96O4jcbdnTx9oEnv7yBNCH6gj5BDPrbfemmIlQLIIfaCOkO6b7s7ePnKD0AcScPvttxP8yAVCH2hgcHAwcttz587pwIEDKVYDJIPQBxoYHx8P6sVz8uTJFKsBktGddQFAu7pwXD/qODs9PT1plgMkgj19YA0hJ3SXlpZSrARIRqTQN7O9ZvaymZ0ws3satPkTMztuZsfM7CvJlgnkAydz0e7WDX0z65L0oKTrJQ1LusXMhmva7JR0r6RPuvsuSZ9LoVYgEyEndMvlsjZs2KChoSF+ANCWouzpXyvphLu/4u7vSnpC0k01bf5c0oPu/qYkufvryZYJZGd8fDyovbtrbm5Oo6OjBD/aTpTQH5B0qmr5dGVdtY9I+oiZ/beZPW9me5MqEMha6Dj7FywvL9ONE20nqRO53ZJ2SrpO0i2SHjazLbWNzGzUzGbNbHZhYSGhlwbSNzY2Fms7unGi3UQJ/XlJO6qWt1fWVTstacbdz7n7/0r6kVZ/BC7i7lPuPuLuI/39/XFrBlpuYmIi1nZXXnllwpUAzYkS+i9I2mlmV5nZRkk3S5qpafPvWt3Ll5n1afVwzysJ1glkLmQQNmm1337o+QAgbeuGvrufl/QZSc9IeknSU+5+zMzuM7MbK82ekbRoZsclfVPSX7v7YlpFA1kIPbY/NTUV+3wAkBYmUQECMLkKssYkKkAL7d69O3JbumuiHbGnDwRibx9ZYk8faLG43TeBdkDoA4FCum9yiAfthtAHUhR1WGagVQh9IIZt27ZFbsvePtoJoQ/EMD9fe1F6Y3feeWeKlQBhCH0gZUtLSwy3jLZB6AMxhRziYbhltAtCH4gp5BDPBQy3jKwR+kCLMdwyskToA03YvHlz8DYMt4wsEfpAEw4ePBjUnuGWkTVCH2gCwy0jbwh9oIUIfGSN0AeatHXr1sht9+/fn2IlwPoIfaBJ999/f+S2k5OTKVYCrI/QB5rEIRvkCaEPJIAZtZAXzJwFJCTqjFq9vb166623Uq4GnYqZs4A2sWXLlkjtfv7zn6dcCdAYoQ8k5M0334zclkM8yAqhD2Tg1ltvzboEFBShDyQo6nH9rM6lAYQ+kKBDhw5lXQKwJkIfSFBIn32uzkUWCH0gI1ydiywQ+kDCNm3alHUJQEOEPpCwd955J+sSgIYIfSBhITNj7dq1K8VKgEtFCn0z22tmL5vZCTO7Z412nzIzN7PYlwgDeRcyM9bx48dTrAS41Lqhb2Zdkh6UdL2kYUm3mNlwnXYflHS3pO8kXSSQJ6VSKXJ/fYlePGitKHv610o64e6vuPu7kp6QdFOddn8r6e8lvZ1gfUAu3XXXXZHb0osHrRQl9AcknapaPl1Z9z4zu0bSDnf/eoK1Abk1MTGRdQlAXU2fyDWzDZL+WdJfRWg7amazZja7sLDQ7EsDba2rqyvrEoBLRAn9eUk7qpa3V9Zd8EFJvybpOTN7VdInJM3UO5nr7lPuPuLuI/39/fGrBnJgdHQ0ctu+vj5G3kRLRAn9FyTtNLOrzGyjpJslzVx40N1/5u597j7k7kOSnpd0o7szQwoKLeQQz+LiokZHRwl+pG7d0Hf385I+I+kZSS9Jesrdj5nZfWZ2Y9oFAkWxvLysAwcOZF0GOlykY/rufsTdP+Luv+ru45V1X3D3mTptr2MvH1g1NjYW1P7kyZMpVQKs4opcIEWhvXhCruYF4iD0gZRt3rw5ctuQq3mBOAh9IGUHDx6M3DZkPH4gDkIfSFlIkA8MDKzfCGgCoQ+0kddeey3rEtDhCH2gBbZu3Zp1CYAkQh9oifvvvz9y2z179qRYCYrO3D2TFx4ZGfHZWbrzozhChlvO6u8S7c/MXnT32HOWsKcPtEjIIR7G2EdaCH2gRUIO8YR08wRCEPpAi4R03XzvvfdSrARFRugDLUQvHmSN0AdaKOQQz65du1KsBEVF6AMtFHKI5/jx4ylWgqIi9AGgQAh9oMUGBwcjt2UmLSSN0AdaLGT45HK5rKGhIcIfiSH0gRYLHT55bm6O+XORGEIfyEDoNIrMn4ukEPpABkKnUZSYPxfJIPSBnGD+XCSB0AdygvlzkQRCH8hIb29vUHvmz0USCH0gIw899FBQe3rvIAlMogJkaGBgIPK8uJs2bdLZs2dTrgjtjklUgBybn5+P3Pbtt99OsRIUBaEPAAVC6AM5wjSKaBahD2Qs5OrcqampFCtBEUQKfTPba2Yvm9kJM7unzuN/aWbHzewHZvasmUUfRhAouJCrc1dWVujFg6asG/pm1iXpQUnXSxqWdIuZDdc0+66kEXf/DUlPS/qHpAsFOlnINIoMvoZmRNnTv1bSCXd/xd3flfSEpJuqG7j7N919ubL4vKTtyZYJdLaQaRQZfA3NiBL6A5JOVS2frqxr5A5J32imKKBo4gy3DMTRneSTmVlZ0oik32nw+KikUYnBowAgC1H29Ocl7aha3l5ZdxEz2yPpgKQb3f2dek/k7lPuPuLuI/39/XHqBQA0IUrovyBpp5ldZWYbJd0saaa6gZldLemgVgP/9eTLBDpfyMlcSdqzZ09KlaCTrRv67n5e0mckPSPpJUlPufsxM7vPzG6sNPtHSb2S/tXMvmdmMw2eDkADISdzJenZZ5+lFw+CMeAa0Ea6u7u1srISuf3g4KBeffXV9ApC22HANaCDPP7440Ht6cWDUIQ+0EaYKAVpI/SBNhMyFg8QitAH2kzIWDyStHHjxpQqQSci9IE2FNJ989y5cylWgk5D6ANtKLT7JuPsIypCH2hDoSd0Jycn6bOPSAh9oE3t3r07qP3dd9+dUiXoJIQ+0KaOHj0a1H5xcTGlStBJCH2gjYV23+QQD9ZD6ANtLLT7ZrlcJvixJkIfaHOho2/u27eP4EdDhD7Q5kK7b66srHBSFw0R+kCbK5VK6urqCtpmcXGRvX3URegDORA6+qYkJk9HXYQ+kANxRt9k2GXUQ+gDOXHZZZcFb8MhHtQi9IGceOSRR4K3ufPOO1OoBHlG6AM5USqV1N3dHbTN0tISe/u4CKEP5Mhjjz0WvA0XbKEaE6MDOdPb26ulpaXg7bL6W0eymBgdKJiDBw/KzIK3u+KKK1KoBnlD6AM5UyqVdOjQoeDtzpw5IzPT0NAQh3sKjNAHcqhUKgWPt3/B3NycRkdHCf6CIvSBnAodb7/a8vIyV+wWFKEP5Nj09HTsbblit5gIfSDHSqWStm3bFnt7M+M4f8EQ+kDOzc/PN/0cc3NzKpfL2r9/fwIVoZ0R+kAHCJ1WsZHJyUmZmfr6+tjz71CEPtABJiYmNDw8nNjzLS4uqlwuv3/4x8z4X0CHiBT6ZrbXzF42sxNmdk+dxy8zsycrj3/HzIaSLhTA2o4dO5Zo8Ne68L+ATrtdfvnlFy339vaqt7f3onV9fX3as2ePuru7ZWbq7u6u+yN4+PBhDQ0NacOGDe17nsTd17xJ6pL0Y0kflrRR0vclDde02S/pocr9myU9ud7zfvzjH3cAyRsbG3NJ3FpwGxsbe/9zn56e9p6enose7+np8enp6UT/fSXN+jr5utZt3bF3zOy3JH3R3X+/snyvJLn7l6raPFNp820z65b0f5L6fY0nZ+wdID2HDx9WuVzOuoyO19XVpfPnz0uShoaG6naDHRwc1KuvvprYa7Zi7J0BSaeqlk9X1tVt4+7nJf1M0tbaJzKzUTObNbPZhYWFeBUDWFepVGqqDz+iWVlZef/+yZMn67ZptD4rLT2R6+5T7j7i7iP9/f2tfGmgcEqlktw99nANWF/1hPVXXnll3TaN1mclSujPS9pRtby9sq5um8rhnQ9JWkyiQADNOXr0qKanp7V16yX/+UaTRkdH378/Pj6unp6eix7v6enR+Ph4q8taU5TQf0HSTjO7ysw2avVE7UxNmxlJ+yr3Py3pv9Y6ng+gtUqlkt544w25Oz8ANTZt2nTR8ubNm7V58+aL1m3dulW7d+9+f8++q6tLY2NjmpiYeL9NqVTS1NSUBgcHZWYaHBzU1NRUrEnt0xRpEhUzu0HSl7Xak+dRdx83s/u0ehZ5xsw2STok6WpJP5V0s7u/stZzciIXAMI1eyI30oSb7n5E0pGadV+ouv+2pD+OWwQAoDW4IhcACoTQB4ACIfQBoEAIfQAoEEIfAAqE0AeAAiH0AaBAIl2clcoLm70l6eUYm35IqwO6xWlTb331upD7fZLeiFz1+rVFaVO7fq3lTqy/el1a9cf97tQur1V/3NrXqm+9xzu9/kbvpVO++9X3P+ruH4xedo1mxmVu5qaYY0JLmorbpt766nUh99uh/rWWO7H+mnWp1B/3uxNSf9zaqb/xukbvpVO++0l9/u6ey8M7X2uiTb31X2vifhxJ1r/WcifW32ztUZ4j7nendpn6o9cR5fGQv93q5U757kd9/XVleXhn1psYPyJr1J+tPNef59ol6s9as/Vnuac/leFrJ4H6s5Xn+vNcu0T9WWuq/sz29AEArZfHY/oAgJgIfQAoEEIfAAqkbUPfzDab2ayZ/UHWtYQys4+Z2UNm9rSZjWVdTygz+yMze9jMnjSz38u6nhBm9mEze8TMns66lqgq3/XHK595e82tF0EeP/NqOf++h2dNM538G1yE8Kik1yX9sGb9Xq1egXtC0j0Rnuc+SZ+X9AdJ19iK+ivbbJA0neP6r5D0SE5rf7qVn3sz70XSrZL+sHL/ySzrbubfIuvPPIH6W/p9T7j2yFmTRsG/Lema6oK1OrfujyV9WNJGSd+XNCzp1yX9R83tlyX9rlYnYL8tg9Bvuv7KNjdK+oakP81j/ZXt/knSNTmtPevQD3kv90r6zUqbr2RZd5z62+UzT6D+ln7fk6o9NGsizZEbwt2/ZWZDNauvlXTCK5Olm9kTkm5y9y9JuuTwjZldJ2mzVv8gzprZEXd/L+la60mi/srzzEiaMbOvS/pKehVf8rpJfP4m6e8kfcPd/yfdin8hqc++HYS8F0mnJW2X9D21ySHXwPqPt7a69YXUb2YvKYPveyOhn31o1rTqCzYg6VTV8unKurrc/YC7f06rb+DhVgX+GoLqN7PrzOwBMzuomgnlMxJUv6TPStoj6dNmdleahUUQ+tlvNbOHJF1tZvemXVygRu/l3yR9yswmldCl9impW3+bf+bVGn3+7fR9b6TRZx+cNYnv6SfJ3R/LuoY43P05Sc9lXEZs7v6ApAeyriMOd1+U1K5/uHW5+5Kk27OuI648fubVcv59f06BWdOqPf15STuqlrdX1uUF9Wcnz7XXyvt7of7sJFZ7q0L/BUk7zewqM9uo1ZO0My167SRQf3byXHutvL8X6s9OcrWncOb5q5J+IumcVo873VFZf4OkH2n1DPSBrM+QU3/2tXZS7Z32Xqi/c2tnwDUAKJC26B4GAGgNQh8ACoTQB4ACIfQBoEAIfQAoEEIfAAqE0AeAAiH0AaBACH0AKJD/B1Jj/RACtXAoAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.scatter(gene_info.mean_, gene_info.drop_rate, c='k')\n",
"plt.xscale('log')\n",
"plt.xlim(1e-4, 1e3);"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [],
"source": [
"NB_drop_prob = lambda mu, phi: ((1. / phi) / (mu + (1. / phi))) ** (1. / phi)"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEACAYAAABfxaZOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvFvnyVgAAH+pJREFUeJzt3X2cXVV97/HPLxOGkAklNAliEpKJbbiSSBWYYm/ptdgEGx4EUargGUQKHZhIBfUW5carlsoFHwsUZshQEJoJT1LqjQhSk4oiL+EyKKgJYkNgIAlXwpjwkAnkafWPtWeyz559zpyHfc4+D9/363VeZ+9z1p75ZTLznT1rr72WOecQEZHmMCHtAkREpHoU+iIiTUShLyLSRBT6IiJNRKEvItJEFPoiIk1EoS8i0kQU+iIiTUShLyLSRBT6IiJNZGJan3j69Omuvb09rU8vIlKXHn/88ZedczNKPT610G9vb2dgYCCtTy8iUpfMbLCc49W9IyLSRBT6IiJNRKEvItJEFPoiIk1EoS8i0kQU+iIiTWTc0Dezm83sJTP7VY73zcyuNbP1ZvYLMzs6+TJFRCQJhZzp3wIsyfP+icD84NEF9JZfVm0zszEPEZF6MG7oO+d+DPwuT5PTgH9x3iPAVDN7a1IF1ppcAV9I8BfyyyJXm7jXC33MmjWrvH+0iDSMJO7InQW8ENrfGLz2YgIfu+ZMBW4D9gJ7gsfINmeeCRMn7nu0tIxuf+Paa7kS2BU8dgM7gU+acW1vL7S2kjnvPM4A3gTeCB47gCPN+MNgezvwenB8oTZv3lz0XyPd3d309PQUdYyI1D5zzo3fyKwduNc5946Y9+4FrnLO/STYXwN81jk3Zo4FM+vCdwExZ86cYwYHy7qbOBUzzdicdhH4Xxjbg8drwCuhx6vB8++AIeDl0PPLwBaCX1Il6u/vJ5PJlPERRKRUZva4c66j5OMTCP3lwIPOuduD/aeB451zec/0Ozo6XD3OvTPJjPfi+8VaQo8JwF233QZ79sDu3f4R2v7MJZewH2Q9WoPnT3Z1wc6drLzlFvYH9gcmhR4HBM+TgbbgUc6faHuBl4DNkcdzwIbgsRkY/ztjH/0iEKmOWgj9k4GLgJOAdwPXOueOHe9j1mvo5+smyfe1LOS4YrpgWvHhPwU4EPg94KDI4/eBacD00PN0YAbjX8x5E3gWeAZYC/wK+CXwVPBeIdra2li+fLl+GYgkqOKhb2a3A8fj8+K3wBfxJ6g4524wn1TX4Uf4DAPnxnXtRNVr6EN8OBfwdcz5XimhX44W4C3ATOCtwfNsoB14W/A4NMexu4H/BJ4Efgo8DDxBYd1FhZxgiEh+VTnTr4R6Dv1SFfLLIlebag8LbQPm4cfhvgM4MnjMx//SCHsdeBT/C+BHwEP4i9W56CKxSOkU+jJGJX9BTAKOAI4GjgP+DP+LIOwV4PvAKuB+YOs4H1N/AYgUrtzQ1zQMDcg5V9SjGG8APwduAv4aOBzfVfRB4Bv4fv+DgI8AK/EXjP8Df9Fnao6POXI/wdKlS0v414pIMXSmL3mV8ldDO/B+4FTgzwkuAOHvM/g20IfvCspl0aJFrF69uujPK9IMdKYvFVXKXwbPAf8EnIAfKXQW8AB+6OnHgJ/gRwRdAhwcc/yaNWswM1auXJnMP0JERulMX8pW6F8D7cD5+G6hkXk6XsV3C/0j/iazOOrzF9lHZ/qSukL/EngO+DwwBzgdWI2/v+Dv8TeE/U/8XwNRZsbChQuTLlukKSn0JXHj/QLYDXwH3/3zHvwQz+nA14D1wFL2XQcYsW7dOnX5iCRAoS8VNd4vgIfwwb8EeBx/o9j1+BFCcQszdHZ2cvDBcVcCRKQQCn2pmnzh/wDQAXwIf8fvQvwNX19i7Fn/tm3btIaBSIkU+lJ1+cL/HuCdwNX4SeW+iA//I2Paamy/SPEU+pKaXOG/A/gUfoz/BuAoYAD4X4ydAqK3t1dn/SJFUOhL6nKF/4+BP8Kvv9kKXAHch59VNErBL1IYhb7UDOccCxYsyHptO340zwn4KV7fh7/4G7cApIJfZHwKfakpa9eujT3rXw38CfBrfJ//I/i/AqIU/CL5KfSlJjnn2G+/7HE7zwF/ip++eTb+jP+EmGPNjPb2do3pF4mh0JeatXPnzjFn/VvxXTy34+/mvQ84N+bYwcFBurq6FPwiEQp9qXnR4N8JZIAr8cM6bwY+GXPc8PAwy5Ytq3h9IvVEoS91IRr8Dj+E88Jg/xr8HP5Rg4ODlS1MpM4o9KVuOOc44IDsKdmWA38XbP8L8N6Y43RxV2Qfhb7UleHhYaZOzV6D6+v4O3hb8RO5vTPmuJaW6G1dIs1JoS91Z+vWsavufhq4E39x935gbuT9vXv3aqI2ERT6Uqfi+vg/BvwQv0DL94Hfjxyzbds2zdUjTU+hL3UrblTPB4AngbcD9wKTIsf09vZWpTaRWqXQl7oWDf5XgROBQeC/4/v7o3RhV5qZQl/qXjT4XwROw5/5fwI4OeYYdfNIs1LoS0OIBv+T+HH84G/eOiTSvre3V8EvTUmhLw2jv78/a/+bwBp84N8c017BL81IoS8NI5PJsGjRotF9B5wD/A7fxRMX77qwK81GoS8NZfXq1Vn7m4CuYPvrwBExx2hSNmkmCn1pONH+/X8FvgUcANyGv3M3rLOzszqFidSAgkLfzJaY2dNmtt7MPhfz/hwz+6GZ/dzMfmFmJyVfqkjhosH/SeAZ4F3Al2PaaxinNItxQ9/MWoDr8cOfFwBnmdmCSLPPA3c5544CzgR6ki5UpFjhpRdfx0/HvBs/Qdsfx7TXRV1pBoWc6R8LrHfObXDO7QTuwA+DDnP4aU8ADgI2J1eiSGnWrl2btf8o+27WujqmvS7qSjMoJPRnAS+E9jcydl3qLwGdZrYRv5jR3yZSnUiZot08/wf4//hlF+Pm39fZvjS6pC7kngXc4pybDZwErDCzMR/bzLrMbMDMBrZs2ZLQpxbJLzwV82v4vkiAr6K5eaT5FBL6m4DDQvuzg9fCzgPuAnDO/RT/szQ9+oGcc33OuQ7nXMeMGTNKq1ikSNGpmL8FPAHMAT4T037y5MlVqEokHYWE/mPAfDObZ2at+Au1qyJtngcWAZjZEfjQ16m81IzwTVt7gU8F25/DT8UctmPHDnXzSMOyaJ9nbCM/BPNqoAW42Tl3hZldDgw451YFo3luBKbgL+pe6pz793wfs6Ojww0MDJT9DxApVHRY5j3A6fgz/7+OaV/Iz4ZItZnZ4865jpKPT+sbW6Ev1bZy5cqsG7H+AFiHv1nrGOBnkfYKfalF5Ya+7siVppHJZLIWVn8GuCbYjhvCqb59aUQKfWkqw8PDWftfBl4C/gdwRqTtjh07NC+PNByFvjSdcLfNq8AXgu2vAvtF2mpeHmk0Cn1pev8MPAXMw99wItLIFPrSlLq7u0e39wBfCbYvBaJTr2kyNmkkCn1pSj092XMC3oafa2QhcEoaBYlUiUJfmlb4bH8XfnlFgM/GtNXZvjQKhb40rejZ/o34pRWPCx5RuktXGoFCX5paeHqG7cB1wfaYlYLQZGzSGBT60tSia+r+EzCM79dfGNNe4/al3in0pemFV9h6Gbgp2L40pu3HP/7xKlQkUjkKfWl60RW2voFfVvGj+OmXw3bv3l2lqkQqQ6EvAsycOXN0exC4E5gIfDqmrbp4pJ5plk2RQHhY5pHAL/AXd+cCQ5G2moFT0qJZNkUSEg79X+IXe24DLoppq7N9qVcKfZHAihUrsvavCp4/gZ9zP0wTsUm9UuiLBDKZTNb+Q8CTwAzg1Jj2OtuXeqTQFwlpaWnJ2v/n4Pn8mLbnnHNOxesRSZpCXyTk1ltvzdrvB94ATsBf0A3bs2dPlaoSSY5CXyQkk8nQ398/ur8NuBv/gxK3eLq6eKTeaMimSIzwSJ4/Bx4ENgLt+Pn3wzR8U6pJQzZFKuxHwH8Cs4G/jHlfZ/tSTxT6IjHCXTyQ/4Lu2WefXfF6RJKi0BeJER2+eSt+oZX3A4dG2jrndLYvdUOhL1KA3wLfxc/HEzdQ89xzz1XwS11Q6IvkMHdu9iDNcBdPdPHEXbt2sWzZsmqUJVIWhb5IDldccUXWKJ4H8Iun/yF+RE/U888/X6XKREqn0BfJIZPJZM3Hsxe4OdiOu6A7efLkapQlUhaFvkge0Qu638KH/4eAgyNtt2/fXqWqREpXUOib2RIze9rM1ptZ3JrRmNmHzWydma01s9uSLVOkNgwCPwAmAXHzbOpirtS6cUPfzFqA64ETgQXAWWa2INJmPnAZcJxzbiFwSQVqFUlF9ILujcHzuTFtOzs7mTBhAu3t7foFIDWpkDP9Y4H1zrkNzrmdwB3AaZE2fwNc75zbCuCceynZMkXSc8UVV2Tt3wu8AhwFzI9p75xjcHCQrq4uBb/UnEJCfxZ+0MKIjcFrYYcDh5vZw2b2iJktSapAkbRF+/XfBL4TbH8kz3HDw8Maxik1J6kLuRPxJz3HA2cBN5rZ1GgjM+syswEzG9iyZUtCn1qk8rq7u7P27wye84U+aBin1J5CQn8TcFhof3bwWthGYJVzbpdz7lngN8T85euc63POdTjnOmbMmFFqzSJV19PTk7W/Gr9Y+juAhXmOmzNnTgWrEileIaH/GDDfzOaZWStwJrAq0uY7+LN8zGw6vrtnQ4J1iqQuPAnbLuCeYDvX2f7kyZPHXA8QSdu4oe+c2w1chL8h8SngLufcWjO73MxGlg59ABgys3XAD4G/c84NVapokTRE+/bvCJ7PzNG+r69vzDEiadMiKiJFCE/L0ILv53wLcDTw80hbLa4ilaBFVESqaNGiRaPbe/BLKUJ8F4+Ga0ot0pm+SJHCZ/t/BjwEPAfMi2mrs31Jms70RaosPHzzYfatnfvulOoRKYZCX6RI4eGbDrgr2I67oKsuHqk1Cn2RMo3cqPVhxv5AdXbGTcsmkh6FvkgJZs6cObr9/4BngZn4Pv4one1LLVHoi5Rg06bsm9LzTctwwQUXVLwekUIp9EUSMBL6Z+DH74dt375d0y1LzVDoi5Qo3MXzBPA0cAjw3pi2mm5ZaoVCX6REubp4PpznGE23LGlT6IskZOTu3FPJ/4Ol6ZYlTQp9kTK0tbWNbv8SP4rnLeS/UUvTLUuaFPoiZVi+fHnW/siKWh/I0V7TLUvaFPoiZYhOnTxe6Gu6ZUmbQl8kQQ/jV9Q6HHh7zPsKfEmbQl+kTNOmTRvd3gN8N9iOO9tfunRpNUoSyUmhL1Kma665Jmt/pIvntJi2vb29Fa9HJB+FvkiZol02PwB2AH8CvDWNgkTyUOiLJCC8otYw8O/B9qkxbXVHrqRJoS+SgNWrV2ft5+viufDCCytej0guCn2RhEydOnV0+178Rd1FwIGRdq+//noVqxLJptAXScjWrVtHt1/GD99sBU6MaasuHkmLQl+kQvJ18Zx99tnVLEVklEJfJEFmNrr9f4Pnk4H9Iu2cc9UqSSSLQl8kQStWrBjd3oCfhO0g4PiU6hGJUuiLJCg6Zn/kbD+ui0d350oaFPoiFRTu17fIe7o7V9Kg0BdJ2KRJk0a3Hwc2ArOBY9IqSCREoS+SsDfffDNrP18Xj0i1KfRFEhZdGWtV8Pz+mLYLFy6seD0iYQWFvpktMbOnzWy9mX0uT7sPmZkzs47kShSpL9GVsR4EXgPeCUQXSly3bl11ihIJjBv6ZtYCXI+/sXABcJaZLYhpdyBwMfBo0kWK1JNMJpM1Xn8n+yZgOyWmvUbxSDUVcqZ/LLDeObfBObcTuIP47sl/AL4CvJFgfSJ1KTqp2sjCKnFdPBrFI9VUSOjPAl4I7W8MXhtlZkcDhznnvpdgbSJ1q6enJ2v/PmAv8F5gShoFiQTKvpBrZhOAbwKfKaBtl5kNmNnAli1byv3UIjWtpaVldHsL8AiwP3BCWgWJUFjobwIOC+3PDl4bcSDwDuBBM3sOv2DQqriLuc65Pudch3OuY8aMGaVXLVIHurq6svbzdfFMnz5dM29KVRQS+o8B881snpm1AmeybxQazrlXnHPTnXPtzrl2/AnNqc65gYpULFInol08I6F/MmN/8IaGhujq6lLwS8WNG/rOud3ARcADwFPAXc65tWZ2uZnFrQYnIjHWAs8Ch+BHR0QNDw+zbNmy6hYlTaegPn3n3H3OucOdc3/gnLsieO0LzrlVMW2P11m+iNfd3Z21n6+LB+D555+vaD0iuiNXpIJydfHkCv3o3bwiSVPoi1RYW1vb6PaPgFeBI4G5MW2jd/OKJE2hL1Jhy5cvH93ehb84BvFn+9H5+EWSptAXqbBokOfr4pk1a1bMqyLJUeiLVNl9wB78EooHRt7bvHlz1euR5qLQF6mCadOmjW4PAT8FWoH3pVWQNC2FvkgVXHPNNVn7+bp4Fi9eXPF6pHmZcy6VT9zR0eEGBjScX5pHeLrlI4B1+Dl5DsVPxhaW1s+l1D4ze9w5V/KaJTrTF6mScBfPU8AzwAz8ZFVRmmNfKkWhL1Il0S6efGvnhod5iiRJoS9SJdGhm98Jnk+Pabt3b7TDRyQZCn2RKgp38TyM79Ofj1+HVKQaFPoiVRTu4tnLvjnKPxDTduHChdUoSZqMQl+kinJ18cSF/rp16ypejzQfhb5IilYD24E/xi9JJ1JpCn2RKps7d9/8mm8A3w+240bxaCUtSZpCX6TKotMn/1vwHNfF09nZSXt7u8JfEqPQF6myaL/+94Dd+AnYpsa0Hxwc1Pq5khiFvkgKwssobgMeBCYCp+Ror/VzJSkKfZEURJdRzNfFM0Lr50oSFPoiNWBkvP4SYFKONlo/V5Kg0BepARuBx4A24IQcbbR+riRBoS+SkilTpmTt57tRC7R+riRDoS+SkhtuuCFrf6Rf/1SgJaa9Ru9IEhT6IinJZDLMnDlzdP8p4DfAdOBPY9qff/75VapMGplCXyRFmzZtytrPN93yG2+8UfF6pPEp9EVqSCFDN0XKodAXqSGPAi8C84CjY97XMopSLoW+SMrCd+c64NvBdtxYnb6+vmqUJA2soNA3syVm9rSZrTezz8W8/2kzW2dmvzCzNWY2N+7jiMhY0btz+4PnjzJ2FM+ePXs0ikfKMm7om1kLcD1wIn5Vt7PMLLq628+BDufcHwF3A19NulCRRhZeRvEx4GngUGBRTFtNviblKORM/1hgvXNug3NuJ3AHkam/nXM/dM4NB7uPoPUgRIoSXkYR9p3tnx3TVpOvSTkKCf1ZwAuh/Y3Ba7mcB9xfTlEizSZ6t+3Iefzp+KkZogYHBytdkjSoRC/kmlkn0AF8Lcf7XWY2YGYDW7ZsSfJTizSUZ4Gf4AM/bsy+SKkKCf1NwGGh/dnBa1nMbDGwDDjVOfdm3AdyzvU55zqccx0zZswopV6RpjHSxdOZahXSaAoJ/ceA+WY2z8xagTPZNxMsAGZ2FLAcH/gvJV+mSOMLX8wFuAvYCSzGX9SNWrx4cRWqkkYzbug753YDFwEP4KcHucs5t9bMLjezU4NmXwOmAN82syfMbFWODyciOUQv5m7FL6XYApwV037NmjUaxSNFM+dcKp+4o6PDDQwMpPK5RWrVxIkT2bNnz+j+6cA9wM+AY2Laz507l+eee646xUlNMLPHnXMdpR6vO3JFasitt96atX8f/oz/aPxNMlEaxSPFUuiL1JDo0M032Tctgy7oShIU+iI1JjwXD8CK4DkDWNWrkUaj0BepMdG5eB4GngPmAO+Jad/a2lr5oqRhKPRFalB4+KZj3x26cdMy7Nq1qxolSYNQ6IvUoFxz8fwVfmx0lObZl0Ip9EVqUPSC7q+Bh4DfA7pi2vf29mrMvhREoS9SoxYtyp5Y+SvB86eBuF78iy++uNIlSQNQ6IvUqNWrV2ftfw/4JX6K27i+/aGhoSpUJfVOoS9Sw6LDN68Kni8l/odXXTwyHoW+SA2LDt+8Ez/t8uHAB2Pad3Z2KvglL4W+SI0LD9/cw77FKsYsVh0455xzFPySk0JfpMZFh29+C/gtfgK2E2La79mzRxd1JSeFvkiNy2QytLS0jO6/AVwdbOc62x8aGtLZvsRS6IvUgejsmz3AK8BfAMfmOEaLp0schb5IHYjerPUq0Bts5zrb17TLEkehL1In9t9//6z9q/FdPacDb89xjLp4JEqhL1Inbrrppqz93+Iv6gJcluOYCy64oJIlSR1S6IvUiUwmw8SJE7Ne+zqwC/gY8Jcxx2zfvl1n+5JFoS9SR2655Zas/Q3AF4LtW4FDYo7RDVsSpoXRRerMlClT2L59++j+BGANcDx+fp5TchyX1s+6JEsLo4s0meXLl2O2b+HEvfgJ2H4HnAxclOO4gw8+uPLFSc1T6IvUmUwmw4oVK7Je2wj8TbD9NeDImOO2bduGmdHe3q7uniam0BepQ5lMZsx8+/cAfcAk4PbgOc7g4CBdXV0K/ial0BepU9H59gE+hV9layF+ZE8uw8PDumO3SSn0RepYf39/1v4wcBawE/gE8dMvj9Adu81JoS9SxzKZDDNnzsx67Qn23ax1F/C/yf2Dbmbq528yCn2ROrdp06Yxr/0j8PeAAZcDDxA/hn/E4OAgnZ2dLF26tCI1Su1Q6Is0gOiyig74Ev4u3ZeAxfi/AI4f5+P09vZiZkyfPl1n/g1KoS/SAHp6eliwYMGY11cD7wIeBN4a7H8eaBnTMtvQ0BCdnZ2j3T9mpr8CGkRBoW9mS8zsaTNbb2ZjZnI1s/3N7M7g/UfNrD3pQkUkv7Vr18YG/4v4M/3L8d09/wBsxg/vPBFoLfDjj/wV0GiPAw44IGt/ypQpTJkyJeu16dOns3jxYiZOnIiZMXHixNhfgitXrqS9vZ0JEybU7nUS51zeB/6k4BngbfjvjyeBBZE2S4Ebgu0zgTvH+7jHHHOME5HkdXd3O3wPz5jHCeB+Dc6FHq+AWwnuI+COAvcWcJbjeD2yH93d3aNf9/7+fjd58uSs9ydPnuz6+/sT/f8FBtw4+ZrvkT1lX7xjgfXOuQ0AZnYHcBqwLtTmNHwXIsDdwHVmZkGBIlJFPT09HHfccXR2do557wf4ufcX4odzfhDf/fPR4DFiN/4vhM34awJv4oeBjjx2BW3yJWIzmHDDDTBlCgCvLV/OF4eHsxsMD/Pa0qXw5JMpVBdv3AnXzOwMYIlz7vxg/2zg3c65i0JtfhW02RjsPxO0eTnysbqALoA5c+Yco3HCIpWzcuXK2OCPmodfiOUvgFnBY0ZlS5MyGJQ14VohZ/qJcc714bsS6ejoaJaTAZFUZDIZMpkMixcvZs2aNTnbPQt8M3iMaAUOBWbifwG0xjwm4q8RxD2axQQzrrrqKgCuvPJKtm7bNqbNwVOnctlluZa5KcFnP1vW4YWE/ibgsND+7OC1uDYbzWwicBAwVFZlIpKI1atXs3LlSi6++GKGhgr7sdwJPB88JLfuCy+ESy8FYM6sWXy5q4vhUBfP5MmT6bvuOoiscVyWMkO/kNE7jwHzzWyembXiL9SuirRZBZwTbJ8B/If680VqRyaT4eWXX8Y5R39/P9OmTUu7pJoxaVL21HRtbW20tbVlvTZt2jQWLVpES4sf7NrS0kJ3dzc9PT2jbTKZDH19fcydOxczY+7cufT19Y1Z1D5tBS2iYmYn4ddhbgFuds5dYWaX468irzKzScAK4Cj8tN5njlz4zUWLqIiIFK/cRVQK6tN3zt0H3Bd57Quh7TeAvyq1CBERqQ7dkSsi0kQU+iIiTUShLyLSRBT6IiJNRKEvItJEFPoiIk1EoS8i0kQKujmrIp/Y7DXg6RIOPQh4pcQ2ca+HXytmezqQNaFcgZKsP99+I9Yffq1S9Zf6vRPdz1d/qbXnq2+89xu9/lz/lkb53g9v/zfn3IGFlx1RzrzM5TwocU5ooK/UNnGvh18rZrsW6s+334j1R16rSP2lfu8UU3+ptav+3K/l+rc0yvd+Ul9/51xddu98t4w2ca9/t4ztUiRZf779Rqy/3NoL+Rilfu9E91V/4XUU8n4xP7vh/Ub53i/0848rze6dAVfG/BFpU/3pquf667l2UP1pK7f+NM/0+1L83ElQ/emq5/rruXZQ/Wkrq/7UzvRFRKT66rFPX0RESqTQFxFpIgp9EZEmUrOhb2ZtZjZgZqekXUuxzOwIM7vBzO42s+606ymWmX3AzG40szvN7H1p11MMM3ubmd1kZnenXUuhgu/1W4OveW2trVeAevyah9X593vxWVPOIP8cNyHcDLwE/Cry+hL8Hbjrgc8V8HEuBy4FTkm6xmrUHxwzAeiv4/oPBm6q09rvrubXvZx/C3A28P5g+8406y7n/yLtr3kC9Vf1+z3h2gvOmkoU/B7g6HDB+LV1nwHeBrQCTwILgCOBeyOPQ4AT8AuwfzyF0C+7/uCYU4H7gY/WY/3Bcd8Ajq7T2tMO/WL+LZcB7wra3JZm3aXUXytf8wTqr+r3e1K1F5s1Ba2RWwzn3I/NrD3y8rHAehcslm5mdwCnOeeuBMZ035jZ8UAb/gdih5nd55zbm3StcZKoP/g4q4BVZvY94LbKVTzm8ybx9TfgKuB+59zPKlvxPkl97WtBMf8WYCMwG3iCGulyLbL+ddWtbnzF1G9mT5HC93suxX7ti82aan2DzQJeCO1vDF6L5Zxb5py7BP8PuLFagZ9HUfWb2fFmdq2ZLSeyoHxKiqof+FtgMXCGmV1YycIKUOzXfpqZ3QAcZWaXVbq4IuX6t9wDfMjMeknoVvsKia2/xr/mYbm+/rX0/Z5Lrq990VmT+Jl+kpxzt6RdQymccw8CD6ZcRsmcc9cC16ZdRymcc0NArf7gxnLObQfOTbuOUtXj1zyszr/fH6TIrKnWmf4m4LDQ/uzgtXqh+tNTz7VH1fu/RfWnJ7HaqxX6jwHzzWyembXiL9KuqtLnToLqT0891x5V7/8W1Z+e5GqvwJXn24EXgV34fqfzgtdPAn6DvwK9LO0r5Ko//VobqfZG+7eo/satXROuiYg0kZoYHiYiItWh0BcRaSIKfRGRJqLQFxFpIgp9EZEmotAXEWkiCn0RkSai0BcRaSIKfRGRJvJfZnr1mfeHtYYAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.scatter(gene_info.mean_, gene_info.drop_rate, c='k')\n",
"plt.plot(mm, NB_drop_prob(mm, phi), c='r', lw=2)\n",
"plt.xscale('log')\n",
"plt.xlim(1e-4, 1e3);"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment