{ "metadata": { "name": "", "signature": "sha256:54b0ca29b543517953a1bf9cff1ed8eae17d232857488587a9dd3b3a3dab14ca" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "from __future__ import print_function, division\n", "import thinkbayes2\n", "import thinkplot\n", "%matplotlib inline" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 40 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Ignore the first few cells for now -- they are experiments I am working on related to the prior." ] }, { "cell_type": "code", "collapsed": false, "input": [ "mu = 1\n", "pmf = thinkbayes2.MakeExponentialPmf(mu, high=1.0)\n", "thinkplot.Pdf(pmf)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Warning: Brewer ran out of colors.\n" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEACAYAAABRQBpkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VOXZx/FvSBAQkLAIaMAGN1QuAaVAhCrDJgRs4lK0\naEtBrbSKW21V9H1rqrVqVQRFkSpYrEW0SMUFEBTiWhRkKSgg8BqEgIqIEAICCfP+cU88M8mQOQmZ\nObP8Ptc118w855zM7Smde54dREREREREREREREREREREREREREREXBsMrAXWA7cd5pxHA8dXAme5\nuLYH8BGwHFgCdK/bkEVEpC6lAxuAbKA+sAI4vdI5Q4A5gdc9gcUuri0EBgVe5wKL6jpwERFxr16E\n4z2wL/Qi4CAwA8ivdE4eMC3w+kMgE2gb4dptQLPA60yguJbxi4hIHciIcDwL2Bz0fgv26z/SOVnA\n8dVcezvwHvAQlpDOqVHUIiJSpyLVDPwu/05aDT93CnADcAJwMzC1hteLiEgdilQzKAbaB71vj/3C\nr+6cdoFz6ldzbQ9gQOD1TODpcB9+0kkn+Tdu3BghRBERqWQjcHJd/sGMwB/NBo4icgdyDk4HcnXX\nLgP6BF73x0YUheMXc9ddd3kdQtzQvXDoXjh0Lxy4b9UJ+bKvThkwBngDGx00BVgDjA4cn4wlgiFY\nZ3EpMCrCtQDXAI8DDYB9gfdhzZkDQ4bU5D9JRERqKlIyAJgbeASbXOn9mBpcC7CUqh3RYf3xj7B1\nK1x1FaTVtGdCRERcidSBHBeefBLuvhsOHvQ6Eu/4fD6vQ4gbuhcO3QuH7sWRifff2v5u3Zymr+7d\n4a9/haZNPYxIRCTOpVkzSo2+3+O+ZpAfNMVtyRJrLtq61bt4RESSUdzXDA4d8vPMM/DEE05hy5bw\nyCNwxhneBSYiEq9qUzOI+2Rgo6Rg3jz405+cfoOGDeHee6FPn2quFhFJQUmdDACWLYPf/x5277b3\naWnwu9/B8OEeRSciEoeSPhkAbNoEN94IW4LmQV92mSWF9PQYRyciEodSIhkA7NxpX/6rVjllvXvD\nX/4CjRvHMDoRkTiUlKOJwmne3OYe9O/vlL3/Plx5JRRrMWwRkRpLyJpBhUOHLClMDVrztHlzeOgh\n6NIlBtGJiMShlGkmqmzOHLjnHmekUf368L//qzWNRCQ1pWwyAFi50kYa7dzplI0aBb/9LdRLyMYw\nEZHaSelkADYz+eabIXgLhH79bH5Co0ZRiE5EJA6lfDIAKC2FO+6wDuUKp50G48ZB69Z1HJ2ISBxS\nMggoL4cJE2D6dKfs2GMtIZxeeWseEZEko2RQyaxZtsppWZm9b9jQ9kc4//w6ik5EJA4pGYSxZAnc\neiuUlDhl6lgWkWSmZHAYmzbBLbdAUZFTdu65Nhy1SZMj/vMiInElZWYg19SPfgTPPGNLVlR4910Y\nORK++MKzsERE4kZK1AwqlJfbvgjTpjllTZvamkbnnFNnHyMi4ik1E7k0b541Ee3fb+/r1YMbboAr\nrrBlsUVEEpmSQQ2sWWMzlr/6yikbOtTmKDRoEJWPFBGJCSWDGtqxw0YarVzplHXqBA8+qAlqIpK4\n1IFcQy1bwqRJkJ/vlH3yCYwYEZogRESSXUrXDJwPgX/9Cx5+2DqZATIybAOdYcPUjyAiiUXNREdo\n6VK4/Xb47junTP0IIpJolAzqwLZt1o+wZo1Tdtpp8MADkJUV01BERGpFfQZ14Ljj4OmnIS/PKVu7\n1voRFi/2Li4RkWhykwwGA2uB9cBthznn0cDxlcBZLq6dASwPPD4PPMeNBg1sp7SxY63vAGDXLrj+\nepvJHOPKiohI1EWqRqQD64ABQDGwBBgOBDWiMAQYE3juCUwAclxeC/AQ8B3w5zCfH/NmospWrbJm\no+3bnbK+faGgABo39iwsEZHDikYzUQ9gA1AEHMR+0edXOicPqFjg4UMgE2jr8to04FLg+ZoEHUtn\nngnPPQdnn+2ULVoEv/oVfP65d3GJiNSlSMkgC9gc9H5LoMzNOce7uPZc4CtgI3GsZUtb0+jyy52y\noiJLCAsXehaWiEidiZQM3LbR1HZU0nBgesSz4kDFvIM//9k2yQHYu9eakB57zJmfICKSiDIiHC8G\n2ge9b4/9wq/unHaBc+pHuDYDuAgIaoCpqqCg4IfXPp8Pn88XIeToGjwYTjoJ/vAH2BL4r5k2DVav\nhnvvhVatPA1PRFJQYWEhhYWFR/Q3Iv2iz8A6gfsDW4GPqL4DOQcYH3iOdO1gbIRR32o+3/MO5MPZ\nvdtGHL3/vlPWsqUlhB//2Lu4RESi0YFchn3RvwF8CryAfZmPDjwA5gD/h3UWTwaujXBthcuI447j\nSI45Bh55BH7zG2e5ih074NprYepUOHTI2/hERGpCM5DrwEcfwf/8D3z7rVPWqxfcfTdkZnoXl4ik\nJi1H4aHt220No+VB0+fatoX77rPhqSIisaJk4LFw22pmZNguasOHa/VTEYkNJYM48c47cNddUFLi\nlPXvbx3OTZp4F5eIpAYlgzhSXGxrG336qVPWvj3cfz907OhdXCKS/JQM4syBAzB+PLz4olN21FE2\nR+HCC9VsJCLRoWQQp+bPt5nLe/c6Zeefbx3OajYSkbqmZBDHNm2C226DDRucsqws+MtfoFMn7+IS\nkeSjZBDnvv8exo2DWbOcsowMuO46uOIKqKethkSkDigZJIj5861GsGePU9arl+2R0KKFZ2GJSJJQ\nMkggxcVw5522wF2FVq3gnnuge3fv4hKRxKdkkGAOHoRJk+DZZ52ytDS48kq45hpIT/cuNhFJXEoG\nCeqDD2yS2s6dTlnXrjYCqW1b7+ISkcQUjVVLJQZ69YLnnw9tHlqxwnZWO8IlykVEXFHNII6Ul9u6\nRpMnh+6cNmwY3HQTNGjgXWwikjjUTJQkVq60zuUvv3TKOnSwZiMtZSEikSgZJJFdu+zLf9Eipywj\nA8aMseYjzUkQkcNRMkgyfj/Mng0PPWQT1ir06GFzElq39iw0EYljSgZJatMmW/46eAXUZs2sKalf\nP+/iEpH4pGSQxMrK4G9/g2eesRpDhfx8uOUWOPpo72ITkfiiZJACli2DP/4xtHO5fXvrX9CCdyIC\nSgYpo6TENsl54w2nLD0dRo+GX/1KM5dFUp2SQQrx+2HuXHjgASgtdcq7doU//cmWxxaR1KRkkIKK\ni20pixUrnLKjj4abb9ZuaiKpSskgRZWXW8fyU0+Fzlzu3dtGIbVq5V1sIhJ7SgYp7pNPrJZQVOSU\nNWsGt98OAwd6FpaIxJiSgbB/P0ycaAvfBRs0CG691ZKDiCQ3JQP5wdKlNks5eAhqq1Y2LLVXL8/C\nEpEYUDKQEKWl8PDD8MoroeUXX2yroGqimkhyUjKQsN55xyalffutU5aVZTWHs87yLCwRiZJobW4z\nGFgLrAduO8w5jwaOrwSCv16qu/Z6YA2wGnjAfchSU+edBy+8AP37O2XFxba15oQJ1s8gIqktUuZI\nB9YBA4BiYAkwHPsSrzAEGBN47glMAHIiXNsXuCNwzUHgWGB7mM9XzaAO+f0wb55NVNuzxynPzrYh\nqF26eBaaiNShaNQMegAbgCLsS3sGkF/pnDxgWuD1h0Am0DbCtb8F7guUQ/hEIHUsLQ1yc62W0LOn\nU15UBFdfDY88ErpUtoikjkjJIAvYHPR+S6DMzTnHV3PtKcB5wGKgEPhxTYKWI9OmjQ0/veMOaNzY\nyvx++Oc/YfhwWL7c2/hEJPYyIhx320ZT047oDKA51pzUHXgRODHciQUFBT+89vl8+Hy+Gn6UhJOW\nZqOKevWyzuXFi61882brS7j0UttVrVEjb+MUkcgKCwspLCw8or8R6Us8ByjAOoIBxgKHCO3wfRL7\ndT8j8H4t0AfoUM21c4H7gbcDxzZg/Q07Kn2++gxiwO+HV1+FceNC+xLatbO+hG7dvItNRGouGn0G\nS7EmnWzgKOAyoNKodV4BRgRe5wDfAV9FuPZloGKPrlMDxysnAomRtDTIy7O+hN69nfItW2xZ7Pvv\nD10ZVUSSj5vMkQuMx0YHTcE6fkcHjk0OPE/EagClwChgWTXXAtQHpgJdgQPALVjtojLVDGLM74c5\nc2zf5ZISp/y446yW0KOHd7GJiDuadCZ1Zvt2qxG8/XZo+YUXwo03QtOm3sQlIpEpGUid8vttN7UH\nH4Rdu5zyVq1s0bu+fbVfgkg8UjKQqNixwyaqLVwYWt6nD9x2G7Ru7U1cIhKekoFE1cKF8Ne/wjff\nOGVNmtgQ1IsvhnpuFjcRkahTMpCoKymBxx6DWbNCy7t2hTvvhA4dvIlLRBxKBhIzy5bBvffCpk1O\nWf36MGoUjBwJRx3lWWgiKU/JQGJq/36YOhWmTYOyMqf8xBOtlqCF70S8oWQgntiwwZa0WL3aKUtL\ng2HD4LrrnPWPRCQ2lAzEM+Xl8K9/weOPw759Tvmxx8Lvfw/9+mkYqkisKBmI57Zts8lq778fWt67\nt81NyKq85q2I1DklA4kLfj/Mn28L3+0IWnGqQQO46ir4xS/UwSwSTUoGEldKSuCJJ2DmTEsQFbKz\nYexYrYYqEi1KBhKXPvkE7rsP1q4NLR861NY5atHCm7hEkpWSgcStig7mSZNCl8M+5hi4/nrIz9cM\nZpG6omQgce/rr60v4c03Q8s7d7amo1NO8SYukWSiZCAJ44MPbJ2jLVucsvR0uOwy23azSRPvYhNJ\ndEoGklAON4O5ZUu44QbIzVXTkUhtKBlIQioqsrkJS5eGlnfpYnMTOnb0JCyRhKVkIAnL77d+hEce\nsX6FCvXqwSWXwG9+A82aeRefSCJRMpCEt3evNR0991xo01Fmpu2bkJenpiORSJQMJGls2mTbbS5e\nHFreqZM1HXXq5E1cIolAyUCSit8PhYU2FHXbNqc8Lc3mJVx3HTRv7ll4InFLyUCS0vff24ijadPg\nwAGnvGlTG4Y6bBhkZHgXn0i8UTKQpLZli9US3nkntDw7G373O+jVy5OwROKOkoGkhPfeg4cfhs2b\nQ8t794abb7bkIJLKlAwkZRw4ADNmwJQpoWsdZWRYs9Gvf23rHomkIiUDSTk7dtjid7Nnhy6TnZlp\ncxMuusiWuRBJJUoGkrLWrbOmo2XLQstPPtn6E3r08CYuES8oGUhK8/th4UKYMAG2bg091qcP3HQT\ntG/vTWwisVSbZOBmLudgYC2wHrjtMOc8Gji+EjjLxbUFwBZgeeAxuCZBi4STlgb9+9u+CddeC40a\nOcfeftv6EsaNg127vItRJF5FyhzpwDpgAFAMLAGGA2uCzhkCjAk89wQmADkRrr0LKAHGRfh81Qyk\n1rZvh8cfh9deCy1v2hSuvtqSg/ZilmQUjZpBD2ADUAQcBGYA+ZXOyQOmBV5/CGQCbV1cG+9NVJLg\njj0WCgpsslrXrk55SYktiDdsGCxYENrxLJKqIiWDLCB4NPeWQJmbc46PcO31WLPSFCyBiERFp07w\n1FO2mU5wn0Fxse2uNmoUrFzpXXwi8SDSJH63v5lq+it/EnB34PU9wMPAVeFOLCgo+OG1z+fD5/PV\n8KNErD+hXz8491x46SVLDhV9B6tXw1VXWX/DmDHqZJbEU1hYSGFh4RH9jUhf4jlYZ29FB+9Y4BDw\nQNA5TwKFWDMQWIdxH6CDi2sBsoFXgTPDfL76DCQqSkrgmWfg+efh4EGnvGLS2tVXa/8ESVzR6DNY\nCpyCfWEfBVwGvFLpnFeAEYHXOcB3wFcRrj0u6PqLgFU1CVrkSDVtaltrvvQSDBrklJeVWYLIz4dn\nn7WtOUVSgZvMkQuMx0YHTQHuA0YHjk0OPE/EagClwChgWTXXAjwLdMWaoT4P/L2vwny2agYSE598\nAuPHw/LloeWtW8Po0XDBBZrJLIlDk85EjoDfb/MRHn0Uvvgi9Fh2tu2f4PNZ/4NIPFMyEKkDZWW2\n1tFTT8E334QeO/NM62Tu1s2b2ETcUDIQqUP79ln/wbPPwp49ocd697aawqmnehObSHWUDESi4Lvv\nbOTRiy+GjjxKS4PBg2111KzKs29EPKRkIBJF27bB5Mnw+uuhs5YzMuCSS+DKK6FlS+/iE6mgZCAS\nAxs32ppHlbffbNgQfv5z+OUvNUdBvKVkIBJDK1bAxIn2HKxJE7j8crjiCmjc2JvYJLUpGYjEmN8P\n775ru62tXx96rFkzGDECLr00dDltkWhTMhDxyKFD8NZb1qdQVBR6rEULWwzv4ouhQQNPwpMUo2Qg\n4rHycpg3D/72N1sVNVjr1rbm0U9/CvXrexOfpAYlA5E4UVYGr7wCU6bAV5UWWsnKgmuusWGpWuJC\nokHJQCTO7N8P//63zVPYsSP0WHa2JYX+/ZUUpG4pGYjEqX37bG/madOq7sHcoYPtpzBwoJKC1A0l\nA5E4V1oK06fDP/9ZdYmL7GybuDZokJKCHBklA5EEsWuXJYUZMyxBBDvhBKspqE9BakvJQCTB7Npl\ni+HNmFG1ptCundUUhgyxJS9E3FIyEElQJSWWFJ5/3l4Ha9fO5ikMHaqkIO4oGYgkuD17rJYwfTrs\n3h167PjjLSlccIHmKUj1lAxEkkRpKbzwgnU0Vx591LatLXORn68ZzRKekoFIkikttSGpzz1n+yoE\na9HCFsP72c+0IJ6EUjIQSVJ798LMmfCPf8DOnaHHmja1xfCGD4fMTG/ik/iiZCCS5L7/Hl5+2ZJC\n5WUuGja0xfCuuALatPEmPokPSgYiKeLgQZg7F/7+d/jii9BjGRnWyTxihM1ZkNSjZCCSYsrLYeFC\nW/vos89Cj9WrBwMGwMiRcOqpnoQnHlEyEElRfj/85z8wdWrVndcAfvITG5bapUvsY5PYUzIQEZYv\nt5rCBx9UPda5s+3RfN55WuoimSkZiMgP1q61pLBwodUcgp1wAvziFzarWXMVko+SgYhUUVRko4/m\nzLGO52DNm8Nll8GwYbZnsyQHJQMROaxvvrGlLmbOrLooXsOGkJdnw1KzsryJT+pOtJLBYGA8kA48\nDTwQ5pxHgVxgLzASWO7y2luAB4FWwLdh/q6SgUgdKy2F2bNt/aMvvww9Vq8e9Otn/QqdOnkTnxy5\naCSDdGAdMAAoBpYAw4E1QecMAcYEnnsCE4AcF9e2B54COgLdUDIQiamyMliwwJqQKg9LBejWzZJC\nr16WJCRx1CYZRPqfuAewASgCDgIzgPxK5+QB0wKvPwQygbYurh0H3FqTYEWk7mRkQG6uLYb3+OPQ\ns2fo8Y8/hptusqUuZs60rTsleUVKBlnA5qD3WwJlbs45vppr8wPv/1vDeEWkjqWlWSJ4/HFLDLm5\nocNOi4rg/vtt5NFjj1VdBkOSQ6StMty20dSkOtIIuAMY6Ob6goKCH177fD58Pl8NPkpEaqJjR7jn\nHrjuOttoZ/Zsp7N5926YNs1WUB0wAC6/XP0K8aKwsJDCwsIj+huRvsRzgAKsIxhgLHCI0I7gJ4FC\nrBkIYC3QB+hwmGtfB97COpsB2mF9Cj2Aryt9vvoMRDxUWgqvvmqjkLZsqXq8c2dLCn37ahJbPIlG\nB3IG1gncH9gKfET1Hcg52OihHJfXAnyOOpBF4lp5Obz3no1A+vjjqsfbtrW+hYsusiW1xVvRGlqa\nizM8dApwHzA6cGxy4HkiVgMoBUYBy6q5trL/A36MkoFIQli3zpLCG2/YiKRgjRrBT39qE9l+9CNv\n4hNNOhORGNqxw3Zhmzmz6i5sAOecY7WFXr3UhBRrSgYiEnP798O8edbhvGFD1eNZWbY1Z16elryI\nFSUDEfGM3w9LllhSeO+9qovjNWgAgwdbbaFjR29iTBVKBiISF4qL4aWXbIvO3burHu/a1ZJC375Q\nv37s40t2SgYiEle+/x7mz4cXX7QltStr1cr2bb7oIjj22NjHl6yUDEQkLvn9sGqVJYU336w6Cikj\nw2oJP/sZnH22zYqW2lMyEJG4t2MH/Pvf1oy0fXvV49nZVlsYOlQdzrWlZCAiCaOsDAoLrbawbFnV\n4w0a2LIXl1wCZ56p2kJNKBmISEJavx5mzbLd2EpLqx4/+WRLCkOGQOPGsY8v0SgZiEhC27vXZja/\n9FL4DudGjWx46sUXw+mnxz6+RKFkICJJ49NPLSm88YaNSqrsjDOstnD++ZYkxKFkICJJZ88eaz6a\nNSv8DOfGjS0hXHihJQj1LSgZiEgS8/vhv/+12sKbb8KBA1XPOflkW/ZiyBDIzIx9jPFCyUBEUsKu\nXbbPwssv205sldWvDz4f5OdDjx6pt4ezkoGIpJSK2sLLL8OCBeH7Fo47zpbVzsuzfRdSgZKBiKSs\n0lJb+mL2bFi9uurxtDTIybHawnnnwVFHxT7GWFEyEBHBOppnz7aO5127qh7PzLQZznl5cNJJsY8v\n2pQMRESCHDgAb79tzUgffVR1WW2A006DCy6AQYOgefPYxxgNSgYiIoexdat1Or/6Knz5ZdXjGRnQ\nu7clhp/8JLGX1lYyEBGJoLwcPvzQksLbb4cfopqZaTWFoUNtpnOizV1QMhARqYHdu23Owmuv2aik\ncE480WoLubmJs+eCkoGISC1t2mQdzq+/Hr4ZqV496NnTEoPPZ6uqxislAxGRI3ToEHz8sdUWFi6E\nffuqntOkCfTvb7WFs8+Ov0ltSgYiInVo715LCK+9BkuXhj+ndWtbGyk3F049NT76F5QMRESiZNs2\na0Z67TXYvDn8OR06WFIYNAiysmIbXzAlAxGRKKvYz3nePJvx/N134c/r3Nn2Xhg4MPbzF5QMRERi\nqKzMJrPNnQuLFoVfGyk93ZbByM2FPn1is/eCkoGIiEf27bN5C3PnwuLFNp+hsoYNbSRSbq6tphqt\niW1KBiIicWDnTpu/MG8erFwZ/pxjjoG+fa0ZqXt3q0HUlWglg8HAeCAdeBp4IMw5jwK5wF5gJLA8\nwrX3AHmAH9gRuCZcl4ySgYgktOJi27pz7lz4/PPw52Rm2lDVgQPhrLOOPDFEIxmkA+uAAUAxsAQY\nDqwJOmcIMCbw3BOYAOREuLYpUBK4/nqgC3B1mM9XMhCRpOD3w2efWW1hwYLwE9sAWrZ0EkOXLrWb\nwxCNZHAOcBf2Cx/g9sDz/UHnPAksAl4IvF8L+IAOLq4FGAs0CzoeTMlARJLOoUO258KCBdactH17\n+PNat4YBA+xx5pnu5zDUJhlkRDieRWjzzRbs13+kc7KA4yNcey/wS6xpKcd9yCIiia1ePRt62rkz\n3Hyz9SssWABvvQU7djjnff01TJ9uj7ZtrbYwcGB0Fs+LVAFx+7O8NmHdCZwA/B14pBbXi4gkvHr1\nrJ/g1lttUtukSXDJJdaPEOzLL+Ef/4ARI2xTnvHjbb7DoUN1E0ekmkEx0D7ofXvsF35157QLnFPf\nxbUA04E5hwugoKDgh9c+nw+fzxchZBGRxJSebiOLuneHP/zBlsBYsMCWxCgpcc7btg2ee84ebdrA\nCScUcuhQIW3a1H6dpEi/6DOwTuD+wFbgI6rvQM7BRg/lRLj2FGB94PrrgR5Yk1Fl6jMQkZR38KBN\nbps/3+Yy7NkT/rxWrWy46u23R2doaS7O8NApwH3A6MCxyYHniVhHcSkwClhWzbUAM4GOQDmwEfgt\n8HWYz1YyEBEJcvAgLFli/QuFheH3eP74Y006ExFJGWVlsGyZJYZFi+Dbb61cyUBEJEWVl8OKFda/\ncOutSgYiIimvNvMM4mx/HhER8YKSgYiIKBmIiIiSgYiIoGQgIiIoGYiICEoGIiKCkoGIiKBkICIi\nKBmIiAhKBiIigpKBiIigZCAiIigZiIgISgYiIoKSgYiIoGQgIiIoGYiICEoGIiKCkoGIiKBkICIi\nKBmIiAhKBiIigpKBiIigZCAiIigZiIgI7pPBYGAtsB647TDnPBo4vhI4y8W1DwJrAufPApq5jlpE\nROqUm2SQDkzEvtTPAIYDp1c6ZwhwMnAKcA0wycW184FOQBfgM2Bsbf8jUkFhYaHXIcQN3QuH7oVD\n9+LIuEkGPYANQBFwEJgB5Fc6Jw+YFnj9IZAJtI1w7QLgUNA17WoRf8rQP3SH7oVD98Khe3Fk3CSD\nLGBz0PstgTI35xzv4lqAK4E5LmIREZEocJMM/C7/VlotY7gTOABMr+X1IiISAznAvKD3Y6naifwk\n8POg92uBNi6uHQm8DzQ8zGdvwJKRHnrooYce7h8biIIMYCOQDRwFrCB8B3JFM08OsNjFtYOBT4BW\n0QhaRETqXi6wDss2FaN+RgceFSYGjq8Ezo5wLdhQ003A8sDjiWgELiIiIiIiCeJIJrQlm0j34grs\nHvwX62vpHLvQYs7NvwuA7kAZcHEsgvKIm3vhw2rYq4HCmETljUj3ohXWT7kCuxcjYxZZbE0FvgJW\nVXNOQn1vpmPNR9lAfSL3R/TE6Y9INm7uxTk4M7UHk9r3ouK8hcBrwCWxCi7G3NyLTKz/rWKuTrL2\nw7m5FwXAfYHXrYAdWN9lsjkX+4I/XDKo8fem12sT1XZCW5sYxRdLbu7Ff4BdgdfJPFHPzb0AuB6Y\nCWyPWWSx5+ZeXA68hM3jAfgmVsHFmJt7sQ04JvD6GCwZlMUovlh6F9hZzfEaf296nQxqO6EtGb8E\n3dyLYFeRvBP13P67yMdZ+sQfg7i84OZenAK0ABYBS4Ffxia0mHNzL57ClrnZijWP3Bib0OJOjb83\nva4+uf0/cOUJbcn4f/ya/Df1xWZt945SLF5zcy/GA7cHzk2j9pMe452be1EfG8HXHzgaq0EuxtqL\nk4mbe3EH1nzkA07Clr3pApREL6y4VaPvTa+TQTHQPuh9e5yq7uHOaRcoSzZu7gVYp/FTWJ9BddXE\nRObmXnTDmgnA2oZzsaaDV6IeXWy5uRebsaahfYHHO9gXYLIlAzf3ohdwb+D1RuBzoCNWY0olCfe9\neSQT2pKNm3txAtZmmhPTyGLPzb0I9gzJO5rIzb04DXgT62A9GutUPCN2IcaMm3sxDrgr8LoNlixa\nxCi+WMvGXQdywnxvHsmEtmQT6V48jXWIVUzU+yjWAcaQm38XFZI5GYC7e/F7bETRKuCGmEYXW5Hu\nRSvgVeyvdR5LAAAANklEQVS7YhXWuZ6Mnsf6RQ5gNcMrSd3vTRERERERERERERERERERERERERER\nERERERERqWv/D+3ChglzYzYwAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 41 }, { "cell_type": "raw", "metadata": {}, "source": [ "Ignore" ] }, { "cell_type": "code", "collapsed": false, "input": [ "mu = 5\n", "pmf = thinkbayes2.MakeExponentialPmf(mu, high=1.0)\n", "thinkplot.Pdf(pmf)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEACAYAAABRQBpkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VPW9//FXdrJvhCSQICA7IovsiEShCNiCS92ta6td\n0NpbvdYuinqrtVWr/rRqFRWXKiq2RRFE0QgIRJAt7GuAsIQQSAhZJ8n8/jiHyRCSzIQkc2Z5Px+P\necz5nvM9Mx8OcN5zdhAREREREREREREREREREREREREREXHbZGArsAN4oIk+z5vT1wNDzHGZwNfA\nJmAjcI9T/5lAPrDWfE1u66JFRKTthAA7gW5AGLAO6Negz1TgM3N4JLDSHE4DBpvDMcA2oK/Zfhj4\nn3apWEREWizYxfQRGGGQB9iA94HpDfpMA2abwzlAApAKHMYID4CTwBagi9N8QWdbtIiItC1XYdAF\n2O/Uzuf0FXpTfTIa9OmGsfsox2nc3Ri7lWZhBIiIiFjEVRjY3fychr/yneeLAT4Cfo2xhQDwEtAd\nYzfSIeBpN79HRETaQaiL6QcwDgSfkonxy7+5PhnmODCOM8wF3gH+49TniNPwa8AnjX15RGSCvaqi\n2EWJIiLiZBfQs6UzudoyWA30wtjNEw5cC8xr0GcecLM5PAooBgowthZmAZuBZxvMk+40fAWQ29iX\nV1UU86e/vIfdbg/418MPP2x5Dd7w0nLQstCyaP4FnOtivd4oV1sGNcAM4HOMM4tmYRwIvsuc/grG\nmURTMQ40lwG3mdPGAjcBGzBOHwV4EFgIPImxi8gO7HH6vDOs/H47NlsNYWGuShURkbPlzhp2gfly\n9kqD9oxG5ltG01seNzcx/gxlZZWs25TH8MEt3uoRERE3udpN5BWWrtxidQmWy8rKsroEr6DlUE/L\nop6WRet5+7n+9jFTHyQ1JYG5b9xPUJC3lysiYi1zPdnilaVPbBkUFBazY/chq8sQEfFbPhEGoF1F\nIiLtyWfCYMnKzVaXICLit7w+DELNU0p37j7EgcPHLK5GRMQ/eX0YDBtUf/2EdhWJiLQPrw+D8WP6\nO4a/Wb7JwkpERPyX14fBuJH9CAo2ytyweS9Hj52wuCIREf/j9WGQmBDD4AHdjIbdzpIV2lUkItLW\nvD4MALLGDnAMf7N8o4WViIj4J58Ig4tG1x83WJO7h5IT5RZWIyLif3wiDDp1jOe8fl0BqKutY6mu\nORARaVM+EQYAWWPqdxV9/a3OKhIRaUs+EwYXOYXB6vW7OFlWaWE1IiL+xWfCoEtaEr17dgagxlbD\n8lVbLa5IRMR/+EwYwOm7irK1q0hEpM34VhiMPc8xvGL1dsrKtatIRKQt+FQYnJORwrnd0wCorrax\nLEe7ikRE2oJPhQHAhHEDHcOLl+ZaWImIiP/wvTC46HzH8Mrvt1N6ssLCakRE/IPPhUFGejJ9e2UA\nUFtTy5IVugBNRKS1fC4MQLuKRETamk+GwSVOYbBq3U6KS8osrEZExPf5ZBikdUo47V5F2XrojYhI\nq/hkGMDpWwdfaVeRiEir+G4YXDgQgoIA47bWRcdLLa5IRMR3+WwYpCTHMeS87gDY6+r44psNFlck\nIuK7fDYMAH6QNcgx/EX2egsrERHxbT4dBhePHUBoWCgAW3fksze/0OKKRER8k0+HQVxsFGOG93G0\nF329zsJqRER8l0+HAcClTruKFmWvx263W1iNiIhv8vkwGD2sDzExkQAcPHyM3C37LK5IRMT3+HwY\nRESEcfHY+ofeLNKBZBGRFvP5MACYlDXYMfzV0lxsthoLqxER8T1+EQaDz+tGakoCACUnyshZs8Pi\nikREfIs7YTAZ2ArsAB5oos/z5vT1wBBzXCbwNbAJ2Ajc49Q/CfgC2A4sAhJaWriz4OBgJo6vf87B\ngsVrW/NxIiIBx1UYhAAvYARCf+B6oF+DPlOBnkAv4E7gJXO8DfgNMAAYBfwK6GtO+x1GGPQGFpvt\nVplyyRDH8LLvtlJyory1HykiEjBchcEIYCeQh7Fyfx+Y3qDPNGC2OZyD8Ss/FTgMnDrx/ySwBejS\nyDyzgcvPqnon3c9JpX+fTABqbDV88Y0OJIuIuMtVGHQB9ju186lfoTfXJ6NBn24Yu49yzHYqUGAO\nF5jtVpsyYahj+LMv17TFR4qIBIRQF9PdvYIrqJn5YoCPgF9jbCE09h1Nfs/MmTMdw1lZWWRlZTVZ\nxMSLBvL8a/OxVdewbecBduUd5txuaS6LFxHxVdnZ2WRnZ7f6cxquxBsaBczEOGYA8CBQBzzp1Odl\nIBtjFxIYB5vHY/ziDwM+BRYAzzrNsxXIwtiVlI5xoLkvZ7K39Irih558n8VLjDuYXnvFhdzz06kt\nml9ExJcFGbf2d7VuP4Or3USrMQ4MdwPCgWuBeQ36zANuNodHAcUYQRAEzAI2c3oQnJrnFnP4FuA/\nLS28Kc67ihZ9vY6amtq2+mgREb/lKgxqgBnA5xgr9TkYB4LvMl8AnwG7MQ40vwL80hw/FrgJuBhY\na75ObWH8BfgBxqmll5jtNjFiSE86JscBcLz4JCu/395WHy0i4rdavCnhYS3eTQTw0puf886H3wBw\n0ZgBPPGHG9u6LhERr9Reu4l80tSJ9buKvv1uqx6JKSLigl+GwTkZKZw/oBsAtTW1LNBppiIizfLL\nMACYdukwx/C8Raupq6uzsBoREe/mt2FwyYUDHc85OHCwiDW5eyyuSETEe/ltGEREhHHpxfW3tv7k\n89UWViMi4t38Ngzg9F1F2cs3UVxSZmE1IiLey6/DoGf39NNuXrfgK93aWkSkMX4dBgDTJg93DM9b\nuIqzuW5BRMTf+X0YTBg3kMjICAD25ReybmOetQWJiHghvw+DqMgILs0a5Gj/+7OcZnqLiAQmvw8D\ngCsuG+kYzl6+SVcki4g0EBBh0LN7+mlXJM9buMragkREvExAhAHAVT8c5Rj+78JV1NbqimQRkVMC\nJgzGj+5PUmIsAIVHS1i6crPFFYmIeI+ACYOwsFB+5HQR2sfzdSBZROSUgAkDgOmThxMcYvyRv1+/\ni7x9RyyuSETEOwRUGKSmJHDhyH6O9tz5Ky2sRkTEewRUGABc6XSa6WdfrqH0ZIWF1YiIeIeAC4Nh\ng86lR7c0ACorq/l0ke5mKiIScGEQFBTENdPHONoffbJSp5mKSMALuDAAmDR+EPFx0QAcPnJcp5mK\nSMALyDCIiAjj8qkjHO05/11uYTUiItYLyDAAuGLqSEJCQwDYsCmPLdvzLa5IRMQ6ARsGKclxTLzo\nfEf7g3naOhCRwBWwYQBwzbT6A8mLl+ZSUFhsYTUiItYJ6DDo26sLg5zuZvqBjh2ISIAK6DAAuOGq\nixzD/124ShehiUhACvgwGDO8N926dgKgoqKK/yz4zuKKREQ8L+DDIDg4mBuuHOdofzhvOdXVNRZW\nJCLieQEfBgCTsgbRMTkOgKJjpXyevc7iikREPEthgPGsA+czi/41dyl1dbpFhYgEDoWBafqUEURH\ndwBgX34hS1dusbgiERHPURiYYqI7cPmU+ltUzJ6Tjd1ut7AiERHPURg4ufbysYSHhwGwbecBVn6/\nw+KKREQ8Q2HgJDkx9rTnJM+e87W2DkQkILgTBpOBrcAO4IEm+jxvTl8PDHEa/zpQAOQ26D8TyAfW\nmq/Jblfczm68apzjBna5m/eyNnePxRWJiLQ/V2EQAryAsbLuD1wP9GvQZyrQE+gF3Am85DTtDRpf\n0duBZzCCYwiwsKWFt5fUlAQumzjU0Z4952sLqxER8QxXYTAC2AnkATbgfWB6gz7TgNnmcA6QAKSZ\n7aXA8SY+O6iFtXrMTVePJzjEWDSr1+1i49Z9FlckItK+XIVBF2C/UzvfHNfSPo25G2O30iyMAPEa\nXdKSmDR+kKP9+r++srAaEZH25yoM3D162vBXvqv5XgK6A4OBQ8DTbn6Px9x8TRZBwcbiyfl+Oxs2\n77W4IhGR9hPqYvoBINOpnYnxy7+5PhnmuOYccRp+DfikqY4zZ850DGdlZZGVleXio9vGOZkpXJo1\niIVfrQXgtXe+5PnH7/DId4uIuCs7O5vs7OxWf46r/fahwDZgAnAQ+A7jILLz5blTgRnm+yjgWfP9\nlG4YK/uBTuPSMbYIAH4DDAduaOT77Vae2pl/qIjr7/o7dbXGrSmee/wOhg0617J6RERcCQoKgrM4\nJutqN1ENxor+c2AzMAcjCO4yXwCfAbsxDjS/AvzSaf73gOVAb4zjCreZ458ENmAcMxiPEQheJyM9\n+bQzi15750tddyAifslrz+gxWbplAHD4SDHX3vkMNTbjttZPPXIro4f1trQmEZGmtNeWQcBL65TA\n9MnDHe1/vr1IWwci4ncUBm64+ZrxjnsWbd95kCUrNltckYhI21IYuKFjUhxX/bD+mPir73xJba2e\ndyAi/kNh4KYbrxpHZGQEAHv2FrB46QaLKxIRaTsKAzclJsRwzfT6p6G9+s5ibDY9K1lE/IPCoAWu\nv+JCYmMjATh4qIi583MsrkhEpG0oDFogNiaSW6+7xNF+8/2vOFFabmFFIiJtQ2HQQlddNpLO6ckA\nlJZW8OacbGsLEhFpAwqDFgoLC+VXt13qaH/0yQoOHD5mYUUiIq2nMDgL48cMYGD/cwCoranl5Tc/\nt7giEZHWURichaCgIGbcMcXR/mppLrlbdItrEfFdCoOzdF7frlwyrv5GrC/MWqDbVIiIz1IYtMLP\nb72U0DDjkRAbt+wj+9tNFlckInJ2FAat0CUtiR873abixdcXUFVls7AiEZGzozBopVuuzSI+LhqA\nQwXHefvDbyyuSESk5RQGrRQXG8Vdt0xytN+Zu4T8Q0UWViQi0nIKgzbwo0kX0K93BgC26hqefeVT\niysSEWkZhUEbCA4O5n9+/iMwnjDEilXbWJazxcVcIiLeQ2HQRvr3yTztiWjP/XO+DiaLiM9QGLSh\nu26eRFxsFAAHDx/jnY+WWFyRiIh7FAZtKD7u9IPJb3/0je5bJCI+QWHQxn40aRh9e9UfTH7mpU90\nZbKIeD2FQRsLCQnmt7+oP5i8cvU2vlyiR2SKiHdTGLSD/n0yufKykY7231/+lOKSMgsrEhFpnsKg\nnfz8lkmkpiQAUHKijOdenW9xRSIiTVMYtJPoqA7c96vpjvair9exfNU2CysSEWmawqAdjRneh0kX\nD3a0n3rxv5SVV1pYkYhI4xQG7ezXP7vMcSO7gsJiXp69yOKKRETOpDBoZwnx0fzm5z90tD+en8P6\nTXnWFSQi0giFgQdMvOh8xozoazTsdh5/7mMqKqutLUpExInCwAOCgoK475fTiIqKACD/wFFefH2h\nxVWJiNRTGHhIakoC995Zv7vo3/NXsmL1dgsrEhGppzDwoKkTh3LRmAGO9hPPzdXFaCLiFRQGHhQU\nFMT//mo6SYmxABQdK+WvL/xH9y4SEcspDDwsMSGGB399paP9zfJNLFi81sKKREQUBpYYM7wPl08Z\n4Wj//ZVPOKhbXYuIhdwJg8nAVmAH8EATfZ43p68HhjiNfx0oAHIb9E8CvgC2A4uABPdL9g8zfjqV\nLp2TASgvr+KxZz6itrbO4qpEJFC5CoMQ4AWMQOgPXA/0a9BnKtAT6AXcCbzkNO0Nc96GfocRBr2B\nxWY7oER2COfh+64hOMT4K9iwKY9X3/nS4qpEJFC5CoMRwE4gD7AB7wPTG/SZBsw2h3MwfuWnme2l\nwPFGPtd5ntnA5S0p2l8M6JPJT2+c4Gi//UG2TjcVEUu4CoMuwH6ndr45rqV9GkrF2H2E+Z7qor/f\n+snV4xk+pKej/djTH1JQWGxhRSISiFyFgbvnPAad5Xyn+gbsuZXBwcE8dN81dEyOA4xnH8z82wc6\nfiAiHhXqYvoBINOpnYnxy7+5PhnmuOYUYOxKOgykA0ea6jhz5kzHcFZWFllZWS4+2vckJcQw8/5r\nufv3s7DX1bFhUx7/fPsLfnHrpVaXJiJeLjs7m+zs7FZ/TsNf9A2FAtuACcBB4DuMg8hbnPpMBWaY\n76OAZ833U7oBnwADncb9FSgCnsQ4eJxA4weR7YF0QdbsOV/zz7e+cLSfeuRWRg/rbWFFIuJrgozn\nr7tat5/B1W6iGowV/efAZmAORhDcZb4APgN2YxxofgX4pdP87wHLMc4a2g/cZo7/C/ADjFNLLzHb\nAe8nV49n+NBejvYjT83hgK4/EBEPaHF6eFhAbRkAHC8+ya33vMDRohMA9OiWxitP3UVUZITFlYmI\nL2ivLQPxsMSEGP78+xsIDTMO5+zOO8yf/z5X9y8SkXalMPBC5/Xtyv2/nOZoZ3+7kdlzsq0rSET8\nnsLAS/1w0jB+PG2Mo/3q21+wLGdLM3OIiJw9hYEXu/uOKQw9v4ej/chTH5C3r8mzcEVEzprCwIuF\nhobw2O+uJz01ETBuaHf/o2/rgTgi0uYUBl4uIT6aJ/54ExER4QAcPFTEA4+9TVWVzeLKRMSfKAx8\nQK8e6Tx839VgnDLGxi37ePTpD6mr0y0rRKRtKAx8xPgxA7jnZ1Md7exvN/LC6wstrEhE/InCwIdc\nO30s11w+1tGe8+9lfDhvhYUViYi/UBj4mBm3T+GiMQMc7Wf/+SlLV262sCIR8QcKAx8TEhLMw7+9\nmgF9uxoj7HYeenIOazbstrYwEfFpCgMf1KFDOE/+6SY6pxvPUK6utvHAY2+zZXvDu4uLiLhHYeCj\nEhNi+Ptjt5KcFAsY1yD8z0NvsjuvoPkZRUQaoTDwYRnpyTz72O3ExUYBcKK0nHv/9Dr5h4osrkxE\nfI3CwMf16JbK04/cQlSUcYvromOl3PuH1/UcZRFpEYWBH+jfJ5Mn//QTwsPDADhUcJx7//gGR4+d\nsLgyEfEVCgM/MfT8Hvz59zcQEhoCwL78QmY8OIsjR0ssrkxEfIHCwI+MGd6HR+6/luAQ4691f34h\nMx58TbuMRMQlhYGfufjC83jsgescWwgHDhYx48FZHD6iQBCRpukZyH5q6crN/OGJ96itqQUgPTWR\n5x+/g85pSRZXJiLtSc9AltOMG9Wfx39/o+NZyocKjjPjd6+xd3+hxZWJiDfSloGfW7F6Ow/++R1s\n1TUAxMdF89TMm+nfJ9PiykSkPZztloHCIAB8t3YnD/7fO1RWVgPG7Swe/8ONjBzay+LKRKStKQyk\nWZu27ef+mW9RcsJ4ZGZIaAgP/fZqJl50vsWViUhbUhiIS3v3F/KbP71Rf6ppUBD33vlDrp422trC\nRKTNKAzELQWFxfzmT2+yd/8Rx7hrLh/LjNunEBKi8wlEfJ3CQNxWcqKc+x95i01b9znGjRnRl5n3\nX0N0VAcLKxOR1lIYSItUVlbz6NMf8s3yTY5x53ZP428P30xqSoKFlYlIaygMpMXq6up4efYi3v1o\niWNcclIsf/njTTr1VMRHKQzkrH26aDV/ffG/jquVw8PDuP9X05k6cajFlYlISykMpFXWbNjN7x9/\nl9LSCse4Ky4bxa9/NpUw8ypmEfF+CgNptX0HjvLg/71D3r76M43O69eVx353PZ06xltYmYi4S2Eg\nbaK8oorHn/2Yr5flOsYlJsTw6APXMfT8HhZWJiLuUBhIm7Hb7bz/72X8483PqautAyA4JJg7bpjA\nT64er+sRRLyYwkDa3JoNu3noyfc5XnzSMW7IwB48dN/V2m0k4qUUBtIuCgqLeeRvH7B+U55jXFxs\nFH/4zVVcOLKfdYWJSKPa83kGk4GtwA7ggSb6PG9OXw8McWPemUA+sNZ8TW5J0eI5qSkJ/L8nfsrt\nN0wgKNj453KitJwHHn2bZ17+xHEnVBHxba7SIwTYBkwEDgCrgOuBLU59pgIzzPeRwHPAKBfzPgyU\nAs+4+H5tGXiRtbl7ePTpDzhSWOIYl5mRwh/uvZKB/c6xsDIROaW9tgxGADuBPMAGvA9Mb9BnGjDb\nHM4BEoA0N+b19l1U0sCQgd158/m7uXBU/e6h/fmF/OJ/X+UfbyykqspmYXUi0hquwqALsN+pnW+O\nc6dPZxfz3o2xW2kWRoCID4iPi+Ivf7yJB+65gqioCADsdXW8+9ESbr/3RbZsz7e4QhE5G64uLXV3\nH01Lf+W/BDxqDj8GPA3c0VjHmTNnOoazsrLIyspq4VdJWwsKCmLapcMZPrgnTzz3Md+v3wVA3r4j\n3Hnfy1w7fSx33DiByA7hFlcq4v+ys7PJzs5u9ee4WomPwjjYe+oA74NAHfCkU5+XgWyM3UBgHDAe\nD3R3Y16AbsAnwMBGvl/HDLxcXV0dH8/P4R9vfE5VVf3B5LROifz2l9MYM7yPhdWJBJ72OmawGuiF\nscIOB64F5jXoMw+42RweBRQDBS7mTXea/wogF/FJwcHB/PhHo3nrxbsZMrD+CuXDR45z/8zZ/PGJ\n9zh67ISFFYqIO9xJjynAsxhnB80CngDuMqe9Yr6/gLEFUAbcBqxpZl6At4DBGLuh9pifV9DId2vL\nwIfY7XY++3INL8xawInScsf46OgO3Hb9Jfz4h6N00zuRdqaLzsRrHC8+yYuvL2TB4jWnje+akcI9\nP7uM0cN6W1SZiP9TGIjXWb1+F0/9Yx778wtPGz96eB/u+dlldO3S0aLKRPyXwkC8ks1Ww0efrOD1\n976ivLzKMT4kNITLp4zg1usuJikhxsIKRfyLwkC8WtHxUv751hd8+sX34PR3GhkZwQ1XXsh1V1xI\nVGSEhRWK+AeFgfiELdvzef61z9jgdOM7MJ6ZcOt1FzN98nAdZBZpBYWB+Ay73c63323lpTc/P+2p\namDcGO/ma7O4bOJQhYLIWVAYiM+pra1jweI1zPrX4tNufgcKBZGzpTAQn1VVZWPupyt5d+4SikvK\nTpuWmpLANZePZdqlw3RMQcQNCgPxeRWV1fx7fk6joRATE8kVU0dy9bTRJCfGWlShiPdTGIjfqKis\n5j+f5fDu3KWnPXITIDQslCmXDOb6K8ZxTmaKRRWKeC+FgfidyspqPlu8lvf+vYyDh4rOmD52ZD+u\nvGwkI4b0JDjYnYf2ifg/hYH4rdraOpau3My7c5eyedv+M6Z3Tk/myqkjmDrxAuLjoiyoUMR7KAzE\n79ntdtZvyuNfHy/j25wtZ0wPDw9j4kUDuWLqSPr3ybSgQhHrKQwkoOzNL+S/C1cx/4vvOXmy4ozp\nvc7tzJQJQ5g0fhCJut2FBBCFgQSkyspqvlyygbnzV7J958EzpoeEhjB6WG8um3gBo4f11jUL4vcU\nBhLQ7HY7W7bn8/H8HBYvzaW62nZGn/i4aH6QNYhJWYPo3zvj1H8aEb+iMBAxnSyrZPHSXBYsXkPu\n5r2N9knrlMgl4wYyYdxA+vTsrGAQv6EwEGnE3vxCFn61loVfrT3jlhendE5PZsK4gVxy4Xn06pGu\nYBCfpjAQaUZtbR3fb9jNl9+s55sVmxs96AzGFsPYkX0ZN6ofQ87rTmhoiIcrFWkdhYGIm2y2GnLW\n7OTrZbksWbn5tIfuOIuO7sDoYX0YN6ofI4f2IjYm0sOVirScwkDkLFRV2fhu7Q6+XJLLitXbKCur\nbLRfUHAwA/pkMGJIL0Ze0It+vTIICdFVz+J9FAYirWSz1bBuUx5LV25h2cotFBQWN9k3NjaSYYN6\nMvKCXlxwfg86pyV5sFKRpikMRNqQ3W5n557DLMvZwrKcrWzdeeC0x3U2lJqSwJCB3Rk8sDtDBnan\nS1qSDkSLJRQGIu2ouKSMVet2kvP9Dr5bu4OiY6XN9k/pGM/g87oxZGAPBg/oRmaXZN1MTzxCYSDi\nIXa7nd17C8j5fger1u9kw6a9VFZWNztPTEwk/ftkcF6frvTvk0n/3hm6qZ60C4WBiEVqamrZtusg\n63L3sCZ3Dxs25zV5hpKzrhkpDOiTSf8+mfTp2Zme3dKIiAjzQMXizxQGIl6itraO7bsOsm7jHtZu\nzGPjln2UnChzOV9QcDDdMlPofW5nevVIp3ePdHqf21mntEqLKAxEvJTdbudgwXE2bd3Hxq372bxt\nP9t3H6K2ptat+TunJdGjWyrdu6bSvWsnunftxDkZKdqKkEYpDER8SFWVjW27DrJp23627jjA9l0H\n2XfgaLNnLDkLCg6mS3oS3bt2osc59SGR2bmjQiLAKQxEfFx5RRU79xxm+66D7Nh9iO27D7Irr8Dt\nLQgAgoJI7RhPZpdkMtKTyejSkczOyWR27kh6aiLh4bqFt79TGIj4IZuthj37jji9Cti99wgHDx9z\neyvilKDgYNJS4snonEyXtCTSUxNJS00kNSWe9NREkhJidPqrH1AYiASQyspq9uYXsmffEXbvLWDP\nviPk7TvCoSPF2OvqzuozQ8NCSeuUQHqnBFJTEoyw6JRAWqcEUjrG0zExVrugfIDCQESorq7hUMFx\n9h88yv6DReQfOEr+oSL2Hyii4GhJi7cmGoqPiyY5KZaOybF0So43hpPiSEk23jsmx5GUEKP7NllI\nYSAizaqqsnHg0DHyDxVxsOA4BYXFHDLfDxcUc6K0vE2+Jyg4mIS4KBLio0mMjyYxIea04VPv8XFR\nJCXEEBPdQbfuaEMKAxFplbLySgqOlHDoyHEzJEocYXH02AmOHiulrvbsdkE1JyQ0hMT4aOLjooiL\niSIuNpLYmEjiYp2GYyKJjTXfYyKJj4siKjJCIdIIhYGItKva2jqKT5RxtOgEhUWlRkAUneDoMaNd\ndLyUwqMlFJe4vsCuLQSHBBMbXR8S0VERREVGEB0VQWRkBDHRHYiOjCAqKsKc1oHoU8NREURHRhAd\n1cHvzrBqzzCYDDwLhACvAU820ud5YApQDtwKrHUxbxIwBzgHyAOuARq7X7DCQMTH1NTUUnyijOPF\nZRwvKaO45CTHissoLinjePFJc5wx7XjxSSoqXN+6oz2FhoUaoREZTkREGJEdwukQUT8c2SGMiIhw\nc3yY8erg1HYaH9nBmC8iIoyI8DDCQkMIDw/16BZMe4VBCLANmAgcAFYB1wNbnPpMBWaY7yOB54BR\nLub9K3DUfH8ASAR+18j3KwxM2dnZZGVlWV2G5bQc6vnLsqisrOZ4SRknSss5cbKC0tIK4/1kBSdK\nKyg9Wc5bsAVCAAAEzElEQVSJUnO4rH66c4iUFO0lPvkcC/8UzQsLDyU8LJSI8DDCw0IJDzdf5riw\n8BBjvPmKiAgjLCyUCLNPeHgoYWGhhIWGOL2HEBYaQmhoqNNwCIMGdIOzCANX20cjgJ0Yv94B3gem\nc3oYTANmm8M5QAKQBnRvZt5pwHhz/Gwgm8bDQEz+8h+/tbQc6vnLsujQIZz0DuGkpya2aD6brcYR\nHn/96xPccvttVFRUcbKskrKKKsrLqygrr6K8wng/WVZJeYXRPjXtZHllyy7qO0u26hps1TVNPknP\nG7gKgy7Afqd2Psavf1d9ugCdm5k3FSgwhwvMtoiI28LCQklOjCU5MZZOHeMZObRXiz/DbrdTXV3D\nyfJKKqtsVFZWU1llo6Ky2jFcWVlNRaWNqqpqKiqrqaiyUeXUr6KqmqpKGxVV1VRW2qiqslFVbcNW\nU0tVdQ01tpp2+NO3PVdh4O4+Gnc2SYKa+Dx7C75HRKTNBAUFOfbxt5e6ujpstlqqqm1U22rMgKhv\nV1cb42y2Gkfb6GfDZk6rNqfZbLXU1JjvtXXmuBpqamqpqTHayz9rnz/HKGChU/tBjH38zl4GrnNq\nb8X4pd/cvFsxdiUBpJvtxuykPiz00ksvvfRy/dpJOwgFdgHdgHBgHdCvQZ+pwKksGgWsdGPeUweO\nwThW8Jc2r1xERNrUFIyzgnZi/LoHuMt8nfKCOX09MNTFvGCcWvolsB1YhHHQWUREREREAtlkjGMG\nOzjzeMQpz5vT1wNDPFSXFVwtixsxlsEG4FvgfM+V5nHu/LsAGA7UAFd6oiiLuLMssjAu9tyIcaq2\nv3K1LDpiHKtch7EsbvVYZZ71OsaZmLnN9PGp9WYIxi6kbkAYro9JjKT+mIS/cWdZjAbizeHJBPay\nONXvK+BT4CpPFedh7iyLBGATkGG2O3qqOA9zZ1nMBJ4whzsCRbg+a9IXjcNYwTcVBi1eb1p9n1nn\ni9ps1F+Y5qyxi9r88boEd5bFCqDEHM6h/j+/v3FnWQDcDXwEFHqsMs9zZ1ncAMzFuJYHjKv7/ZE7\ny+IQEGcOx2GEgW+c6N8yS4HjzUxv8XrT6jBo6oI1V338cSXozrJwdgf1ye9v3P13MR14yWzbPVCX\nFdxZFr0wTsr4GlgN/MQzpXmcO8viVWAAcBBj98ivPVOa12nxetPqzSd3/wM3vKjNH//jt+TPdDFw\nOzC2nWqxmjvL4lmM05LtGP8+vP0OvGfLnWURhnEW3wQgCmMLciXG/mJ/4s6y+D3G7qMs4FzgC2AQ\nUNp+ZXmtFq03rQ6DA0CmUzuT+k3dpvpkmOP8jTvLAoyDxq9iHDNobjPRl7mzLC7A2E0Axr7hKRi7\nDua1e3We5c6y2I+xa6jCfC3BWAH6Wxi4syzGAH82h3cBe4A+GFtMgcTn1putuajN37izLLpi7DMd\n5dHKPM+dZeHsDfz3bCJ3lkVfjOt2QjC2DHKB/p4r0WPcWRbPAA+bw6kYYZHkofo8rRvuHUD2mfVm\nay5q8zeulsVrGAfE1pqv7zxdoAe58+/iFH8OA3BvWdyHcUZRLnCPR6vzLFfLoiPwCca6Ihfj4Lo/\neg/juEg1xpbh7QTuelNERERERERERERERERERERERERERERERERERNra/wcK/e6aP2qyxwAAAABJ\nRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 42 }, { "cell_type": "raw", "metadata": {}, "source": [ "Ignore" ] }, { "cell_type": "code", "collapsed": false, "input": [ "metapmf = thinkbayes2.Pmf()\n", "for lam, prob in pmf.Items():\n", " if lam==0: continue\n", " pmf = thinkbayes2.MakeExponentialPmf(lam, high=30)\n", " metapmf[pmf] = prob\n", " \n", "interarrival = thinkbayes2.MakeMixture(metapmf)\n", "thinkplot.Pdf(interarrival)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEACAYAAAC+gnFaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4HOWB5/Fvqw/dlizJSLYsWfJ9YHwwvuIYRCDE+CFA\nQhJCDjzkgEkgYQg7Azy7Ozi7sxNgJ0yWMDFk4hCSDIEsGRiz5iYITAAb4/uQL9mWZVuyLUuyjm61\nWur9o0qtVltSt2S1qo/f53n66beq3up622X1r+t9q7pARERERERERERERERERERERERERESkj5VA\nFXAQuH+AOo+by3cAC8x5acAmYDuwF/hJUP084E3gAPAGkDvirRYRkRFhBw4BZYAT40N9VkidVcAr\nZnkJ8FHQsgzz2WHOX25OPwr8vVm+H3h4JBstIiIjZxnwWtD0A+Yj2JPALUHTVUBhSJ0M4GNgdj91\nisxpERGxQEqY5cXA8aDpWnNeuDoTzbId4yiiHngHo4sIjBCoN8v1XBgcIiIySsIFgT/C17ENsF4X\nMB8jGK4AKgbYRqTbERGREeYIs/wEUBI0XYLxjX+wOhPNecGagQ3A5UAlxlFAEVAHjAdO97fxKVOm\n+A8fPhymiSIiEuQwMHUoK4Q7ItgCTMMYLHZhjAWsD6mzHrjNLC8FmjA+6AvoPRsoHfgsRjdRzzqr\nzfJq4KX+Nn748GH8fn9CPh566CHL26D3p/en95d4D2BK2E/+EOGOCHzA3cDrGP3964B9wJ3m8qcw\nzhhahXF2URtwu7lsPPAMRtikAL8D3jaXPQz8Efg2cBT4ylAbLiIiIyNcEAC8aj6CPRUyfXc/6+0C\nFg7wmueAayLYtoiIRFm4riGJkoqKCqubEFV6f/FN7y+5hJ7tE2v8Zp+XiIhEwGazwRA/22P+iMDj\n8VndBBGRhBbzQXCyrsXqJoiIJDQFgYhIkov5IKg9qSAQEYmmmA8CHRGIiERX7AfBKQWBiEg0xX4Q\n1LegU0hFRKIn5oPA4/FxrtFtdTNERBJWzAcBQK26h0REoiYugkDjBCIi0aMgEBFJcnERBOoaEhGJ\nnrgIghOnzlvdBBGRhBUXQXC2oZ12d6fVzRARSUhxEQQAJ07qqEBEJBriJgiOn1AQiIhEg4JARCTJ\nKQhERJJc3ARBrcYIRESiIuaDICXFuPXm6bNtum2liEgUxHwQjC/MDpR1VCAiMvJiPghKiscEyhon\nEBEZeXEWBM0WtkREJDHFfBCUTswJlI8eVxCIiIy0mA+CstLeIDimIBARGXExHwRFl2ThctoBaGxy\n09LaYXGLREQSSyRBsBKoAg4C9w9Q53Fz+Q5ggTmvBHgH2APsBn4YVH8NUAtsMx8rB9q43Z7CxKBx\nAh0ViIiMrHBBYAeewPigng3cCswKqbMKmApMA+4A1przO4F7gTnAUuAuYKa5zA88hhEaC4DXBmvE\npBKNE4iIREu4IFgMHAKOYnywPwfcGFLnBuAZs7wJyAUKgTpguzm/FdgHFAetZ4u0kWUluYFyjYJA\nRGREhQuCYuB40HQtfT/MB6ozMaROGcY3/01B836A0ZW0DiM8BlRaogFjEZFoCRcE/ghfJ/TbffB6\nWcALwD0YRwZgdB+VA/OBU8BPB3vxSUGnkNacaKa7O9JmiYhIOI4wy09gDPr2KMH4xj9YnYnmPAAn\n8Cfg98BLQXVOB5V/Bbw8UAPWrFkDQHXVfrJzpzOuaA51p1uZUJQ90CoiIkmjsrKSysrKi3qNcP30\nDmA/cDVwEtiMMWC8L6jOKuBu83kp8DPz2YYxdtCAMWgcbDzGkQDmskXA1/rZvt/vN779//jR99ix\nuw6A++5axvIlJf1UFxFJbjabDYYwBgvhu4Z8GB/yrwN7gecxQuBO8wHwClCNMaj8FPB9c/5y4BvA\nVVx4mugjwE6MMYIruTAoLlA+qXcY4cixpnDVRUQkQuG6hgBeNR/BngqZvruf9d5n4KC5LYLt9jE5\nKAiO1igIRERGSsxfWdxDRwQiItERN0EwvjCb1FTjAKax2U1jk8fiFomIJIa4CYKUFBvlpcFHBY0W\ntkZEJHHETRAAlJWqe0hEZKTFVRAEjxNUKwhEREZEHAeBuoZEREZCXAXBpJIcnA7j3gT1p1s536J7\nE4iIXKy4CgKnw97nqOBQ9TkLWyMikhjiKggApk3OC5QPHFYQiIhcrLgLgqlBQaAjAhGRixd3QTBt\nSm8QHKw+R8+P0omIyPDEXRCML8wiM8MFQEtrB/Vn2ixukYhIfIu7ILDZbOoeEhEZQXEXBNB3wPig\ngkBE5KLEZxBMURCIiIyUuAyCqeW9QVB9pBGfr9vC1oiIxLe4DIKxuWmMy88EwNvZRU1ts8UtEhGJ\nX3EZBND3egJ1D4mIDF/cBkHwOIHOHBIRGb74DQIdEYiIjIi4DYLJZbnYbDYAamrP0+7utLhFIiLx\nKW6DID3NSenEHHPKr6MCEZFhitsgAJg9oyBQ3lt1xsKWiIjErzgPgnGB8t79CgIRkeGI6yCYFXRE\ncODQObydXRa2RkQkPsV1EOTlpjO+MBuATl+XTiMVERmGuA4CCBknUPeQiMiQxX8QzAweJzhrYUtE\nROJT/AdB0IBx1cGzdHXpB+hERIYikiBYCVQBB4H7B6jzuLl8B7DAnFcCvAPsAXYDPwyqnwe8CRwA\n3gByh9rwHpcUZFCQnwGAx+PjyLGm4b6UiEhSChcEduAJjDCYDdwKzAqpswqYCkwD7gDWmvM7gXuB\nOcBS4C5gprnsAYwgmA68bU4Pi81mY/b03qOCPbqeQERkSMIFwWLgEHAU44P9OeDGkDo3AM+Y5U0Y\n3+4LgTpguzm/FdgHFPezzjPATcNqvWmWBoxFRIYtXBAUA8eDpmvp/TAfrM7EkDplGF1Gm8zpQqDe\nLNeb08M2J2jAeN+Bs3R3+y/m5UREkoojzPJIP1Ftg6yXBbwA3INxZNDfNgbczpo1awLliooKKioq\nLqhTPD6bMdmpnG/poLXNy/ETzUwqGfawg4hI3KisrKSysvKiXiNcEJzAGPTtUYLxjX+wOhPNeQBO\n4E/A74GXgurUA0UY3UfjgdMDNSA4CAZis9mYM3McH35sNG1P1RkFgYgkhdAvyD/+8Y+H/Brhuoa2\nYAwClwEu4BZgfUid9cBtZnkp0ITxQW8D1gF7gZ/1s85qs7yaviExLMGnke7cM2CuiIhIiHBHBD7g\nbuB1jDOI1mEM+t5pLn8KeAXjzKFDQBtwu7lsOfANYCewzZz3IPAa8DDwR+DbGAPRX7nYNzLv0t5h\nhl17T+PzdeNwxP1lEiIiURfatx9r/H5/ZMMUfr+fO3+0gbMN7QD843+9qs9RgohIMjBv2DWkz/aE\n+cpss9lYMLcoML19V/0gtUVEpEfCBAHA/KAg2LarzsKWiIjEj4QKgrmzLwncx/jwkUbOt3RY3CIR\nkdiXUEGQleli+pQ8c8rPzj3qHhIRCSehggBg3qXqHhIRGYqEC4IFl/UdMI70rCMRkWSVcEEwtXws\nWZkuABqb3NTUNlvcIhGR2JZwQWC3pzB3du/FZTqNVERkcAkXBADz5/YGgcYJREQGl5BBEHxh2d79\nZ2h3d1rYGhGR2JaQQVCQn0F5qfHroz5fN1t3nLK4RSIisSshgwBg8eW998/ZvPWkhS0REYltCRsE\nS4KC4JPtp/B2dlnYGhGR2JWwQTCpJIfCcZkAuD2d7N6rexSIiPQnYYPAZrOpe0hEJAIJGwTQt3vo\n460ndVN7EZF+JHQQzJiaz5jsVAAam90cONxgcYtERGJPQgeB3Z7CogUTAtObPzlhYWtERGJTQgcB\n9D2N9KNPTuhH6EREQiR8EMybU0hqqgOAuvpWjp84b3GLRERiS8IHgctlZ+Fl4wPTGz+ssbA1IiKx\nJ+GDAGDFspJAeeOHNeoeEhEJkhRBsHDeeDIzjHsUnD7bRtVBnT0kItIjKYLA5bSzbNHEwLS6h0RE\neiVFEABc8anSQPkvm47j83Vb2BoRkdiRNEEwe8Y48vMyAGhp7WD7bt2wRkQEkigIUlJsrFjaO2j8\n3gfqHhIRgSQKAoAVyyYFypu3nsTt0Z3LREQiCYKVQBVwELh/gDqPm8t3AAuC5v8aqAd2hdRfA9QC\n28zHyohbfBHKSnMoKc4BwOv1sUk/OSEiEjYI7MATGB/Us4FbgVkhdVYBU4FpwB3A2qBlT9P/h7wf\neAwjNBYArw214cNhs9m4YlnvoHHlxmOjsVkRkZgWLggWA4eAo0An8BxwY0idG4BnzPImIBfouXv8\nRqBxgNe2DbGtI+LK5ZOw2YxN79xbz6m6FiuaISISM8IFQTFwPGi61pw31Dr9+QFGV9I6jPAYFQX5\nGVw+v/cnJ95898hobVpEJCY5wiyP9LcYQr/dh1tvLfA/zPL/BH4KfLu/imvWrAmUKyoqqKioiLBJ\nA7v2qsls2WbcsezP7x3h1pvn4HTYL/p1RURGW2VlJZWVlRf1GuG6Z5ZiDOz29PM/CHQDjwTVeRKo\nxOg2AmNg+UqMQWKAMuBlYO4A2xhsuT8avwvU1dXN39z3Cg3n2gG4765lLF9SEmYtEZHYZ3Z9D6nr\nPVzX0BaMQeAywAXcAqwPqbMeuM0sLwWa6A2BgYwPKn+BC88qiiq7PYWrrygPTL/xzuHR3LyISEwJ\nFwQ+4G7gdWAv8DywD7jTfAC8AlRjDCo/BXw/aP0/AB8A0zHGEW435z8C7MQYI7gSuPci38eQXXNl\neWDQeNfe05zUoLGIJClLztwZgqh0DfX4p395PzBWcNOqGdz21XlR25aIyGiIRtdQQrv2qsmB8tvv\nHcXr7bKwNSIi1kjqIFgwt4iC/N4foqv8iy4wE5Hkk9RBYLencP210wLTL792gO5u3b1MRJJLUgcB\nwNVXlpOe5gTgxKnzbN15yuIWiYiMrqQPgswMF9dU9J5Kuv7VAxa2RkRk9CV9EABcf+00UlKMQfbd\n+05TfXSgn0cSEUk8CgJgXEEmn1rce2WxjgpEJJkoCEw3rJweKP9l83HONrRb2BoRkdGjIDBNnZzH\n7BnjAOO3iF7cUGVxi0RERoeCIMgXr58ZKL9ZWa2jAhFJCgqCIAsuK2La5HwAfD4dFYhIclAQBLHZ\nbHz1i3MC0zoqEJFkoCAIMX9uIdOn9B4V/Mf/01GBiCQ2BUGI0KOCt97VUYGIJDYFQT/mXVrIjKm9\nRwX/d/1ei1skIhI9CoJ+XHBUUHmEmtpmC1skIhI9CoIBXDankMvmFALg9/v57XM7LW6RiEh0KAgG\nYLPZWP3VefTc6GfrzlPs2BPuVswiIvFHQTCI8km5XH1FWWD6N8/uoKur27oGiYhEgYIgjFtvvhSX\nywHAseNNVL6vu5iJSGJREISRNzadm1bNCEw/+8Ju2t2dFrZIRGRkKQgicNOqGYzNTQegsdnNsy/s\nsrhFIiIjR0EQgbQ0B7d/bV5g+pU3D3Oo+pyFLRIRGTkKgggtX1LCvEuLzCk/Tz79iQaORSQhKAgi\nZLPZuHP1QlxOOwDVxxp59a1DFrdKROTiKQiGoKgwiy/dOCsw/eyfdut3iEQk7ikIhuiG62ZQPH4M\nAB6Pj1/8egt+v9/iVomIDJ+CYIhcTjvf+9bl9FxxvH1XHW+8U21to0RELkIkQbASqAIOAvcPUOdx\nc/kOYEHQ/F8D9UDo+ZZ5wJvAAeANIDfyJltv9oxxfW52//SzOzhV12Jhi0REhi9cENiBJzDCYDZw\nKzArpM4qYCowDbgDWBu07Glz3VAPYATBdOBtczqufO1Ll1JSbHQReb0+fv5vH+ssIhGJS+GCYDFw\nCDgKdALPATeG1LkBeMYsb8L4dt9znuVGoLGf1w1e5xngpqE0Oha4XHZ+cMdi7Hbjn7Dq4Fle2rDf\n4laJiAxduCAoBo4HTdea84ZaJ1QhRpcR5nNhmPoxaWp5Hl++cXZg+g//sZu9+89Y2CIRkaFzhFke\n6ekwtmGu11N3wPpr1qwJlCsqKqioqBjCS0ffzZ+fybadp9h/qIHubj8//dePeOwfP0vOmDSrmyYi\nSaCyspLKysqLeo3QD/BQS4E19PbzPwh0A48E1XkSqMToNgJjYPlKer/xlwEvA3OD1qkCKoA6YDzw\nDjCzn+374+HUzLMN7dz339+kpbUDgHmXFvHf7vt0oNtIRGS02Gw2CP/Z3ke4T6otGIPAZYALuAVY\nH1JnPXCbWV4KNNEbAgNZD6w2y6uBlyJrbmwqyM/gnjsXB6Z37K7jhfX7LGyRiEjkwgWBD7gbeB3Y\nCzwP7APuNB8ArwDVGIPKTwHfD1r/D8AHGGcHHQduN+c/DHwW4/TRz5jTcW3hvPHc/PneE6qef3Ev\nH287aWGLREQiM6TDBwvERddQj66ubtY88i57qowB47Q0Bw//w9WUTsyxuGUikiyi0TUkQ2C3p3Df\nXcsYl58JGD9B8U+PvU/zeY/FLRMRGZiCYITl5qTx4L3LSU01Tsg6fbaNRx//kE5fl8UtExHpn4Ig\nCspKc/nR95bSc3S278AZfv7Lj+nujp9uLhFJHgqCKFm0cALfvKX3jNn3P6rh6X/frl8qFZGYoyCI\noptWzeBzn5kSmN7w5kH+9LJOKxWR2KIgiCKbzcZ3vrmAZYtKAvOefWE3r//5sIWtEhHpS0EQZXZ7\nCvf8zWIum937c0pP/WYrb797xMJWiYj0UhCMApfTzt/f8ymmlOeZc/z867otCgMRiQkKglGSke7k\nH/5uBZMnjTXnKAxEJDboyuJR1tLawY8feY/qY723afjWNxZw/bXTLGyViCSK4VxZrCCwQH9h8OUb\nZ/PVL87p2YkiIsOiIIgjra1e/tdjG9l/qCEw73OfmcJ3vrlAP18tIsOmIIgzHo+Pf37iQ7buPBWY\nt2jBBP72e0tIT3Na2DIRiVcKgjjk83Xz83/bzMYPawLzyktzefDeT1OQn2Fhy0QkHikI4lR3t5/f\nPb+T/3x1f2De2Jx0Hrx3OVMn5w2ypohIXwqCOPdmZTW/fGYrXV3dADgddr67egHXXDnZ4paJSLxQ\nECSAnXvq+d8//5C2dm9g3tVXlPPd2xbictktbJmIxAMFQYI4WdfCo49/QE1tc2De5Elj+dFdS5lQ\nlG1hy0Qk1ikIEojH4+PJpz/hvQ+PBealpjr47jcXcNWKMl1vICL9UhAkGL/fz6tvHeLpZ3cExg0A\nli8p4c7Vl5OV5bKwdSISixQECar6aCP/snYTJ06dD8wbm5POHX+9kCWXF1vYMhGJNQqCBObx+Pj1\nv2/nrXer+8z/1OISvvPNBeTmpFnUMhGJJQqCJLDpkxP88jdbaWx2B+ZlZ6Xy7W/MZ8WyUo0diCQ5\nBUGSaG318vQftvPOxqN95l82p5BvfX0+pRNzrGmYiFhOQZBktu2sY+3TWzjb0B6Yl5JiY9Vnp/GV\nm2aTlanBZJFkoyBIQu3uTp59YRevvnWY4H+rMdmpfP3Lc/nMijL9mqlIElEQJLGjNU2s+/029lSd\n6TO/dGIOX7v5UhYtnKDxA5EkoCBIcn6/nw821/LMczv6dBcBTJucz9e/fCmXzSm0qHUiMhqiFQQr\ngZ8BduBXwCP91HkcuA5oB/4a2BZm3TXAd4Cer68PAq/187oKgmHo6PDx4ob9/Oer++no8PVZdtns\nQr7yhdnMnjHOotaJSDRFIwjswH7gGuAE8DFwK7AvqM4q4G7zeQnwf4ClYdZ9CGgBHguzfQXBRWhq\n9vDihipee+swnb6uPstmTi/g5s/PYuFlReoyEkkg0QiCZRgf2ivN6QfM54eD6jwJvAM8b05XARVA\n+SDrPgS0Aj8Ns30FwQg429DOH1/ay583HqG7u++/Z1lpLl+8fibLFk3UoLJIAhhOEIT7yy8GjgdN\n15rzIqkzIcy6PwB2AOuA3MibLENVkJ/B97/9Vzz+8EquvqK8zwf+0ZomHvvFR3zvv7zCixuqaGnt\nsLClImKFcEEQ6dfxofYtrMU4YpgPnCL8kYGMgAlF2dz1nUWs/edVXP+56bhcjsCysw3t/O75nXz3\nbzfwi3VbOFrTZGFLRWQ0OcIsPwGUBE2XYHyzH6zORLOOc5B1TwfN/xXw8kANWLNmTaBcUVFBRUVF\nmCZLOAX5GXzr6/P50g2z2PDGQV57+3DgSMDr9fHWu9W89W41M6cVcPUV5XxqyUTS05wWt1pE+lNZ\nWUllZeVFvUa4b/IOjAHfq4GTwGYGHyxeinGW0NIw647HOBIAuBdYBHytn+1rjGAUeL1dvP9RDa+8\neYjqY40XLE9Lc7BiaSnXXFnO1Ml5GlwWiWHROn30OnpPAV0H/AS401z2lPn8BMagcBtwO7B1kHUB\nfovRLeQHjpivV9/PthUEo8jv97PvwFk2vHGIzVtP9LkHQo/i8WP49NISViwr1d3SRGKQLiiTEdPY\n5OHdD47x9rtH+twHIdiU8jxWLC1h+ZIS8vMyRrmFItIfBYGMuJ6jhLffO8IHm2svuEDNYGPOzAKW\n/NVEFi+cwCUFmaPeThExKAgkqjweH1u2n2TjhzVs21WHz3dh1xFAeWkuixYWs3jhBMon5WpMQWQU\nKQhk1LS2etn0yQk2fljDrn2nGWg/FeRnsHjhBOZdWsSls8bp7CORKFMQiCUamzxs3nqCj7eeZOfe\n+gGPFOz2FGZOy2fB3CLmzy2irDSXlJRY/y8oEl8UBGK5dncn23bWsXnrCT7Zfop2d+eAdXPGpDF3\n9iXMmTmO2TPGMXFCtrqRRC6SgkBiis/XzZ79Z9ixq45tu+o5dnzwq5XHZKcGQmHOzHGUTszREYPI\nECkIJKada3KzY1c923fXsX1XfdjfNcrKdDFzegHTp+QxbUo+U8vHkpmh22+KDEZBIHGju9vP0Zom\n9lSdYU/VGfbuP0NrmzfMWjYmTshm2pQ8pk/JZ9rkPCaV5OhXU0WCKAgkbnV3+6k9eZ7d+06zp+os\n+w6coanZE3Y9l8tBWUkOZaW5lE/KZfKkXCaV5OJy2Ueh1SKxR0EgCcPv93OyroX9Bxs4cPgch6rP\ncfR40wX3U+hPSoqN4vFjmDzJCIdJpbmUFI9hbE6aBqMl4SkIJKF1dPg4fLSRA4fPcfBwAwerz11w\nb+bBZGW6mFg8hpLiMZQU55jPCghJLAoCSTqNTR6O1jRx5Fgj1ceaOFrTxMm6ViK/lUZQQEwYQ1Fh\nFsVF2RQVZlFUmIXLqS4miS8KAhGMaxmOHW/myLFGjhxroqb2PLUnz+P2DHxNQ/9sXFKQQVFhFhOK\nsplQlMX4omyKLslkXEGmQkJikoJAZAB+v5+Gc25qaps5fsIIhuEHBICNvLFpXDIuk8ICIxgKx2Uy\nriCDS8ZlUpCXgcOhs5lk9CkIRIYoOCBO1rVwqr6Vk3WtnKpr4fTZdobSxRQsJcVG/tgMLhmXQUF+\nBgV5GeTnpZM/1nzOS2dMdqrGJmTEKQhERpC3s4v6022cqmvhZH0rp+pbOHmqhfozbTSccw/4Q3uR\ncjhSyB+bTp4ZDgV56eTlZZA/Np38semMzU0jJydNXVAyJAoCkVHS6eui4Zyb+jNtnDnbxukz7Zw+\n08bphjbqT7fR2OQesW1lZrjIzUkjNyeVsbnpRnlMqjnPeIzNTWdMtksX14mCQCRWeL1dnGkwAqKh\nsZ2zDW7ONbppaGyn4ZybhnNu2trDXUk9VDZyzIDIyU4lO9vFmOxUsrNSGdNTzk5ljDmdnZWqC+8S\nkIJAJI64PZ2ca3Rz9pybc+eCQqLRmG5s9tB83hPRRXTDlZbmYEyWGRDZrj7lzAwXWZkuMjOcxnOm\nMzBPA+GxS0EgkmC6u/20tHbQ1Oyh6XwHjU1umpo7aGrymPOM58YmT9gf8RtJqamOviGR4SQz88Lg\nyMp0kZHuJDPDSXq6k4w0J+npDnVhRZGCQCSJ+XzdNJ83AqOlpYPzLR2cb/XS0tJB8/kOWlqNeS2t\nXlpavJxv7aCrq/+bCEWby+UgM90IhfR0JxnpDtLNkMjMcJKR7iI93UFGupP0tJ46RjkjwyinpTlw\nOlJ05lUIBYGIRMzv99Pu7qSlxUuzGRwtZnicb+mgvb2T1jYvrW1e2sxym1mOlb/LlBQbaakOUlMd\npKc5AuW0NDtpqQ7SgueZdVwue2/dNAfpZj2jjp30NCdOZ/wGjIJARKKuu9uP29NJa1tnIBiMsPAG\n5gWHh9vdSbvbh9vTSXu7UR7u9Rmjx4bLZSfVZSc11U6qy0Gqy44r6JEaeHYElXuWO3C5Ui5Ylprq\nuGB9h31kQ2c4QeAYsa2LSFJISbGRmWEMJjMuc8jrd3f76fD6AqHgdncaIeE25vWU3e5O2j2duM35\nwWW3pxO3xxfFri0/Xq8Pr9dHS2uUNmFKSbHhdNpxOe04nSl9yi6XHZejt+wMKbtcwesYz8OhIBCR\nUZWSYjPGA9Kc5F/ka3X6uujwdOHp8BkPjy9Q7ujowuPx4e7w0RE03+Px0eE1l3k6A+Xg9X2+0Rs7\n6e7209Hho6PDN2rbDKUgEJG45XTYcWbZycoa2VuYdnV10+HtoqOjC29nFx0dPrydXXi9xqMj8OwL\nTPfU9Xp9Qev1ffYG1e95tmrAPpiCQEQkhN2eQkZ6Chnpzqhvq6urG29nF52d3Xi9XXR2dtHpM8re\noHJnZ1egnlE2ngPr+YxwefF3Q2+DBotFRBLIcAaLI7mqYyVQBRwE7h+gzuPm8h3AggjWzQPeBA4A\nbwC5Q2m0iIiMnHBBYAeewPhAnw3cCswKqbMKmApMA+4A1kaw7gMYQTAdeNucTiqVlZVWNyGq9P7i\nm95fcgkXBIuBQ8BRoBN4DrgxpM4NwDNmeRPGt/uiMOsGr/MMcNMw2x+3Ev0/ot5ffNP7Sy7hgqAY\nOB40XWvOi6TOhEHWLQTqzXK9OS0iIhYIFwSRjtRGMjBhG+D1/EPYjoiIjLKlwGtB0w9y4YDxk8BX\ng6arML7hD7ZuFUb3EcB4c7o/h+gNCj300EMPPcI/DjHCHMBhoAxwAdvpf7D4FbO8FPgognUfpTcU\nHgAeHumGi4jIyLkO2I+RMg+a8+40Hz2eMJfvABaGWReM00ffQqePioiIiIhIfyK5iC2eHQV2AtuA\nzdY2ZUSwMhxyAAACY0lEQVT8GuPsr11B8xLposH+3t8ajDPhtpmPlaPfrBFRArwD7AF2Az805yfK\n/hvo/a0hMfZfGsZp+9uBvcBPzPlxv//sGF1JZYCT/scl4t0RjB2VKFZgXFEe/EH5KPD3Zvl+4nsc\nqL/39xDwI2uaM6KKgPlmOQujK3cWibP/Bnp/ibL/ADLMZwfGGO2nGeL+i8Ubh0ZyEVsiiPXfeRqK\njUBjyLxEumiwv/cHibEP6zC+bAG0AvswrvdJlP030PuDxNh/AO3mswvji3QjQ9x/sRgEkVzEFu/8\nGIPlW4DvWtyWaEmGiwZ/gHGCxDri8NC7H2UYRz6bSMz9V4bx/nrObEyU/ZeCEXb19HaDDWn/xWIQ\n+K1uwChYjvEf8jrgLoyuh0TWc35zIlkLlGN0O5wCfmptcy5aFvAn4B6gJWRZIuy/LOAFjPfXSmLt\nv26M9zERuAK4KmR52P0Xi0FwAmOAp0cJxlFBIjllPp8BXsToDks09fS9aPC0hW2JhtP0/oH9ivje\nh06MEPgd8JI5L5H2X8/7+z297y+R9l+PZmADcDlD3H+xGARbMH7JtAyjz+sWYL2VDRphGUC2Wc4E\nrqXvIGSiWA+sNsur6f0DTBTjg8pfIH73oQ2ja2Qv8LOg+Ymy/wZ6f4my/wro7dZKBz6LcRZUQuy/\ngS5ESwTlGP152zFOZ0uE9/cH4CTgxRjfuZ3Eumgw9P19C/gtxinAOzD+yOK1D/3TGF0L2+l7KmWi\n7L/+3t91JM7+mwtsxXh/O4G/M+cnyv4TEREREREREREREREREREREREREREREREREZGR9P8BnE+U\n9lmQAogAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 43 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Ok, let's start here. Suppose we know $\\lambda$. We can compute the distribution of interarrival times (times between logins)." ] }, { "cell_type": "code", "collapsed": false, "input": [ "lam = 0.1 # average arrival rate in logins per day\n", "interarrival = pmf = thinkbayes2.MakeExponentialPmf(lam, high=90)\n", "thinkplot.Pdf(interarrival)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEACAYAAAC+gnFaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt8VNW99/HPZCYXEpJAuCUkwYQ7KiBVEetttGqRttK7\ntRetp6f6qtW2nnNa9DznOU3b56mlz2lrreeg9dIXvam9aaEVL7WM2loRFAIYAsg1ARIgJJB7Jsk8\nf6w9MztDkpmQDHsm832/XuOsvffaMz9isn97r7X2XiAiIiIiIiIiIiIiIiIiIiIiIiIiItLHUqAa\n2A2sGKDOg9b2SmBRxDY3sBlYa1tXAdRa6zdb3yEiIgnIDbwLlAHpwBZgXkSdZcBzVvkS4I2I7f8C\n/ApYY1v3TWu9iIg4LC3K9sWYRLAf8ANPAcsj6twIrLbKG4BxwBRruQSTKB4DXBH7RS6LiIgDoiWC\nYqDGtlxrrYu1zo+ArwO9/Xz23ZimpMcxyUNERBwQLREEYvyc/s72PwgcxfQBRG5fBZQDFwBHgB/E\n+D0iIjLCPFG2HwJKbculmDP+weqUWOs+hmk2WgZkAXnAz4FbMAki6DH6diSHzJgxI7Bnz54oIYqI\niM0eYOZIfqDH+tAyIIPoncVLOL2zGOAq+h7si2zle4BfD/D9gUTzzW9+0+kQTqOYYpOIMQUCiRmX\nYopNIsZE7C05IdGuCLqBu4AXMCOIHgd2AHdY2x/BJIFlmE7lVuC2gQ7qtvJKTLNQANhn+7zTdwoE\ncLnUrywiEi/REgHAOutl90jE8l1RPuMV6xV0SwzfC0BDcycT87JirS4iIkMUrbPYcQePtjgdQh9e\nr9fpEE6jmGKTiDFBYsalmGKTiDGdiURvcwn84W/7+MhlZU7HISKSFKym9CEd23VFICKS4hI+ERyo\nVyIQEYmnhE8ENcda6Ont78ZkEREZCQmfCLq6e6lvbHc6DBGRUSvhEwHAfjUPiYjETVIkAnUYi4jE\nT1IkAnUYi4jET1IkAl0RiIjET1IkgsMNbXR29TgdhojIqJQUiaA3EOCArgpEROIiKRIBwP76ZqdD\nEBEZlZIoEeiKQEQkHpInEdTpikBEJB6SJxHUt2Am3xERkZEUSyJYClQDu4EVA9R50NpeCSyK2ObG\nTGBvn6qyAHgJ2AW8CIwb6MtzstIBaO3w03CqM4ZwRURkKKIlAjfwECYZnAvcTP9zFs8EZgG3A6si\ntn8VqKLvVJX3YhLBbOBla7lfZVPGhsr71DwkIjLioiWCxZi5iPcDfuApYHlEnRuB1VZ5A+bsfoq1\nXIJJFI/Rd6IE+z6rgQ8PFEBZYW6orEQgIjLyoiWCYqDGtlxrrYu1zo+ArwORz5GeAtRb5XrCieM0\n022JYK8SgYjIiIs2eX2svbOR06K5gA8CRzH9A94o3zHg96z91UNUb6oFoOf4xXDTwhhDEhEZ/Xw+\nHz6fb1ifEW1eyyVABaaPAOA+zNn9SludhwEfptkITMeyF/gK8DmgG8gC8oDfA7fY6tQBRcB6YG4/\n3x/o8vdw83f/ir/HXFT8coWX3OyM2P51IiIpJh5zFm/CdAKXARnATcCaiDprMAd3MImjCXOA/3eg\nFCgHPgX81VZvDXCrVb4VeHagANI9aZyjDmMRkbiJlgi6gbuAFzAjf54GdgB3WC+A54C9mE7lR4A7\nB/gse/PP94DrMMNHr7GWB1Ru6yfYc0SJQERkJA3p8sEBgUAgwLqNNTz8px0AXDG/kH/7+AKHwxIR\nSUzxaBpKCNOLbENIdUUgIjKikiIRlE3JxZ1mEtyhhlbaO7sdjkhEZPRIikSQme6meGIOAIGAHkkt\nIjKSkiIRAMwoyguV96p5SERkxCRNIrD3EygRiIiMnKRJBH2HkJ5yMBIRkdElaRKB/Yqg5lgLXd2a\nzF5EZCQkTSLIyUqnqCAbgO6eAAePtjockYjI6JA0iQD6dhjvOazmIRGRkZBUiaC8SHMTiIiMtKRK\nBDOK1GEsIjLSkioRzJwabhrad6RZHcYiIiMgqRJBbnYGhVaHsb+nl/11LQ5HJCKS/JIqEQDMLs4P\nlXcfOulgJCIio0PSJYJZxeHmod2H1E8gIjJcSZcI5pSErwh26YpARGTYYkkESzFzDO8GVgxQ50Fr\neyWwyFqXBWwAtmBmN7vfVr8CqMVMbL+Z8JzIUZUV5uJxW4+kPt5Kc7s/1l1FRKQf0RKBG3gIc6A+\nF7gZmBdRZxkwEzO38e3AKmt9B3A1cAGwwCpfZm0LAD/EJI1FwPOxBpyZ7qZsim0YqW4sExEZlmiJ\nYDFmLuL9gB94ClgeUedGYLVV3gCMA6ZYy23WewYmqTTa9jvjaTJn25uHatU8JCIyHNESQTFQY1uu\ntdZFq1Nild2YpqF6YD2miSjobkxT0uOY5BGzWRo5JCIyYjxRtgdi/JzIs/vgfj2YpqF84AXAC/gw\nzUfftup8B/gB8IX+PriioiJU9nq9eL3ePkNId9WeJBAIBCdsFhFJKT6fD5/PN6zPiHb0XILp2A12\n5t4H9AIrbXUexhzcn7KWq4GrMFcBdv8baAf+K2J9GbAWmN/P9wcCgdNzUW9vgM+uXE9rh5m7+NF7\nrmDyuDFR/ikiIqOfdVI8pDPjaE1DmzCdwGWYdv6bgDURddYAt1jlJUATJglMJNzkMwa4DjNCCKDI\ntv9HgG1DCjrNpeYhEZEREi0RdAN3YZp1qoCngR3AHdYL4DlgL6ZT+RHgTmt9EfBXTB/BBsxZ/8vW\ntpXAVkwfwVXAPUMN3H5jmTqMRUTOXLQ+AoB11svukYjlu/rZbxvwngE+85YB1sfMfkWwS3cYi4ic\nsaS7szjIngj2HD5FT2+vg9GIiCSvpE0EBbmZTMrPAqDT38OBej2JVETkTCRtIgCYWxq+/WDHwSYH\nIxERSV7JnQimKRGIiAxXUieCc22JoOpgE/3dcyAiIoNL6kRwzpSxjMkwA58aTnVw7GSHwxGJiCSf\npE4E7rQ05pSGRw+peUhEZOiSOhGAOoxFRIYr6RPBueowFhEZlqRPBLNL8nGnmecrHTjaTItmLBMR\nGZKkTwRjMj2UF5oZywIB2KnnDomIDEnSJwKAebbmoWo1D4mIDMmoSwRVBxsHqSkiIpFGRSKwjxza\nfegU/m49gE5EJFajIhFMyMuicLyZoazT38PeOj2WWkQkVqMiEQDMmzY+VK46oH4CEZFYxZIIlmLm\nId4NrBigzoPW9kpgkbUuCzMz2RbM7Gb32+oXAC8Bu4AXCU9pecbs/QTv7Fc/gYhIrKIlAjfwECYZ\nnAvcDMyLqLMMmImZ2/h2YJW1vgO4GrgAWGCVL7O23YtJBLMx01feO5x/BMD8ctsVwcFGTVQjIhKj\naIlgMWYu4v2AH3gKWB5R50ZgtVXegDm7n2Itt1nvGZik0tjPPquBDw899L6KCrKZkGcmqmnt6GbP\n4ebhfqSISEqIlgiKgRrbcq21LlqdEqvsxjQN1QPrMU1EYBJFvVWuJ5w4zpjL5epzVbB9/4nhfqSI\nSEqINnl9rA/4dw2wXw+maSgfeAHwAr5+6g74PRUVFaGy1+vF6/UOGMT8sgJ8lUcA2LrvBB+9vDx6\n5CIiSczn8+Hz+Yb1GdESwSGg1LZcijnjH6xOibXO7iTwZ+BCTCKoBwqBOqAIODpQAPZEEM388oJQ\necfBJrp7evG4R83AKBGR00SeIH/rW98a8mdEO0puwnQCl2Ha+W8C1kTUWQPcYpWXAE2YA/1EwqOB\nxgDXYZqJgvvcapVvBZ4dcuT9mDJ+DFPGmfsJOrp62H1Izx0SEYkmWiLoBu7CNOtUAU8DO4A7rBfA\nc8BeTKfyI8Cd1voi4K+Yg/8GYC1mhBDA9zCJYRdwjbU8IuxXBds1jFREJKrItv1EExjqPMTrKw/z\nwB+2A7BwegHfvvWieMQlIpKQXC4XDPHYPuoa0OeX9e0n6OrucTAaEZHEN+oSwcT8LKZOyAagq7uX\nnTXqJxARGcyoSwSgfgIRkaEYnYnA1jy0bZ9uLBMRGcyoTATn2+4w3lnbRHtnt4PRiIgktlGZCMaP\nzQzNY9zdE9BVgYjIIEZlIgBYNHNCqLx5T4ODkYiIJLZRnAgmhspvv6tEICIykFGbCOaVjiMrww1A\n3Yk2Dje0RdlDRCQ1jdpEkO5J6zN6aPO7xx2MRkQkcY3aRAB9+wm2qJ9ARKRfozwRhPsJtu47gb9b\n01eKiEQa1Ylg6oRsCgvM4yY6unrYUdPkcEQiIolnVCcCgPfYh5Gqn0BE5DSjPhEsmmFPBOonEBGJ\nFEsiWApUA7uBFQPUedDaXgksstaVYiasfwfYDnzFVr8CM+XlZuu1dIhxx+z88gI8bvNo7n11zZxo\n7ozXV4mIJKVoicANPIQ5UJ8L3AzMi6izDJiJmdLydmCVtd4P3AOch5nC8svAXGtbAPghJmksAp4f\nzj9iMNmZHuaWjgstv7VbzUMiInbREsFizBSU+zEH9qeA5RF1bgRWW+UNmHmKp2Ampg/OUdyCmeKy\n2LbfWZsdbfGcSaHyxp1Hz9bXiogkhWiJoBiosS3X0vdgPlCdkog6ZZgz/w22dXdjmpIeJzzJfVxc\nbEsEW/acoLNLs5aJiARFSwSxThgceXZv328s8Dvgq5grAzDNR+XABcAR4Acxfs8ZmTohh9JJOQB0\n+nuo3KtOYxGRIE+U7Ycwnb5BpZgz/sHqlFjrANKB3wO/BJ611bG3zzwGrB0ogIqKilDZ6/Xi9Xqj\nhNy/i+dMouZYKwBv7jzG4rmTz+hzREQSic/nw+fzDeszorXTe4CdwPuAw8CbmA7jHbY6y4C7rPcl\nwAPWuwvTd9CA6TS2K8JcCWBtuxj4dD/fHwgEYr0oGdyOg43c+/hGwMxX8MS/Xkla2lnrphAROStc\nLhcMsQ822hVBN+Yg/wJmBNHjmCRwh7X9EeA5TBJ4F2gFbrO2XQZ8FtiKGSIKcB9mhNBKTLNQANhn\n+7y4mV2ST35OBidbu2hs6WT3oZPMKY1r14SISFJI9FPiEbsiAPjJs+/wl82m1erjV5TzuWtnjdhn\ni4gkgjO5Ihj1dxbbLZ4bHj305s5jDkYiIpI4UioRLJxeQIbH/JMPHm3hyAlNViMiklKJICvDw8Lp\n4WcPbdRVgYhIaiUC6Htz2T+q6h2MREQkMaRcIrhk7iTSTGcKVQebOH6qw+GIRESclXKJYNzYTOaX\njw8t/6NKzx4SkdSWcokA4LLzCkPl19+pczASERHnpWQiWDJvspqHREQsKZkI8nMy+jQPvf6OOo1F\nJHWlZCKAvs1DGj0kIqksZRPBknmTcafZmodOqnlIRFJTyiYC0zxUEFp+XVcFIpKiUjYRALz33Cmh\nshKBiKSqlE4El9qah3YcbKK+sd3hiEREzr6UTgR5ORlcMCP87KFXth4ZpLaIyOiU0okA4OqFU0Pl\n9ZWHGcn5D0REkkEsiWApUA3sBlYMUOdBa3slsMhaVwqsB94BtgNfsdUvAF4CdgEvAo5NFbZ47iRy\nssxEbYcb2thVe9KpUEREHBEtEbiBhzDJ4FzMfMXzIuosA2YCs4DbgVXWej9mPuLzMHMYfxmYa227\nF5MIZgMvW8uOyEx3c6mt09in5iERSTHREsFizFzE+zEH9qeA5RF1bsRMUg+wAXN2PwWoA7ZY61sw\ncx0X97PPauDDZxT9CPEuKAqVX9teh7+718FoRETOrmiJoBiosS3XEj6YD1anJKJOGabJaIO1PAUI\njtest5Ydc94545mUnwVAc5uft3YfdzIcEZGzKloiiLXnNHKiZPt+Y4HfAV/FXBn09x2O9tCmpblO\n6zQWEUkVnijbD2E6fYNKMWf8g9UpsdYBpAO/B34JPGurUw8UYpqPioABJwWoqKgIlb1eL16vN0rI\nZ8a7sIjfvLoXgE27jtHc1kVudkZcvktEZKT4fD58Pt+wPiPyTD6SB9gJvA84DLyJ6TDeYauzDLjL\nel8CPGC9uzDt/w2YTmO771vrV2I6isfRf4dx4GwO5/z6oxtCo4b++Ya5fGjJtLP23SIiI8FlHrEf\n7djeR7SmoW7MQf4FoAp4GpME7rBeAM8BezGdyo8Ad1rrLwM+C1wNbLZeS61t3wOuwwwfvcZadtz7\nLgg3D730dq3uKRCRlDCkrOGAs3pF0NbZzW3/9QodXT0AfO8LFzNv2vgoe4mIJI54XBGklOxMD1ec\nH56n4MW3Dg1SW0RkdFAiiHD9heGRr3/bXkdzu9/BaERE4k+JIMKs4jxmFOUB0NXdqwfRiciop0QQ\nweVycf2F4XvmXtikTmMRGd2UCPpx5YIisjLcABw82kJ1TZPDEYmIxI8SQT8iO42f3xR5D52IyOih\nRDCApReHb5b+2/Y6Gls6HYxGRCR+lAgGMHNqHvOmmWkSunsCrHuzJsoeIiLJSYlgEPZHTKzbWEun\nv8fBaERE4kOJYBBL5k0OPZ76VFsXr22rczgiEZGRp0QwCHdaGh+4JHxVsOaNAxpKKiKjjhJBFNe+\npzg0lPRAfQtb955wOCIRkZGlRBBF7ph0rrE9lXTtGwcdjEZEZOQpEcTgg7bmoY27jnHwaH8TrYmI\nJCclghgUT8zhkrmTQsu//9s+B6MRERlZSgQx+tjl5aHya9vqqDvR5mA0IiIjJ5ZEsBSoBnYDKwao\n86C1vRJYZFv/BGZ+4m0R9Sswcx9HzlyWsOaUjmNBeQEAPb0Bnvn7fmcDEhEZIdESgRt4CHOgPhcz\nX/G8iDrLgJnALOB2YJVt28/o/yAfAH6ISRqLgOeHGrgTPnFl+Krg5c2HOX6qw8FoRERGRrREsBgz\nF/F+wA88BSyPqHMjZpJ6gA2YieiDT2x7DWgc4LMTfZrM08wvL2BOST4A/p5e1rx+wOGIRESGL1oi\nKAbsD9mptdYNtU5/7sY0JT2OSR4Jz+Vy8Ykrp4eWn99Uy6nWLgcjEhEZPk+U7bHeRht5dh9tv1XA\nt63yd4AfAF/or2JFRUWo7PV68Xq9MYYUHxfNnkh5YS776prp9PfwzOv7ufW62Y7GJCKpy+fz4fP5\nhvUZ0ZpnlmA6doPt/PcBvcBKW52HAR+m2QhMx/JVmE5igDJgLTB/gO8YbHsgER/p8Pd36vn+byoB\nyEx3s+orlzEhL8vhqERETMsFQ2x6j9Y0tAnTCVwGZAA3AWsi6qwBbrHKS4AmwklgIEW28kc4fVRR\nQrt03uTQvMad/h5++6ruKxCR5BUtEXQDdwEvAFXA08AO4A7rBfAcsBfTqfwIcKdt/yeB14HZmH6E\n26z1K4GtmD6Cq4B7hvnvOKvS0lx85n0zQ8svvV2r+wpEJGkl+sidhGwaAggEAvyvn23inQNmUNTV\nC6fytY+e73BUIpLq4tE0JANwuVx87trwVYFv62EO1OsZRCKSfJQIhmHetPFcNHsiAIEA/PLl3Q5H\nJCIydEoEw/SZa8JXBW/uPEblngYHoxERGTolgmGaXpTH+2zzFTz2/E56ensdjEhEZGiUCEbAZ6+d\nxZgMc2/ewaMtPL+x1uGIRERip0QwAgpyM/n4FWWh5Sd9e2hu06MnRCQ5KBGMkA9deg6F48cA0Nzm\n5ynfXocjEhGJjRLBCMlMd3Pr9eFnDq3bWMO+umYHIxIRiY0SwQi6dN5k5tsmr1m1tore3sS8IU5E\nJEiJYAS5XC7u+MBcPG5zU9/O2pM8v0kdxyKS2JQIRljppLF87IrwTGa/+MtuzWQmIglNiSAOPnZ5\nOcUTcwBo6+zmsXU7HY5IRGRgSgRxkJnu5ksfDE/t/I+qet7YcdTBiEREBqZEECfzywu4dlF4xs5V\na6s4qWktRSQBKRHE0eevn8X43EwAmlq7+J+1VSTqY7VFJHUpEcRRbnYGdy8/L7T8xo6j+CqPOBiR\niMjpYkkESzHzEO8GVgxQ50FreyWwyLb+Ccy0lZFTURYALwG7gBeBcbGHnFwunDWRGy4uDS0/uq6a\nY03tDkYkItJXtETgBh7CJINzgZuBeRF1lgEzMXMb3w6ssm37GeGJ7+3uxSSC2cDL1vKo9fnrZ1FU\nkA1Aa0c3DzyzXU8oFZGEES0RLMbMRbwf8ANPAcsj6twIrLbKGzBn94XW8mtAYz+fa99nNfDhoQSd\nbLIyPHztI+eTZqaQY/v+Rp5ar2cRiUhiiJYIijGTzgfVWuuGWifSFEyTEdb7lCj1k97caeO4yTs9\ntPzb1/ay+d3jDkYkImJ4omyPdYhL5ETJQxkaExisfkVFRajs9Xrxer1D+OjE8okry9lxsIktexoI\nBOCHv9/Oj760hIl5WU6HJiJJyufz4fP5hvUZ0Wa6XwJUEG7nvw/oBVba6jwM+DDNRmA6lq8ifMZf\nBqwF5tv2qQa8QB1QBKwH5vbz/YHRNtyysaWTex5+g8bmTgDOnTaO73z+IjxuDeASkeFzmSboaMf2\nPqIdfTZhOoHLgAzgJmBNRJ01wC1WeQnQRDgJDGQNcKtVvhV4NrZwk9/4sZn868fmh/oLqg428fjz\negSFiDgnWiLoBu4CXgCqgKeBHcAd1gvgOWAvplP5EeBO2/5PAq9jRgfVALdZ678HXIcZPnqNtZwy\n5pcX8JlrZoSWn3uzhnUbawbZQ0QkfoZ0+eCAUdc0FBQIBPjB77bx2vY6ANxpLipuuZAF1nwGIiJn\nIh5NQxInLpeLu5efx8ypeYCZyOb7v6nkcEObw5GJSKpRInBQZoab+26+IPQ8ouY2P9/65ds0tnQ6\nHJmIpBIlAodNzMvi3z91ARke87+i7kQb3/nlZto6ux2OTERShRJBAphdks/XP7EgNJJoz5FTrHy6\nEn+3HkMhIvGnRJAgFs+dzJc+FH6M05Y9DfxYzyQSkbNAiSCBXH9hCZ+5ZmZo+bXtdfz4mXeUDEQk\nrpQIEswnriznA5dMCy2/svUIP3m2SslAROJGiSDBuFwuvnjDHJZeVBJat77yMP/9RyUDEYkP3VCW\noHp7A6z60w5efKs2tO7y8wr52kfPJ92j/C0i/TuTG8qUCBJYb2+A/15bxV/ePhRat2jGBFbctJAx\nmdEeHCsiqUiJYBTq7Q3w+PM7+dOGg6F1s0vy+c/PLCI3O8PByEQkESkRjFKBQICnX9nLk+v3hNZN\nnZDNf3x6EcUTcxyMTEQSjRLBKPfnDQf56XPVoeWxY9JZ8ckFLJg+wcGoRCSRKBGkgNer6nngD9vp\n9PcA4HG7uH3ZPN5vG2UkIqlLiSBF7D50kv/75JbQLGcA1y4q5vYPzCUz3e1gZCLiNCWCFHL8ZAff\nfXILe46cCq0rL8zlG59cyNQJ2Q5GJiJOitd8BEsxcwzvBlYMUOdBa3slsCiGfSuAWmCz9VqKDMnE\n/Czu/6eLuXrh1NC6fXXN/NtP3+DVbUccjExEkk20rOEGdgLXAoeAjcDNmOkqg5ZhprNcBlwC/Bgz\nd/Fg+34TaAZ+GOX7dUUQRSAQ4MW3DvHoc9X4e8J3Hl85v4g7PjCXsWPSHYxORM62eFwRLMbMRbwf\n8ANPAcsj6twIrLbKG4BxQGEM+yZ6s1RScLlcvP+iEu7/wsUUjh8TWv/qtiN8bdU/qNzT4GB0IpIM\noiWCYsyk80G11rpY6kyNsu/dmKakxzHJQ4ZhVnE+P/rSpVy7KPwjPnayg//8+Vv85I/v0NzudzA6\nEUlk0Z5TEGu7zFDP7lcB37bK3wF+AHyhv4oVFRWhstfrxev1DvGrUkd2poe7P3we75k1kVV/qqK5\nzRz8//L2Id7afZwv3jCH9547JXjpKCKjgM/nw+fzDeszoh0RlmA6doOdufcBvcBKW52HAR+m6QdM\n5/BVQHkM+wKUAWuB+f18v/oIzlBjSyc//XM1r1fV91m/oLyAf75hLudMGetQZCIST/EYPurBdPi+\nDzgMvMngncVLgAes98H2LQKCQ1vuAS4GPt3P9ysRDNMbO47y8J939LnnwJ3m4oaLS7npqunk5eh5\nRSKjSbzuI7gBc3B3Y9rz7wfusLY9Yr0/hDnzbwVuA94eZF+AnwMXYJqe9lmf1/fU1VAiGAEt7X6e\nXL+HdRtr6OkN/zyzMz3ceOk5LH/vOWTraaYio4JuKJNBHahv4dF11Wzbd6LP+rzsDD52eRk3LC7V\nnckiSU6JQKIKBAK8UX2UX7+8h4PHWvpsK8jNZPml53D9RSW6QhBJUkoEErOe3l5e3VrHk+v3UN/U\n3mdbTpaH919UwgcvmcaEvCyHIhSRM6FEIEPm7+7lxbdq+e1r+/p0KIN5sql3wVRuWFzKzKl5DkUo\nIkOhRCBnrNPfg6/yCH/8xwEOHW89bfuMojzef1EJV8wvVLORSAJTIpBh6+0N8ObOYzzz9/1U1zSd\ntn1MhocrFxRy1YIi5pWOIy0t0X+FRFKLEoGMqB0HG3lh0yH+9k4d/u7e07ZPys/i8vMLuXJ+IeWF\nubpjWSQBKBFIXDS3dbG+8ggvvlVLzbHTm40ASibmsGTeZBbPmcSs4nxdKYg4RIlA4ioQCLDjYBOv\nbD3C36vqQ88yijR+bCYXzp7I4jmTmF9eoD4FkbNIiUDOGn93L5V7G3h1Wx0bqo/S0dXTbz2P28Xs\n4nwWTJ/ABTMKmFWcj8cdy3xIInImlAjEEZ1dPWzZ28CbO4+xcecxTrZ2DVh3TIaHOaX5zCnNZ27p\nOOaU5JOTpclzREaKEoE4rqe3l921p9iw8yhb9jSw90jzoPVdLpg2eSzzSscxuySf6UV5lE7K0VWD\nyBlSIpCEc7K1i617T1C5t4HKvSc4GnEXc3/S3WlMmzyW8qJcphfmMr0ol3Om5KqvQSQGSgSS0AKB\nAMdOdlBd00R1zUmqDzaxr66Z3hj/HxfkZlIyMYeSSTnmfWIOxZNymJCbqaGrIhYlAkk67Z3d7D50\nkh0Hm9hz5BT76lpiumqwy/CkMXncGKaMN69gudAqjx2jPghJHUoEMio0t3Wxr66ZfXXN7LXeDx1v\npbvnzH4XsjM9FORmUpCbyYS8LMbnZjAhL8ssW+vzcjL0CG4ZFeKVCJYSnlzmMU6fahLgQcwkNG3A\n54HNUfYtAJ4GzgH2A58ETn+egRKBWHp6e6lvbKf2WCu1x1s5dLyNmmMtHGpoo6W9//sZhiorw01e\ndgZ52enbE9bZAAAHG0lEQVTk5WSQHyxby3nZGeRkecjJ8pCdab1neXCnqWNbEkc8EoEbM93ktcAh\nYCODT1V5CfBjzFSVg+37feC49b4CGA/c28/3J1wi8Pl8eL1ep8PoI9Vjau3wU9/YTn1TO0cb26lr\nbKe+sZ2jTe0cbeqg02/ucTi+v5KJZQtH/PvHZJiEYE8S2VkecjI9ZKa7ycpwk5nhJivdHVrOyjDl\nzHQ3mze+ztVXe822dDfpnjTH78xO9d+pWCViTGeSCKINw1gMvIs5awczQf1y+iaCG4HVVnkDMA4o\nxExeP9C+N2ImuMfa10f/iSDhJOL/+FSPKScrnelF6UwvOv1R2YFAgJZ2Pw3NnXz3//yFTy4/j4ZT\nHZxo7qSxpZOGU52caO7kVFvXGTc9tXd1097VTcOpM4u/2vdr5m7p2yzlcbtId5ukEHxluM27x+0K\nr/OkheqF1rvTcLtdeNJMQvG4XbjTrJc7LVxOc5ntacH6Ztmd5uLpZ9YxqXxhn3XutDTS0iDN5cLl\nCr67SHOZg09aWnh9mgtwmf1cmHW4zHzZ9n2HItV/z+MpWiIoBmpsy7WYs/5odYqBqYPsO4XwHMX1\n1rLIiHO5XORmZ5CbncHUCdlc+57ifusFAgHau3o42drFqbYuTrX6OdXWZS37aW4zyy0d3bR1dNPW\n2U1rh5+2zm7icdHa3ROgu6eb9oHvzYur6i2HOfzExrh/TzAxmIRivdIAzLvLOrF1uaDy7/vZ//9e\nCZ3rBtOIy0pMfZZD5eB/zPbgd/Xd3/oeV6iqtT0YWzgOwh+Hy+Xi9bdqaX10w4CxnMb2BZF1BsqL\nrgH26RPrMEfNRUsEsf6KxxKFa4DPCwzhe0TiwuVymSadTA9FBdkx79fbG6C9yySH1s5uWq1E0dLh\np72zh46ubjr9vXT4e+js6qGjq4dOfw8dfqvc1cPh7HQm5meFlv09pz/pdbTq6Q3+6Uc/BHR09dDY\n0hm13tnUcKqTXbUnnQ4j7pYAz9uW78O06ds9DHzKtlyNOcMfbN9qTPMRQJG13J93CScKvfTSSy+9\nor/eZYR5gD1AGZABbAHmRdRZBjxnlZcAb8Swb7CTGEzfwPdGOnARERk5N2BG/7yLOasHuMN6BT1k\nba8E3hNlXzDDR/8C7AJexHQwi4iIiIiIGEsx/Qa7Ob1P4mx5AjOiaZttXQHwEs5dyZQC64F3gO3A\nVxIkrizM0OEtQBVwf4LEBeZ+ls3A2gSJaT+w1YrpzQSJaRzwO8zQ7irM6D4nY5qD+fkEXycxv+tO\n/5zuw/ztbQN+DWQmQEwAX7Vi2m6VSZC4hsWNaUoqA9Lpv1/ibLgCWETfRPB94BtWeQVnv2+jELjA\nKo/FNLvNS4C4AIJDbTyYfqLLEySufwF+Bayxlp2OaR/mj9TO6ZhWA/9klT1AfgLEFJQGHMGcBDkZ\nUxmwF3PwB/NkhFsdjgngfMwxKgtz7HwJmJEAcQ3bpfQdbXQvzt1sVkbfRBAcEQXmoDzQaKez5VnM\nnduJFFc25i7y83A+rhJMX9TVhK8InI5pHzAhYp2TMeVjDnCRnP45BV0PvGaVnYypAHPiNR6TLNcC\n1zkcE8DHMY/vCfoPTAJwOq5h+zjwqG35s8BPHIqljL6JoNFWdkUsn21lwAEgl8SIKw1z9daMORsB\n5+P6Leaq7irCicDpmPZimjs2AV9MgJguwDTr/Qx4G/O3l+NwTHZPAHdaZadjuh3z+30U+EWCxDQX\nk6AKMCdhr2Oe/TakuBLxaVkBpwOIUXDMrhPGAr/HtAdGTgHmVFy9mINKCXAl5izc7mzH9UHMH+xm\nBr7h0Ymf1WWY5HQD8GVME6STMXkwI/3+x3pv5fQrcKd+pzKAD2ESeqSzHdMM4GuYE7CpmL/Bzzoc\nE5gz/ZWYfoB1mJOxyAnEo8aViIngEKY9MKgU83iKRFBP3xvhjjoQQzomCfwC0zSUKHEFnQT+DFyI\ns3G9F/NMq33Ak8A1mJ+Z0z+rI9b7MeAZzPO8nIyp1noFnyfxO0xCqHMwpqAbgLcwPytw9ud0EeZs\nuwHoBv6AacZOhJ/TE1Z8V2HO/HcxxJ9VIiaCTcAswjei3US4o89pazAdRFjvzw5SNx5cwOOYkR0P\nJFBcEwmPShiDaTvd7HBc/445iSjH3Pn+V+BzDseUjWnKA9P8cj2m6dHJmOowzwSbbS1fixkZs9bB\nmIJuxiTxICd/TtWYG2bHYP4Or8X8HSbCz2my9T4N+ChmRJPTx4QRMdCNaGfTk8BhoAvzh3Ibzt8I\ndzmmCWYL4aF1SxMgrvmY9uUtmKGRX7fWOx1X0FWETyacjKkc8zPaghnqF/zddvrntBBzRVCJOdPN\nT4CYcjCPqs+1rXM6pm8QHj66GnN17nRMAK9acW0h3CSbCHGJiIiIiIiIiIiIiIiIiIiIiIiIiIiI\niIiIiEii+f/zR8D3x3Q8cAAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 44 }, { "cell_type": "markdown", "metadata": {}, "source": [ "If we observe someone, we are more likely to land during a longer interval." ] }, { "cell_type": "code", "collapsed": false, "input": [ "observed = interarrival.Copy()\n", "for val, prob in observed.Items():\n", " observed[val] *= val\n", "observed.Normalize()\n", "\n", "print(interarrival.Mean(), observed.Mean())\n", "thinkplot.Pdf(observed)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "9.76490343621 19.9050636319\n" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEACAYAAAC+gnFaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt8XHWd//HXXHJPmqRJ2jSXNr3TtFCKNBQQCAhaqsJe\nXKGrsuKusItV3HUVdNe1xf3pqrv8XBaX5bcKoisXZZGtCgLKDiBIoVBK77f0lrRNkzZp7pfJzO+P\n78nMyTTJTNIkZy7v5+Mxj5xz5nuSD6GZz/neQUREREREREREREREREREREREREREZIjVwG5gH3DX\nCGXus97fCqywXX8IaAS2RZSvAd4AtgBvAisnMF4REZlAHmA/UAWkAe8ASyLKrAGesY4vAV63vXcF\nJjFEJgIf8AHr+HrgfycqYBERGRt3lPdrMIngENAPPA7cGFHmBuAR63gTUACUWuevAC3DfN/jQL51\nXAA0jCVoERGZON4o75cDR23n9Zin/mhlyoETo3zfu4HfAf+MSUaXxhKsiIhMvGg1gmCM38c1xvt+\nAHwOmA38NaYvQUREHBCtRtAAVNrOKzFP/KOVqSB6U08NcK11/CTw/eEKzZ8/P3jgwIEo30pERCIc\nABbEWjhajWAzsBDTWZwO3ARsjCizEbjFOl4FtGJGCo1mP3CVdXwNsHe4QgcOHCAYDMbV62tf+5rj\nMSim5IpLMSmmiX4B86N8Bg8RrUbgB9YBz2FGEP0A2AXcbr3/IGbE0Brrw70TuNV2/2PWB34Rph/h\nH4CHgduA7wEZQLd1LiIiDoiWCACetV52D0acrxvh3rUjXN/M2Z3OIiLigGhNQxKhtrbW6RDOophi\nF49xKabYKKbJEznaJ94ErfYuERGJkcvlgjF8vqtGICKS4pQIRERSnBKBiEiKUyIQEUlxSgQiIilO\niUBEJMUpEYiIpDglAhGRFBfLEhMpLxAM0tTVx/H2XtI9boqz0ynJThuctCEiktCUCEbR2tPPM/ub\neXZ/E6d7+oe8V56XyRWVhdyweAb5Gfo1ikjiivdHWseWmNjU0Mp33zhMW69/1HLZXg8fXVrKHy2e\niccd779OEUkFY11iIt4/uRxJBE/sOM6Pth0bci03zcPs/Cz8gSBH2rrp8QeGvL+sJI+/vbSKkuz0\nqQxVROQsSgTn6Ok9jfznlvAmbCXZ6dy6vJzLKwvxWk/8vf4Amxpa+cn249S394TKTs9M457aBcwt\nyJ7SmEVE7JQIzsGrR1v4xqt1ofPlM/O4+7J5TBuhD2AgEOSnO0/w6I7jBKw4c9I8bLhqAUuKc6ck\nZhGRSJOx+uhqYDewD7hrhDL3We9vBVbYrj+E2bZy2zD3fBaz29l24FsxxjtpTnX18W9vHgmdLy3O\n5R+umD9iEgDwuF2sXTaLf6xdQE6aB4DO/gG+9tJ+DrR0TXrMIiITIVoi8AD3Y5JBNWbHsSURZdZg\nNkleiNly8gHbew9b90a6GrgBuABYBvzzWAOfSMFgkH994zDtfaZjuCQ7na9eMZ9Mryem+5fPnMY3\nr1kUGj3U2T/AV337ONHRO2kxi4hMlGiJoAazF/EhoB94HLgxoswNwCPW8SagACi1zl8BWob5vn8F\nfNP6ngBNYwl6or1a38pbJ9oAcOHib1ZVkTfGIaHzC7P5eu3CUM3gTK+fe145QFf/wITHKyIykaIl\ngnLMpvOD6q1rYy0TaSFwJfA64AMujhboZOkfCPDDrQ2h8w8vKuGCGXnj+l7zC7NZf+UC0qxO5cNn\nuvnO7w+G+g9EROJRtMfeWD/BIjslot3nBQqBVcBK4KfAvOEKrl+/PnRcW1s74XuEPnugmeNWE05u\nmoe1S2ed0/erLsnlczVz+JfXDwHwxrEzPLW7kY8sKR39RhGRcfL5fPh8vnHfHy0RNACVtvNKzBP/\naGUqrGujqQeeso7fBAJAEXAqsqA9EUy0voEAP915InR+09JZo3YOx+qaqiLqWrr5+Z5GAH707jGq\ni3OpLtFIIhGZeJEPyRs2bBjT/dGahjZjmnGqgHTgJmBjRJmNwC3W8SqgFTNSaDRPA9dYx4us731W\nEphsvsOnabGWjijKSuNDC0sm7Ht/cnl5aAjpQDDIvZsO0eNXf4GIxJ9oicAPrAOeA3YCT2CGfN5u\nvQCeAeowncoPAnfY7n8MeA3zYX8UuNW6/hCmKWibVeYWplgwGOR/9pwMnX940QzSPRO3GKvX7eJL\nl84NdR4f7+jlh1uPRblLRGTqpeyEsi0n2vh73z4AMr1ufvjh88c8UigWvzl4iv+76VDo/BtXL2L5\nzPF1RouIxGIyJpQlpWcPNIeOr5tbPClJAOB9VdOpKcsPnf/rG4fp1pBSEYkjKZkI2nr9vNHQGjpf\ns6B40n6Wy+Vi3cWzybWaiBo7e3loa7S+dBGRqZOSieDlIy30B0yT0+KiHGbnZ03qzyvKTuf294QH\nVj2zv4ndzR2T+jNFRGKVkongxUPhAUpXz5k+JT/z6jlDm4geeOsoAwFNNBMR56VcIqhv62HPqU4A\nvC4XV01RInC5XNx+UWVoZNL+li6eq2uOcpeIyORLuUTwWn24b2BlWf6ETCCLVWluBh9ZMjN0/qN3\nj3Emyg5oIiKTLeUSweu2TuLLKgum/Od/5LxSSnMyAGjv8/OIOo5FxGEplQiau/pCzUIel4uVs/Kj\n3DHxMrxubruoInT+fN2pUEwiIk5IqUSwqeFM6Pj8GXmTNncgmkvKC0Idx0GC/OCdepzYm1lEBFIs\nEdibhVaVT31twO4vVlTgNbP/2NHUweu2JCUiMpVSJhH0+AfYdrI9dH5J+dT3D9iV52WyxrbI3Q+3\nNuDXcFIRcUDKJILtTR2hSWRz8rOYkZPucERw89JZoUXp6tt7eHa/oxu1iUiKSplE8M6JcG1gRZws\n+paf4eWj1eENax7dcZzOPq1DJCJTK3USQWNb6Hh56TQHIxnqwwtnUJJtaidtvX5+tutElDtERCZW\nSiSClu5+DrZ2A2Y28flxtFNYhtfNJy8Ib/G8ce9JWrr7HYxIRFJNSiSCrY3hZqHzinPJstrl48WV\ncwqZX5gNQO9AgCd2qlYgIlMnlkSwGtgN7APuGqHMfdb7W4EVtusPYbat3DbCfV/A7Fc8qQv+bLWN\nFrowTvoH7NwuFx9bNit0/usDTZzs7HMwIhFJJdESgQe4H5MMqoG1wJKIMmuABZi9jW8DHrC997B1\n73AqgeuAw2MLeex2NoWXfD5/Rvw0C9nVlOVzXlEOAP2BII/tOO5wRCKSKqIlghrMXsSHgH7gceDG\niDI3AI9Yx5uAAmBwKMwrQMsI3/te4EtjC3fsWnv6qW/vASDN7WLh9JzJ/pHj4nK5uOWCstD5bw+e\nosGKW0RkMkVLBOWYTecH1VvXxlom0o1WuXdjiPGc7G4Or+OzYHoOGd747RZZPnNaaD/jgWCQn2xT\nrUBEJl+0xXZineoauUnyaPdlA1/BNAuNdH/I+vXrQ8e1tbXU1tbGGJKxw7YTWHVxfNYG7G45v4wv\nNO4BzE5qH60upapgcndQE5HE5vP58Pl8474/WiJowLTlD6rEPMmPVqbCujaS+UAVpmN5sPxbmGao\nk5GF7YlgPHbZagTVxfHZP2B3XnEuNWX5vHHsDEGC/HjbMb56xXynwxKROBb5kLxhw4Yx3R+tnWQz\nphO4CkgHbgI2RpTZCNxiHa8CWjEjhUayDZgJzLVe9cBFDJMEzlWvP8A+2xLP5yVAjQDg4+eH+wpe\nb2jVMtUiMqmiJQI/sA54DtgJPAHsAm63XgDPAHWYTuUHgTts9z8GvAYswvQj3DrMz5i0ldb2t3Th\nt5Z3rsjLpCAzbbJ+1ISaX5jNFZWFofP/2nbMwWhEJNnFsiD/s9bL7sGI83Uj3Ls2hu8/L4Yy47Lv\ndOLVBgZ9/PwyXq1vJRAM8vaJNnY3d3BeAjRtiUjiid8hNBNg3+mu0HG8DhsdScW0TGrnhGsFj+7Q\nbGMRmRxJnQj2D0kE2Q5GMj43Vc/CbW1e89bxM+orEJFJkbSJoLNvIDSRzOtyUZWfeEMwK6ZlcuVs\nW61gu+YViMjES9pEsL8lXBuoKsiK64lko7mpuhSXNc1is2oFIjIJEvPTMQb2juIFCdgsNGh2ftaQ\nWsHjWoNIRCZY0iYCe//AgsLETQQANy0N1wreOHZmSJITETlXyZsIWhJ3xFCkOflZvLeyIHT+uEYQ\nicgESspE0Nk3wPGOXsB0FM/Jz3Q4onN389LwfgWvN7RywJboRETORVImgsNnukPHlfmZpHkS/z+z\nqiCLy22zjR/TCCIRmSCJ/wk5jEO2RDA3iVbu/FNbreD3Da3UqVYgIhMgKRPB4Eb1YNrXk0VVQRaX\nVYT7CrSLmYhMhKRMBIdtiSDZ1vJfa6sVvFbfysFW1QpE5NwkXSIIBoNDm4aSqEYAMK8wm0vLNYJI\nRCZO0iWCk119dPYPAJCX7mV6VmIsPT0Wa5eFawWvHm3lkK0GJCIyVkmXCA5FNAu5XCPugpmw5hdm\ns8qqFQQJ8tOdqhWIyPgldyJIsmYhu5uqS0PHLx9p4Whbj4PRiEgiizURrAZ2A/uAu0Yoc5/1/lZg\nhe36Q5itK7dFlP8OZrezrcBTQH6MsYzq8JnwB2KydRTbLSrK4eJZ5lcWJMgTGkEkIuMUSyLwAPdj\nkkE1ZtexJRFl1gALMPsb3wY8YHvvYeveSM8DS4HlwF7gy2MJfCT2J+NkmFE8mj+19RW8dKSFetUK\nRGQcYkkENZj9iA8B/cDjwI0RZW4AHrGONwEFwGDbxStAyzDf9wUgYLunItagRzIQCIb2IAAoz0vu\nRLC4KIf3lE4DIBAM8oT6CkRkHGJJBOWYjecH1VvXxlpmNJ8CnhlD+WE1dfXRN2ByS0FGGtMyYtmS\nObHZRxC9dPg0x9pVKxCRsYnlkzIY4/eKHJ4T631/B/QBjw735vr160PHtbW11NbWjviN7M1ClUne\nLDRoSXEuK2ZOY0tjGwPBID/b1cidNXOcDktEppDP58Pn8437/lgSQQNQaTuvxDzxj1amwroWzScx\n/QvvG6mAPRFEY28jr0jyZiG7m5eWsqWxDYAXD57ipupSSnMzHI5KRKZK5EPyhg0bxnR/LE1DmzGd\nwFVAOnATsDGizEbgFut4FdCKGSk0mtXAFzH9DRPSnmHvH6icljqJYNmMPC6YkQeAPxjkZ7vUVyAi\nsYslEfiBdcBzwE7gCcywz9utF5j2/TpMp/KDwB22+x8DXgMWYfoRbrWu/xuQi+k03gL8+zn8dwBw\n1DZ0tGJaaj0R2/sKflN3isbOXgejEZFEEu/TboPBYKxdDbD251tp6/UD8NCHlzEzJ3WSQTAY5K4X\n97KjqQOA6+eXsG7lbIejEhEnWCsqxPz5njQzi8/0+kNJIMPjpiQ73eGIppbL5RqyX8ELdc00dfU5\nGJGIJIqkSQT2juLKaZm4k3CNoWiWz8yjujgXMH0FT6qvQERikDSJoKEtdSaSjcTlcnHz0vAaRM8d\naKZZtQIRiSJpEsGxjnDnaFle6vQNRLqodBqLi3IA6A8EeWp3tMFbIpLqlAiSjMvlGrKL2bMHmjnd\n3e9gRCIS75ImEZywJYJZKT6Z6uJZ01g43dQK+gYCqhWIyKiSIhEEg0GOtSsRDDK1gnBfwTP7m2hR\nrUBERpAUiaC1x0+332xPmZPmIT8FFpuLpqYsn/mF2QD0DgT4+R7VCkRkeEmRCI5HNAsl4/aUY+Vy\nubjZtovZL/c10dqjWoGInC0pEsEx9Q8Ma1VFAXOtXdp6BwI8veekwxGJSDxKikRwvF0jhobjdrm4\n2TaC6Jf7mkKzr0VEBiVHIrDVCLT88lCXVRQwJ9/UCrr9AxpBJCJnSYpEMGQOgRLBEO6IeQUb957U\nCCIRGSLhE4EZOhpeXqIsRZeXGM3llQXMKwiPINIaRCJil/CJoL1vgM5+M3Q0w+OmMFNDRyO5XS4+\nfn64VvDMfq1MKiJhsSSC1cBuYB9w1whl7rPe3wqssF1/CLNT2baI8tMxG9LsBZ4HCmIPeaiTneEP\ntJk5Gjo6kpqy/NAaRH2BAE/sUK1ARIxoicAD3I9JBtXAWmBJRJk1wALMdpa3AQ/Y3nvYujfS3ZhE\nsAj4rXU+LvaduGbkpNYeBGPhcrn4+LKy0PkLdc1DluUQkdQVLRHUYLafPAT0A49j9hi2uwF4xDre\nhHm6H5zJ9ArQMsz3td/zCPAHYwnabmiNQIlgNCtK81haEt6v4LEdxx2OSETiQbREUI7ZZ3hQvXVt\nrGUizSS8uX2jdT4ujRFNQzIyl8vFJ84P1wpePHR6yIY+IpKaoiWCWDcMjmyYj32jYVN2LOWHsDcN\nzcxVjSCa82fksWLmNAACwSCPbletQCTVRRti0wBU2s4rMU/8o5WpsK6NphHTfHQCmAWMuPbB+vXr\nQ8e1tbXU1tYOed/eNDQjxfYpHq9PXFDGlhfaAHj5SAsfrS6lylqKQkQSj8/nw+fzjfv+aENsvMAe\n4H3AMeANTIfxLluZNcA66+sq4LvW10FVwC+A823Xvg2cAr6F6SguYPgO42AwOHJlIRgM8if/vTW0\n8uijf7hcK4/G6J6X97Pp2BnAzD7+u/fOdzgiEZko1ujJmIdQRmsa8mM+5J8DdgJPYJLA7dYL4Bmg\nDtOp/CBwh+3+x4DXMKODjgK3Wtf/CbgOM3z0Gut8zNr7BkJJINPrZlq6ZzzfJiV9zNZX8Fp9K/tO\ndzoYjYg4Kd4H3Y9aI9h3upPPP78bgKr8LL53ffVUxZUUvvFqHa8eNYO63jMrn3uuWuBwRCIyESa6\nRhDX7COGNIdg7D62bBYu69/KW8fP8O7JdocjEhEnJHYi6NDQ0XMxJz+Lq6umh85/uLWB0WpgIpKc\nEjoRnOzSrOJz9fHzZ5HmNrWCPac6ea2+1eGIRGSqJXQiaNSs4nM2MyeDDy2cETr/0bvH8AdUKxBJ\nJQmdCJptK2gWaw7BuH20upScNDPiqr69hxfqmh2OSESmUoIngvAGK5pMNn7TMrz8yZLwRvePbj9O\njzUsV0SSX8Imgh7/AO19Zv9dr8tFvvYhOCcfXjSDoqw0AE739PM/2uheJGUkbCJostUGirPTcWsf\ngnOS6XXzMdsy1U/uauSMNroXSQkJmwjUPzDxrp1bRIW11WeXf4AntEy1SEpI2ERg32qxJDvNwUiS\nh8ft4tbl4RXEf7Wvach+0CKSnJIiEahGMHEuKc9naXF485qHt0ZbSFZEEl3CJoLmiD4CmRgul4s/\nX1EROn+tvlVLT4gkuQROBGoamiyLi3K4ek546YkfbKlnQJPMRJJWwiYCNQ1Nrj9bXk6Gx/zz2N/S\nxW8PnXI4IhGZLAmZCILB4JCmoRIlgglXkp3OH50X3kr6x+8eo6tfk8xEklFCJoKu/kBoQ5oMj5s8\nbUgzKT6yZOaQSWb/vavR4YhEZDLEkghWA7uBfcBdI5S5z3p/K7AihntrMNtebgHeBFaOJeihQ0fT\nBzdhkAmW6fVwywXh4aRP7WmksbN3lDtEJBFFSwQe4H7MB3o1Zr/iJRFl1gALgIXAbcADMdz7beCr\nmKTxD9Z5zJq71T8wVa6pms6CwmwA+gYCPLL1mMMRichEi5YIajB7ER8C+oHHgRsjytwAPGIdb8Js\nRF8a5d7jQL51XACMabD60KGjGjE0mdwuF5+2DSd96chptms4qUhSiZYIyjGbzg+qt67FUqZslHvv\nBv4FOAJ8B/jyWII+Za8RZKlGMNmWzcjjvZWFofMH3jqq4aQiSSRaIoj1r32sjfQ/AD4HzAb+Gnho\nLDef7g7XCAqzVCOYCn+xoiI0nPTQmW5+ua/J4YhEZKJEW7u5Aai0nVdinuxHK1NhlUkb5d4a4Frr\n+Eng+yMFsH79+tBxbW0ttbW1nLIlgiIlgilRkp3O2qWz+OG7phXvJ9uPceXsQiVikTjg8/nw+Xzj\nvj/ak7wX2AO8DziGGemzFthlK7MGWGd9XQV81/o62r1vY2oCL1nv/xPDjxwKDreZ+p3P7WJ/SxcA\n9153HouLcqL+h8q56x8IsO7Xu6i3FqK7pqqIL6yqcjYoETmLNZIy5paaaE1DfsyH/HPATuAJzAf5\n7dYL4BmgDtMx/CBwR5R7wYwu+jbwDvCP1nnM7DWC6XoinTJpHjd/+Z5wJe/FQ6fYpo5jkYQX7wPw\nz6oR+ANB/uCnWwgSxIWLpz+6Aq873v8zkss3X63jd0dbAKjKz+K+DyzBo/8HInFjomsEcaelp5+g\n1YddkOlVEnCAOo5FkkvCJYLT6ih2XEl2OmuXzQqd/9e2Y0NWgxWRxJLQiUD9A875g0Uzhmxr+eDb\nR6PcISLxKuESgYaOxoc0j5vPrJwdOn+tvpVXrX4DEUksCZcIWjSZLG5cMCOPD8wrDp3/x1tH6ejz\nOxiRiIxHwiUC+/ISRVpewnG3Li9nemZ4qeofalE6kYSTgIlAfQTxJC/Dy+22uQXPHmjS3AKRBJNw\niUCjhuLP5RUFrCovCJ3f/+YRev0BByMSkbFIuEQwpEaQqUQQD1wuF3/1nkqyvWanuPr2Hn6684TD\nUYlIrBIqEfT6A7RbnZEel4v8zGhr5slUKc5O55PLwyuUP7nrBAes9aBEJL4lVCJo7bWNGMpMw60t\nKuPK9QuKWVqcC4A/GOTe1w/RN6AmIpF4l1iJoCc8NLFAtYG443a5uPOSOUOWn/jJ9uMORyUi0SRU\nIrDPIShQ/0BcKs/L5FZbE9FTuxvZ2dThYEQiEk1iJYIeDR1NBB9cWMLymXkABIJB7t10iB7/gMNR\nichIEiwR2JqGMtQ0FK/cLhefr6kiJ82MIjre0ctD7zQ4HJWIjCShEkFrj5aXSBQzctL59IqK0Pmv\n9jex5USbgxGJyEhiSQSrgd3APuCuEcrcZ72/FVgR472fxexYth34VizB2puG1Fkc/66dWzRkotl3\nNx2mvVdrEYnEm2iJwAPcj/lAr8bsObwkoswaYAGwELPl5AMx3Hs1cANwAbAM+OdYgh06akg1gnjn\ncrlYd/FsplnNeM3dfdz35mGG24daRJwTLRHUYPYiPgT0A48DN0aUuQF4xDreBBQApVHu/Svgm9Z1\ngJi2uBrSWaxEkBAKs9L43Mo5ofPX6lt5Zn+zgxGJSKRoiaAcsO84Um9di6VM2Sj3LgSuBF4HfMDF\nsQSrpqHEdGlFAR9cUBI6//6Weg62ataxSLyI9mkaax1+rFN8vUAhsApYCfwUmDdcwfXr1wPQPxCg\nwVVOUfXFpLvdoREpkhj+/MIKdjZ3cLC1m75AgG+9dpDvvv88Mr36/yhyrnw+Hz6fb9z3R0sEDUCl\n7bwS82Q/WpkKq0zaKPfWA09Zx28CAaAIOBUZwGAiON7eyzu/2g5AfqYXl5aXSCgZXjdfunQun39+\nN70DAY629fDg2/XcWTMn+s0iMqra2lpqa2tD5xs2bBjT/dGahjZjmnGqgHTgJmBjRJmNwC3W8Sqg\nFWiMcu/TwDXW8SLr/bOSgJ0mkyW+2flZ/KVt74Ln65p56fBpByMSEYieCPzAOuA5YCfwBGbI5+3W\nC+AZoA7TMfwgcEeUewEewjQFbQMeI5xIRmRPBIXqKE5Y180tonbO9ND5/W8eoaG9x8GIRCTe21eC\ng0MNf7mviQfeOgLA6vnFfHalmhQSVVf/AJ97bhfHO3oBmD0ti3uvW0yW+n1EJoTVdB7z53vCzCxu\n7dGCc8kiO83D3ZfNI91t/vkdaevmu29ofoGIUxIzEWidoYS3YHo2n1k5O3T+u6Mt/PfuRgcjEkld\nCZMI7AvOqbM4OVw7t2jI/IJH3j3GO1qPSGTKJUwiUNNQcvr0iorQrmaBYJBv/f4gjZ29DkclkloS\nJhGcsS1Wlq+moaSR5nFz9+XzQkuGtPX6+cbv6uj1a4tLkamSOImgR4kgWU3PSuPLl8/Da00S3N/S\nxb9sOkRAncciUyIhEkGvP0CXtcOV1+UiN13DDJNNdUkut10Unmz26tEWfvTuMQcjEkkdCZEI2vrC\ntYFpGVpeIll9cGEJNyyaETr/2a4TPF+nlUpFJltCJAJ7R3G+Vh1Nan9xYQU1Zfmh8++9eUQjiUQm\nWUIkAntHcUGGRgwlM4/bxZcuncv8wmwA/MEg33i1jsNnuh2OTCR5JUQiaOsd2jQkyS0rzcNXr5hP\nkTVfpLN/gA0vH6Cluz/KnSIyHgmRCIZuUalEkApKstP52pULyPSaf6KNnb38w0v76ejTnsciEy0h\nEoHmEKSm+YXZfOnSuXiswQF1rV3c8/IBejTHQGRCJUYiGNJZrD6CVHJJeQGfs21es6O5g2++Wkf/\ngJKByERJjESgGkFKu3ZuEbetCM8x2Hz8DPduOsRAQBPORCaCEoEkhBsXz2Dt0lmh85ePtPDAW0e1\ndLXIBIglEawGdgP7gLtGKHOf9f5WYMUY7v0CZr/i6cO8FzJkeQl1Fqesjy2bxYcXhiecPXugif/c\nUq9kIHKOoiUCD3A/5gO9GlgLLIkoswZYgNmf+DbggRjvrQSuAw5HC1I1AgGz69JtF1VwtW2ry//Z\ne5IH31YyEDkX0RJBDWYv4kNAP/A4cGNEmRuAR6zjTUABUBrDvfcCX4oWYK8/QLdtnaEcbWeY0twu\nF5+/pIr3VhaGrv1i30n+/a2jWqROZJyiJYJy4KjtvN66FkuZslHuvdE6fzdagENqA5laZ0jAa80+\nvnJ2OBk8s7+Jf998RMlAZByitbPE+lc1lk/nLOArmGahqPd//Z717D1gFh47v+Zy4IIx/ChJVh63\ni79dNRe3y4Xv8GkAnj3QTCAI61bOxq0HBkkhPp8Pn8837vujJYIGTFv+oErMk/xoZSqsMmkj3Dsf\nqMJ0LA+WfwvTlHQyMoA/+/xd1L20H4DqmdOihCupxON28TeXVOF2wYuHTDJ4rq6ZHn+Az18yh3RP\nQgyKEzlntbW11NbWhs43bNgwpvuj/aVsxnQCVwHpwE3AxogyG4FbrONVQCvQOMq924GZwFzrVQ9c\nxDBJADRiSEbncbv4fE0V76sqCl176chp7nn5AF39Aw5GJpI4oiUCP7AOeA7YCTwB7AJut14AzwB1\nmI7hB4H4bWkEAAAQQklEQVQ7otwbadTmJ40Ykmg8bhd31szhQwtLQte2NLbx5Rf3DlnCXESGF+8N\nqcGH36nnZ7tOAPCJ88u42TapSMQuGAzyxM4T/HhbeGezstwMvl67kNLcDAcjE5la1qCamD/f474R\ntb1PS1BLbFwuFzcvncVnV84JdRYf6+jlb3+zh72nOh2OTiR+xX0isO9FkJeuRCDRrZ5fzFcunxfq\nLG7p6eeuF/fysjW6SESGSoBEEO7wU41AYnVpRQFfv2oBudYExL6BAN/6/UF+vO2Y5hqIRIj7RDC0\naUiziiV2y2bkce/7z6MiLzN07fEdx/nmq3X0+DWiSGRQ3CeCIdtUqmlIxqg8L5N/uW4x7ykNz0F5\nrb6VL/5mLyc6eh2MTCR+xH0iaLf3EahpSMYhN93L165cwB8sCq9cWtfaxZ3P7WJTQ6uDkYnEh7hP\nBH6rPTfT69ZMURk3j9vFpy+q5M6Vc/BaI4o6+ge455UDfH9LPX5tciMpLGE+WdUsJBPh/fOL+fa1\niynJTg9d+/meRu5+cS9NXX0ORibinMRJBGoWkgmyuCiH+z6whJqy/NC1Xc0dfPbXu/h9vZqKJPUo\nEUhKmpbh5atXzOfW5eV4rKai9j4///i7A9z7+iE6+zSqSFJHwiQCTSaTieZ2ufjIklK+ec0iirPC\nTUW/PXSKz/x6J1sb2x2MTmTqJEwiUI1AJsvSklzuX72EWtsWmE1dfXzlf/fyn28fpdcfcDA6kcmX\nMIkgT5PJZBLlZXj54qVzueuyuUNqn0/vPclnfr2Tt4+3ORidyORKmESgUUMyFa6cPZ3vrV7CxbPC\nHcnHO3r56kv7+M7vD9LSrWWtJfkkTCLQXgQyVYqy01l/5Xw+t3JOaK0iAN/h0/zVszt57kCz1iuS\npJIwiUCzimUquVwuPjC/mP9Ys3RI30F7n5/73jzMF17Yw86mDgcjFJk4sSaC1cBuYB9w1whl7rPe\n3wqsiOHe72B2LNsKPAXkMwp1FosTCrPS+OKlc7nnqgWU5oQ3t9l7upMv/nYP33qtjsZOrVkkiS2W\nROAB7sd8oFcDa4ElEWXWAAswexTfBjwQw73PA0uB5cBe4MujBaHho+Kk98zK53vXV3Nz9SzS3eE/\nm5ePtPCXv9rJj95t0B7JkrBiSQQ1mP2IDwH9wOPAjRFlbgAesY43AQVAaZR7XwACtnsqRgtCNQJx\nWqbXzScuKOOBNdVcUVkYut4XCPDEzhN86hfbeXLXCS1xLQknlkRQDhy1nddb12IpUxbDvQCfAp4Z\nKYB0j5tMb8J0Z0iSK83N4O7L5/Ht9y1m4fSc0PX2Pj8Pb23gL365g417T2r+gSSMWB6zYx0eEfNG\nyRH+DugDHh3uzb1PPki218P6HTOpra2ltrZ2nD9GZGItLcnl3usW4zt8mp9sO84Jq6+gpaefB98+\nylO7G/mj82by/nnFepCRSeXz+fD5fOO+P5YP71XAekw7P5i2/ADwLVuZ/wB8mKYfMJ3DVwFzo9z7\nSeDTwPuAnmF+dnDNY5uZW5DF/aurYwhVxBn9AwF+c/AUj+84QXP30FVMp2V4uWHhDD64sERNnDIl\nXGb9rJgfzmN5TNmM6QSuAtKBm4CNEWU2ArdYx6uAVqAxyr2rgS9i+gyGSwIhueooljiX5nFz/YIS\n/t8Hl3LbikoKMtJC77X1+vmv7ce49Rfb+P6Weu2MJnEnlk9YP7AOeA4zCugHmGGft1vvP4hp31+D\n6RjuBG6Nci/Av2GSwwvW+e+BO4YLIC9dy0tIYsjwurlx8Qw+ML+YF+qaeWp3IyetfQ56/AF+vqeR\n/9l7kpqyfD60sIQLZ+YNPr2JOCbe/wUG1zy2mffPK+bOmjlOxyIyZv5AkFeOnObJXY0cOtN91vsV\neZl8cGEJ11RNV81XJsxYm4YSIhH88Xkz+dSFo44uFYlrwWCQzcfb2Lj3JG+fOHsBu3SPm8sqCrhu\nbhEXzMzDrVqCnIOxJoKEeATRk5IkOpfLxcqyfFaW5VPf1sOv9jfx24On6LQmofUNBPAdPo3v8Glm\nZKdz7bwirp4znbK8TIcjl1QQ748dwTWPbWbdxbO5fkGJ07GITKju/gFePHSa5+ua2d/SNWyZhdNz\nuHJ2IVfMLhyyz7LIaJKyaejuy+ZxxezC6KVFEtSBli5eqDuF7/Bp2vv8w5ZZWpzLZZUFrCovoDQ3\nY9gyIpCkieD/1C7kwtJpTsciMun6BgJsajiD7/BpNh87g3+E5a7n5GdxSVk+NeX5LC7KUZ+CDJGU\nieBf37+EBdOznY5FZEq19/r5fUMrrxxpYWtjOwMjJIWCjDRWlk3jotJpXDAzj4LMtGHLSepIykTw\ngw8tU1VYUlprTz+vN5zhjYZWtjS20zcw8jpG8wqyWT4zjwtL81hWkkumV/NwUk1SJoKf/fGFZKfp\nH7MImIlpWxvbrMRwhtbekbfP9LpcLC7KYUlJLtXFOSwpztUyFykg6YaPelwusrRgl0hIptfNJeUF\nXFJeQCAYZM+pTt4+0cbWE+3sOdU5pF/BHwyyo7mDHc3h3dQq8jKpLsllSXEOi6bnUDktE4873p8J\nZTLF+//94M1PvcNjf7jc6ThEEkJX/wDbmzrYeqKNdxrbh53NHCnD42Z+YTbzC7NZMD2bBYXZSg4J\nLulqBNM0mUwkZtlpHmrK8qkpMzu/tnT3s7O5g13Nnexs7uDA6a6zRiL1DgTY2dzBTlutIcPjpnJa\nJnPys5idH/46IztdayMlobj/lM3VgnMi41aYlcbllYVcbu2o1usPsO90J7uaO9nV3MH+li5OdZ/d\nx9A7EGB/S9dZE92yvR4qrcRQkZdJWV4Gs3IzKM3N0J4LCSwBEkHchyiSMDK8bpbNyGPZjLzQtZbu\nfva1dLH/dBcHrK+ReyoM6vIPsOdUJ3tOdZ71XnFWOrOsxFBmJYcZOemUZKdTkOnVXIc4FvefslqC\nWmRyFWalUZMVbk4CM1z1SFsPh1u7OdLWw5EzPRw+0z3irGeA5u4+mrv72Hay/az3vC4XxdnpFGen\nMyMnjeIsc1ySk05JdhqFmWnkpXvVL+GQuE8EqhGITL2CzDQKMtO4wFZzCAaDtPb4OXymm8Nnejje\n0cuxjl6Ot/fQ2Nk34oQ3MKOXTnT2mu08m4Yv43a5KMjwUpCZRmGml/zMNKZnpVnXvCZZZHjJS/cy\nLcNDhset/ooJEsun7Grgu5iNZb7P0C0qB90HXA90Ybaf3BLl3unAE8Ac4BDwUcyuZmfJy1CNQCQe\nuFwuCrPSKMxKO2vJF38gSFNnH8c6ejjW3suJjl5OdPbR3NXHya4+2npHrkkMCgSDnO7p53TPyPMi\n7NLdbnLTPaHkkBc6Nl9z0zzkpHnISvOQneYhO81Nlnfw2INXtY+QaL8JD7AHuBZoAN4E1hLeZQzM\nzmTrrK+XAP+K2a5ytHu/DTRbX+8CCoG7h/n5wad3N3Lj4hnj+E+bHD6fj9raWqfDGEIxxS4e40qF\nmHr8AZq7womhqauP5q5+mjrN8Zle/6jNTgCndm6mqPriCYsp3eMm22sSRPZgwvC6yfR6yPC6yfC4\nw19Dx67QtUyvm62vv8oVV11FulUm3eMmze0izePG48KxGstEDx+twWw/ecg6fxyzx7A9EdwAPGId\nbwIKgFLMxvUj3XsDZnN7rHt9DJ8I4q5GkAp/tBMhHmOC+IwrFWLK9LqpmJZJxbSR91foGwhwptdP\nS3c/LT39nOnx09LTT4v19VfPbqfq0ito7/PT3jtAX2DkZTZi0TcQoG8gQOs5bCG998mfs6h7+AdV\nFy4rKbjwul2ke9x43S7S3O7QtTS3e0iZwffS3G48btO34na7hnz1ul24Q1/B63bhcZtzj8s1rn6W\naImgHDhqO6/HPPVHK1MOlI1y70zM5vZYX2eOFID6CERSQ7rHTUl2+oj7LvS+UMT666sB01/ROxC0\nkoKf9r4BOvr8tPX66egboK3PfO3uH6DLHzBf+wfo7A9fC4zSpzERggTpCwTpO7d8NSWifcrG+puK\nJQW5Rvh+wdF+Tq7WGBKRCC6Xi0yvi0zvyIljNIOJpNtKDl39A3T7B+jqD9DrD9A7YL38AXoGbNci\n3judlU5VflbovG8gQH8giD8QHLXzPNGsAn5tO/8ypk3f7j+Am23nuzFP+KPduxvTfAQwyzofzn7C\niUIvvfTSS6/YXvuZQF7gAFAFpAPvAEsiyqwBnrGOVwGvx3DvYCcxmL6Bf5rIoEVEZGJdjxn9sx/z\nVA9wu/UadL/1/lbgoij3ghk++htgL/A8poNZRERERETEWI3pN9jH2X0SU+UhzIimbbZr04EXcK4m\nUwn8L7AD2A58Lk7iysQMHX4H2Al8M07iAjOfZQvwiziJ6RDwrhXTG3ESUwHwJGZo907M6D6nY1qM\n+R0Nvs5g/r07HdeXMX9/24BHgYw4iOlOK57t1jFxENM582CakqqANIbvl5gKVwArGJoIvg18yTq+\ni6nv2ygFLrSOczHNbkviIC6AwU2lvZh+ovfGSVx/A/wE2GidOx3TQcwfqZ3TMT0CfMo69gL5cRCT\nnRs4jnkQcjKuKqAO8+EPZnWEP3M4pmWYz6hMzGfnC8B8h2OaEJcydLTR3Yww2WwKVDE0EQyOiALz\noTzSaKep8jRm5nY8xZWNmUW+FOfjqsD0RV1NuEbgdEwHgaKIa07GlI/5cIvk9O/J7v3AK9axk3FN\nxzx8FWIS5i+A6xyO6SOY5XsG/T0mAcTT/79x+Qjwn7bzjwP/5lAsVQxNBC22Y1fE+VSrAg4DecRH\nXG5M7a0d8zQCzsf1M0yt7irCicDpmOowTR2bgU/HQUwXYpr1Hgbexvzt5TgcU6SHgDusY6fjug3z\nb/wk8OM4iOk8THKajnkIew2z9tuYYorHnSSCTgcQo8Hxuk7IBf4b0x4YueavU3EFMB8qFcCVmKdw\nu6mO60OYP9YtjDzh0Ynf1eWY5HQ98BlME6STMXkxI/3+3fraydk1cCf/racDH8Yk9UhTHdd84POY\nh7AyzN/hxx2OaTdmMc/ngWcxD2MDY40pHhNBA6YtcFAlZnmKeNDI0IlwJx2IIQ2TBH6MaRqKl7gG\nnQF+BbwHZ+O6DLOm1UHgMeAazO/M6d/VcetrE/BzzHpeTsZUb73etM6fxCSEEw7GZHc98Bbhxaud\n/F1djHniPgX4gacwTdlO/64esmK7CvPkv5cx/p7iMRFsBhYSnoh2E+GOPqdtxHQOYX19epSyk8EF\n/AAzsuO7cRRXMeFRCVmYdtMtDsf1FcxDxFzMzPcXgU84HFM2pikPTPPL+zFNj07GdAKzJtgi6/xa\nzKiYXzgYk91aTCIf5OTvajdm0mwW5m/xWszfotO/q8FV72YDf4QZzeT0Z8KEGGki2lR6DDgG9GH+\nUG7F+Ylw78U0wbxDeFjd6jiI63xM+/I7mKGRX7SuOx3XoKsIP0w4GdNczO/oHcxQv8F/207/npZj\nagRbMU+5+XEQE5hk2Uw4eRIHcX2J8PDRRzA1dKdjetmK6R3CTbJOxyQiIiIiIiIiIiIiIiIiIiIi\nIiIiIiIiIiIiIvHo/wOHXTbs8FNk2wAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 45 }, { "cell_type": "markdown", "metadata": {}, "source": [ "If we land during an intererval of duration $x$, the time since last login is uniform between 0 and $x$. So the distribution of time since last login (`timesince`) is a mixture of uniform distributions." ] }, { "cell_type": "code", "collapsed": false, "input": [ "metapmf = thinkbayes2.Pmf()\n", "for time, prob in observed.Items():\n", " if time == 0:\n", " continue\n", " pmf = thinkbayes2.MakeUniformPmf(0, time, 101)\n", " metapmf[pmf] = prob\n", " \n", "timesince = thinkbayes2.MakeMixture(metapmf)\n", "print(timesince.Mean())\n", "thinkplot.Pdf(timesince)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "9.95253181595\n" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEACAYAAAC+gnFaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFz1JREFUeJzt3XuQXOV9p/Gn5yoJCY0kLrKEsGQEtrBjB3BAJDaMDXhl\nVRaqUsl62Yrt4CpD4iiOs7sxkN3ErtrFNjiVYEJiSBZixRWM1w7Boixu3jDl+AJIIAlkSZZG19Ft\nRjPSaKSRRjPT3fvHe0bquXb3qDXdmvf5VHV1n3Pet/s3LfX5nvOe06dBkiRJkiRJkiRJkiRJkiRJ\nGmQZsAXYBtw7SptHkuUbgGty5j8JtAJvD2k/G3gZ2Aq8BDSUsF5JUglVA83AQqAWWA8sGdJmObA6\neXwD8GrOsg8TgmFoEDwEfDF5fC/wtZJVLEkqqRuBF3Km70tuuR4DPpEzvQWYmzO9kOFBsAW4NHk8\nN5mWJJVBVZ7l84GWnOm9ybxi2wx1KWHIiOT+0jHaSpLOoXxBkC3weVLj7DfQtpj2kqQSqsmzfB+w\nIGd6AWGLf6w2lyXzxtJKGBI6CLwDaBup0UULLs+2t+zJ81SSpCG2A4sLbZxvj2AtcCVhnL+OcCxg\n1ZA2q4BPJY+XAp2cGfYZzSrg08njTwPPjtSovWUPu451k81mK+b2pS99qew1WNPkqsuarKnUN+CK\nPOvgQfIFQT+wAngR2AR8F9gM3JPcIJwxtINwdtHjwOdy+n8H+BlwFeE4wl3J/K8BtxFOH/0onjUk\nSWWTb2gI4PnkluvxIdMrRul75yjzDwO3FvDakqRzLN8egYZobGwsdwnDWFPhKrEuayqMNZ07Q8/2\nqTTZXce6eef0aeWuQ5LOG6lUCopYv7tHIEmRMwgkKXIGgSRFziCQpMgZBJIUOYNAkiJnEEhS5AwC\nSYqcQSBJkTMIJClyBoEkRc4gkKTIGQSSFDmDQJIiZxBIUuQMAkmKnEEgSZEzCCQpcgaBJEXOIJCk\nyBkEkhQ5g0CSImcQSFLkDAJJipxBIEmRMwgkKXIGgSRFziCQpMgZBJIUOYNAkiJnEEhS5AwCSYqc\nQSBJkTMIJClyhQTBMmALsA24d5Q2jyTLNwDXFND3euB1YB2wBvi1oqqWJJVMviCoBh4lrNCvBu4E\nlgxpsxxYDFwJ3A18s4C+DwF/TgiNv0imJUllkC8IrgeagV1AH/A0cMeQNrcDK5PHrwENwNw8fQ8A\nM5PHDcC+cdYvSTpLNXmWzwdacqb3AjcU0GY+MG+MvvcBPwH+khBGNxZVtSSpZPIFQbbA50kV+bpP\nAJ8H/hX4HeBJ4LaRGv71Aw/QUF8LQGNjI42NjUW+lCRNbk1NTTQ1NY27f74V+FLgy4RxfoD7gQzw\nYE6bx4AmwtAPhIPDNwOLxujbBVyYU0MnZ4aKcmV3dnWzcMa0Qv4WSRKQSqWgiA30fMcI1hIOAi8E\n6oBPAKuGtFkFfCp5vJSwUm/N07eZEBYAHwW2FlqwJKm08g0N9QMrgBcJZwE9AWwG7kmWPw6sJpw5\n1Ax0A3fl6Qvh7KK/BeqBk8m0JKkMih3bn2gODUlSkUo9NCRJmuQMAkmKnEEgSZEzCCQpcgaBJEXO\nIJCkyBkEkhQ5g0CSImcQSFLkDAJJipxBIEmRMwgkKXIGgSRFziCQpMgZBJIUOYNAkiJnEEhS5AwC\nSYqcQSBJkTMIJClyBoEkRc4gkKTIGQSSFDmDQJIiZxBIUuQMAkmKnEEgSZEzCCQpcgaBJEXOIJCk\nyBkEkhQ5g0CSImcQSFLkDAJJipxBIEmRKyQIlgFbgG3AvaO0eSRZvgG4psC+fwRsBjYCDxZesiSp\nlGryLK8GHgVuBfYBa4BVhBX4gOXAYuBK4Abgm8DSPH0/AtwOvB/oAy4uyV8jSSpavj2C64FmYBdh\nhf00cMeQNrcDK5PHrwENwNw8ff8A+GoyH+DQOOuXJJ2lfEEwH2jJmd6bzCukzbwx+l4J3AS8CjQB\nHyymaElS6eQbGsoW+DypcbzuLMIQ0q8B/xd4V5HPIUkqgXxBsA9YkDO9gLBlP1aby5I2tWP03Qs8\nkzxeA2SAOUDH0AL++isPMKu+FoDGxkYaGxvzlCxJcWlqaqKpqWnc/fNtydcAvwRuAfYDrwN3Mvxg\n8YrkfinwcHI/Vt97CENHXwKuAn4EXD7C62d3dHWzaMa0cfxpkhSnVCoFRYzU5Nsj6Ces5F8knAX0\nBGdW5ACPA6sJIdAMdAN35ekL8GRyexvoBT5VaMGSpNIqdmx/orlHIElFKnaPwG8WS1LkDAJJipxB\nIEmRMwgkKXIGgSRFziCQpMgZBJIUOYNAkiJnEEhS5AwCSYqcQSBJkTMIJClyBoEkRc4gkKTIGQSS\nFDmDQJIiZxBIUuQMAkmKnEEgSZEzCCQpcgaBJEXOIJCkyBkEkhQ5g0CSImcQSFLkDAJJipxBIEmR\nMwgkKXIGgSRFziCQpMgZBJIUOYNAkiJnEEhS5AwCSYqcQSBJkTMIJClyhQTBMmALsA24d5Q2jyTL\nNwDXFNH3vwEZYHaB9UqSSixfEFQDjxJW6FcDdwJLhrRZDiwGrgTuBr5ZYN8FwG3A7vGXL0k6W/mC\n4HqgGdgF9AFPA3cMaXM7sDJ5/BrQAMwtoO9fAV8cd+WSpJLIFwTzgZac6b3JvELazBuj7x3J9FtF\n1itJKrGaPMuzBT5PqojXnAr8GWFYaDz9JUkllC8I9hHG8gcsIGzJj9XmsqRN7Sh9rwAWEg4sD7R/\ngzCU1Da0gIe/8gCz6msBaGxspLGxMU/JkhSXpqYmmpqaxt0/35Z4DfBL4BZgP/A64aDv5pw2y4EV\nyf1S4OHkvpC+ADuB64DDI7x+dkdXN4tmTCv8L5KkyKVSKShipCXfHkE/YSX/IuEsoCcIK/J7kuWP\nA6sJIdAMdAN35ek7VKHDT5Kkc6DSx+bdI5CkIhW7R+A3iyUpcgaBJEXOIJCkyBkEkhQ5g0CSImcQ\nSFLkDAJJipxBIEmRMwgkKXIGgSRFziCQpMgZBJIUOYNAkiJnEEhS5AwCSYqcQSBJkTMIJClyBoEk\nRc4gkKTIGQSSFDmDQJIiZxBIUuQMAkmKnEEgSZEzCCQpcgaBJEXOIJCkyBkEkhQ5g0CSImcQSFLk\nDAJJipxBIEmRMwgkKXIGgSRFziCQpMgZBJIUuUKDYBmwBdgG3DtKm0eS5RuAawro+3Vgc9L+GWBm\nwVVLkkqmkCCoBh4lrNCvBu4ElgxpsxxYDFwJ3A18s4C+LwHvBT4AbAXuH+8fIUkav0KC4HqgGdgF\n9AFPA3cMaXM7sDJ5/BrQAMzN0/dlIJPT57Jx1C9JOkuFBMF8oCVnem8yr5A28wroC/AZYPXIL58t\noERJ0njVFNCm0DVxapw1/A+gF3hqpIUPf+UBZtXXAdDY2EhjYyPpbJatR49z6dR6ZifLJClWTU1N\nNDU1jbt/ISvvpcCXCeP8EMbyM8CDOW0eA5oIQz8QDg7fDCzK0/f3gM8CtwA9I7x2dkfXcRbNuIB0\nNkt1KpT7fEsrP9hzgLrqKr563dVcUFtInklSHFJhXVnwxnkhQ0NrCQeBFwJ1wCeAVUParAI+lTxe\nCnQCrXn6LgP+lHDMYKQQOG19x1H+++sbeWTTDjLZLD/YcwCA3nSGl/cfKuBPkCSNppBN6X5gBfAi\n4SygJwinfd6TLH+cML6/nHBguBu4K09fgL8hhMPLyfTPgc+NVMBjW3YCsOlIF292HB20rC+TGamL\nJKlAhY6pPJ/ccj0+ZHpFEX0h7CkU7VDPqTGXHznVy3N7DjJnSh3LL7t0YBdJkjSKSTe4/k/NLWzu\nPAbAvhM93P3uheUtSJIq3Hl/iYnUkOMhAyEA8GZ7J81dxye6JEk6r5z3QQCQyY5+huuL+9omsBJJ\nOv9U/NDQg29tG3P5j/a38aP9bVw1czp/fPUVZ/VavekML+5rI5WC/zD/EmqrJkVOStKYKj4ICrX1\n6HGaDrYPmz+ws3DwZA81qRQbDndxQU01N1w8a9iB5Jf3t/HDloMA1FdVcdv8S8553ZJUbpMmCAD2\nnxj56wjrOjp5fMuuQfOm19bwvlkXDpr33J6Dpx+v2nPQIJAUhUkVBCPZeKSLjUe6hs1/dveBYUEg\nSTE67wbBf7D7QLlLkKRJZVLtEWTHOHtoWNs8y9PZLP8vuXzFTXPneOBY0qQ1uYKghM+VyWb53s59\np6dvmXdxCZ9dkirHpNrMPVe/XJAbCJI02UyqIJAkFW9SBcH2Y93lLkGSzjuTKggOnRz7yqSSpOEm\nVRAU42hvHyf6+wtuv2rPAbp6+4BwIPmVA4dY3dJKTzp9rkqUpAkRbRAc7+vnvrWb8v6+wYDVLa18\nZ0c4aPxGeyff3bGPVXsO8PzeNrLZLK+2HealfW0Gg6TzTrRBAOEic09t31tw+3UdnQD8sKX19LwX\n97ayufM439q2h2d27Wd1zjJJOh9EHQQQfr+g7eQpXthb+Ao8M+RE1R/uPXONopfGednrbDbLyX73\nJiRNvEn1hbLx+os3N+dvlPiXXftpG3JQenvX2Z2tlM5meeitbeztPsknFy9g6SWzz+r5JKkY0e8R\nFOvlc/BDN2sOHWH38ROks1m+tW1PyZ9fksZiEFSAI8nZSCPp6OllbXsnvenMBFYkKSYODVWA1Cjz\ne9MZvvrWVo739XPT3Iv4L1dcNqF1SYqDewQTIJvNcuBEz6i/rZwaJQo2HD7K8b7wXYcfj/Dra5JU\nCu4RnCO96Qwn02lm1tXy5LY9rDl0hA/MmckfvGcR2WyW/mz29KWtU6PsEjgYJGkiGATnwMYjXXxr\n2x5O9Ke5+90LWXPoCAAbOo7Sk07zjV/s4MCJHj5z1eW8f/bMMlcrKXYODZ0Dj27awfG+fjLZLI9t\n2Tlo2RdefZudx7rpSaf5u807R3mGYLRjB5JUSgZBmfWk067wJZWVQ0Nldv/aTdxw8axB89KZLN/e\n3sKrbYfLVJWkmBgEZXayP03TgcFnBP3hzzeM2PbgiR6+vG4LAEsaZvCRd1zEs7sPsKRhBr+5YC5T\na6pPt/33gx38tK2D2+ZdwnUXNZy7P0DSea/SRyWy9/xkXblrOG/cfvk7WL7gUnrTGT7/6lun5z/2\nG786YvvDp3qZVVdLarTTliSdl5LPdMEfbPcIJpFVew6QSo184buedJq3D3eRSqX4lVkz+N7O/fyk\ntYPrLmrgs+9eOPHFSqoYBsEk84PdB0ac/+imnTR3HQfgqpnT2Xo0PH6jvZNPLk4zpbp6xH6SJj/P\nGorA7/90/ekQAE6HwICh33je3tXN32/Zxdr2zjGftzed4Y32zoJ/3Geo/Sd6Tn9zWlL5uEcgvv52\nMzdeMpsDJ3qoSqX4aWsHAG92dNJQt5jFF06no6eXn7R2cHXDDK6cOR2A7+/af/rSF48sfT911WG7\nYv+JHubU11I/xl7Gz9sOs3LbHuqrq/nf1y1hRq3/FaVyqfSjhB4srgA3zb1o0LWOvvDeK/hZ22Fe\nT74xPWB2fR2ZbJbO3j5m1tXyv65dwsbOLo719nPjJbNPBwWEvZTc5893Qb10Jstrh45QX13FtXNm\neoBbGkOxB4sLabgMeBioBv4P8OAIbR4BPg6cAH4PWJen72zgu8A7gV3AfwJGGocwCM5js+vrOHyq\nF4C506bw6cWXs2jGNA71nOLP3zjzY0DXzmlg8YUXcCqd4ZZ5F7Pu8FH+cetuAFZc/S7eN+tCfnyw\ng6e2t5ye956Z00mlUlQbCNIwpQ6CauCXwK3APmANcCeQ+5Ney4EVyf0NwDeApXn6PgS0J/f3ArOA\n+0Z4/YoLgv1vrmXetR8sdxmDxFzTJxcvoPXkqdNnStVUVfH5q9/F5dOnsuZQJy3dJ1nTfoST/Wmu\nmdPA5Xu3874bf53jff0suGAqU6qr6MtkmVJdNehCgBCOnfRlMmMOcZVCU1MTjY2N5/Q1imVNhanE\nmqD0p49eDzQTttoBngbuYHAQ3A6sTB6/BjQAc4FFY/S9Hbg5mb8SaGLkIKg4+9dV4Eo34pq+3dwy\naLo/k+GvNjaP2HZdRyf/8NT3+ODUOSV57brqKpZePJt0NstPWzuYXlvDx+ZfQm86w6VT69l+rJsp\n1dUsaZjO9JoassDbh7u44sILmFVXe/pj+tTqF7jhQx+mOpUiS5b2nl66+9PMqq+lOpViVl3toCvR\nVqdSHO3tY3ptDQdP9DB36hSqq8b+zJ/sTzOluoosUJWzF5XJZgdND3jllVcqbgVXiSvdSqxpPPIF\nwXwg95O2l7DVn6/NfGDeGH0vBQZ+Lb41mZbOK73pzKBjJ8f7+nlm1/5h7V7Y2zpsXq432zv5k9fe\nLnl9Z2NtSyutOcdxivXO6dPYd6KH/kzhF1OfWlPN7Po6Dp7o4eIp9dRWpaiuSlGTSnHw5CnWHGhn\n5sbtZMmSyQJkae7qDmFZX8fFU+rIZKEqBcf6+jnW18/cqVPIMvLvgMBom8zD5440ApkCXm07zDd+\nsb2gvy+VGu2XR0rjbJ47XxCM/g4WX0NqlOfLFvE6ks4Du4+fKLrPyf40+/pPAnDwZM+w5d19/fzy\n6LFh89PZLO09p2gf4TTmo2P8DGwptPf0srlzeE2TzVLghZzp+wlj+rkeA/5zzvQWwhb+WH23EIaP\nAN6RTI+kmTNB4c2bN2/eCruNPD46TjXAdmAhUAesB5YMabMcWJ08Xgq8WkDfgYPEEI4NfK2URUuS\nSuvjhLN/mglb9QD3JLcBjybLNwDX5ukL4fTRHwFbgZcIB5glSZIkKVhGOG6wjeHHJCbKk4QzmnJP\n55gNvEz59mQWAK8AvwA2Ap+vkLqmEE4dXg9sAr5aIXVB+D7LOuC5CqlpF/BWUtPrFVJTA/B9wqnd\nmwhn95W7pncT3qOB21HC//dy13U/4fP3NvAUUF8BNf1xUs/G5DEVUNNZqyYMJS0Eahn5uMRE+DBw\nDYOD4CHgi8nje5n4YxtzgYEfF5hOGHZbUgF1AUxL7msIx4k+VCF1/Vfgn4FVyXS5a9pJ+JDmKndN\nK4HPJI9rgJkVUFOuKuAAYUOonHUtBHYQVv4Qro7w6TLX9D7COmoKYd35MnBFmWsqiRsZfLbRfZTv\ny2YLGRwEA2dEQVgpj3a200R5lvDN7UqqaxrhW+Tvpfx1XUY4FvURzuwRlLumncDQb7SVs6aZhJXb\nUOV+n3J9DPj35HE565pN2PiaRQjM54DbylzTbxMu3zPgfxICoJL+/cblt4F/yJn+XeBvylTLQgYH\nQe5V1lJDpifaQmA3MIPKqKuKsPd2jLA1AuWv63uEvbqbORME5a5pB2GoYy3w2Qqo6VcJw3r/CLxJ\n+OxdUOaahnoS+FzyuNx13U34P94GfLsCanoPIZxmEzbCfka49ltRNVXi7xFky11AgQbO1y2H6cC/\nEMYDh36bpVx1ZQgrlcuAmwhb4bkmuq7fJHxY1zH6Fx7L8V79BiGcPg78IWEIspw11RDO9Pu75L6b\n4Xvg5fy/Xgf8R0KoDzXRdV0BfIGwETaP8Dn83TLXtIVwMc+XgOcJG2PpYmuqxCDYRxgLHLCAcHmK\nStDK4C/CDf9NyHOvlhAC3yYMDVVKXQOOAj8ErqO8df064ZpWO4HvAB8lvGflfq8GfkLuEPCvhOt5\nlbOmvcltTTL9fUIgHCxjTbk+DrxBeL+gvO/VBwlb3B1AP/AMYSi73O/Vk0ltNxO2/LdS5PtUiUGw\nFriSM19E+wRnDvSV2yrCwSGS+2fHaHsupIAnCGd2PFxBdV3EmbMSphLGTdeVua4/I2xELCJ88/3f\ngE+WuaZphKE8CMMvHyMMPZazpoOEa4JdlUzfSjgr5rky1pTrTkKQDyjne7WF8KXZqYTP4q2Ez2K5\n36tLkvvLgd8inM1U7nVCSYz2RbSJ9B1gP9BL+KDcRfm/CPchwhDMes6cVresAur6FcL48nrCqZF/\nmswvd10DbubMxkQ5a1pEeI/WE071G/i/Xe736QOEPYINhK3cmRVQE4SwbOdMeFIBdX2RM6ePriTs\noZe7ph8nNa3nzJBsuWuSJEmSJEmSJEmSJEmSJEmSJEmSJFWi/w86A6JqvrSuhwAAAABJRU5ErkJg\ngg==\n", "text": [ "" ] } ], "prompt_number": 46 }, { "cell_type": "markdown", "metadata": {}, "source": [ "The data is in the form of \"time since last login\", so we need to be able to look up a time, $t$, and get the probability density at $t$. But we have a PMF with lots of discrete times in it, so we can't just look it up. One option: Compute the CDF, generate a sample, and estimate the PDF by KDE:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "cdf = thinkbayes2.Cdf(timesince)\n", "thinkplot.Cdf(cdf)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 47, "text": [ "{'xscale': 'linear', 'yscale': 'linear'}" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAEACAYAAABMEua6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHBZJREFUeJzt3XtwXOd93vHv3rC47OIOgiRAEryJFCWRlGxJdBRHSCU7\nkkeJ0kwmsmp72iTjaKZVL5OZSJGT1JhpU9dtM+NJ1DiKK3lSt7XaOoojjS3Jci0ojmxRUixKoiSQ\nBO8ACAIEiesC2NvpH2d59gKAuwAX+57dfT4jDN737MHZn5bAg4N33/MeEBERERERERERERERERER\nEREREREpC88AF4H3r7HPnwIngHeBW0tRlIiIrM0nsYN6pVD/DPD9VPtO4I1SFCUiImvXw8qh/hfA\nQxn9AaBzvQsSEZGlvEU4RhdwPqM/BHQX4bgiIrJKxQh1AE9O3yrScUVEZBX8RTjGMLAlo9+d2pZl\n586d1smTJ4vwdCIiVeUksKvQnYsR6s8DjwLPAoeASezZMtlVnTyJZbnrBL6vr4++vj7TZWRxY03g\nzrq+/OUv89gf/QGReIyJhQjRZJyFRJz5eIz51OeFRJyFRIxoMlGSml548hv88qNfXNPXevDg83jw\nejz4PF6n7fV47c9c7af382Rsy3zck+p7PPDMf/4aX3zsd/EAXo8XD+DxeJiOLtASrMOT2s/+7LEf\nz9yW2v/qZ7tWsvYlo7+knTpW+v8T/ssff4Xf+4MvZf//Zxw78zVJ/Ze9DXKO6clo2w/mDh8s/5pn\n+8q/+/c88Ud/mGev7Ocu9NiF7pX1/+qBhkBwZ0GHSikk1L8N3A20Y4+dfxkIpB57Cnvmy2eAQWAO\n+M3VFCCyklgywVR0gbl4lJnoItOxRebiUU7NXMaLh59NjPA/Bt9Z9zo8eAh4vdR4fQRSH/5U3+/x\nEreSNPhrCPr8HK4N8YkN2/B7PPhTj/u8Xvuzx4Pf63UC25dqOyGdCtVieyXUzKENW4t+3OsRDgTZ\nWB82XcYSQZ+fxppa02Vcl0JC/eEC9nn0eguR6hRPJhlbmGViIeKE+ExskanYIpF4dMWvS2JhreKt\nGw8egj4/tTkfQZ+fGq+PhGXREqwj6PMT8Hid7TU+O5gLDdvvhZrZ37qx4LpEiq0Ywy9lq7e313QJ\nS7ixJrj+upKWxXR0gaHItH32HYtyevbydR3zhjtuo9YXoN4fIJQ6U/Z6PHTUNlDnD1DvC9jh7Q9Q\n4/XhXYez4Fxu/PdTTYVza12rsf7f5WmW28bUZX0kLYvLixEuL84zNj/LxGKE0fmZVR/Hi4fGmiAN\n/iCNgSChQA2h1OdwIEi9P4DPU6wJXCLulPorseCsruozdSmOxUScsflZRudnuTg/w3BketXHaAvW\nU+8PsKEuRGddmMZAkHAgWJKza5FKolCXVVtMxDk7O8nEwhwjkRkmFiMFjW/X+QLMJ2Lsb91E0Ouj\nu6GJ5po6any+ElQtUh0U6pKXZVlcWoxwbnaSobmpgoZS6n0Bkljsa+6ksy5EW209Df6aElQrUt0U\n6rKseDLJublJDo+dJ5ZMMJ+IrbivBw8twTo21YXprAuxsT5MOBAsYbUicpVCXRyxZIJjU5cYmZvm\n3NwkCSu54r4tNXV01DawNdRMd0MTQZ++lUTcQD+JVS5pWYxEpjk2Oc7p2SvXDPId4Va2hVrY0tBE\nnT+w4n4iYo5CvUpNRhc4NjnG8emJFS/yCfmDbA+30BNuYWNdWDNRRMqAQr2KJC2Lc7OTvH9llJEV\nph2GA0F2htvY2dhKW7B+XS5bF5H1o1CvAguJOMenLnH08igz8cUlj/s9XnY2trG3qYPOupCCXKSM\nKdQr2FwsyntXRhmYHFt2lcINtSFuaulkR7gVv1dXZopUAoV6BZqOLvLm+HlOzVxeclGQ3+NlT1MH\n+1s3lv1qdCKylEK9gkTiUfovnOL83NSSx+p8AW5r7+KGxnZdwSlSwRTqFWAhEeedS8N8MDm2ZEpi\nyF/DnRu2siPcqtkrIlVAoV7GkpbFR5NjvDl+fsmYea3Pz12dPQpzkSqjUC9Tw3PT/Hj0NFOxhazt\nzTV1fLy9i+0Kc5GqpFAvM5F4lB8Mn+Di/GzW9qDXz8fau7ippVNhLlLFFOplZGBynNdGTy3Zvr9l\nI7e1d2n9FRFRqJeDSDzK3579iOmcoRaA39i+n5ZgnYGqRMSNFOoud252kh+NnGQxGc/afs/mXewM\nt+rqTxHJolB3qYSV5I2xcxy9cjFruxcPn915QOuVi8iyFOoudHkxwvfOH8taPdHv8XLP5l30hFsM\nViYibqdQd5ljU+P0X8h+M3RDbYh7Nu/UZf0ikpdC3SUsy+InY2eXDLfc3NLJJzZs0zRFESmIQt0F\nElaSF88fZziSvWbLJzu3s69lg6GqRKQcKdQNm40t8tLQcSYWI862kL+GB7beSJOGW0RklRTqBl1e\njPDdMx8Ss9Lrtmyqa+SBrXs13CIia6JQN+Tc7CQvDh3L2nZr22Zub+/W3HMRWTOFugHDc9O8Mnwi\na9ununazI9xqqCIRqRQK9RI7MXWJH104mbXt/u49bA01G6pIRCqJQr2Ezsxc4dWMOeh+j5f7t+xh\nc32jwapEpJIo1Evk+NQlXs05Q/+VbfvoqG0wVJGIVCKFegnkLpnr93h5cNs+2hXoIlJkCvV1ttwa\n6P+45yZag/WGKhKRSqZQX0cX52eXBPpndxzQRUUism68BexzHzAAnAAeX+bxduAl4AhwFPhnxSqu\nnE1FF3jxfPY89Ae27FWgi8i6yhfqPuBJ7GDfBzwM3Jizz6PAO8BBoBf4E6r8L4BoIsFLQ8ezbmxx\nZ8cWuhqaDFYlItUgX6jfAQwCZ4AY8CzwYM4+F4Crc/IagQkgTpVKWhb/ffBnTEbnAfDg4dNduznY\nttlwZSJSDfKdUXcB5zP6Q8CdOft8A/gRMAKEgd8oWnVl6G/PfkjCSjr92zu62a4rRUWkRPKFulXA\nMb6EPZ7eC+wEXgEOADO5O/b19Tnt3t5eent7C6uyTLw8dJyxhVmnvz3Uyq06QxeRVejv76e/v3/N\nX59v5ahDQB/2mDrAE0AS+GrGPt8H/hh4PdX/f9hvqL6dcyzLsgr5HVGezsxc4eXh406/zhfgc7sO\n4vMU8l60iMjyUgv8FbzKX77EeRvYDfQANcBDwPM5+wwA96bancAe4BRVZCq6kBXoAL++/WYFuoiU\nXL7hlzj27JaXsWfCPA18BDySevwp4D8A3wTexf4l8RhweT2KdSPLsnj21LtZ2x7Yspd6f42hikSk\nmpVy4e6KHH55/eKZrPuK6hZ0IlJMxR5+kWs4MjGSFeg9oRYFuogYpVBfoyuL8xweT8/23FAb4t6u\nXQYrEhFRqK/Z/zn9Xlb/rs5temNURIxTCq3BuxMXsvqf2rybDXUhQ9WIiKQp1Ffp0sIcb4yfc/o+\nj5cdjbpiVETcQaG+CpZl8cPhwaxtn92x31A1IiJLKdRX4bXR00zFFpz+vZt3EQoEDVYkIpJNoV6g\ny4sRjk2NO/1toRZ2NrYZrEhEZCmFeoG+l3PDi95NOwxVIiKyMoV6AY5NjhOJR53+nR1bqPVV9X1A\nRMSlFOp5zMdj9GfcZzTkr9ENL0TEtRTqefzk4tms/q/23GSoEhGR/BTq13BxfpbBmQmnf1NzJw1a\nfVFEXEyhvgLLsvju2Q+cfnuwgZ/r3GawIhGR/BTqKzg9eyWrf1fnNryeUq5ULCKyegr1ZcSTSX48\netrptwcb2FgfNliRiEhhFOrLODx+noVEHAAvHj7dvdtwRSIihVGo54jEo3yQceOL29q7CGspABEp\nE7qCJse3Bt9x2vW+AAfbNhmsRkRkdXSmnmEqupDVP9C2WTe+EJGyosTK8NOxc1n9W1o6DVUiIrI2\nCvWU0cgMZzOmMX5q8+6rd/EWESkbCvWUf5gYdtohf5Dt4RaD1YiIrI1CHRiJTDM0N+X0P97RpbN0\nESlLCnXg70fPOO2mQC17mjrMFSMich2qPtRHItNcic47fa3vIiLlrOpD/cWMOxo1BmrZGmo2WI2I\nyPWp6lAfn58jbiWd/l06SxeRMlfVof7c2aNOuy1Yr7N0ESl7VRvqwxmzXQBubtloqBIRkeKp2lB/\nLWNpXYC9zZrxIiLlrypDfT4eYya26PQPtGrRLhGpDFUZ6kcmRrL6d3ZsMVSJiEhxVV2ox5NJjk9f\ncvofa9fVoyJSOaou1AenJ5y7GgHcojdIRaSCFBLq9wEDwAng8RX26QXeAY4C/cUobD1YlsVro6ec\n/sHWTQR9uk+IiFSOfInmA54E7gWGgbeA54GPMvZpBv4r8EvAENBe/DKL48PJMaft93i5pVVn6SJS\nWfKdqd8BDAJngBjwLPBgzj7/BPhr7EAHuIRLnZ657LQ31oWp99cYrEZEpPjyhXoXcD6jP5Talmk3\n0Aq8CrwNfKFo1RXRxEKE4ci007+hybV/UIiIrFm+4RergGMEgNuAe4B64KfAG9hj8K7x6oWTTrs9\n2MBuhbqIVKB8oT4MZE7i3kJ6mOWq89hDLvOpj78DDrBMqPf19Tnt3t5eent7V1vvmsSTSSYWI05/\nS6ipJM8rIrJa/f399Pf3r/nr803Q9gPHsM/CR4A3gYfJfqN0L/abqb8EBIHDwEPAhznHsiyrkBP/\n4vvZpWHeupT+XfTbN9yO31t1szlFpAylrqMp+GKafGfqceBR4GXsmTBPYwf6I6nHn8Ke7vgS8B6Q\nBL7B0kA36thU+r3bXeE2BbqIVKxSXkpp5Ex9fGGO586kl9j93M6DhALBktchIrIWqz1Tr/hT1qOX\nR512e7BBgS4iFa2iQz1hJRmcnnD628MtBqsREVl/FR3q718eJZmalen3eDnQpiV2RaSyVXSoHx5P\nXze1PdyKz1PR/7siIpUb6tPRhaz+vuYNhioRESmdig31o1cuZvU31ocNVSIiUjoVG+ofZazIeKhj\nq8FKRERKpyJDfTQyQ9xKOn3dVFpEqkVFhvpPx8457e76Jt0IQ0SqRsWFeiyZYGxh1unrLF1EqknF\nhfq52cms/raQLjgSkepR0aHeWRvS4l0iUlUqKvHiySQnMpYF+HhHt8FqRERKr6JC/e8vnsFKLQvQ\n4K+hq77RcEUiIqVVUaF+bGrcae9p6ri6ZKWISNWomFDPXRZga6jZUCUiIuZUTKgPzU1n9TfUNhiq\nRETEnIoJ9R9fPO209zVv0NCLiFSligj1hUQ8q3+gVeumi0h1qohQH8y4sTRAY02toUpERMyqiFB/\nfeys0765pdNgJSIiZpV9qEcTiaz+jboZhohUsbIP9fcuX8jqtwbrDVUiImJe2Yf6cCQ9lXF7qNVg\nJSIi5pV1qMeTSUbnZ5z+/taNBqsRETGvrEN9MGPxrpA/qPuQikjVK+tQPzGdnspYq7sbiYiUb6gn\nLYuRjPH0O7TMrohI+Yb60NxUVn9zg5bZFREp21A/NXPZaXc3NOHzlO3/iohI0ZRtEmaunb4z3Gaw\nEhER9yjLUJ/MWTt9R1jz00VEoExD/UTOAl41Pp+hSkRE3KUsQ/39y6NOe09Th8FKRETcpexCPZZM\nELPSi3gdbNtssBoREXcpu1B/L+MsHaBZa6eLiDgKCfX7gAHgBPD4Nfa7HYgDv1aEulZ0dvaK09as\nFxGRbPlC3Qc8iR3s+4CHgRtX2O+rwEvAut0c1LIsxhfmnP7mel1wJCKSKV+o3wEMAmeAGPAs8OAy\n+/1L4DvA+DKPFU1moHvxcENT+3o+nYhI2ckX6l3A+Yz+UGpb7j4PAl9P9a3ilLZU7tIAfm/ZvSUg\nIrKu8qViIQH9NeD3U/t6WMfhl3Ozk077rs6e9XoaEZGylW+92mFgS0Z/C/bZeqaPYQ/LALQD92MP\n1Tyfe7C+vj6n3dvbS29vb8GFRpMJLi7MOv3uhqaCv1ZEpFz09/fT39+/5q/Pd1btB44B9wAjwJvY\nb5Z+tML+3wReAJ5b5jHLstY+MjMwOc5ro6cAqPfX8IVdt675WCIi5cLj8cAqRkDynanHgUeBl7Fn\nuDyNHeiPpB5/avUlrs3YfPosfVOd7nAkIrKcdRv/XsZ1nak/NXDYaR/q2MqBtk3FqElExNVWe6Ze\nFtNHIvFYVn9Ps9Z7ERFZTlmEeuYNpkH3IxURWUlZhHrm0gB7tSqjiMiKyiLUpzJuirE11GywEhER\nd3N9qE9HF5iLR53+lgaFuojISlwf6q9fPOu0N9c3amkAEZFrKKuE1NrpIiLX5upQT1oW5+bS671s\n1w2mRUSuydWhPhKZdtp+j1frp4uI5OHqUL8QmXHacSuJ11PKC2BFRMqPq0P9yMSI076lZaPBSkRE\nyoNrQ92yLJIZy7nrLkciIvm5NtTP5dzlqDVYb6gSEZHy4dpQPz6Vvt1pvb9G4+kiIgVwbahPRxed\n9jZdRSoiUhBXhrplWVxanHP6u5raDFYjIlI+XBnqkxkLeAFs1J2OREQK4spQH5gcc9rNNXUaTxcR\nKZArQ/3DrFDXei8iIoVyZagHvD6nfWPzBoOViIiUF9eF+nR0kflE+p6kWu9FRKRwrgv1D66MOu3m\nmjqtny4isgquS8xIxln6dM4sGBERuTbXhfrg9ITT/kTnNoOViIiUH1eFejyZxEN6+uLuRl10JCKy\nGq4K9aG5KayMlRmDPr/BakREyo+rQj3zTkfbQ7p1nYjIarkq1DMvOmoJ1hmsRESkPLkm1C3LImEl\nnf62kFZmFBFZLdeE+lTO9MWO2gZDlYiIlC/XhHrmeDqAR4t4iYismmtC/dJCxGnf3t5tsBIRkfLl\nnlDPuClGm4ZeRETWxBWhHksmss7UNyjURUTWxBWhfiEy41x01Biopc4fMFyRiEh5ckWon5uddNqN\ngaDBSkREyluhoX4fMACcAB5f5vHPAe8C7wGvA/tXU8T5uXSot2voRURkzQpZXMUHPAncCwwDbwHP\nAx9l7HMK+AVgCvsXwF8ChwotImml13vpbtBNMURE1qqQM/U7gEHgDBADngUezNnnp9iBDnAYKHhO\n4lwsymw86vQ768KFfqmIiOQoJNS7gPMZ/aHUtpX8NvD9QgvIGnoJNuhORyIi16GQ4Rcr/y6OXwR+\nC7hruQf7+vqcdm9vL729vVxZnHe21WvWi4hUuf7+fvr7+9f89YVci38I6MMeKwd4AkgCX83Zbz/w\nXGq/wWWOY1nW0t8PTw0cTj9Rx1YOtG0qoCQRkeqQWjKl4HVTChnreBvYDfQANcBD2G+UZtqKHeif\nZ/lAL0iX3iQVEbkuhQy/xIFHgZexZ8I8jT3z5ZHU408B/xZoAb6e2hbDfoP1mjKHXgBag/WF1Cwi\nIiso5VKIS4ZfBibHeW30lNN/ZO+dJSxHRMT91mP4Zd0cn7rktLV+uojI9TMa6hfm02uo39q22WAl\nIiKVwVio5w7FtGk8XUTkuhkL9ZnYYlY/rIW8RESum7FQP5OxMiPo9nUiIsVgLNQvLaTvdLRJ672I\niBSFsVDPXPOlJ9xqqgwRkYpiJNQty2IhEXf6Wm5XRKQ4jIT6TCya1W+uqTNRhohIxTES6gNTY057\nY10Yr94kFREpCiOhHk8mnXbCSl5jTxERWQ0joT4dW3DauxrbTJQgIlKRjIT62Yw56pvq9SapiEix\nlDzUIxn3I/V5vFoeQESkiEoe6pcWIk47YSX1JqmISBGVPNSPXhl12ttCLaV+ehGRilbyUF9MJJy2\nlgcQESkuA6GevpI0FKgp9dOLiFS0koZ6NJFgKmM6Y3dDUymfXkSk4pU01Eci6TsdhQNBgr5C7nst\nIiKFKmmoD0yNO+0ar6+UTy0iUhVKGuq1GWfmOksXESm+kob6VDQ9nn5DU3spn1pEpCqUNNQnFtMX\nHm3UdEYRkaIraajHkvYc9YDXR6NuNC0iUnRGFvQK+Wt0o2kRkXVgJNQj8ZiJpxURqXhGQv1g2yYT\nTysiUvGMhHp7bYOJpxURqXhGQn1DbcjE04qIVDwjoV7j09WkIiLroeSh3lWvRbxERNZLyUP96lx1\nEREpvpKHek9YdzsSEVkvJQ/1zfWNpX5KEZGqUUio3wcMACeAx1fY509Tj78L3Hqtg7UG61ZTn4iI\nrEK+UPcBT2IH+z7gYeDGnH0+A+wCdgO/A3z9WgcMuGgd9f7+ftMlLOHGmsCddammwqimwrm1rtXI\nF+p3AIPAGSAGPAs8mLPPrwB/lWofBpqBzuUOtj3UutY614Ub/wHdWBO4sy7VVBjVVDi31rUa+UK9\nCzif0R9Kbcu3T/dyBwvrRtMiIusqX6hbBR4nd8nFZb+u3q9QFxFZT/nWvz0E9GGPqQM8ASSBr2bs\n8xdAP/bQDNhvqt4NXMw51iCwc+2liohUpZPY71sWhT91wB6gBjjC8m+Ufj/VPgS8UawnFxGR4rsf\nOIZ9pv1EatsjqY+rnkw9/i5wW0mrExERERGRtSnk4qVSeAZ7nP/9jG2twCvAceAH2NMxS2kL8Crw\nAXAU+FcuqKsWe2rqEeBD4CsuqOkqH/AO8IJLajoDvJeq6U2X1ETqOb8DfIT9b3in4br2YL9GVz+m\nsL/XTb9WT2D/7L0P/C8g6IKa/nWqnqOpNi6oKYsPe1imBwiw/Jh8qXwS+2rXzFD/T8BjqfbjwH8s\ncU0bgYOpdgh7mOtGF9RVn/rsx36P5OddUBPA7wL/E3g+1Tdd02nsH7hMpmsC+7qR30q1/UCTS+oC\ne8bdBewTGpM19QCnsIMc4H8D/9RwTTdj51Mtdna+gj25xC3/dgB8Angpo//7qQ9TesgO9QHSF0pt\nTPVN+i5wL+6pqx54C7jJBTV1Az8EfpH0mbrpmk4DbTnbTNfUhB1WuUzXddWngR+n2iZrasU+iWrB\n/sX3AvApwzX9OvDfMvp/iB3mbvm3A+wiv5HR/zzwZ4ZqgaWhfiWj7cnpl1oPcBYIY74uL/ZfVTPY\nZwm4oKb/i/2X1t2kQ910TaewhxPeBr7okpoOYg+ffRP4GfbPX4ML6rrqGeCfp9qma/od7O/xMeBb\nLqhpL/YvmlbsE6qfYK+rtaqa1nuVxkIvXnIDC3P1hoC/xh5Dm8l5zERdSexw6AZ+Afvs2GRND2D/\n4L3DytdWmHid7sL+RXM/8C+wh/hM1+THnoH256nPcyz969jU93oN8MvYv6BzlbqmncC/wT6Z2oz9\nM/h5wzUNYF8D9APgRewTq9wbUOStab1DfRh77OyqLdjLCLjFRew/ZwA2YQdHqQWwA/1b2MMvbqkL\n7De0vgd8zHBNP4e9xtBp4NvAP8J+vUy/ThdSn8eBv8FeK8l0TUOpj7dS/e9gh/uo4brA/uX3D9iv\nF5h9rT6OfSY8AcSB57CHi02/Ts+karsb+4z8OKt8ndY71N/GXr2xB/u39EOk3+Ryg+ex3xwh9fm7\n19h3PXiAp7FnKHzNJXW1k353vQ57nPEdwzV9CfuEYDvwWeBHwBcM11SPPVQG9vDGp7GH9kx/T41i\nr8V0Q6p/L/YMjxcM1wX2Kq/fzuibfK0GsC+WrMP+ObwX++fQ9Ou0IfV5K/Br2LNyTH9PLbHcxUsm\nfBsYAaLY3/S/iT129UPMTRX6eeyhjiOkp3vdZ7iuW7DHYo9gT9f7vdR206/VVXeTPjEwWdN27Nfo\nCPb0s6vf2254nQ5gn6m/i30G2uSCuhqAS6R/EeKCmh4jPaXxr7D/ajZd09+lajpCetjTdE0iIiIi\nIiIiIiIiIiIiIiIiIiIiIiIiIiIiUgn+P0bYPOBv1R50AAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 47 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Get a sample:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "sample = cdf.Sample(10000)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 48 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Estimate the PDF:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "pdf = thinkbayes2.EstimatedPdf(sample)\n", "thinkplot.Pdf(pdf)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEACAYAAACwB81wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmQHOd53/FvT/dce2AX900BAsBLJ6SYgizKXNm0CkJJ\nZCrlKooVmQ5TZbJKZiQnVTZNVxLtfw5dcYVmWCYZm1QhKklQWZZdUEyFkmKNrJQiiJRAiCIACgAJ\n4r6xiz1m5+ju/PH2zvYMdndmdmene3d+H9bU9PH2zMMB8Mw7T7/9NoiIiIiIiIiIiIiIiIiIiIiI\niIiISJvsBo4Cx4DHZ2jzdLD/ELAztP0J4A3gdeBrQHrhwhQRkfmwgePAFiAJvAbcUdNmD/BSsPwR\n4CfB8hbgLaaS/DeA31u4UEVEZDaJOvvvwiT8k0AJ2AfcX9PmPmBvsHwA6AfWAjeCY7oAJ3g+24qg\nRUSkefUS/kbgdGj9TLCtkTbXgL8ATgHngCHg+/MJVkRE5q5ewvcbfB1rmm3bgD/ElHY2AD3Av244\nMhERaSmnzv6zwObQ+mZMD362NpuCbQPAj4GrwfZvAb8OfDV88LZt2/wTJ040FbSIiHAC2N7MAfV6\n+K8COzC99BTwALC/ps1+4KFgeRemdHMReDNYz2J+AdwLHL4p4hMn8H0/do8vfelLkcegmBRTJ8al\nmBp7YKooTanXwy8DjwEvY0bsvAAcAR4N9j+PGaGzB3Nydwx4ONj3GvA/MV8aHvBz4H80G6CIiLRG\nvYQP8J3gEfZ8zfpjMxz758FDREQiVq+k07EGBgaiDuEmiqkxiqlxcYxLMS2c6UbXtJsf1KNERKRB\nlmVBkzlcPXwRkQ6hhC8i0iGU8EVEOoQSvohIh1DCFxHpEEr4IiIdQglfRKRDKOGLiHQIJXwRkQ6h\nhC8i0iGU8EVEOoQSvohIh1DCFxHpEEr4IiIdQglfRKRDLIqEX3DzlL1S1GGIiCxqjdziMFLXChc4\nO3Yc23LY2vtesk5P1CGJiCxKjfTwdwNHgWPA4zO0eTrYfwjYGWy7DTgYegwDX2gmON/3uZw/DYDr\nlzk1elQ9fRGROaqX8G3gGUzSvxN4ELijps0eYDuwA3gEeDbY/iYm+e8EPgyMA3/fTHDj5RGKXqGy\nXvQmODN2DN0SUUSkefUS/l3AceAkUAL2AffXtLkP2BssHwD6gbU1be4FTgCnmwluqHjppm0jpWtc\nnjjTzMuIiAj1E/5GqpP0mWBbvTabatp8FvhaM4F5vsdw8UplvSfZX1m+mH+H8fKNZl5ORKTj1Uv4\njdZOau+cHj4uBXwG+NtGgwIYKV3H9cvmBRJp3tVzJ93Ossr+a4WLzbyciEjHqzdK5yywObS+GdOD\nn63NpmDbpE8BPwMuz/Qmg4ODleWBgQEGBgaqyjl9qTUkrARrs1t4a+QXgCnt+L6PZdV+14iILD25\nXI5cLjev16iXLR3MydffAs4BP8WcuD0SarMHeCx43gU8FTxP2gd8h6k6fy2/9iSs65U5MnQAP/ih\nsKPvQ2TsLnzf58jQgUrPf9uyD9Dl9Nb7fxQRWXKCzm5TPd56JZ0yJpm/DBwGvoFJ9o8GD4CXgLcw\nJ3efBz4fOr4bc8L2W80ENVy6Ukn2WbuHjN0FmP/B3uSKSruR0rVmXlZEpKPFoR5yUw//7ZHXGS0N\nA7C+ayurMlPniYeLVzg1ehSAjN3Njr6diIh0moXo4UdiojxeWV6WXFm1ryfZjxX8P064Y1Xj9EVE\nZGaxS/ie71L2zdW0FhbJRLpqv205dCf7KusjRZV1REQaEbuEXwr12JOJ1LSjcFTHFxFpXgwTfrGy\nXNu7n7QslPBHS0N4vrvgcYmILHaxS/jFqh7+9Ak/ZWdIByN3fHxGS0NtiU1EZDGLXcIvNZDwobqX\nr7KOiEh9sU74qVkSfm9yeWV5tDy8oDGJiCwFsUv4RW+ispy0Z074WaenMjyz6E5U1f5FRORmsUv4\nJTdc0snM2C5h2XSFJlMbK6mXLyIym1glfN/3bxqWOZvwePwxlXVERGYVq4Rf9kuVOXRsy8G2Zp/M\nMzxdshK+iMjsYpXwGx2hM6nL6a3U8QtunrLq+CIiM1rUCT9h2WRD0yOP6S5YIiIzim3Cn21IZli3\nE6rj68StiMiMYpXwG7nKtlZ3UnV8EZFGxCrhN1vSAXPidmq65HHKXmlBYhMRWezilfDDY/Bnuegq\nzNTxeyrr6uWLiEwvXgl/DjV8qB2eqRO3IiLTiU3Cr73xiWPNftFVWNUFWDpxKyIyrUYS/m7gKHAM\neHyGNk8H+w8B4ZvM9gPfxNz4/DCwa6Y3aeTGJzMJT7Ew4Y7heuWGjxUR6RT1Er4NPINJ+ncCDwJ3\n1LTZA2wHdgCPAM+G9v0l8FJwzPsxiX9acxmhUwnScsjaU3X8cZV1RERuUi/h3wUcB04CJWAfcH9N\nm/uAvcHyAUyvfi3QB3wceDHYVwZmrLdU9/BnnjRtJtXDM5XwRURq1Uv4G4HTofUzwbZ6bTYBW4HL\nwJeBnwN/DXTN9EZzGZIZ1qV5dUREZjX77GQEM5nVV1tw94PX/hDwGPAK8BTwJ8B/rj14cHCQ4eJV\nJtwxfu1jO/mXn9zW4NtOCV9xmy+P4vkuCctu+nVEROIol8uRy+Xm9Rr1zozuAgYxNXyAJwAPeDLU\n5jkghyn3gDnBe0/w2v8P09MHuBuT8D9d8x6+7/u8deP1Ss98S++d9IZuYdioXw3/jIKbB2Br7/vo\nCY3eERFZSoKBLY2PbqF+SedVzMnYLUAKeADYX9NmP/BQsLwLGAIuAhcwpZ5bg333Am/M9EbzreFD\nzbw6KuuIiFSpV9IpY0oyL2NG7LyAGWnzaLD/ecwonD2Yk7tjwMOh4/8d8FXMl8WJmn0Vzd74ZCZd\nzjKuFS4AMK6ELyJSpamfAwvEL7slDg/9BDBTJbxn+Ufn9EJFr8CbQ68Er5Pgzv5dWFZsri0TEWmZ\nhSjptMXkFbZA3btczSaVSFemZPB8j7w7Nu/YRESWilgkfNefujLWmUfCB+jS/PgiItOKScIP9fAT\nyXm9li7AEhGZXiwSfjk09818SjpQPXPmeHkY32/0UgIRkaUtFgk/3MN3rPn18FOJLE7wK8H1XSZU\nxxcRAWKT8EM9/MT8eviWZWk8vojINGKR8KtLOvPr4YNubC4iMp1YJHy3RcMyJ1XdEKV8Q3V8ERFi\nk/BDwzLnWdIBSFfV8cuq44uIEJeE74V7+PMv6aiOLyJys1gk/HLVhVfzT/igOr6ISK1YJPxW1/BB\ndXwRkVqxSPie7wFgYbXspiXpRLbya0F1fBGRmCT8SbblTM4AN2+WZdX08lXWEZHOFq+E34IROmGq\n44uITIlXwm/RCdtJquOLiEyJVcKf79TItVTHFxGZEquEP9+pkWupji8iMiVeCb/FPXyAnlAdf7Q0\n1PLXFxFZLBpJ+LuBo8Ax4PEZ2jwd7D8E7AxtPwn8AjgI/LTeG7Xqoquw7mR/ZdnU8b2Wv4eIyGJQ\nr0ttA88A9wJngVeA/cCRUJs9wHZgB/AR4FlgV7DPBwaAa40E0+pROgCpRIZkIk3JK+D5LuPl0aq7\nYomIdIp6Pfy7gOOYnnoJ2AfcX9PmPmBvsHwA6AfWhvY3PLB+IUo6lmXRU9XLV1lHRDpTvYS/ETgd\nWj8TbGu0jQ98H3gV+P16wSxESQeqx+OPajy+iHSoel3qRgeuz9SLvxs4B6wGvoc5F/Cj2kZ/9eSL\nAKxMr+Pe3/wkAwMDDb5tY3pCI3XGyzfwfLdlUziIiLRDLpcjl8vN6zXqlVt2AYOYE7cATwAe8GSo\nzXNADlPuAZPU7wEu1rzWl4BR4C9qtvu/uGq+A27vv4tkItVw8M341fDPKLh5ALb2vreqzCMistgE\n09A0NRdNvZLOq5iTsVuAFPAA5qRt2H7goWB5FzCESfZdQG+wvRv4JPD6bG+2EDX8ST3OVILX8EwR\n6UT1MmwZeAx4GTNi5wXMCJ1Hg/3PAy9hRuocB8aAh4N964Bvhd7nq8B3Z3qjhGWTsBbusoDuZB9X\nC+cBGNUFWCLSgVozNeX8+L+4+iOSiTS39//agr1J2StxZOhAZf3O/l0LMgxURKQdFqKk0zatnkfn\nptdPJMnaPZV1TbMgIp0mNgm/1fPoTCc8Wkd1fBHpNPFJ+AvcwwfoSS6vLI/qAiwR6TCxSfgLddFV\nWJfTixWcGC64eYpeYcHfU0QkLmKT8NtxAjVh2XQ7U/PoqKwjIp0kPgm/DT18qB2Pf70t7ykiEgcx\nSvjtGSJZNZFaaVi3PRSRjhGbhO+0aUx8xu6unC8o+yXd9lBEOkZsEn67Sjq1tz1UHV9EOkVsEn47\nRulMCpd1NDxTRDpFbBJ+u2r4UH3idiyYLllEZKmLRcK3sNo6P33KzpC2swD4vsd4eaRt7y0iEpVY\nJHzbciYnAmqbcC9/RMMzRaQDxCPhRzBrZVUdXyduRaQDxCPht/GE7aTuZB9WMLPohDtGSdMsiMgS\nF4uEf0vPbW1/T9ty6NI0CyLSQWKR8JOJdCTv2xuaPVN1fBFZ6mKR8KNSNV1yaUjTLIjIktbRCT9j\nd+EkUgC4fpm8OxpxRCIiC6eRhL8bOAocAx6foc3Twf5DwM6afTZwEPj2HGNcMJZl0ZvU8EwR6Qz1\nEr4NPINJ+ncCDwJ31LTZA2wHdgCPAM/W7P8icBiIZb2kxwmXdZTwRWTpqpfw7wKOAyeBErAPuL+m\nzX3A3mD5ANAPrA3WN2G+EP6GJu+u3i7h8fjj5RHKXinCaEREFk69hL8ROB1aPxNsa7TNfwP+CPDm\nEeOCchJJupzeyromUxORpareJa6NlmFqe+8W8GngEqZ+PzDbwYODg5XlgYEBBgZmbd5yPcnllfl0\nRorX6U+tbuv7i4jUk8vlyOVy83qNemWWXcAgpoYP8ASmt/5kqM1zQA5T7gFzgncA+ALwu0AZyADL\ngL8DHqp5Dz/q4ZDj5RFO3DgEmGmab++/q+1z+4iINCPIUU0lqnolnVcxJ2O3ACngAWB/TZv9TCXx\nXcAQcAH4U2AzsBX4LPBP3JzsYyFr91SGZ5b9kmbPFJElqV5Jpww8BryMGbHzAnAEeDTY/zzwEubE\n7HFgDHh4hteK5SgdmByeuZzrhYsAjJSu0p1cVucoEZHFJQ51i8hLOgA3ild5Z/QIAGm7i1v7PhRx\nRCIiM1uIkk7H6En2Y1nm4yi44xTcfMQRiYi0lhJ+IGHZNTdFuRZhNCIiraeEH7IsNJnajaISvogs\nLUr4Ib2pFZXlsfINXXUrIkuKEn5IMpEm6/QEa75uiiIiS4oSfo1lyale/o3S1QgjERFpLSX8Gr2h\nhD9Suo7nuxFGIyLSOkr4NTJ2Nyk7A4Dnu5ojX0SWDCX8GpZl0ZecmjxtuHglwmhERFpHCX8a/alV\nleWR0jVclXVEZAlQwp9GxukmY3cB4PkeIxqTLyJLgBL+DPpCvfzh4uUIIxERaQ0l/Bn0hW6CMlK6\njuuVI4xGRGT+lPBnkLazZG1zEZaPrzH5IrLoKeHPIlzWGVJZR0QWOSX8WYQT/mhpmJJXiDAaEZH5\nUcKfRcrO0O30BWs+14I7YomILEZK+HWsyKyrLF8vXCAOd+cSEZkLJfw6liVX4lhJAEpeUTdGEZFF\nq5GEvxs4ChwDHp+hzdPB/kPAzmBbBjgAvAYcBv5sXpFGJGElWJ5eW1m/VrgQYTQiInNXL+HbwDOY\npH8n8CBwR02bPcB2YAfwCPBssH0C+ATwQeD9wfLdLYm6zVakp8o6I6XrFN2JCKMREZmbegn/LuA4\ncBIoAfuA+2va3AfsDZYPAP3AZJd4PHhOYb48FmU9JGVn6A3d/lC9fBFZjOol/I3A6dD6mWBbvTab\ngmUbU9K5CPwAU9pZlJanwydvL+L5XoTRiIg0z6mzv9EhKdYMx7mYkk4f8DIwAORqDx4cHKwsDwwM\nMDAw0ODbts+y5AqSiRQlr0jZLzFUuMiKzPqowxKRDpHL5cjlcvN6jdpEXWsXMIip4QM8AXjAk6E2\nz2GS+L5g/ShwD6ZXH/afgDzwX2u2+4tlqOPl/Bku5E8C5v63t/Z9mISlgU4i0n6WZUH9HF6lXrZ6\nFXMydgumDv8AsL+mzX7goWB5FzCESfarMPV8gCzw28DBZoKLm5WZ9aEhmgWGipcijkhEpHH1Sjpl\n4DFMOcYGXgCOAI8G+58HXsKM1DkOjAEPB/vWY07mJoLHV4D/08LY2y5h2azKbKz08i/nT7M8tQZL\nvXwRWQSa+jmwQBZNSQfA9V3eHHoF1zfTJW/s3sGK0Dh9EZF2WIiSjtSwg17+pMv50/gasSMii4AS\n/hyszKzHtkw1rOhNaFy+iCwKSvhzYFsOqzObKusX86coe6UIIxIRqU8Jf45WZtaTSmQAcP0yF/Pv\nRByRiMjslPDnKGHZrO/aWlm/VrhAvjwaYUQiIrNTwp+H3uQKepL9lfXz429pvnwRiS0l/HmwLIv1\nXe/GCkZGjZVvMFy8EnFUIiLTU8Kfp4zdxcrQnDrn82/jeuUIIxIRmZ4Sfgusyd6Ck0gBUPaKXNAJ\nXBGJISX8FrAthw1d766sXyucZ7w8EmFEIiI3U8JvkWXJlVU3STk7dlwncEUkVpTwW8SyLDZ0batM\npDbhjnG1cC7iqEREpijht1DKzrAms7myfjF/ipJXiDAiEZEpSvgttiqzkbTdBYDnu5wffzviiERE\nDCX8FktYCTZ0bausDxevMFK6HmFEIiKGEv4C6En20Z9aU1k/N3YCz3cjjEhERAl/wazv2oJt2YCZ\nQvnyxJmIIxKRTqeEv0CcRIq12S2V9cv5M0y449EFJCIdTwl/Aa1IryPr9ALg43NOY/NFJEKNJvzd\nwFHgGPD4DG2eDvYfAnYG2zYDPwDeAH4JfGHOkS5ClmWxsWt71eRq14sXI45KRDpVIwnfBp7BJP07\ngQeBO2ra7AG2AzuAR4Bng+0l4N8D7wF2AX8wzbFLWtbprroH7vnxtyl5xQgjEpFO1UjCvws4DpzE\nJPB9wP01be4D9gbLB4B+YC1wAXgt2D4KHAE2zCviRWhNdnPl7lie73Ju/IRKOyLSdo0k/I3A6dD6\nmWBbvTabatpswZR6DjQX4uKXsGw2dm+vrN8oXmWoeDnCiESkEzkNtGm0K2rNclwP8E3gi5iefpXB\nwcHK8sDAAAMDAw2+5eLRk+xneXot1wumhn9u/ATdzjJSdibiyERkMcjlcuRyuXm9Rm2Sns4uYBBT\nwwd4AvCAJ0NtngNymHIPmBO89wAXgSTwv4DvAE9N8/p+p5Q3XN/l+PBBit4EAF1OL+/ufT+W1cgf\ng4jIlCBvNJU8GinpvIo5GbsFSAEPAPtr2uwHHgqWdwFDmGRvAS8Ah5k+2XcU27LZ3HNbZdTOeHmE\nSxOn6xwlItIajST8MvAY8DImcX8Dc/L10eAB8BLwFubk7vPA54PtHwM+B3wCOBg8Jn8pdKQup5c1\n2Vsq65fypxgtDUUYkYh0ijjUEjqmpDPJ933eHnmdsfINwNwxa3vfTlKJdMSRichisVAlHWkxy7LY\n3HMbTiIJgOuXOTV6RBOsiciCUsKPSDKR5pbu2yv1/Hx5lHPjb0UclYgsZUr4EepO9rGua2tl/Xrh\nItcKFyKMSESWMiX8iK1Mr6c/tbqyfm7sBOPlkQgjEpGlSgk/YpZlsbF7Oxm7GzCzap4aPaL5dkSk\n5ZTwYyBh2byr5w5sy1z4XPKKnB49iu97EUcmIkuJEn5MpOwMm3turayPlW9wIf9OhBGJyFKjhB8j\nvckVrM2+q7J+ZeIsw8UrEUYkIkuJEn7MrM5soje5orJ+ZuwYBTcfYUQislQo4ceMZVls6t5RuerW\n811OjR7VRVkiMm9K+DHkJJLc0jN1UdaEO8bZMd00RUTmRwk/prJOLxu63l1ZHype4mrhfIQRichi\np4QfY8vT61ieXltZPz/+tmbWFJE5U8KPMcuy2ND1brJOb7DF59ToUYruRKRxicjipIQfc+airNtx\nEinAzKx5cvQNyl4p4shEZLFRwl8EamfWLLh5To4extXIHRFpghL+ItGdXMam7qkrcfPlkWAOfU2/\nICKNUcJfRPrTq9nQta2yPloa4vTYm0r6ItKQRhP+buAocAx4fIY2Twf7DwE7Q9tfxNzQ/PU5xigh\nKzPrq+6Je6N4VXfLEpGGNJLwbeAZTNK/E3gQuKOmzR5gO7ADeAR4NrTvy3T4jctbbU1mM6syGyrr\nI6XrqumLSF2NJPy7gOPASaAE7APur2lzH7A3WD4A9APrgvUfAdfnG6hMsSyLddmtrM5sqmwbKw3z\n9sgvKWsefRGZQSMJfyNwOrR+JtjWbBtpIcuyWNe1pWp2zXx5hBM3DjHhjkcYmYjEVSMJv9EJXKw5\nHifzsCa7uWoKhqJX4MSNQ4yU9KNKRKo5DbQ5C2wOrW/G9OBna7Mp2NaQwcHByvLAwAADAwONHirA\nyswGkol0ZcSO57ucHHmDNdnNrMncgmXVfheLyGKTy+XI5XLzeo1GMoEDvAn8FnAO+CnmxO2RUJs9\nwGPB8y7gqeB50hbg28D7pnl9X7NAtka+PMo7o4er7ofb5fSyufs2UnYmwshEpNWCjlxTvblGSjpl\nTDJ/GTgMfAOT7B8NHgAvAW9hTu4+D3w+dPzXgR8Dt2Lq/A83E6A0Luv0sG3ZB+h2+irbxssjHLvx\ncy7m38H1yhFGJyJRi8NvffXwW8z3fS5PnOFS/hR+6FSKbTmsymxkeXoNyeAGKyKyOM2lh6+Ev4SN\nl0c4O3Zs2lE7Xc4y+lKr6HH6SNtdqvOLLDJK+HIT3/cZKl7mUv4URW/6aZVty6HbWUZ3sp+eZD/p\nRFZfACIxp4QvM/J9j+vFywwVLzFWGp61rZNIsSy5kuXpNWTtHiV/kRhSwpeGlL0iN4pXGSkPMVYa\nxvVnPpmbtrOsSK9jeXodtmW3MUoRmY0SvjTN930KXp6x0jCj5SFGS0PTTsRmTvhuYGV6A3aikcs3\nRGQhKeHLvPm+z1j5BkPFiwwXr9w09XLCslmd2cSqzAYS6vGLREYJX1rK9V2GCpe4MnH2phO+TiLF\nmsxmVqTXYlm6rYJIuynhy4KYHOlzeeI0BTdftS+VyLAmewv9qdU6uSvSRkr4sqB83+d68SIX86du\nmoY5bWdZndlMf2qVevwibaCEL23h+S5XJ85zeeLMTSN8kokUKzMbWJFap5O7IgtICV/ayvXKXCmc\n48rE2ZtG9iSsBMuSq1iRXkeX06tyj0iLKeFLJMpeiWuF81ydOE/ZL920P5XI0JdaTX9qFRmnO4II\nRZYeJXyJlOe7DBUvc3XiPBPu2LRt0naWHqefnuRyupPLsC2VfUTmQglfYsH3ffLuKNcLFxkqXp72\nQq5JqUSGjNNN1u4mlciStrOk7Kyu6hWpQwlfYsfzXUZK1xkuXmGkdO2mC7lmkrazZO0esk4P3U4f\nGbtb5wFEQpTwJdY832WsfIPRkpnCwUzb3NifvZNI0eP005tcTk+yHyeRXNhgRWJOCV8WFc93Kbjj\n5N0xCu44BXeCopen4E5Q74ugy1lGt9NH1jG/ApJWSr8ApKMo4cuS4PkuE+4Y+fJo5RfBbDN6ghkG\n6lhJ7EQyqP9bWFjm2bJIkMAK2jiJJI6VIplIk7IzOFZSXxay6Cjhy5I0eRJ4pHSdkdJ18uWRlr5+\nwkqQTKRJJlI4Vhon+NJIBA+mvjqwrETl2bacqoe+NKSdFirh7waeAmzgb4Anp2nzNPApYBz4N8DB\nJo5VwpemlL0So6Uh8u4o+fIoeXd01pFA7ZCwbFKJDKlEmrTdRdruIhM8JzTVhCyAhUj4NvAmcC9w\nFngFeBA4EmqzB3gseP4I8JfArgaPhZgm/Fwux8DAQNRhVFFM0/N9H893Kfslyn6JH+Z+yMfvuRt8\nn8n/PN/Dn2zjlSj5RUruBEWvULdcNF+pRIaDP/4l9wz8Bk4iRdJK4yRS2Fai8ivCIoFlWeY5KEO1\nQxz+/GoppsbMJeHXu+rlLuA4cDJY3wfcT3XSvg/YGywfAPqBdcDWBo6NrTj+ASum6VmWZcoqOKTJ\n8tP/+zP23PuZho93vTIlr0DJL1DyirheGdd38fwyHh7g44e+PMyyF7QrU/aLsw43LXoT/PCH/8x7\ndm1v/P+ppnxkTZ6TsCYLTFRiIhRX9cluc+zk+YuEZWNbdlUZ6h+/t5/3f/R2bMsJ9juh9omq8yBT\nkc0WtXkG5vylFYe/U7XiGNNc1Ev4G4HTofUzmF58vTYbgQ0NHCsSOTvhYCccMsxt2gff93H9MkVv\ngoKbD0YcjTPhjs944/i6r4mPP1mmms8P4DrHjpaHODt2fB5vUJ9V9QVghb4WLJjmS+HyxFmODB2Y\nOXZr6vjKklV9kj501iX03tXxzBTtzXyGCpc4OfLGjEf5VcH6TBUtwnumllZnNrMstWKWOBZGvYTf\n6F81na2SjmVZVmX0T5fTW7XPDD3N05dcydrsuyh7RUp+EdcrBb8iXDzcyq+GqV770jH5/zNt6Xaa\nTZ7vUvZunpNptmMWWsGbYKR0vWWvV/aL9RtFYBfwv0PrTwCP17R5DvhsaP0osLbBY8GUfXw99NBD\nDz2aerT8p5kDnAC2ACngNeCOmjZ7gJeC5V3AT5o4VkREYuRTmNE2xzG9dIBHg8ekZ4L9h4AP1TlW\nRERERESWqt2Ymv8xpq/vt8OLwEXg9dC2FcD3gF8B38UMNW2nzcAPgDeAXwJfiEFcGcyw29eAw8Cf\nxSCmSTbmYr9vxyimk8Avgrh+GpO4+oFvYoZGH8aMmosyptswn8/kYxjzdz3qz+kJzL+914GvAekY\nxPTFIJ5fBsvEIKam2JhSzxYgSXQ1/o8DO6lO+H8O/HGw/DjwX9oc0zrgg8FyD6YsdkcM4uoKnh3M\nuZq7YxATwH8AvgrsD9bjENPbmH+QYVHHtRf4t8GyA/TFIKZJCeA8prMTZUxbgLcwSR7gG8DvRRzT\nezH5KYPJm98DtkUcU9M+SvUonj8JHlHYQnXCnxxpBCb5Hm13QDX+AXPFclzi6sJcOf2eGMS0Cfg+\n8AmmevhV1fV2AAACrklEQVRRxwQm4a+s2RZlXH2YRFYrDp8VwCeBHwXLUca0AtPBWo75Uvw28NsR\nx/Q7mKlpJv1HTKKPy59dQ34H+OvQ+ueA/x5RLFuoTvjhAbdWzXq7bQHeAXqJPq4E5pfYCKZ3QQxi\n+lvML7R7mEr4UccEJrkeBF4Ffj/YFmVcH8SU5L4M/Bzzb6874pjCXgQ+HyxHHdMjmL/jl4CvxCCm\n2zFfQiswna0fY+YvazqmKGd18iN872ZMjnmNQg/wd5iaXe0UkVHE5WESxybgNzC96ihj+jTmH+VB\nZr74L6o/v49hvog+BfwBpnQY1u64HMwIur8Knse4+Rd1VJ9VCvgM5su7Vrtj2gb8IaajtQHzb/Bz\nEcd0FDPx5HeB72A6XbWzBTYUU5QJ/yymXjdpM2b6hTi4iPmJBLAek1TaLYlJ9l/BlHTiEheYk2v/\nCHw44ph+HTOX09vA14HfxHxecficzgfPl4G/x8xLFWVcZ4LHK8H6NzGJ/0KEMU36FPAzzGcF0X5O\n/wLTg74KlIFvYcrPUX9OLwax3YPpyf+KOXxOUSb8V4EdTF2Y9QBTJ92ith9zoobg+R9mabsQLOAF\nzEiKp2IS1yqmRgFkMXXNgxHH9KeYjsJWzNXe/wT8bsQxgfnZPTnHQjemPv16xHFdwMxtdWuwfi9m\nJMq3I4xp0oOYL+xJUX5ORzEXkGYx/w7vxfw7jPpzWhM83wL8K8zooaj/njctDhdmfR04BxQx/yAe\nxtTKvk90w53uxpRPXmNqyNruiON6H6b2+xpmuOEfBduj/qwm3cNUhyHqmLZiPqfXMMPoJv9uRx3X\nBzA9/EOYnmtfDGLqBq4w9QVJDGL6Y6aGZe7F/NqOOqZ/DmJ6jalSatQxiYiIiIiIiIiIiIiIiIiI\niIiIiIiIiIiIiIgsfv8fZLH2Cj251wYAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 49 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Second option: use numerical differentiation to compute the derivative of the CDF, which is the PDF:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "import scipy\n", "import numpy\n", "xs = numpy.linspace(0, 90, 101)\n", "ys = [scipy.misc.derivative(cdf.Prob, x) for x in xs]" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 50 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Numerical differentiation is more accurate, especially near zero. The value at zero is wrong: there are ways we could fix it, but it's not necessary because we won't get zero as data." ] }, { "cell_type": "code", "collapsed": false, "input": [ "thinkplot.plot(xs, ys)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Warning: Brewer ran out of colors.\n" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEACAYAAACwB81wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucFNWd9/FPMwPDnVFQUBgzyM1BgyJZRINrGzUBNJKY\nbJToYxI3AYM80TybDWs2TzK+NjdNsmsIPoob9YWYCGqCwQvBxNgbNXEUucllcACRmxC8zMidwenn\nj1+1VdN0T1dPd0/VdH/fr1e9+pyq09U/Bvh1zTmn6oCIiIiIiIiIiIiIiIiIiIiIiIiIiHSQSUA9\n0ADMTtNmjnN8NTDWs/9m4DVgrVMWEZGQKgM2AdVAV2AVUJPUZgrwtFM+D3jJKZ+FJfvuznn+CAwr\nbLgiIpJOlwzHx2MJfyvQDCwEpia1uRKY75TrgEpgEPbFUAccBj4A/ge4Kh9Bi4hI9jIl/MHAdk99\nh7MvU5tTsav7C4ETgZ7A5cCQXIIVEZH2K89wPO7zPJEU++qB24FngAPASqDFf2giIpJPmRL+TqDK\nU6/CruDbajPE2Qdwv7MB/AjYlvwBw4YNi2/evNlvvCIiYjYDw/N5wnLnpNVANzIP2k7AHbQFONl5\nPQ3YAPRN8RnxeDweX7YsHr/kknj8e9+Lx1ta4oH7/ve/H3QIx1FM/igm/8IYl2LyB/89MK0SeluO\nAbOAZdhMm/ucxD3DOT7PSfZTsMHdA8BXPO9/DOiPDfjOBN5P90ELFkBjIzz1FHzhC3Dmmdn+UURE\npC2ZEj7AUmfzmpdUn5Xmvf/oN5B9+9zy6tVK+CIi+ZZplk6HOXzYLa9ZE1wcCdFoNOgQjqOY/FFM\n/oUxLsVUOKlm13S0eDwe5+KL3av8gQOta0dERFKLRCKQZQ4PzRX+kSNuec8e20REJH9CkfBbWuDo\n0db7wtCtIyJSTEKR8JOTPSjhi4jkWygSvrc7J+G11zo+DhGRYhbahF9fn3q/iIi0T2gT/rFjsGFD\nx8ciIlKsQpvwQf34IiL5pIQvIlIiQpHwvXfZVla65TVrIJ7144FERCSVUCR877TMkSOhd28rv/su\n7NoVTEwiIsUmFAnf26XTowd89KNuffXqjo9HRKQYhS7hV1S0Tvhr13Z8PCIixSiUCf+ss9y6pmaK\niORHKBK+d9C2ogLOOMOtv/66zckXEZHc+En4k7AFyRuA2WnazHGOrwbGevbfCqwDXgN+A1SkenPy\nFf6JJ8KgQe4xLXkrIpK7TAm/DJiLJf3RwDRSr2k7HBgBTAfudvZXA18DzgU+6pzrmlQf4p2lU+F8\nJXhXvFK3johI7jIl/PHYWrVbsXVpFwJTk9pcCcx3ynVAJTAQW7+2GeiJLaXYE9iZ6kOSr/ABajxf\nK0r4IiK5y5TwBwPbPfUdzj4/bd4Ffg5sA3YBjcCfUn1IqoQ/erS7b/36DFGKiEhGmRYx93ufa6pl\ntoYBt2BdO03Ao8C1wK+TGz75ZO2HN1ht2hQFoq0GbhsarNunWzef0YiIFJlYLEYsFsvpHJkS/k6g\nylOvwq7g22ozxNkXBf4KvOPs/x1wASkS/oUX1tLUZOVzz7XXvn2hqgq2b7dZOg0Nrfv1RURKSTQa\nbbWY+m233Zb1OTJ16SzHBmOrgW7A1cCSpDZLgOud8gSs62YPsNGp98B+A7gUSNk54+3S6d7dLXu7\nddSPLyKSm0wJ/xgwC1iGJetFwAZghrMBPA1swQZ35wEznf2rgAexL43Ecy/vTfUh3lk63m4bDdyK\niORPpi4dgKXO5jUvqT4rzXvvcLY2Jd94leBN+OvWZTqLiIi0JRR32qaapQN2x23EGQ7esqX1F4OI\niGQn1Am/Vy+orrZySwts3NihYYmIFJXQJXzvoC1oPr6ISL6EIuGnG7QFDdyKiORLKBJ+ukFb0BW+\niEi+hCLhp+vDB1vysKzMym++CQcOdFxcIiLFJHQJP7kPv3t3GDrUyvE41Nd3XFwiIsUkFAm/rT58\naN2to4QvItI+oUj4CRUV7rx7L+/ArfrxRUTaJ1QJP93TMPVMHRGR3IUq4ScP2CYMHw7lzkMgtm2D\n/fs7LiYRkWIRqoSfPGCbUFEBw4a5dfXji4hkL1QJP90VPqgfX0QkV50m4asfX0QkN6FK+G0tYahH\nLIiI5CZUCb+tK/xhw9yB2x07+HBJRBER8cdPwp8E1AMNwOw0beY4x1cDY519o4CVnq0J+EZbH5Ru\n0Bbs6n/ECLeuRyWLiGQnU8IvA+ZiSX80MA2oSWozBRiOrX07Hbjb2b8RS/5jgXHAQWBxWx/W1hU+\n6EFqIiK5yJTwx2Nr1W4FmoGFwNSkNlcC851yHVAJDExqcymwGdje1odlSvjqxxcRab9MCX8wrZP0\nDmdfpjZDktpcA/wmUzDZJHxd4YuIZCfTIuZxn+dJfgKO933dgE+Tvv+fXbtqAXjuOZgwIUo0Gk3Z\nbtgw68s/ehTeegsaG6Gy0meEIiKdWCwWIxaL5XSOTAl/J1DlqVdhV/BttRni7EuYDLwK7E33Iaee\nWgvA1KmQJtdbsOX2fPy1a62+di1MnNhm/CIiRSEabX0xfNttt2V9jkxdOsuxwdhq7Er9amBJUpsl\nwPVOeQLQCOzxHJ8GPOwnmExdOgBjxrjl1av9nFVERCBzwj8GzAKWAeuBRcAGYIazATwNbMEGd+cB\nMz3v74UN2P7OTzB+Ev7ZZ7tlJXwREf8ydekALHU2r3lJ9Vlp3nsAGOA3mGwT/rp10NwMXbv6/QQR\nkdLVae60TRgwAAY784SOHNGTM0VE/Op0CR/gnHPcsrp1RET86ZQJX/34IiLZK4qEH/d7t4CISAnr\nlAl/6FDo08fK775rT88UEZG2dcqE36WL5uOLiGSrUyZ8aN2ts2pV/mMRESk2RZHwdYUvIpJZp034\nZ57proD1xhtaAUtEJJNOm/C7d4czznDra9bkPx4RkWISqoTf1hKHqagfX0TEv1Al/G7dsmuvfnwR\nEf9Ck/C7dYNI8jIqGXgT/vr1tjCKiIikFpqEn03/fUL//nDaaVY+elTr3IqItKVTJ3zQg9RERPwK\nTcLPdsA2QQO3IiL++En4k4B6oIH0C5HPcY6vBsZ69lcCj2GrZK3HlkBMKV9X+C0t7TuPiEixy5Tw\ny4C5WNIfja1PW5PUZgowHFv7djpwt+fYL7AlEGuAMVjiTynbGToJp50GlZVWbmqCrVvbdx4RkWKX\nKeGPx9aq3Qo0AwuBqUltrgTmO+U67Kp+INAPuBC43zl2DEh7P2x7r/AjkdZX+erWERFJLVPCHwxs\n99R3OPsytRkCDAX2Ag8AK4D/Bnqm+6D2JnzQwK2IiB+ZFjH3u7RI8gz6uHPuc7EFzl8B7gT+Dfhe\n8pt37aplxQqorYVoNEo0GvX5sUYJX0SKXSwWIxaL5XSOTLc6TQBqsT58gFuBFuB2T5t7gBjW3QM2\nwHuRc+6/YVf6ABOxhH9F0mfEx42L88lPwo9+lP0fAKC5GaJRW9QcYOlSOOmk9p1LRKQziNidqlnd\nrpqpS2c5NhhbDXQDrgaWJLVZAlzvlCcAjcAeYDfW1TPSOXYpsC7dB+XSpdO1K5x1llvXVb6IyPEy\nJfxjWJfMMmxa5SJsps0MZwObhbMFG9ydB8z0vP9/A7/GpmuOAdJew7d3lk6CBm5FRNqWqQ8fYKmz\nec1Lqs9K897VwD/4CaS9N14l6AYsEZG2heZO21y6dMDWuO3i/Gk2boTGxtxjEhEpJkWT8Hv3dvvx\n43F4+eXcYxIRKSZFk/ABzj/fLf/tb7mfT0SkmIQm4ec6aAutE/5LL9mVvoiImNAk/FwHbQFqaqBf\nPyvv3QubN+d+ThGRYhGahJ+PLp2yMhg/3q3/9a+5n1NEpFgUVcKH47t1RETEFF3CP+88t7xyJRw6\nlJ/zioh0dkWX8AcOhOHDrdzcDCtW5Oe8IiKdXdElfIAJnnW1ND1TRMQUZcJXP76IyPGKMuGfc447\nzXPrVti1K3/nFhHprIoy4VdUwLhxbv3FF/N3bhGRzqooEz7Axz/ull94Ib/nFhHpjEoi4b/yiqZn\nioj4SfiTsGULG4DZadrMcY6vBsZ69m8F1gArgTafX5mPZ+l4DR4Mp59u5aNHLemLiJSyTAm/DJiL\nJf3RwDSgJqnNFGA4thTidOBuz7E4EMW+BMaTRrdu7rPs8+nCC92y+vFFpNRlSrPjsaULtwLN2ELl\nU5PaXAnMd8p1QCUw0HM84yK7+e7OSUjux9fTM0WklGVK+IOxhcgTdjj7/LaJA3/CFkP/WroPKVTC\nP/ts6NPHynv2QENDYT5HRKQzyJTw/V4Tp7uKn4h150wGbgIuTNUoH49GTqWsDC64wK1rto6IlLJM\ni5jvBKo89SrsCr6tNkOcfQCJW572AouxLqLnkz+koaGW2lorR6NRotFoxsD9mjgRli2z8vPPww03\n5O3UIiIdJhaLEYvFcjpHpv71cmAjcAmWvF/GBm43eNpMAWY5rxOAO53Xntig7z6gF/AMcJvz6hW/\n7ro4Cxbk9OdIq6kJLrsMWlogEoFnnoETTijMZ4mIdJRIJAI+xki9MnXpHMOS+TJgPbAIS/YznA3g\naWALNrg7D5jp7B+EXc2vwgZzn+T4ZA/AT3+aTcjZ6dcPxoyxcjyuRVFEpHRl9e1QIPF4gafPPPAA\n3HWXlS+7DH7844J+nIhIwRXiCr8oTJzoll96CT74ILhYRESCUhIJf/hwOPlkK+/bB6+9Fmw8IiJB\nKImEH4m0np6pfnwRKUUlkfCh9V23esyCiJSikkn448dDuXPXwcaNsHdvsPGIiHS0kkn4vXrZSlgJ\nWutWREpNySR8ULeOiJS2kkr43oHbujo4diy4WEREOlpJJfzTT4dBg6y8fz+sWRNsPCIiHamkEn4k\nom4dESldJZXwQYuiiEjpKrmE/7GPuQuubN6sbh0RKR0ll/B79oRJk9z6b34TXCwiIh2p5BI+wLRp\nbvm55+Ctt4KLRUSko5Rkwh8+3O68BVsY5ZFHgo1HRKQjlGTCB/jiF93y4sVw8GBwsYiIdAQ/CX8S\nUA80ALPTtJnjHF+NLVruVQasBJ5oZ4wFccEFcNppVt6/H558Mth4REQKLVPCLwPmYkl/NLaebU1S\nmynAcGAEMB24O+n4zdjyiKGaANmlC1xzjVtfuNC6d0REilWmhD8eW6t2K9AMLASmJrW5EpjvlOuA\nSmCgUx+CfSH8inAsp9jKFVdA795W3rZNz8kXkeKWKeEPBrZ76jucfX7b/Bfwr0Aor5179oTPftat\nP/xwcLGIiBRapoTvtxsm+eo9AlwB/B3rvw/d1X3CF75g3TtgD1TbvDnYeERECqU8w/GdQJWnXoVd\nwbfVZoiz73NYd88UoDvQF3gQuD75Q2praz8sR6NRotGon9jz4pRTIBqFP//Z6osWwXe+02EfLyLi\nSywWIxaL5XSOTFfe5cBG4BJgF/AyNnC7wdNmCjDLeZ0A3Om8el0EfAv4dIrPiMcDfqDNihUwfbqV\nu3eHp56Cfv0CDUlEpE2RSASy7D3J1KVzDEvmy7CZNouwZD/D2QCeBrZgg7vzgJlpzhWqWTpeY8fC\nqFFWPnwYHn882HhERAohDH3rgV/hAzzxBNx2m5UHDYLf/x7KyoKNSUQknUJc4ZeMT34STjjByrt3\n2zN2RESKiRK+o6ICrrrKrS9eHFwsIiKFoITvcdVVtioWwCuvwDvvBBuPiEg+KeF7DBwI555r5ZYW\nWLYs2HhERPJJCT+Jd3GUP/whuDhERPJNCT/JJZdA165WXr8e3nwz2HhERPJFCT9J376tFzrXVb6I\nFAsl/BQmT3bLS5dCCG4TEBHJmRJ+ChMnuo9N3rED1q0LNh4RkXxQwk+hogI+8Qm3vnRpcLGIiOSL\nEn4a3tk6y5bBgQPBxSIikg9K+GmMG2fP1AFobIQFC4KNR0QkV0r4aZSVwY03uvWHHoK9e4OLR0Qk\nV0r4bZg8GUaOtPLhwzBvXrDxiIjkQgm/DWVlcPPNbn3JEi2BKCKdlxJ+BuedB+efb+WWFpgzJ9h4\nRETay0/CnwTUAw3A7DRt5jjHVwNjnX3dgTpgFbZa1o9zijRA3/iG+xTNF1+ElSuDjUdEpD0yJfwy\nYC6W9Edj69nWJLWZAgwHRgDTgbud/YeBi4FzgDFOeWJeou5gI0bA5Ze79YceCi4WEZH2ypTwx2Nr\n1W4FmoGFwNSkNlcC851yHVAJDHTqB53XbtiXx7u5hRucL3/ZLf/lL7BtW2ChiIi0S6aEPxjY7qnv\ncPZlajPEKZdhXTp7gOewrp1OqbraHrkA9mydhx8ONBwRkayVZzju97FhyQvpJt73Adal0w9YBkSB\nWPKba2trPyxHo1Gi0ajPj+1Y114LL7xg5SeesHn6/foFG5OIlIZYLEYsFsvpHJlWPJ8A1GJ9+AC3\nAi3A7Z4292BJfKFTrwcuwq7qvf4vcAj4WdL+eLyTPI4yHofrroONG61+003wla8EG5OIlKaIzSTJ\nlMNbydSlsxwbjK3G+uGvBpYktVkCXO+UJwCNWLIfgPXnA/QALgM69fyWSMSu8hMWLYLm5uDiERHJ\nRqaEfwyYhXXHrAcWARuAGc4G8DSwBRvcnQfMdPafAvwZ68OvA54Ans1j7IG47DIYMMDKb7+tdW9F\npPPI6teBAuk0XToJDzwAd91l5UGD4LHHoHv3YGMSkdJSiC4dSeHzn4dKp7Nq9264//5g4xER8UMJ\nvx369LG7bxMefFCLnYtI+Cnht9MVV8CYMVY+dgzuuENr34pIuCnht1OXLjB7tr0C1NXBs51+SFpE\nipkSfg5GjYJ/+ie3/vOfw/79wcUjItIWJfwcff3r0L+/lffuhV/+Mth4RETSUcLPUe/e8K1vufXf\n/hZefTW4eERE0lHCz4NLL4WLLnLrP/whHDkSXDwiIqko4edBJGIDuL17W33bNrj33mBjEhFJpoSf\nJyef3Hr924cegvr64OIREUmmhJ9Hn/kMjBtn5Q8+gP/4D5ujLyISBkr4eRSJwHe/CxUVVt+4Ucsh\nikh4KOHnWVUVzJjh1u+9V49dEJFwUMIvgGuvhRpnqfejR+EHP4CWlmBjEhFRwi+AsjLr2ikrs/rK\nlbB4cbAxiYgo4RfIqFHwpS+59blz4b33gotHRMRvwp+ErVXbAMxO02aOc3w1MNbZVwU8B6wD1gLf\nSP3W4vTP/wxDhlh53z5L+iIiQfGT8MuAuVjSHw1MA2qS2kwBhmPr304H7nb2NwPfBM7E1ru9KcV7\ni1ZFBXz7227997+HNWuCi0dESpufhD8eW692K5bAFwJTk9pcCcx3ynXY4uUDgd3YmrYA+7H1cE/N\nKeJO5oILWj924Y47bI6+iEhH85PwBwPbPfUdzr5MbYYktanGunrqsgux8/uXf3Hn5tfXw+9+F2w8\nIlKayn208buOU/Jiut739QYeA27GrvRbqa2t/bAcjUaJRqM+P7JzOPVUuOEGuNvp6PrlL2HECDjn\nnGDjEpHOIxaLEYvFcjqHnxXPJwC1WB8+wK1AC3C7p809QAzr7gEb4L0I2AN0BZ4ElgJ3pjh/PF4C\nawMeOQLXXAPbnd+DevSAX/wCzj032LhEpHOKRCLgL4d/yE+XznJsMLYa6AZcDSxJarMEuN4pTwAa\nsWQfAe4D1pM62ZeMigpbESuxWMqhQ7YQ+vLlwcYlIqXDT8I/BswClmGJexE2+DrD2QCeBrZgg7vz\ngJnO/o8D1wEXAyudLfGbQsk5/XSYNw8GDLD64cNwyy12Y5aISKFl9etAgZREl47Xtm1w443w979b\nvU8f+NWvYNiwYOMSkc6jPV06SvgB2b4dvvpVeOcdqw8cCPffb68iIpkUqg9fCqCqCubMgZ49rb5n\nj/Xpv/9+sHGJSPFSwg/QqFHw059CuTM5dvNmuzNXi6aISCEo4QfsvPPg+99368uXw89+Flw8IlK8\nlPBDYPJk+PrX3fpjj8GjjwYXj4gUJyX8kLjhBvjUp9z6z34GL78cXDwiUnw0SydEjhyBr30N1q+3\neu/e9jiGmpJ5vqiI+KVpmUVg7164/np7BejbF+65B0aODDYuEQkXTcssAiedZM/Y6dvX6u+/DzNn\n2gweEZFcKOGH0MiRcNdd1qUD0Nhog7pK+iKSCyX8kKqpaZ30330Xpk93+/dFRLKlhB9iZ55pd+P2\n6mX1pia70l+xIti4RKRzUsIPuTFjbKZOZaXVDxywRzDkuA6CiJQgzdLpJDZvhptugrffdvddfbUl\n/8TyiSJSOjQts8ht325Jf9cud9+IEfDDH9qz9kWkdGhaZpGrqoIFC+Dii919DQ02b//xx0HfmyLS\nFr8JfxK2Tm0DMDtNmznO8dXAWM/++7HlDl9rZ4zi0a8f3HEH3Hqr25Vz+DD84Afw7/8O+49bIl5E\nxPhJ+GXAXCzpjwamAck3+08BhmNr304H7vYce4ASXtawECIR+Nzn4MEHW6+S9cwz8MUvwquvBheb\niISXn4Q/HlurdivQDCwEpia1uRKY75TrgEpgkFN/Hngv10DleMOGwfz5cNVV7r5du2DGDHvO/sGD\nwcUmIuHjJ+EPBrZ76jucfdm2kQLo3h2+8x34yU9sbdyERYtg2jQ9cVNEXOU+2vgdCkweLfY9hFhb\nW/thORqNEo1G/b5VHJdeanP2f/QjeOEF27dzpz2H54or4JZb3Ln8ItL5xGIxYjnegONnSs8EoBa3\nH/5WoAW43dPmHiCGdfeADfBehA3WAlQDTwAfTXF+TcvMo3gcli61Lp19+9z9lZXwzW/ClCk2BiAi\nnVuhpmUuxwZjq4FuwNXAkqQ2S4DrnfIEoBE32UsHikQsqT/6KFx2mbu/sdGWUvzqV6G+Prj4RCQ4\nfr8dJgN3YjN27gN+DMxwjs1zXhMzeQ4AXwEST3x5GLva7w/8HfgeNnMnQVf4BfT883D77bB7t7sv\nEoHPfAZuvBH69w8uNhFpP91pKykdPAj33Qe//jUcO+bu79kTvvQluPZaG/wVkc5DCV/a9Oab8J//\nCS++2Hr/SSfZjJ7LL9cVv0hnoYQvvrz4oq2qtWVL6/1lZTBxos3rv+ACDe6KhJkSvvj2wQewZAnM\nm9f6CZwJw4fDl79sA79lZR0enohkoIQvWTt0CP74R0v+q1Ydf/yUU+DTn7bunsG6lU4kNJTwJSdv\nvgm//S0sXmxfBMnGjoWpU+0mLw3yigRLCV/yoqkJHnnEHs/Q2Hj88d69ba7/pz4FZ5yhBVhEgqCE\nL3l15Aj85S/w1FPwt79Zv3+y8nIYNQrOPtu6fUaN6vg4RUqREr4UzDvvWOJfvNhW3krnzDPhs5+F\naFTP7hEpJCV8KbiWFlixAp5+Glavtn7/dE4/3fr9P/YxGD/eFm8RkfxQwpcO19QEa9fCH/4Azz4L\nR4+mbheJ2NX/+PH2WlNjN3xprr9I+yjhS6Camqzb59lnYd261o9xSKV/f+vzHznStlGjbN3eLlpp\nWSQjJXwJjcOH7cr/1Vehrs7KLS2Z39e7N4webb8BDBtm3ULV1ZoGKpJMCV9C6/33Yfly6/ffsMEe\n0ex3CcZIBIYMgREjbBs61G4IO+UUOOEEdQtJaVLCl06jpcVm+7z+um0bN8L69ann/bele3e7A3jI\nENuqquAjH4HTToOTT9aXgRQvJXzp1OJxeOst6/9vaLCHu73xBuzYkfoegEx69LDkP3So+wUwYICN\nHZxwgm3duuX/zyHSEQqV8CfhLn7yK1ovbZgwB1sk5SDwZWBlFu9Vwpc2HTliif/112HTJvsCeOst\n2/bvz+3cvXpZ4u/Xz9369rWxhD59Wm99+9q9Bf36aUxBgleIhF8GbAQuBXYCrwDTgA2eNlOAWc7r\necAvsGUO/bwXQprwY7FY6BZTV0zH27fPvgB27LAuom3boK4uxtGjUZqaCve5FRX2pdCzp31p9OjR\neuvd2/YnXjdtinH++VEqKuy3iq5d7UujRw87R48etq+jZygF/feXimLypz0JvzzD8fHAJmCrU18I\nTKV10r4SmO+U64BKYBAw1Md7QyuMf8GK6Xh9+tiMnpoad19tbYza2ijvvQdbt9pvBzt32mOg337b\n7hp+7z0bL2hPVxHYbx1Hjti5/Ni1K8Yjj0Qztisvt8RfUWFfCImta9fjt8QXR7du7ubdX15uW1mZ\n+5rYEscWLIjRo0f0w3qXLraVlbmvyceSt0Rb7xaJuK/ZCvrfVCphjKk9MiX8wYD3Rvod2FV8pjaD\ngVN9vFekYBL99GPHpj7e0mK/ITQ22tbUZNv+/bY/eXv/fTve2AjNzYWJ+dgx21I9rbQQdu2ymVOF\nVFbmfgEkfxl4y4njmzbZVN7E8fZs0Pa+RDn5Nbl9Yv+rr9rPKtUXWPK+5PemapOunZ9juciU8P32\ntYRh8FckK126uP32H/mI//fF43afwYED7nbokO07dMimm+7f724HDsCyZTBunLVpbrY7ko8csfaJ\n9xTqSyRo2f4WdfAg7N1bmFjaa/due6RIsZsA/MFTvxWYndTmHuAaT70eGOjzvWDdPnFt2rRp05bV\ntok8Kwc2A9VAN2AVUJPUZgrwtFOeALyUxXtFRCREJmOzbTZhV+kAM5wtYa5zfDVwbob3ioiIiIhI\nsZqE9fk3kLp/vyPcD+wBXvPsOxH4I/A68Aw21bQjVQHPAeuAtcA3QhBXd2za7SpgPfDjEMSUUIbd\n7PdEiGLaCqxx4no5JHFVAo9hU6PXY7PmgoxpFPbzSWxN2L/1oH9Ot2L/914DfgNUhCCmm5141jpl\nQhBTVsqwrp5qoCvB9fFfCIyldcK/A/i2U54N/KSDYxoEnOOUe2PdYjUhiKun81qOjdVMDEFMAP8H\n+DWwxKmHIaY3sP+QXkHHNR+4wSmXA/1CEFNCF+At7GInyJiqgS1YkgdYBHwp4JjOwvJTdyxv/hEY\nFnBMWTuf1rN4/s3ZglBN64SfmGkElnzrOzqgJI9jdyyHJa6e2J3TZ4YgpiHAn4CLca/wg44JLOH3\nT9oXZFz9sESWLAw/K4BPAs875SBjOhG7wDoB+1J8Args4Jg+jz2aJuG7WKIPy9+dL58H/ttTvw74\nZUCxVNNhozfZAAACXklEQVQ64b/nKUeS6h2tGngT6EPwcXXBfhPbh11dEIKYHsV+Q7sIN+EHHRNY\ncl0JLAe+5uwLMq5zsC65B4AV2P+9XgHH5HU/MNMpBx3TdOzf+N+BBSGI6QzsS+hE7GLrr9jzy7KO\nKci1heIBfnY2EnNeg9Ab+C3WZ7cv6VgQcbVgiWMI8I/YVXWQMV2B/adcSfqb/4L6+/s49kU0GbgJ\n6zr06ui4yrEZdP/PeT3A8b9RB/Wz6gZ8GvvyTtbRMQ0DbsEutE7F/g9eF3BM9diDJ58BlmIXXcm3\ns/mKKciEvxPrr0uowh6/EAZ7sF+RAE7BkkpH64ol+wVYl05Y4gIbXHsKGBdwTBdgz3J6A3gY+AT2\n8wrDz+kt53UvsBh7LlWQce1wtlec+mNY4t8dYEwJk4FXsZ8VBPtz+hh2Bf0OcAz4Hdb9HPTP6X4n\ntouwK/nXacfPKciEvxwYgXtj1tW4g25BW4IN1OC8Pt5G20KIAPdhMynuDElcA3BnAfTA+jVXBhzT\nd7ALhaHY3d5/Bv5XwDGB/drdxyn3wvqnXws4rt3Ys61GOvVLsZkoTwQYU8I07As7IcifUz12A2kP\n7P/hpdj/w6B/Tic7r6cBV2Gzh4L+d561MNyY9TCwCziK/Yf4CtZX9ieCm+40Ees+WYU7ZW1SwHF9\nFOv7XYVNN/xXZ3/QP6uEi3AvGIKOaSj2c1qFTaNL/NsOOq6zsSv81diVa78QxNQLeBv3C5IQxPRt\n3GmZ87HftoOO6S9OTKtwu1KDjklEREREREREREREREREREREREREREREREREpPP7/45PXlSVMNWf\nAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 51 }, { "cell_type": "markdown", "metadata": {}, "source": [] }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 51 }, { "cell_type": "markdown", "metadata": {}, "source": [] } ], "metadata": {} } ] }