Skip to content

Instantly share code, notes, and snippets.

@spektom
Created February 5, 2018 06:28
Show Gist options
  • Save spektom/05c5ce564b75a3973fcd6214fa321caa to your computer and use it in GitHub Desktop.
Save spektom/05c5ce564b75a3973fcd6214fa321caa 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": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x10eaa3590>"
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAW8AAAD3CAYAAADSftWOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHSlJREFUeJzt3Xt0lNW9PvDnfeeSydxIomPBRiBRKEaOSpsV9DRG+7M0\n7dGWgkCAGvVY69JK29hSQRYkpV6Qn0tam64WWteStYhoqVVWV09Xu2r0GC415aDUNqAekHBPCCGQ\nuSRze/f5YzITBpKZyWQu7rfP559kru93J+F5N3v23q8ihBAgIiKpqPkugIiIxo7hTUQkIYY3EZGE\nGN5ERBJieBMRSciYi4P09LhzcZg4xcVW9PX5cn7cXNBz2wC2T3ZsX+a4XI5RH9Ntz9toNOS7hKzR\nc9sAtk92bF9u6Da8iYj0jOFNRCQhhjcRkYQY3kREEmJ4ExFJiOFNRCQhhjcRkYQY3kREEkq6wjIY\nDGLlypU4ceIEVFXFE088AaPRiJUrV0JRFEybNg1NTU1QVZ4HiIhyJWl4v/322wiFQnjllVewa9cu\n/PSnP0UwGERDQwNmz56NxsZGtLa2Ys6cObmoN2V/+msn3J7BuPtuu/HT+SmGiCjDkoZ3WVkZwuEw\nNE2Dx+OB0WjEvn37UFVVBQCoqanBrl27EoZ3cbE190tKD/bCYbfE3ZVonwDZ6KktI2H75Mb2ZV/S\n8LZarThx4gS+8pWvoK+vDxs3bsSePXugKAoAwGazwe1OvPFUvjapubjnnY8NsrLB5XLopi0jYfvk\nxvZl9lijSRremzdvRnV1NX7wgx/g1KlTuPfeexEMBmOPe71eOJ3OzFRKREQpSfopo9PphMMRSf8J\nEyYgFAqhoqIC7e3tAIC2tjZUVlZmt0oiIoqTtOd93333YdWqVVi6dCmCwSAeffRRzJw5E2vWrMGG\nDRtQXl6O2traXNRKRERDkoa3zWbD888/f8n9LS0tWSmIiIiS4+RsIiIJMbyJiCTE8CYikhDDm4hI\nQgxvIiIJMbyJiCTE8CYikhDDm4hIQgxvIiIJMbyJiCTE8CYikhDDm4hIQgxvIiIJMbyJiCTE8CYi\nkhDDm4hIQgxvIiIJJb2SzmuvvYbXX38dAOD3+3HgwAFs2bIFTz31FAwGA6qrq7Fs2bKsF0pERMOS\nhvf8+fMxf/58AMDatWtx1113oampCc3Nzbjqqqvw4IMPYv/+/aioqMh6sUREFJHysMk//vEPHDx4\nEHfccQcCgQAmT54MRVFQXV2N3bt3Z7NGIiK6SNKed9SmTZvwyCOPwOPxwG63x+632Ww4duxYwtcW\nF1thNBrSrzIdB3vhsFvi7nK5HLmtIYv01JaRsH1yY/uyL6Xw7u/vx+HDh3HTTTfB4/HA6/XGHvN6\nvXA6nQlf39fnG1+VaXJ7BuNu9/S481JHprlcDt20ZSRsn9zYvsweazQpDZvs2bMHN998MwDAbrfD\nZDLh6NGjEEJg586dqKyszEylRESUkpR63ocPH0ZpaWns9tq1a7F8+XKEw2FUV1fjhhtuyFqBRER0\nqZTC+4EHHoi7feONN2Lbtm1ZKYiIiJLjIh0iIgkxvImIJMTwJiKSEMObiEhCDG8iIgkxvImIJMTw\nJiKSEMObiEhCDG8iIgkxvImIJMTwJiKSEMObiEhCDG8iIgkxvImIJMTwJiKSEMObiEhCDG8iIgkx\nvImIJJTSZdA2bdqEN998E8FgEEuWLEFVVRVWrlwJRVEwbdo0NDU1QVV5HiAiypWkidve3o733nsP\nL7/8MrZs2YKuri6sW7cODQ0N2Lp1K4QQaG1tzUWtREQ0JGl479y5E9OnT8cjjzyChx56CLfddhs6\nOjpQVVUFAKipqcHu3buzXigREQ1LOmzS19eHkydPYuPGjTh+/DgefvhhCCGgKAoAwGazwe12J3yP\n4mIrjEZDZipO1cFeOOyWuLtcLkdua8giPbVlJGyf3Ni+7Esa3kVFRSgvL4fZbEZ5eTkKCgrQ1dUV\ne9zr9cLpdCZ8j74+3/grTYPbMxh3u6cn8UlGFi6XQzdtGQnbJze2L7PHGk3SYZPPfe5z2LFjB4QQ\n6O7uxsDAAG6++Wa0t7cDANra2lBZWZm5aomIKKmkPe8vfOEL2LNnDxYsWAAhBBobG1FaWoo1a9Zg\nw4YNKC8vR21tbS5qJSKiISlNFXzssccuua+lpSXjxRARUWo4OZuISEIMbyIiCTG8iYgkxPAmIpIQ\nw5uISEIMbyIiCTG8iYgkxPAmIpIQw5uISEIMbyIiCTG8iYgkxPAmIpIQw5uISEIMbyIiCTG8iYgk\nxPAmIpIQw5uISEIMbyIiCaV0GbR58+bBbrcDAEpLS1FXV4ennnoKBoMB1dXVWLZsWVaLJCKieEnD\n2+/3QwiBLVu2xO6bO3cumpubcdVVV+HBBx/E/v37UVFRkdVCiYhoWNLw/uCDDzAwMID7778foVAI\n3/nOdxAIBDB58mQAQHV1NXbv3p0wvIuLrTAaDZmrOhUHe+GwW+Lucrkcua0hi/TUlpGwfXJj+7Iv\naXhbLBZ885vfxMKFC9HZ2YlvfetbcDqdscdtNhuOHTuW8D36+nzjrzQNbs9g3O2eHnde6sg0l8uh\nm7aMhO2TG9uX2WONJml4l5WVYcqUKVAUBWVlZXA4HDh37lzsca/XGxfmRESUfUlnm7z66qt45pln\nAADd3d0YGBiA1WrF0aNHIYTAzp07UVlZmfVCiYhoWNKe94IFC/D4449jyZIlUBQFTz/9NFRVxfLl\nyxEOh1FdXY0bbrghF7USEdGQpOFtNpvx3HPPXXL/tm3bslIQERElx0U6REQSYngTEUmI4U1EJCGG\nNxGRhBjeREQSYngTEUmI4U1EJCGGNxGRhBjeREQSYngTEUmI4U1EJCGGNxGRhBjeREQSYngTEUmI\n4U1EJCGGNxGRhBjeREQSSim8e3t7ceutt+LQoUM4cuQIlixZgqVLl6KpqQmapmW7RiIiukjS8A4G\ng2hsbITFYgEArFu3Dg0NDdi6dSuEEGhtbc16kUREFC9peK9fvx6LFy/GFVdcAQDo6OhAVVUVAKCm\npga7d+/OboVERHSJhBcgfu2111BSUoJbbrkFv/rVrwAAQggoigIAsNlscLvdSQ9SXGyF0WjIQLmJ\nhTWB7//0bcyYUoyyKyfAYbfEPe5yObJeQ67oqS0jYfvkxvZlX8Lw/t3vfgdFUfDXv/4VBw4cwIoV\nK3D27NnY416vF06nM+lB+vp84680Bec8fnx84jw+PnEed1aXocRujnu8pyf5iUYGLpdDN20ZCdsn\nN7Yvs8caTcJhk5deegktLS3YsmULrr32Wqxfvx41NTVob28HALS1taGysjKz1Y5DvzcQ+/7td08g\nGOKHqUSkT2OeKrhixQo0Nzejrq4OwWAQtbW12agrLf2+SHjbC01w+wJ4/9CZPFdERJQdCYdNLrRl\ny5bY9y0tLVkpZrzc3iAA4Kv/PhX/9c4R7O/sw9RJTlzmtCR5JRGRXHS1SOf80LDJ5RMsuHVWKYQA\n9n3E3jcR6U/KPW8ZuIeGTRw2MyZPVFFYYIgNpRAR6Ymuet7RDyydtsgsE4vZiMFAOJ8lERFlhb7C\n2xcZ83ZaTQAAi9mAYEhDmEv4iUhn9BXe3gDMJhUWc2Q0qMAcWRjE3jcR6Y2+wtsXgNM6vDCncCjE\nGd5EpDe6CW8hBNy+QGy8G4gMmwDAoJ/hTUT6opvwHvCHEAqLuJ53dNjEHwzlqywioqzQTXifj800\nMcXuY8+biPRKN+HtHppp4rBeOGwSGfMe4Jg3EemMbsL74jnewHDP28/wJiKd0U94D62kdFpH+MAy\nwDFvItIX/YT3CD1vk1GFqiicKkhEuqOf8L5odSUAKIoCS4GB4U1EuqOb8HaP0PMGIkMnHDYhIr3R\nTXif9wWgKgpshaa4+y1mA0JhgVCY+5sQkX7oJrzd3gAcVhPUoYsjR1m4RJ6IdEg34d3vC8TN8Y6y\ncHMqItKhpBdjCIfDWL16NQ4fPgxFUbB27VoUFBRg5cqVUBQF06ZNQ1NTE1Q1f+eBYCiMAX84bnVl\nFKcLEpEeJQ3vt956CwDwyiuvoL29HT/5yU8ghEBDQwNmz56NxsZGtLa2Ys6cOVkvdjT9Q9euvPjD\nSgAoiA6bcIk8EelI0u7yF7/4RTzxxBMAgJMnT8LpdKKjowNVVVUAgJqaGuzevTu7VSYx0gKdqMJo\nzzvI8CYi/UjpGpZGoxErVqzAX/7yF/zsZz/Drl27oAx9MGiz2eB2uxO+vrjYCqPRMP5qR9HZ4wUA\nTLrCAZfLEbnzYC8cdguKh3rcmsDwYzqgp7aMhO2TG9uXfSlfgHj9+vVYvnw5Fi1aBL/fH7vf6/XC\n6XQmfG1fny/9ClNw7NR5AIAqNPT0DJ9I3J5BaOHw0Pf+uMdk5nI5dNOWkbB9cmP7Mnus0SQdNtm+\nfTs2bdoEACgsLISiKJg5cyba29sBAG1tbaisrMxQqemJLo2fMMKYN6cKEpEeJe15f+lLX8Ljjz+O\nb3zjGwiFQli1ahWuvvpqrFmzBhs2bEB5eTlqa2tzUeuoRtoONspoUGBQFc42ISJdSRreVqsVzz//\n/CX3t7S0ZKWgdCTqeSuKMrREnj1vItIPXSzSic42cVgvnecNRIZOBgNhCCFyWRYRUdboI7y9ARQW\nGGEaZUaLpcCAsCbg53RBItIJfYS3Lxi3FezFLCZD7HlERHogfXhrmoDbFxhxdWWUpSAS3tFtY4mI\nZCd9eHsGgxBi5NWVUdEl8m72vIlIJ+QP79g0wdGHTaJL5KMfbBIRyU7+8B6IhLc90Zj3UHi7Gd5E\npBPSh3d0KMReyGETIvrXIX14ewYivWl74ejrjSwcNiEindFBeCfveQ8Pm7DnTUT6oJvwTvSBpdGg\nwmhQOFWQiHRD/vAe6k1ffNX4i1nMRpznsAkR6YT84R3teScJ78ICI9zeIDSN+5sQkfx0Ed4GVYmN\na4+msMAATYhY2BMRyUz68HYPBGEvNMUuyzaawoLIbJTzHPcmIh2QPrw9vmDCBTpRsfD2+JM8k4jo\nk0/q8A5rGnz+UNLxbiAybAKw501E+iB1eHsHIpc2s6cS3mYOmxCRfiS8DFowGMSqVatw4sQJBAIB\nPPzww7jmmmuwcuVKKIqCadOmoampCaqan3OAO7ZAJ/Vhk3McNiEiHUgY3r///e9RVFSEZ599FufO\nncPXv/51zJgxAw0NDZg9ezYaGxvR2tqKOXPm5KreON4UNqWKioZ3P3veRKQDCcP7y1/+cuzK8EII\nGAwGdHR0oKqqCgBQU1ODXbt2JQ3v4mIrjKNcomw8Dna5AQATXQ64XI6LHuyFw26J3bRZBRQF8AXC\nlz5XQnpoQyJsn9zYvuxLGN42mw0A4PF48N3vfhcNDQ1Yv359bFqezWaD2+1OepC+Pl8GSr3Uia7+\nyDfhMHp6Lq3D7RmMu+0oNKGnb2DE58rE5XJI34ZE2D65sX2ZPdZokg5Wnzp1Cvfccw/mzp2Lr371\nq3Hj216vF06nMzNVpiGVTaku5LQVoN/LMW8ikl/C8D5z5gzuv/9+/PCHP8SCBQsAABUVFWhvbwcA\ntLW1obKyMvtVjiKVTakuNMFuxoA/zKvIE5H0Eob3xo0b0d/fj1/84heor69HfX09Ghoa0NzcjLq6\nOgSDwdiYeD6kuilVVNHQRYo5XZCIZJdwzHv16tVYvXr1Jfe3tLRkraCxcKe4KVWU0x4J735PAFcU\nFWatLiKibJN8kU5qm1JFTbAVAOBcbyKSn9ThneqmVFFFdg6bEJE+SB3eqW5KFTWBY95EpBPShvdY\nNqWKckbDm8MmRCQ5acM7uilVqjNNAKDIHhnzZs+biGQnbXiPdaYJELmKvNmoMryJSHrShrdn6GLC\nYxnzVhQFTpuZwyZEJD15wzu6l7cl9fAGIkMn/d4gNMELERORvCQO77H3vIHIjBNeiJiIZCdxeI9t\nU6qo6CrL8x6OexORvHQQ3mMcNonN9ea4NxHJS97w9qV+FZ0LTYhOF2TPm4gkJm14pzNVELhgoQ6n\nCxKRxKQNb88YN6WKKuKYNxHpgNThPZZNqaKiOwtyzJuIZCZveI9xU6ooh9UEBex5E5HcpAzvdDal\nijIaVNitJo55E5HUUgrvv//976ivrwcAHDlyBEuWLMHSpUvR1NQETdOyWuBI0tmU6kJF9gL0uf1c\nZUlE0koa3r/+9a+xevVq+P2RMeJ169ahoaEBW7duhRACra2tWS/yYu7oviZphnepyw5/MIzus75M\nlkVElDNJw3vy5Mlobm6O3e7o6EBVVRUAoKamBrt3785edaPoOTcIALjMaUnr9VMnOQAAnV3ujNVE\nRJRLCS9ADAC1tbU4fvx47LYQIjbDw2azwe1OHoDFxVYYjWOb0peIp6MbADC97DK4XI6Rn3SwFw57\nfLhHnztrxkS8/Mb/ovv84Oiv/4STte5UsX1yY/uyL2l4X0xVhzvrXq8XTqcz6Wv6+jI7PHHwaB8A\nwGpU0NMz+snD7RmMux19rqNAhaoo2P9xb8LXf1K5XA4p604V2yc3ti+zxxrNmGebVFRUoL29HQDQ\n1taGysrK9CtLU/dZHxQAVxQVpvX6ApMBV15uw9FuN8J5+MCViGi8xhzeK1asQHNzM+rq6hAMBlFb\nW5uNuhLqOuvDZRMsMJvSH4qZOsmBQFDDqV5+aElE8klp2KS0tBTbtm0DAJSVlaGlpSWrRSUy4A/h\nvDeAmWUl43qfsokO7Hz/FDpPuVHqsmeoOiKi3JBukU7X0PS+T5VYx/U+UydFxuo7u/rHXRMRUa5J\nG94TxxnepS47DKrC6YJEJCX5wrs3M+FtMqooddlxtNuDUJgfWhKRXKQL7+6+zIQ3EPnQMhTWcPKM\nd9zvRUSUS9KFd1evD2ajimJnwbjfa+pErrQkIjlJFd5CCHT1+XBFsRXqGPfxHsnUiUMfWp7ih5ZE\nJBepwrvP7UcgqGHiZeMfMgGAT7tsMBpUHGbPm4gkM+bl8fk03pkm/73vxCX3FdnNOH7ag55zA3Cl\nuWKTiCjXpOp5d8fCO3MhO2NKMcKawEt/+QiC+3sTkSSkCu9TsfC2Zew9yyY5cO2UYrx/qBd7P+zJ\n2PsSEWWTVOHdlYWet6IoqK/9DIwGFS+98RF8g6GMvTcRUbZIFd7dZ31wWk2wWtK7gs5oJpZYcefN\nU3DeE8BrbYcy+t5ERNkgTXgHQxrOnB/MyOKckXzlpimYWGLFm++ewP/f+i4OHj+fleMQEWWCNOF9\nus8HIZCxaYIXMxlVfG/B9fi38svwwdFzeLplL5575T3s+PtJ9PNK80T0CSPNVMEd758CMLwbYCZd\nOIVw1vTL8WmXFe/97xl0dPahozNy1R5XkQVTPuXAov93DS6fwCmFRJRfUoR3z7kBtO49jssnWPD5\nmZOyfrwriq2orZqMfm8Ax057cPy0B6f7BtBzbhD/82EPpk50YNZ0F2Zdczk+7bLFrulJRJQrUoT3\na20fI6wJzL+1HCZj7kZ6nDYzrisrwXVlJRjwh3DstAduXxAHOvvQ2eXG620f4zKnBdNKJ2DqJCem\nTnTgUyVWOK0mBjoRZdUnPrw7u/rRvr8bUyY6UHXtp/JWR2GBEdOvKgIAzCwvwckeL46d9uBkrxfv\n7O/GO/u7Y881GlSUOAvgtJphLzTBVmiEzWKCzWKErdAEh9WMCTYzJtjNKLIVoMCc/uXciOhfU1rh\nrWkafvSjH+HDDz+E2WzGk08+iSlTpmS6NgghsO3NgwCARbddnZHNqDKhwGRA2ZVOlF3phBACbl8Q\nZ84P4mz/IDwDQXgHguj3BtDTN4BU1myaTSpcRYWYYIuEvaPQjEKLASajAWajCpNRhUFVYFBVGAwK\nSor74fP4YTCokcdNKsxGA8wmFQWmyGsMBhVGgxL7mQlEfp5CRL5qAlAAqGrkOYoC/m+BSCJphfcb\nb7yBQCCA3/zmN9i3bx+eeeYZ/PKXv8x0bejscuODo+fwb+WX4dqp47tmZbYoigKnzQynzYzyK+M/\nTBVCIBjS4A+G4Q9qCATD8AfDGAyEMegPwecPwTcY+dp7fhAnevK3r7gCwGBQYTJGThJGgwKjIXLi\nMBrU2Pcmoxo7oUTvNxgUGBQFqjp0EoACKICiAKoSOTkYht7PqCowGlWYLnxv49DxVBWqquCMJ4jz\n/T4oiLzfhYQABASEBmhCxLY00IbOkurQSSh27AtOTkD8CSp6Mou8PvJ99D2jJ7kLX6Nc8N4KlEgd\nAsAFr4ueFOOeO/QzEBCAAHq9QfSe9ULTBMKagCYENE3E2hb3e1GUS9t0YbvUi+5TI8+P/C4UqBh+\nbfTnB8Qf58Kf83D7Io1Qon8cw1+SnuQLvAF4BoIJn5NPF26DEffTFkO3hRjq7MT/jcSeZjSgr38Q\nQPzvOPZTUYZ/VgBgL8zOMGpa4b13717ccsstAIAbb7wR//znPzNaVJSrqBC33ngl/uOmzPfqc0FR\nFJhNBphNBjhSeH4oPBT0gTCCIQ2hsEBY0xDWhoJhKLDMJiN8gwGEw5F//GFNIBTWEA5Hvoa0SBhE\nQ0HBxQEERP68oiGF2HPD2vDXQCgEbSD+PiIam1uun4T//I9rM/6+aYW3x+OB3T58xXWDwYBQKASj\nceS3c7lSia4RXgdgeX16Pe4vp3lMIiIZpDV1w263w+sd/i++pmmjBjcREWVeWuH92c9+Fm1tbQCA\nffv2Yfr06RktioiIElNEGptYR2ebfPRRZA/sp59+GldffXU26iMiohGkFd5ERJRf0mxMRUREwxje\nREQSYngTEUlI6vDWNA2NjY2oq6tDfX09jhw5Evf4tm3bMH/+fCxatAhvvfVWnqpMX7L2AcDZs2dR\nW1sLv9+fhwrHJ1n7Nm/ejIULF2LhwoX4+c9/nqcq05esfS+99BLuuusuLFiwAH/84x/zVGX6Uvn7\n1DQNDzzwAF5++eU8VDg+ydr35JNPYv78+aivr0d9fT3cbnduCxQS+/Of/yxWrFghhBDivffeEw89\n9FDssdOnT4s777xT+P1+0d/fH/teJonaJ4QQbW1tYu7cuWLWrFlicHAwHyWOS6L2HT16VMybN0+E\nQiGhaZqoq6sTBw4cyFepaUnUvt7eXnHHHXeIQCAg3G63qKmpEZqm5avUtCT7+xRCiOeee04sXLhQ\nbN26NdfljVuy9i1evFj09vbmozQhhBBS97wTLdN///33MWvWLJjNZjgcDkyePBkffPBBvkpNS7Jt\nCFRVxYsvvoiioqJ8lDduido3ceJEvPDCCzAYDFAUBaFQCAUFBfkqNS2J2ldSUoLt27fDZDLhzJkz\nKCgokG5jsGR/n3/605+gKErsObJJ1D5N03DkyBE0NjZi8eLFePXVV3Nen9ThPdoy/ehjDsfwEnmb\nzQaPx5PzGscjUfsA4POf/zyKi4vzUVpGJGqfyWRCSUkJhBBYv349KioqUFZWlq9S05Ls92c0GtHS\n0oK6ujp87Wtfy0eJ45KofR999BH+8Ic/4Hvf+16+yhu3RO3z+Xy4++678eyzz+KFF17A1q1bc945\nlDq8Ey3Tv/gxr9cbF+Yy0Ps2BMna5/f7sXz5cni9XjQ1NeWjxHFJ5fd39913Y8eOHdizZw/eeeed\nXJc4Lonat337dnR3d+Pee+/F66+/js2bN8dWZcsiUfsKCwtxzz33oLCwEHa7HTfddBPDeywSLdO/\n/vrrsXfvXvj9frjdbhw6dEi6Zfx634YgUfuEEPj2t7+Nz3zmM/jxj38Mg0G+C1Ykat/HH3+MZcuW\nQQgBk8kEs9kMVZXrn2Oi9j322GP47W9/iy1btmDevHm47777UFNTk69S05KofZ2dnViyZAnC4TCC\nwSDeffddXHfddTmtT+pu3Jw5c7Br1y4sXrw4tkz/xRdfxOTJk3H77bejvr4eS5cuhRACjz76qHRj\npsnaJ7tE7dM0DX/7298QCASwY8cOAMD3v/99zJo1K89Vpy7Z72/GjBmoq6uLjQtXVVXlu+Qx+Vf+\n+7z99tsxd+5cLFq0CCaTCXPnzsW0adNyWh+XxxMRSUiu/6cREREAhjcRkZQY3kREEmJ4ExFJiOFN\nRCQhhjcRkYQY3kREEvo/m+jJztro+eQAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x10ea07dd0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"\n",
"import numpy\n",
"import matplotlib.pyplot as plt\n",
"import seaborn as sns\n",
"\n",
"sns.set(color_codes=True)\n",
"values = numpy.random.beta(0.1, 8, 3000)\n",
"sns.distplot(values)"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x10c20fc10>"
]
},
"execution_count": 25,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD5CAYAAAA6JL6mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl01PW9//HnrMlsmewL2YCwbxJAUamCQQq3QEUiBOIF\nW/HqtdrFe6/Veqz1th4ttb09Vy222qv2571VlrbWpdWqoKAihiVBdgiQkH2dJDNJZjIz398fgbAI\nMyFktsz7cQ4HMt/JzIssr3zymc/381UpiqIghBAiYqlDHUAIIcSVkSIXQogIJ0UuhBARTopcCCEi\nnBS5EEJEOClyIYSIcH6L3OPx8KMf/YgVK1awcuVKjhw5ct7xzZs3U1hYSFFRERs2bAhYUCGEEBen\n9XeHLVu2APD666+zY8cOfv3rX/P8888D0NPTw1NPPcWmTZswGAysXLmSgoICkpOTL/l4jY0dPp8v\nIcFIa2vn5fwfgiacs0F455NsAxPO2SC88w21bCkplkse8zsiv/nmm/nZz34GQE1NDXFxcX3HysvL\nycnJwWq1otfrmT59OiUlJZcV7kJareaK3j+QwjkbhHc+yTYw4ZwNwjtfNGXzOyLvfVItDz30EO+/\n/z7PPPNM3+12ux2L5exPCZPJhN1u9/lYCQlGv/8JXz95Qi2cs0F455NsAxPO2SC880VLtn4VOcDa\ntWv5j//4D5YvX84777yD0WjEbDbjcDj67uNwOM4r9ovx9+tESorF7/RLqIRzNgjvfJJtYMI5G4R3\nvqGW7YqmVt544w1+97vfAWAwGFCpVKjVve+Wl5dHRUUFNpsNl8vFzp07yc/Pv6xwQgghrozfEfnX\nv/51fvSjH3H77bfjdrt55JFHeP/99+ns7KSoqIiHH36YNWvWoCgKhYWFpKWlBSO3EEKI0/wWudFo\n5L//+78vebygoICCgoJBDSWEEKL/5IQgIYSIcFLkQggR4aTIhRAiwkmRCyFEhOv3OnIB724/SYe9\nu1/3nTM1M7BhhBDiNBmRCyFEhJMiF0KICCdFLoQQEU6KXAghIpwUuRBCRDgpciGEiHBS5EIIEeGk\nyIUQIsJJkQshRISTIhdCiAgnRS6EEBFOilwIISKcFLkQQkQ4KXIhhIhwUuRCCBHhpMiFECLCSZEL\nIUSEkyIXQogIJ0UuhBARTopcCCEinBS5EEJEOClyIYSIcFLkQggR4bS+Dvb09PDII49QXV2Ny+Xi\n3nvvZe7cuX3HX3nlFTZu3EhiYiIA//mf/8nIkSMDm1gIIcR5fBb5m2++SXx8PE8//TQ2m40lS5ac\nV+T79u1j7dq1TJo0KeBBhRBCXJzPIl+wYAHz588HQFEUNBrNecf379/PCy+8QGNjI3PmzOGee+4J\nXFIhhBAXpVIURfF3J7vdzr333svy5ctZvHhx3+3PPfccxcXFmM1m7r//flauXMlNN93k87Hcbg9a\nrcbnfcLVu9tP9vu+C64bHqgYQghxHp8jcoDa2lruu+8+iouLzytxRVG44447sFgsAMyePZsDBw74\nLfLW1k6fx1NSLDQ2dvQne0h02Lv7db9Q/B/C+WMn2QYmnLNBeOcbatlSUiyXPOZz1UpTUxN33nkn\nDz74ILfddtt5x+x2O4sWLcLhcKAoCjt27JC5ciGECAGfI/Lf/va3tLe3s27dOtatWwfAsmXL6Orq\noqioiAceeIDVq1ej1+u57rrrmD17dlBCCyGEOMtnkT/66KM8+uijlzy+ZMkSlixZMuihhBBC9J+c\nECSEEBFOilwIISKcFLkQQkQ4v8sPRWB9VFrd7/vOmZoZwCRCiEglI3IhhIhwUuRCCBHhpMiFECLC\nSZELIUSEkyIPoc7uHvqxZ5kQQvgkq1aCrN3hYsfBerbvq+NkXQfGGC0ZyUYyk01kp5nRqOVnqxDi\n8kiRB1HJoQZeeHM/Hq+CWqViVJaVU/V2yqvbKa9uJ8ESw+ypw4gz6UMdVQgRQaTIg6SiroP/efsA\nOq2a5TeM5JoJaVhNerbsqaKl3cnhShvHqtt4+7OTXDcpnREZcaGOLISIEPJ7fBC0O1w8++e99Li9\n3P3Nicy7Ohvr6VG3SqUiyRrL9ZPTuWFKBgDbymrZc7QplJGFEBFEijzA3B4vv/nLl7S0O1k6eyRT\nRyVf8r4jhsWx6PrhWIw6vixv5uDJ1iAmFUJEKinyAHvvi0qOVrVxzfhUvnFtrt/7x5n03DwjC0OM\nhpJDDRyvaQ9CSiFEJJMiDyC3x8vm3dXE6jWsnj8OlUrVr/ezGHvLXKdV8+mXtdQ0OQKcVAgRyeTF\nzgD5qLSak3UdtHY4GZcTzxeH6i/r/RMssRRMy+T9kiq2ldWyeNZwvxtsWcyxdNi7ZXMtIaKMjMgD\n6FBF7xz3uNyEAb1/WqKRGeNScPZ4+GRvLV45eUgIcRFS5AHS3N5NQ2sXw5JNV7QufGxOPFmpZupa\nOtl/vGUQEwohhgop8gA5OxqPv6LHUalUXD8pHWOMltJjTTTaugYjnhBiCJEiD4Bul5sTtR1YjDoy\nk01X/Hixeg1fm5KBosAne2vxeLyDkFIIMVRIkQfAsep2vF6FcTkJ/V6p4k96kpHxuQl0dPawV6ZY\nhBDnkCIPgKoGO9B7gs9gmjo6GWOsln3Hm2ntcA7qYwshIpcU+SBzuT002rpItsYSq9cM6mPrtGpm\nTkhDUeDz/XWyBa4QApAiH3R1zZ0oCgwbhLnxi8lONZObbqHR1s2RU7aAPIcQIrJIkQ+yM2dhDks2\nBuw5rh6Xik6rZveRJrpd7oA9jxAiMkiRD7Kapk50WjXJVkPAnsMYq+WqUUn0uL2Uyi6JQkQ9KfJB\n1O5wYe/qISPJiFo9OKtVLmVcTgJWk56jp9poae8O6HMJIcKbzyLv6enhwQcfpLi4mNtuu40PP/zw\nvOObN2+msLCQoqIiNmzYENCgkaBvWiUpMPPj51KrVcwYl4oClBxskBc+hYhiPjfNevPNN4mPj+fp\np5/GZrOxZMkS5s6dC/SW/FNPPcWmTZswGAysXLmSgoICkpMvvd/2UHd2fjzwRQ6QmWIiK8VEVaOD\nino7ky2Bm84RQoQvnyPyBQsW8P3vfx8ARVHQaM4upysvLycnJwer1Yper2f69OmUlJQENm0Y83gV\n6lo6iTPqMBt1QXveGeNSUatg16EG3HLGpxBRyeeI3GTqHVna7Xa+973v8YMf/KDvmN1ux2KxnHdf\nu93u9wkTEoxotb7XV6ekWHweD5ljzVjMsRc9VN1ox+1RyB1mveR9AsFijmXK6BRKjzSy92gT08al\nhu3HL1xzgWS7EuGcL1qy+d2PvLa2lvvuu4/i4mIWL17cd7vZbMbhOHvBA4fDcV6xX0pra6fP4ykp\nFhobO/w+Tqh02C/+wmL5qd5NslLiYi55n0AZl23lwIlmdh2qJzvVGJYfv3D+vEq2gQvnfEMtm6/i\n9zm10tTUxJ133smDDz7Ibbfddt6xvLw8KioqsNlsuFwudu7cSX5+/mUFG0oaTu9KmJoY/HlqvU7D\nVXnJuNxeyo41B/35hRCh5XNE/tvf/pb29nbWrVvHunXrAFi2bBldXV0UFRXx8MMPs2bNGhRFobCw\nkLS0tKCEDjeKotDS5sRq0qP3M20UKGNy4jlSZePIKRu1zQ4ygrByRggRHnwW+aOPPsqjjz56yeMF\nBQUUFBQMeqhI0+7oocfjJckavLnxC2nUKq6fPIy/bz/Jxi3lfO+2KSHLIoQILjkhaBA0nz4hJyku\ndEUOvbstpiYYKD3W1HdhCyHE0CdFPgia204XeQhH5NB7NaEZY1MAWL/lmFzjU4goIUU+CJraulGp\nIDEuJtRRSI43MHNCGhV1Hew4UB/qOEKIIJAiv0Jer0JrRzfx5hi0mvD4cBbeOBKtRsWfPy7H1eMJ\ndRwhRICFR/NEsDaHC7dHCfn8+LmS4w3cPCOb5nYnH+yqCnUcIUSASZFfobPz46GfVjnXoutyMcVq\neWf7Sdo7XaGOI4QIICnyK9S3YiXEL3ReyBir45avjaDL6eGv206EOo4QIoCkyK9Qc1s3ahUkWMJr\nRA4wJz+T9EQjH5VWU93ofx8cIURkkiK/Ah6vQkuHkwRLDBp1+H0otRo1ywtGoSi9yxGFEENT+LVP\nBLHZnXi9SthNq5zrqrwkxucmsO94C18el31YhBiKpMivQN8LnWG0YuVCKpWKFXNHo1LB+s3H8Hhl\nz3Ihhhop8isQLmd0+pOdauaGKcOoaXLw0Z6aUMcRQgwyKfIr0NrhRK2CeHP4vdB5oaU3jsQQo+GN\nbcexd/WEOo4QYhBJkQ+QoijY7E7iTHrUalWo4/gVZ9LzzVkjcHS7eWPb8VDHEUIMIinyAXJ0u3F7\nFKwRMBo/Y+70LNITjWzZU01VgyxHFGKokCIfoDa7E4B4sz7ESfpPq1GzYu5oFAVe+/AoiuyOKMSQ\nIEU+QDZ772nvkTA/fq4peUlMyUviYEUruw43hjqOEGIQ+L34srg42+kRuTWCRuRnrJg7mv0nWli/\n+Sg2h7NfuzbOmZoZhGRCiIGQIh8gm92FWgVxxvAr8o9Kq/3eZ1xuAvtPtLDveAtTRycHIZUQIlBk\namUAFEWhLYJWrFzMlLwkDDFa9p1ooUN2RxQiokmRD8CZFSuRNj9+Lp1WzYyxKXi9CjsPyVy5EJFM\ninwAbBG4YuVihmdYSEswcKrBTnWjI9RxhBADJEU+AGdWrETSGvKLUalUXDMhFRVQcrAej1eWIwoR\niaTIB6CtY2iMyAESLLGMyYmnvbOHgydbQh1HCDEAUuQD0LtiRYUlDFesDMTU0cnE6DTsLW+ms1v2\nYREi0kiRXyZFUWhzOIkz6SJ2xcqFYnQa8sck4/YocpKQEBFIivwyOboif8XKxYzKspIUF8uJ2g7q\nWzpDHUcIcRmkyC/TUFmxciH16Rc+AUoONeCVfViEiBj9KvKysjJWrVr1ldtfeeUVFi5cyKpVq1i1\nahXHjw/97VHPnpo/tEbkACnxBkYOi6Ol3cmxqrZQxxFC9JPfU/RffPFF3nzzTQwGw1eO7du3j7Vr\n1zJp0qSAhAtHkbpZVn9NG5NCZX0HpUebGJ5uQa/ThDqSEMIPvyPynJwcnn322Yse279/Py+88AIr\nV67kd7/73aCHC0dtfStWdKGOEhDGWC2TRybR7fJQdkwu1ixEJPA7Ip8/fz5VVVUXPbZw4UKKi4sx\nm83cf//9bNmyhZtuusnn4yUkGNFqfY/yUlIs/mKFhHK0iY5OF1aLHmvcV39DCQcW85VfP/SaSRmU\n17RzuLKV/LGpJMTFDsrnJFw/ryDZrkQ454uWbAPe/VBRFO644w4slt4ws2fP5sCBA36LvLXV94qI\nlBQLjY0dA40VUJ1ONy63l3SDjg57d6jjfIXFHDtouaaNSeajPTVs3VNFwfSsK/6chPPnVbINXDjn\nG2rZfBX/gFet2O12Fi1ahMPhQFEUduzYMeTnym2nz+gMx61rB1t2qpm0RANVjQ5qmmQfFiHC2WUX\n+VtvvcX69euxWCw88MADrF69muLiYkaNGsXs2bMDkTFs9BW5aegXuUqlYsbY3uWIuw434pV9WIQI\nW/2aWsnKymLDhg0ALF68uO/2JUuWsGTJksAkC0OtUVTkAEnWWPIy4yivbueTL2u58aphoY4khLgI\nOSHoMpy54HK0FDlA/uhktBoVf9l6nC6nO9RxhBAXIUV+GVo7uonRaYjVR8/aamOsjokjEmlzuPj7\njspQxxFCXIQUeT+5PV7aHS7iTENz/bgvE4YnYjXr+UdJZd+ZrUKI8CFF3k+Nti4UJbqmVc7QadXc\n8rURuHq8/PWTE6GOI4S4gBR5P9Wd3hEwGosc4IYpGWQkGdlaViPLEYUIM1Lk/XSmyK1RWuQatZrb\nZuehKPCnj8tDHUcIcQ4p8n6qa47uETn0XklodJaVPUebOHLKFuo4QojTpMj7qa6lExUM2c2y+kOl\nUrH8plEAbPqoHEX2LBciLEiR91N9SycWkx6NOro/ZHmZVvJHJ3Osuo2yctkdUYhwEN2t1E+d3T20\nd/YQbxmae5BfrqU3jkQF/PnjcrmSkBBhQIq8H2pPv9CZMEQvJnG5MlPMXDcpnapGBzsO1Ic6jhBR\nT4q8H8680Ckj8rOWfG0EGrWKN7Ydx+3xhjqOEFFNirwf6lulyC+UHG9gTn4mjbZutpbVhDqOEFFN\nirwfzo7Ir/zqO0PJouuHE6PT8NanJ3H2eEIdR4ioJUXeD3UtncToNJhiB3xBpSHJatJz84ws2hwu\nNu+++OUAhRCBJ0Xuh1dRqG/tIj3RiEqlCnWcsLNgZg6GGC1/214h29wKESJS5H60tHfT4/aSlhie\nF1sONVOsjgUzc3B0u/lHyalQxxEiKkmR+3Fmj5X0RGOIk4SveTOysBh1vPdFJfaunlDHESLqyKSv\nH/UtXYAU+Uel1T6Pj82JZ+ehRv72eUXfafxCiOCQEbkfZ1aspCdFd5H7MzY7HmOsls27quTiE0IE\nmRS5H3Wn15CnJUiR+6LRqJmSl4TL7eXtz06GOo4QUUWK3I+65k6sZj2GGJmF8mdUppXUeAMfl9bQ\nZOsKdRwhooYUuQ+uHg8t7d1kRPn8eH+p1SpuuWEEHq/CXz+VS8IJESxS5D40tHahAGlS5P02c3wa\nmSkmPttXJ5eEEyJIpMh9kKWHl0+tVnHrDSNRFPjLtuOhjiNEVJAi9+FMkcuI/PLkj05mREYcuw43\ncqK2PdRxhBjypMh9OFPkMkd+eVQqFbfNyQPkknBCBIMUuQ/1LZ1o1CqSrLLr4eUan5vApBGJHKxo\nZf/JllDHEWJI61eRl5WVsWrVqq/cvnnzZgoLCykqKmLDhg2DHi6UFEWhrqWTlHgDWo38vBuIvlH5\nlnK8XhmVCxEofhdHv/jii7z55psYDOdvGtXT08NTTz3Fpk2bMBgMrFy5koKCApKTkwMWNpjsXT04\nut2MzooPdZSIlZNm4doJaXx+oJ5tpdVMyLaGOpIQQ5LfoWZOTg7PPvvsV24vLy8nJycHq9WKXq9n\n+vTplJSUBCRkKMiKlcGx5MaRaNQqXv37QXrcckk4IQLB74h8/vz5VFV99aIBdrsdi8XS97bJZMJu\nt/t9woQEI1qtxud9UlIsPo8HQ9mJ3nndUbkJZ/Mca8ZiDu/58lDn23Ws+Su3TRqZRNmxJl557zDT\nxqb23b7guuFBTOZbOHzNXUo4Z4Pwzhct2QZ83rnZbMbhOHvCh8PhOK/YL6X19N4ll5KSYqGxsWOg\nsQbNkYreIjfp1Ofl6bB3hyqSXxZzbFjmG5cbz6HKVnYerCcr2di33UE4fJ4hfL7mLiacs0F45xtq\n2XwV/4BfxcvLy6OiogKbzYbL5WLnzp3k5+cP9OHCTt+uhzK1csVidBqumZBGj9tL2bGmUMcRYsi5\n7BH5W2+9RWdnJ0VFRTz88MOsWbMGRVEoLCwkLS0tEBlDor61C0OMhjiTPtRRhoSJI5PZe7SJo6fa\nGJuTQIIlJtSRhBgy+lXkWVlZfcsLFy9e3Hd7QUEBBQUFgUkWQh6vl/qWTnLSzHKdzkGiUauYPi6F\nzbuq2XmogZtnZIU6khBDhiyQvoiG1i48XoVhyaZQRxlSMpNNDEs2Udvcycm68Jy7FCISSZFfRE1T\n7/y4FPngUqlUzJyQikatYuehBjq75fqeQgwGKfKLqGnuXY0zLEmKfLBZjHom5yXR5fTwp62yO6IQ\ng0GK/CJqT++jLSPywJg4IgGrSc9Hu6s5XiO7IwpxpaTIL6KmyYFeq5bNsgJEo1Yzc2IaCvCHdw/h\n9sgZn0JcCSnyC3i9CrUtnaQnGVHLipWASU80csOUDE412PnrJ3JZOCGuhBT5BZrauuhxe2VaJQhW\nzB1NsjWWv31ewdEqW6jjCBGxpMgv0LdiRV7oDDhDjJa7Fk0ABV586wBdTneoIwkRkaTIL9C3YkVG\n5EExJjuef7o2l6a2bl778Gio4wgRkaTIL1AjK1aCbskNI8hJM/PJ3lq27a0JdRwhIo4U+QVqmx1o\nNSpS4mXFSrBoNWruXTIJY4yWV987THlNW6gjCRFRBryN7VCkKAo1TZ2kJxrRqOVnXKB9VFp93tvX\nTUpn864q/mt9GYuuz+3b7nbO1MxQxBMiYkhbnaOl3YmzxyPTKiGSmWIif0wyXU43H+2pxiPry4Xo\nFynyc8ip+aE3cUQiwzMsNNq6+bi0Ri7aLEQ/SJGfQ17oDD2VSsWsyelkJBmpanTw6Ze1eBUpcyF8\nkSI/x5kiz5AiDymNWs2c/ExS4mM5UdvB/753GEXKXIhLkiI/R02zA41aRVqCIdRRop5Oq2bu9CwS\nLDF8VFrDS+8clD1ZhLgEKfLTzqxYSU0woNXIhyUc6HUa5l2dxYiMOD7dV8ezf/oSp8sT6lhChB1p\nrNOa27rpcrrJSjGHOoo4R6xeyw9X5jN5ZBJfHm/mF6/txmZ3hjqWEGFFivy0ivreS48NT7eEOIm4\nUIxew3cLJzNrcjonajt4/OUSDla0hjqWEGFDivy0M0WekyZFHo60GjV3fmM8K+eOxtHVwy9f38M7\n20/KihYhkCLvU1FnByAnTaZWwpVKpWLe1dk8dPs04s0x/Onj4zz9xz002LpCHU2IkJIiP62ivoOk\nuBgsRn2oowg/RmVa+cm3ryZ/dDKHT9n4yf98wYe7qmR0LqKWFDlgsztpd7hkWiWCxBn13L90Mncv\nnoBWo+L/3j/C2v/b3XcugBDRRIocqKjrnR/PlRc6I4pKpeLaien87K6ZTBuTwtGqNn7y0he8se04\nPW5Zcy6ihxQ5Z1/ozJUReUSKN8dw/9LJfHfpZOJMet789CSPv/yFXD5ORA3ZxhYZkYe7C7e79WXB\nzBz2HGnkUKWNp/53N2Nz4pk2JgWd9vwxi2yNK4YSGZEDlfUdWE164s0xoY4irpBOq+aaCWksmJmD\n1aTncKWNtz49SUOrrGwRQ5ffEbnX6+Xxxx/n8OHD6PV6nnjiCXJzc/uOP/HEE+zevRuTqXejqXXr\n1mGxRM7ItqPTRXO7kyl5SaGOIgZRaoKBRbNyKTvWzL7jLby3o5LJeUlMyUtCrVaFOp4Qg8pvkX/w\nwQe4XC7Wr19PaWkpP//5z3n++ef7ju/fv5/f//73JCYmBjRooFTWn1k/Hjk/fET/aNRqpo1JITPZ\nxCd7a9lb3kxts4PZMq0ihhi/Uyu7du3ihhtuAGDq1Kns27ev75jX66WiooLHHnuMFStWsGnTpsAl\nDZCzL3TKiUBDVVqikcWzhpOb3nvBirc/O8nhSjnFXwwdfkfkdrsds/lsyWk0GtxuN1qtls7OTv75\nn/+Zb3/723g8HlavXs2kSZMYN27cJR8vIcGIVqvx+ZwpKcEbHdfbugHIn5BBSqLR952PNWMxh/dF\nmcM5X6izLZw1gr3Hmvh0bw1Pv17KmsUTWXzDSCC4X3OXK5yzQXjni5ZsfovcbDbjcJw9ycLr9aLV\n9r6bwWBg9erVGAy9+3dfe+21HDp0yGeRt7Z2+ny+lBQLjY0d/Qo/GI5UtGCK1aJyu/v1vB327iCk\nGhiLOTZs84VLtpEZFkwx2Xy+v54X/7qPwydb+EHxdGyt4XkiUbC/Hy5XOOcbatl8Fb/fqZVp06ax\ndetWAEpLSxkzZkzfsZMnT7Jy5Uo8Hg89PT3s3r2biRMnXla4UOrs7qG+tYucNAsqlbwAFi3SEo38\n+I4Z5KSZ2VpWw09e2I69qyfUsYQYML8j8nnz5vHpp5+yYsUKFEXhySef5OWXXyYnJ4e5c+dyyy23\nsHz5cnQ6HbfccgujR48ORu5BcfhU7wkjo7OsIU4igi0xLpYf3T6d3799gF1HGnniDzv5/rIpZMiF\nt0UE8lvkarWan/70p+fdlpeX1/fvu+66i7vuumvwkwXB4creIh+XkxDiJCIUYvQa7r11Ev/YVc2G\nD47wxP/bxb1LJjJphCxFFZElqs/sPFTRilajJi8zLtRRRJCde7ZoSryBr03J4LN9dfx6fRkzxqUy\nLje+b7pNzgIV4S5qz+y0d/VwqsHOqMw4dH5W0Yihb+SwOOZfk02MXkPJoQa2ldXKxlsiYkRtkR+u\ntKEg0yrirJR4AwuvzyUlPpaTdR38bXuFXB9URIQoLvLeE0LG5UqRi7NMsTrmX5PD+NwE2hwu3vms\ngnd3VOL1ykUrRPiK2jnyQ5Wt6LVqRmTI/Lg4n1qt4urxqaQmGNhxoJ4NW45RcqiBb39jHFkplz4D\n+HJ2aZR5dzGYorLI2ztdVDU6mDA84SvbmwpxRm66hbREA5X1dj7fX8/jL5XwtSnpfHPWCBLjwvcM\nWhF9orLIj5xedjhW5seFH7F6LXcvnsjM8Wls2HKMrWW1fLavnjn5w5g3I5uUeEOoIwoRnUV+6PT8\n+HgpctFPV41KZvLIJLbvr+ONbSf4YGcVH+6s4qpRycydnsX44fK1JEInSovcRoxOw/CM8N1QR4Qf\ntVrFrMkZzJyQRsnBBj7YVUXpsSZKjzWRbI0lM8VEXqYVs0EX6qgiykRdkdvsTmqaHEwakYhWI/Pj\n4vJpNWqum5TOdZPSOV7Tzkel1ZQcbKDsWDNlx5rJTDYxJieezGSTXMRCBEXUFfkXBxsA5IpAYlCM\nHBbHyGFxrJw7mv99/whHT9mobnJQ3eTAGKtlXE48Y7Lj0evkpDMROFFX5J/vr0OtUnHN+LRQRxFD\niCFGy+gsK6OzrLS0d3O0qo3y6jZ2H2niy/IWRmdbmTA8EWNs1H3LiSCIqq+q2mYHJ+s6mJKXRJxJ\nH+o4YohKjItl5oRY8kcnc+SUjYMVrRw42crhShtjsuOZNDIyL4sowldUFfnn++sBuHaCjMZF4Ol1\nGiaNTGL88ATKq9v5sryZgxWtHDllw97Zwz9dmysvjIpBETVFrigK2/fXEaPTkD86JdRxRBTRqNWM\nyY4nLzOOY1W9hf73HZV8VFrN16/OYd6MbJlyEVckar56yqvbaWrr5rqJacTo5YUnEXwatZqxOfGM\nyozD44V3tp/kr5+c4B8lp5g7PZN5M7KxGGXKT1y+qCny7QfqALhuYnqIk4hIczl7qPSHRqNm7vRM\nbrwqgw9oALRxAAAMg0lEQVR3VfGPklO8/VkF/yg5xaxJGcyeOoycNDnHQfRfVBS52+Ol5GADcUad\nnIEnwkasXsvC64Zz84xstpbV8N4XlWzZU82WPdWMyLBw3cR0CmYOj94tSkW/RUWRf7avDntXD/Nm\nZKNRy7eFCC8xOg3zZmRTMC2TL8tb+Li0mr3HmzlR28EfPzhKVoqJiSMSyRtmJS/TSoIlJtSRRZgZ\n8kXudHn4y7bj6LVqFszMCXUcIS5Jo1YzdXQyU0cn09rhpPRYEwcrbZQeaaSq8RRwCoB4s57sVAvZ\nqWayU83kpJlJSzDKWaRRbMgX+T9KKmmzu1h0fa6MZETY6M+8u0oFMyemMzozjua2bhptXTTaumhp\nd/Ll8Wa+PN7cd1+9Vk12qpmRw6zkZcYxKtMqW+1GkSFd5G0OF3/bUYnFqOOfZuaGOo4QA6LTqklP\nMpKeZOy7rdvlobWjm9YOJ63tTlo6nByvbae8pp33d/bex2LUkZ5oJCPZxLAkY982AXJRi6FnSBf5\nm5+cwOnysGxOHoaYIf1fFVEmVq8hI8lERpKp7za3x0tLezeNtm7qW7uob+nkaFUbR6vaUKkgNd5A\nVqqZcTkJpCcafTy6iDRDtt1O1LbzcWkNaYlGbrxqWKjjCBFwWo2a1AQjqQlGJo4Ar1ehub2bmiYH\n1Y2O3nJv7WLX4UbSE41MHZXMVaOSyMu0yk6gEW5IFnlts4NfbyhDURSKbx4tX6QiKqnVKlLiDaTE\nG7hqVDJdTjfVjQ66XR72nWjm3S8qefeLSgwxWiaNSGTiiEQm5CaQLFc9ijhDrshb2rv5r/Wl2Lt6\nWL1gLJNHyna1QkDvDo2jsqwAjM+Np66lk6pGB1UNdkoONVByqHeLZ7NBR1qCga9NyWBUppWMJNlX\nPdwNqSJvsnXx641lNLc7KZw9Ul7UEeISNBo1mSlmMlPMXDM+lTaHi7rmTmqbO6lr6aS8pveFUwC9\nTk12ipnsNAuZySbSE42kJxpJiJNVYOFiSBS5q8fDuzsqeefzCnrcXr5+dTbfuFZWqQjRHyqVinhz\nDPHmGMblJuBVFNrsLhIsMZRXtVFRb+dkXUdfsZ+hUatIjjcQb9JjNeuJM+qxmPSYY7UYYrUYY7TE\n6DTodRp0WjU6rRqNWoVG3fu3StU7/dM71j93xK/gVXr/5fUqeBXlgr97N8HzehUUBUoON/Q9glqt\nQq1SoVariLca6O52oVadfeyhOrjzW+Rer5fHH3+cw4cPo9freeKJJ8jNPVuSGzZs4PXXX0er1XLv\nvfdy0003BTTwGYqiUN3ooPRYE1vLamhq68Zq0rN8wSiunZiGSiW/CgoxEGqVqu+ci7wsK3lZVjwe\nLza7izaHi/bTf+xdPdg7XdS3dIY4sW9ajQq9VkOMXkPJwQYsRh1WUwxWsx6rSX/6h5geqzkGU6w2\nIrvDb5F/8MEHuFwu1q9fT2lpKT//+c95/vnnAWhsbOTVV1/lT3/6E06nk+LiYmbNmoVeP/g7uHV2\nu/niYD31rZ00tHZRWd9Bc7sT6B0ZzL8mm2/OGiHLDIUIAI1GTZI1liTr+ScZWcyx2Nq7cLrcdLk8\ndDs9uHo8uNxeetweejwKHo8Xj1fB4zk7qlYUBQVQFFDOfUBF6StSlap3lK1S9Y7ez/27d8q+9+1z\n3vXsqN2rgEpFV3cPPR4vrh4v9s4eDna0+vx/ajUq4kynf7sw6jEbdFiMOowxWgyn/8ToNei1amJ0\nGrRaNTqNGo1GhUbd+5uARqW6aGatRoUpNjD7z/ttvV27dnHDDTcAMHXqVPbt29d3bO/eveTn56PX\n69Hr9eTk5HDo0CGmTJky6EE/3HWKv2w70fe2KVbLNeNTmToqmUkjk2SDfiFCRKNWYYzVYQxQSQ2U\nxRxLh737vNu+NjmDjs4e2hxObPbe3yxsHU5sdidtDtfp25xUNTpwezoGPdOaheOZNTlj0B/Xb5Hb\n7XbMZnPf2xqNBrfbjVarxW63Y7Gc3W7TZDJht9t9Pl5Kiv/tOS92nzuXTOHOJYP/A+JyLOhHdiFE\neBv8Gh2Y/nRhf/ldYG02m3E4HH1ve71etFrtRY85HI7zil0IIUTg+S3yadOmsXXrVgBKS0sZM2ZM\n37EpU6awa9cunE4nHR0dlJeXn3dcCCFE4KkURVF83eHMqpUjR46gKApPPvkkW7duJScnh7lz57Jh\nwwbWr1+Poijcc889zJ8/P1jZhRBC0I8iF0IIEd5kExIhhIhwUuRCCBHhpMiFECLChbzIu7u7+e53\nv0txcTH/8i//QktLy1fus3btWoqKiigsLGTDhg0Bz+T1ennssccoKipi1apVVFRUnHd8w4YNLF26\nlOXLl7Nly5aA57mcbK+88grLli1j2bJlPPfcc2GV7cx97rrrLl577bWgZutPvo8//pjly5ezbNky\nHn/8cYL58pG/bC+99BJLly6lsLCQ999/P2i5zlVWVsaqVau+cvvmzZspLCykqKgoKN+fF3OpbG+/\n/TbLli1jxYoVPPbYY3i93rDJdsaPf/xjfvnLX17Zkygh9tJLLynPPPOMoiiK8vbbbys/+9nPzju+\nfft25Tvf+Y6iKIridDqVm2++WbHZbAHN9N577ykPPfSQoiiKsmfPHuVf//Vf+441NDQoixYtUpxO\np9Le3t7372Dxla2yslK59dZbFbfbrXi9XqWoqEg5ePBgWGQ741e/+pWybNky5Y9//GPQcp3hK19H\nR4eycOFCpbm5WVEURXnhhRf6/h3qbG1tbcrs2bMVp9Op2Gw2Zc6cOUHLdcYLL7ygLFq0SFm2bNl5\nt7tcrr7vSafTqSxdulRpbGwMi2xdXV3K3Llzlc7OTkVRFOWBBx5QPvjgg7DIdsZrr72mLF++XHn6\n6aev6HlCPiI/dwuAG2+8ke3bt593PD8/nyeffLLvbY/H03dCUjAy+dqWwGKx9G1LECy+sqWnp/P7\n3/8ejUaDSqXC7XYTExO8rUZ9ZQN49913UalUffcJNl/59uzZw5gxY1i7di3FxcUkJyeTmJgYFtkM\nBgPDhg2jq6uLrq6ukGzqlJOTw7PPPvuV28vLy8nJycFqtaLX65k+fTolJSVhkU2v1/P6669jMPRe\nKCPY3w++sgHs3r2bsrIyioqKrvh5grrD1MaNG/nDH/5w3m1JSUl9Z4OaTCY6Os7f3yAmJoaYmBh6\nenp4+OGHKSoqwmQyEUiDvS1BsLLpdDoSExNRFIVf/OIXTJgwgREjRoRFtiNHjvD222/zzDPP8Jvf\n/CZomfqbr7W1lR07dvDGG29gNBq5/fbbmTp1atA+fr6yAWRkZLBw4UI8Hg/33HNPUDKda/78+VRV\nVX3l9lB/P8Cls6nVapKTkwF49dVX6ezsZNasWWGRraGhgd/85jc899xz/P3vf7/i5wlqkZ+Zuz3X\n/fff33eav8PhIC4u7ivv19bWxve+9z2uueaaoHwRh/O2BL6yATidTh555BFMJhM/+clPgpbLX7Y3\n3niD+vp67rjjDqqrq9HpdGRmZnLjjTeGRb74+HgmT55MSkoKADNmzODgwYNBK3Jf2bZu3UpDQwMf\nfvghAGvWrGHatGkB2ZzucoX6+8Efr9fL008/zYkTJ3j22WfDZovad999l9bWVu6++24aGxvp7u5m\n5MiRLF26dECPF/KplWnTpvHxxx8DvV+w06dPP+94d3c33/rWtygsLOS+++4LWqZw3ZbAVzZFUfjO\nd77D2LFj+elPf4pGowlaLn/ZfvjDH7Jx40ZeffVVbr31Vr71rW8FtcT95Zs4cSJHjhyhpaUFt9tN\nWVkZo0aNCotsVquV2NhY9Ho9MTExWCwW2tvbL/VQQZWXl0dFRQU2mw2Xy8XOnTvJz88Pdaw+jz32\nGE6nk3Xr1vVNsYSD1atX8+c//5lXX32Vu+++m0WLFg24xCEMrhC0cuVKHnroIVauXIlOp+NXv/oV\nAL/4xS9YsGABu3fv5tSpU2zcuJGNGzcC8OSTT5KdnR2wTPPmzePTTz9lxYoVfdsSvPzyy33bEqxa\ntYri4mIUReGBBx4I6rybr2xer5cvvvgCl8vFtm3bAPi3f/u3oH1j+fu4hZq/fP/+7//OXXfdBcCC\nBQuC+gPaX7bPPvuM5cuXo1armTZtWtCnCC701ltv0dnZSVFREQ8//DBr1qxBURQKCwtJS0sLi2yT\nJk1i06ZNzJgxgzvuuAPoLdB58+aFPNtgzIufS07RF0KICBfyqRUhhBBXRopcCCEinBS5EEJEOCly\nIYSIcFLkQggR4aTIhRAiwkmRCyFEhPv/kFHHAEh34y0AAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x10c07f3d0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"\n",
"import scipy\n",
"import seaborn as sns\n",
"\n",
"sns.set(color_codes=True)\n",
"values = scipy.stats.gompertz.rvs(3, size=1000)\n",
"sns.distplot(values)\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"ename": "TypeError",
"evalue": "_parse_args() got an unexpected keyword argument 'size'",
"output_type": "error",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[0;31mTypeError\u001b[0m Traceback (most recent call last)",
"\u001b[0;32m<ipython-input-4-de628a3af062>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m 5\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 6\u001b[0m \u001b[0msns\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcolor_codes\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 7\u001b[0;31m \u001b[0mvalues\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mscipy\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mstats\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbeta\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpdf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m2\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m8\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msize\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m100\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 8\u001b[0m \u001b[0msns\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdistplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python2.7/site-packages/scipy/stats/_distn_infrastructure.pyc\u001b[0m in \u001b[0;36mpdf\u001b[0;34m(self, x, *args, **kwds)\u001b[0m\n\u001b[1;32m 1643\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1644\u001b[0m \"\"\"\n\u001b[0;32m-> 1645\u001b[0;31m \u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mloc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mscale\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_parse_args\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwds\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 1646\u001b[0m \u001b[0mx\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mloc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mscale\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmap\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0masarray\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mloc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mscale\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1647\u001b[0m \u001b[0margs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmap\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0masarray\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0margs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;31mTypeError\u001b[0m: _parse_args() got an unexpected keyword argument 'size'"
]
}
],
"source": [
"%matplotlib inline\n",
"\n",
"import scipy\n",
"import seaborn as sns\n",
"\n",
"sns.set(color_codes=True)\n",
"values = scipy.stats.beta.pdf(2, 8, size=100)\n",
"sns.distplot(values)"
]
},
{
"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.13"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment