{ "cells": [ { "cell_type": "markdown", "metadata": { "tags": [ "s1", "content", "l1" ] }, "source": [ "### BINOMIAL DISTRIBUTION - COIN TOSS EXPERIMENT\n", "\n", "\n", "Let us consider an experiment of tossing a biased coin, in which p is the probability of 'heads' and n is the number of tosses. An experiment is defined as a series of trials. In this case, an experiment will consist of n trials/tosses. We will repeat the experiment 1000 times to simulate the binomial distribution. \n", "\n", "We must check if the experiment follows the assumptions of binomial distribution:\n", "\n", "1. Each trial is an independent of the other, meaning outcome of one toss does not affect that of others.\n", "2. There are only two possible outcomes for each trial, heads or tails.\n", "3. The probability of 'success', p is the same across the n trials.\n", "3. The number of trials, n is fixed.\n", "\n", "Consider n tosses of a coin. Since the tosses are independent, we can calculate the probability of an outcome by multiplying the individual probabilities in each toss. Probability of 'heads' is p and that of 'tails' is 1-p. Outcome of one experiment with n tosses, has k heads and n-k tails and has probability **pk(1-p)3-k**.\n", "There are ${n} \\choose {k}$ number of distinct n-toss sequences that contain k heads. This forms the binomial coefficient.\n", "\n", "Binomial probability mass function (pmf) = ${n} \\choose {k}$ **pk(1-p)n-k**\n", "\n", "We can now go ahead and simulate our experiment." ] }, { "cell_type": "code", "execution_count": 35, "metadata": { "collapsed": true, "tags": [ "s1", "ce", "l1" ] }, "outputs": [], "source": [ "%matplotlib inline\n", "import random\n", "from scipy.stats import binom\n", "import matplotlib.pyplot as plt\n", "\n", "def toss(p,n):\n", " heads = 0\n", " for i in range(n):\n", " if random.random() < p:\n", " heads += 1\n", " return heads #This gives the number of heads in a single experiment of n trials\n", "\n", "size = range(1,1000)\n", "n,p = 30,0.6\n", "\n", "k = [] #Holds the number of heads from each experiment\n", "pmf= [] #Holds the binomial pmf of each experiment\n", "\n", "random.seed(12345)\n" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true, "tags": [ "s1", "l1", "hint" ] }, "source": [ "Use toss function to create list of k values, sort them and for each k, calculate pmf using binom.pmf function. Plot k and pmf to get the distribution plot. " ] }, { "cell_type": "code", "execution_count": 36, "metadata": { "tags": [ "s1", "l1", "ans" ] }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 36, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VOXd//H3NxsJJCGEhC17IEAA2QxLZHFXcMPWVkEE\nd0qrVWutD63tc/XXx8dat1ZbqqKigijuSiuKWlFRtoRVwhpCVghJCNnJMsn9+2OGPmlMyABJzpnM\n93VduUjOnJN8MgwfzpzlvsUYg1JKKe/hY3UApZRSXUuLXymlvIwWv1JKeRktfqWU8jJa/Eop5WW0\n+JVSysto8SullJfR4ldKKS+jxa+UUl7Gz+oArYmIiDDx8fFWx1BKKY+xZcuWEmNMpDvr2rL44+Pj\nSU9PtzqGUkp5DBHJcXddPdSjlFJeRotfKaW8jBa/Ukp5GS1+pZTyMlr8SinlZbT4lVLKy2jxK6WU\nl9HiV6oDfZpRyNf7i62OodQp2fIGLqU8TVOT4bE1+3juq4MAXDgskt9dNYLEyGCLkyn1fbrHr9RZ\nOlHfyF2vb+W5rw5y46RYHroimbTs41z+l6/54+o9VNY2WB1Rqf+ge/xKnYWiylruXLaFnfll/PbK\nZG6fmoCIcO24KB5fs5fnv87i3a0FLJo5nB+Oi8LHR6yOrJTu8St1pvYVVvKDxevZX1jJ8zedyx3T\nEhFxFntkSA8e+9EYPrxrCjHhQTzw9g5+8Ox6tueVWZxaKS1+pc7Il/uKuO7Z9Tiamnh7YSqXjRzQ\n6npjYsJ4d+F5PPnjMRwuO8G1i7/lgbd3UFRZ28WJlfo/WvxKnablG7K57ZU0YsJ78sFdUxgV1fuU\n6/v4CNedG83aBy5g4fmD+XB7ARc98RXPf3WQekdT14RWqhkxxlid4XtSUlKMDsus7KaxyfC/H+1h\n6beHuHh4P56ZM45ePU7/NNmhkmoe/udu/rW3iMSIXvzuqhFcOLxfJyRW3kREthhjUtxZ1609fhGZ\nISL7RCRTRBa18vhwEdkgInUi8kArj/uKyDYR+ac7P08pu6muc/CT5eks/fYQt01JYMn8lDMqfYCE\niF68dMsEXr51AgC3vpLGba+kcaikuiMjK9Wmdl+5IuILLAYuBfKBNBFZZYzZ3Wy1UuAe4No2vs29\nwB4g9OziKtX1jpSf4PZX0tlbWMH/zBrJvNT4Dvm+Fw7rx5TBEby6Ppun/3WAy/78FbdNSeDui4YQ\nEujfIT9Dqda4s8c/Ecg0xmQZY+qBlcCs5isYY4qMMWnA9y5YFpFo4ErgxQ7Iq1SX2lVQzrWLvyW3\ntIalt0zosNI/KcDPhzunJ/LFA+dz7dgonv86i4ue/Ip3tuTT1GS/w7Cqe3Cn+KOAvGZf57uWuesv\nwIPAKc9iicgCEUkXkfTiYr3lXVnvs91H+fFzG/Dz8eGdn6ZywbDOOw7fLySQx3/svPwzuo/z8s//\nXb2n036e8m6delWPiFwFFBljtrS3rjFmiTEmxRiTEhnp1nzBSnUKYwwvrstiwfJ0hvYP5v27zmP4\ngK45Snny8s/rU6J5dX22HvdXncKd4i8AYpp9He1a5o4pwDUiko3zENFFIvLaaSVUqgs5Gpv47Qe7\nePijPcwYOYCVC1LpFxLYpRl8fIRfXT6cAD8fnlizr0t/tvIO7hR/GpAkIgkiEgDMBla5882NMb82\nxkQbY+Jd231hjLnpjNMq1Ykqahu49ZU0VmzKZeH5g1l843iCAnwtyRIZ0oMF0xP56LsjbMs9bkkG\n1X21W/zGGAdwN7AG55U5bxljMkRkoYgsBBCRASKSD9wP/FZE8kVEr+BRHiP/eA0/enY9Gw4e40/X\nncOimcMtH1fnzmmJRAT34I8f78WO99soz+XWhcjGmNXA6hbLnmv2eSHOQ0Cn+h5fAl+edkKlOllT\nk2HBsi0cKa/l1dsmMmVIhNWRAOjVw497L0nidx/s4ou9RVyc3N/qSKqb0CEblNf7JKOQ3Ucq+MOs\nkbYp/ZNmT4ghMaIXj368F0ejDu+gOoYWv/JqjU2Gpz7bz5B+wVwz5nSuUu4a/r4+PDhjGAeKqnh3\na77VcVQ3ocWvvNo/dhwms6iKX1wyFF+bjpV/+cgBjI8N46nP9nOivtHqOKob0OJXXsvR2MTT/zrA\n8AEhzBzV+rDKdiAi/PqKZI5W1LH020NWx1HdgBa/8lrvbSvgUEk191861PIreNozIT6cS5L789yX\nBymtrrc6jvJwWvzKK9U7mnj68wOMju7NpSM842qZRTOHUV3v4K9fHLA6ivJwWvzKK72VnkdB2Qnu\nv3Tov6dLtLsh/UK4YUIMr23MIfdYjdVxlAfT4ldep7ahkb99kcm5cX04f6hnjQt1n+sk9BOf6lAO\n6sxp8Suv8/qmXAoravnlZZ6zt39S/9BA7piayKodh9mZrxO3qzOjxa+8yon6Rv7+5UFSE/ty3mB7\n3azlrp+cn0h4rwAe1aEc1BnS4ldeZdmGbEqq6vjlZUOtjnLGQgL9ueeiIaw/eIyv9uvcFer0afEr\nr1FV5+C5rw4yfWgkKfHhVsc5KzdOiiOub08e/XgvjTpTlzpNWvzKa7z8zSGO1zTwy0s9d2//pAA/\nHx64bBh7Cyt5f5u702Mo5aTFr7xCeU0DS9ZlcUlyf8bEhFkdp0Ncec5ARkf35qlP91HboEM5KPdp\n8Suv8OI3WVTWOri/G+ztn+TjIyyaOZzD5bW8uj7b6jjKg2jxq26vtLqepd8c4spzBjJiUPeaH+i8\nwRFcOCySxWszKavRoRyUe7T4Vbf3/NcHqWlo5L5LkqyO0in+a+ZwKuscLF6baXUU5SHcKn4RmSEi\n+0QkU0QWtfL4cBHZICJ1IvJAs+UxIrJWRHaLSIaI3NuR4ZVqT1Gl8zDIrDGDSOofYnWcTjF8QCjX\njY/m1fU55JXqUA6qfe0Wv4j4AouBmcAIYI6IjGixWilwD/BEi+UO4JfGmBHAZOCuVrZVqtM8++VB\nGhoN917SfY7tt8Y55hA89dl+q6MoD+DOHv9EINMYk2WMqQdWArOar2CMKTLGpAENLZYfMcZsdX1e\niXOydvtNc6S6pSPlJ1ixKZfrxkeRENHL6jidalBYELdOSeCD7QVkHC63Oo6yOXeKPwrIa/Z1PmdQ\n3iISD4wDNp3utkqdicVrMzHG8POLuuex/ZZ+esFgegf58+jHe62OomyuS07uikgw8C5wnzGmoo11\nFohIuoikFxfrbejq7OSV1vBmWh7Xp8QQE97T6jhdoneQP3dfOIR1B0pYd0D/Dam2uVP8BUBMs6+j\nXcvcIiL+OEt/hTHmvbbWM8YsMcakGGNSIiM9a6hcZT9//eIAIsLdFw2xOkqXmpcaR3SfIB79eC9N\nOpSDaoM7xZ8GJIlIgogEALOBVe58c3GOefsSsMcY89SZx1TKfYdKqnl3awFzJ8UysHeQ1XG6VA8/\nXx64bBgZhytYteOw1XGUTbVb/MYYB3A3sAbnydm3jDEZIrJQRBYCiMgAEckH7gd+KyL5IhIKTAHm\nAReJyHbXxxWd9tsoBTz9+X78fYWfXjDY6iiWuGbMIEYOCuXxNfuoc+hQDur7/NxZyRizGljdYtlz\nzT4vxHkIqKVvAM+a6UJ5tANHK/lwx2EWTEukX0ig1XEscXIoh3kvbWb5hhzumJZodSRlM3rnrupW\n/vL5AXr6+/KT871zb/+kaUmRTEuK4G9rMyk/0dD+BsqraPGrbiPjcDkffXeE26YmEN4rwOo4lls0\nczjlJxp4+dtDVkdRNqPFr7qNP392gJBAP+6Yqoc2AEYO6s30pEhWbs7D0dhkdRxlI1r8qlvYkVfG\n53uOcue0RHr39Lc6jm3cOCmWwopavthbZHUUZSNa/KpbeOqz/fTp6c+tU+KtjmIrFw/vR//QHry+\nOdfqKMpGtPiVx0vPLuWr/cX85PzBhATq3n5zfr4+3DAhlq/2F+vInerftPiVx3vy0/1EBAcwPzXO\n6ii2NHtCDAKsTNO9fuWkxa882vrMEjZkHeNnFwyhZ4Bbt6V4nUFhQVw0vB9vpuXToCd5FVr8ysP9\n/cuD9A/twY2TYq2OYms3ToqlpKqOz3YftTqKsgEtfuWxDhZX8U1mCTdNiiPQ39fqOLZ2/tB+RIUF\n8fomPdyjtPiVB1u+IQd/X2H2RN3bb4+vjzB7QgzfZJaQXVJtdRxlMS1+5ZGq6xy8uyWfmaMGEhnS\nw+o4HuH6CTH4+ghv6KWdXk+LX3mkD7YXUFnn0Ct5TkP/0EAuTe7P21vyddROL6fFrzyOMYblG3JI\nHhjKuXF9rI7jUW6cFEtpdT2f7Cq0OoqykBa/8jhp2cfZW1jJ/NQ4nHP9KHdNHRJBbHhPPcnr5bT4\nlcdZtiGbkEA/Zo0dZHUUj+PjI8yZGMumQ6VkFlVZHUdZRItfeZSiilo+2VXIj8+N0Ru2ztCPU6Lx\n9xXd6/diWvzKo7yxOQ9Hk2GentQ9YxHBPbhs5ADe3ZpPbYOe5PVGbhW/iMwQkX0ikikii1p5fLiI\nbBCROhF54HS2VcpdDY1NvL45h2lJESRE9LI6jkebOymW8hMNfLTziNVRlAXaLX4R8QUWAzOBEcAc\nERnRYrVS4B7giTPYVim3fLb7KEcr6pifGm91FI+XmtiXxIheOlyzl3Jnj38ikGmMyTLG1AMrgVnN\nVzDGFBlj0oCWk3u2u61S7lq2IZso14Bj6uyICDdOimVLznH2FlZYHUd1MXeKPwrIa/Z1vmuZO85m\nW6X+7cDRSjZmlTJ3ciy+PnoJZ0e4bnw0AX4+epLXC9nm5K6ILBCRdBFJLy4utjqOspnlG3MI8PXh\nhpQYq6N0G316BXDFqAG8v7WAmnqH1XFUF3Kn+AuA5v/aol3L3OH2tsaYJcaYFGNMSmRkpJvfXnmD\nqjoH720t4KrRA+kbrOPydKS5k+OorHPwjx2HrY6iupA7xZ8GJIlIgogEALOBVW5+/7PZVikA3t+a\nT1WdQy/h7AQpcX1I6hesh3u8TLvFb4xxAHcDa4A9wFvGmAwRWSgiCwFEZICI5AP3A78VkXwRCW1r\n2876ZVT3Y4xh2YYczonqzdiYMKvjdDsiwtxJsezIL2dXQbnVcVQXcesYvzFmtTFmqDFmsDHmf13L\nnjPGPOf6vNAYE22MCTXGhLk+r2hrW6XctTGrlANFVczTcXk6zQ/GRxPo78MK3ev3GrY5uatUa5Zv\nzCaspz/XjNFxeTpL7yB/rho9iFXbC6iq05O83kCLX9lWYXktazKOcn1KjE6t2MnmToqlur6RD7e7\ne92G8mRa/Mq2Xt+cS5Mx3DRJT+p2trExYSQPDGXFxlyMMVbHUZ1Mi1/ZUr2jiTc253LB0Ehi+/a0\nOk63d/Ik7+4jFezI15O83Z0Wv7KlNRmFFFfquDxdadbYQfQM8GXFxhyro6hOpsWvbGn5hhxiwoOY\nPlRv5usqIYH+zBo7iH/sPEz5iZbDbqnuRItf2c7ewgo2Z5dy06Q4HZeni904MY7ahiY+2KYnebsz\nLX5lO8s25NDDz4frdVyeLndOdG9GR/dmxaYcPcnbjWnxK1upqG3gg20FXD1mEH16BVgdxyvNnRTL\n/qNVbMk5bnUU1Um0+JWtvLsln5r6RubruDyWuXrMIEJ6+OmdvN2YFr+yDWMMyzfmMCYmjNHROi6P\nVXoG+HHtuCg++u4Ix6vrrY6jOoEWv7KNbzOPkVVczfzJurdvtRsnxVLvaOLdrflWR1GdQItf2cay\nDdmE9wrgytEDrY7i9ZIHhjI+NozXN+udvN2RFr+yhYKyE3y+R8flsZO5k+LIKq5mY1ap1VFUB9Pi\nV7bwxqZcDM4rSpQ9XDl6IL2D/FmxSe/k7W60+JXl6hyNrEzL5eLh/YgJ13F57CLQ35cfjo9iTUYh\nJVV1VsdRHUiLX1nuk12FlFTVM0/H5bGduZNiaWg0vLNFT/J2J1r8ynLLNuQQ37cn04ZEWB1FtTCk\nXwgTE8J5fVMuTU16kre7cKv4RWSGiOwTkUwRWdTK4yIiz7ge3yki45s99gsRyRCRXSLyhogEduQv\noDxbxuFytuQc56bJcfjouDy2NHdSLLmlNXyTWWJ1FNVB2i1+EfEFFgMzgRHAHBEZ0WK1mUCS62MB\n8Kxr2yjgHiDFGDMK8AVmd1h65fGWb8gh0N+HH5+r4/LY1YxRA+jbK4DlOlxzt+HOHv9EINMYk2WM\nqQdWArNarDMLWGacNgJhInLyYmw/IEhE/ICewOEOyq48XHlNAx9sL+DasVH07ulvdRzVhh5+vtww\nIYZ/7TlKQdkJq+OoDuBO8UcBec2+zncta3cdY0wB8ASQCxwByo0xn555XNWdvL0lj9qGJubpuDy2\nN9d1N7VO0tI9dOrJXRHpg/PdQAIwCOglIje1se4CEUkXkfTi4uLOjKVsoKnJ8NrGHM6N68PIQb2t\njqPaERUWxMXJ/XkzLY86R6PVcdRZcqf4C4DmB2CjXcvcWecS4JAxptgY0wC8B5zX2g8xxiwxxqQY\nY1IiI3XWpe5uXWYJ2cdqdBRODzI/NY5j1fWs/u6I1VHUWXKn+NOAJBFJEJEAnCdnV7VYZxUw33V1\nz2Sch3SO4DzEM1lEeoqIABcDezowv/JQL397iIjgAGaMGmB1FOWmKYMjSIzoxbINerjH07Vb/MYY\nB3A3sAZnab9ljMkQkYUistC12mogC8gEXgB+5tp2E/AOsBX4zvXzlnT0L6E8y77CSr7cV8zNqfH0\n8NNxeTyFj49w0+Q4tuWWsaug3Oo46iyIHUfeS0lJMenp6VbHUJ3kV2/v4B87D7Nh0cU6y5aHKT/R\nwORH/sXVYwby2I/GWB1HNSMiW4wxKe6sq3fuqi5VVFHLB9sLuD4lRkvfA/UO8ufacVF8uP0wZTU6\nSYun0uJXXeqV9dk4mgy3T02wOoo6Q/NT46hzNPF2uo7f46m0+FWXqa5zsGJTLjNGDiCuby+r46gz\nlDwwlAnxfXhtU46O3+OhtPhVl3k7PY/yEw3cMS3R6ijqLM1LjSfnWA1fHdB7bjyRFr/qEo7GJl76\n9hDnxvXh3Lg+VsdRZ2nGyAFEBPdguV7a6ZG0+FWXWJNxlLzSE9ype/vdQoCfDzdOjGHtviLySmus\njqNOkxa/6nTGGJZ8fZD4vj25dER/q+OoDnLjpDh8RHhNx+/xOFr8qtOlZR9nR345t09LxFfH3O82\nBvQO5LIR/XkzPY/aBh2/x5No8atO98K6LPr09OdH46OtjqI62LzUOMpqGvjHDh1t3ZNo8atOdbC4\nis/3HGVeajxBATo8Q3eTmtiXpH7BOkmLh9HiV53qpW8O4e/ro6NwdlMiwrzUOHbml7M9r8zqOMpN\nWvyq05RU1fHulnyuGx9FRHAPq+OoTvKDcVH0CvBl2YZsq6MoN2nxq06zfEMOdY4mbp+ql3B2ZyGB\n/vxwfDT/3HmE0modv8cTaPGrTlHb0MjyjTlcktyPIf2CrY6jOtm81DjqHU28mZbX/srKclr8qlO8\nuzWf0up6vWHLSwztH8LkxHBe25hDo47fY3ta/KrDNTUZXlx3iDHRvZmYEG51HNVF5qfGU1B2grV7\ni6yOotqhxa863Od7jnKopJo7piXinHFTeYNLR/Snf2gPlumlnbanxa863AvrsogKC2KmzqfrVfx9\nfZg7KY6v9xeTVVxldRx1Cm4Vv4jMEJF9IpIpIotaeVxE5BnX4ztFZHyzx8JE5B0R2Ssie0QktSN/\nAWUvW3OPk5Z9nNunJuDnq/sV3mb2xBj8fYXXNuZaHUWdQrv/MkXEF1gMzARGAHNEZESL1WYCSa6P\nBcCzzR57GvjEGDMcGINzwnbVTb24LovQQD+unxBjdRRlgX4hgcwYNZC3t+RRU++wOo5qgzu7ZBOB\nTGNMljGmHlgJzGqxzixgmXHaCISJyEAR6Q1MB14CMMbUG2P09r5uKvdYDZ/sKmTu5DiCe/hZHUdZ\nZH5qHJW1Dj7cruP32JU7xR8FNL84N9+1zJ11EoBi4GUR2SYiL4qIzrnXTS399hC+PsIt58VbHUVZ\nKCWuD8MHhLBsQw7G6KWddtTZB2H9gPHAs8aYcUA18L1zBAAiskBE0kUkvbhYp3PzNGU19byZlsc1\nY6LoHxpodRxlIRFhfmo8e45UsCXnuNVxVCvcKf4CoPkB22jXMnfWyQfyjTGbXMvfwfkfwfcYY5YY\nY1KMMSmRkZHuZFc2smJTLicaGrlzeoLVUZQNXDtuECGBfizTqRltyZ3iTwOSRCRBRAKA2cCqFuus\nAua7ru6ZDJQbY44YYwqBPBEZ5lrvYmB3R4VX9lDnaOTlb7OZPjSS4QNCrY6jbKBngB8/Ojeaj3cd\nobiyzuo4qoV2i98Y4wDuBtbgvCLnLWNMhogsFJGFrtVWA1lAJvAC8LNm3+LnwAoR2QmMBR7pwPzK\nBj7cdpiSqjoW6PAMqpl5k+NoaDSs3KyXdtqNW5deGGNW4yz35suea/a5Ae5qY9vtQMpZZFQ2Zozh\nhXVZJA8MZcqQvlbHUTaSGBnMtKQIXt+cy08vGKz3ddiI/k2os/Ll/mIOFFWxYHqCDs+gvmfe5DiO\nlNfy+Z6jVkdRzWjxq7PywtdZDAgN5KrRg6yOomzo4uT+RIUF6Ulem9HiV2dsV0E56w8e49Yp8fjr\n23jVCl8f4cZJsaw/eIzMokqr4ygX/deqztgL67II7uHHnEmxVkdRNjZ7QgwBvj4s171+29DiV2ek\noOwE/9x5hNkTYggN9Lc6jrKxvsE9uHL0QN7dWkBVnY7fYwda/OqMvPzNIQBunao3bKn2zUuNo6rO\nwfvbWt77qaygxa9OW0VtAyvT8rhq9ECiwoKsjqM8wLiYMEZFhbJ8Q7aO32MDWvzqtK3cnEtVnUPn\n01VuExHmT45n/9EqvskssTqO19PiV6el3tHE0m+ySU3sy6io3lbHUR7kmrGDiAoL4k+f7KVJJ2S3\nlBa/Oi0fbC+gsKKWBdN1b1+dnkB/Xx64fCi7CipYtUPH6reSFr9yW/mJBh77ZB9jY8I4f6iOoKpO\n36wxUYyKCuXxNfuobWi0Oo7X0uJXbvvzZ/s5Vl3Hw9eOwsdHh2dQp8/HR/jNzGQKyk7w6vpsq+N4\nLS1+5ZbdhytYtiGbmybF6bF9dVbOGxLBhcMi+dvaTI5X11sdxytp8at2NTUZ/vvDXfTpGcADlw1r\nfwOl2vHrK5KprnPw1y8yrY7ilbT4Vbve21ZAes5xFs0cTu+eepeuOntD+4dwfUoMyzdmk3Os2uo4\nXkeLX51S+YkG/rh6D+fG9eG68dFWx1HdyP2XDsXPx4fH1uyzOorX0eJXp/Tkp/s4XlPPH2aN1BO6\nqkP1Cw3kzumJfLTzCNtydVL2rqTFr9q0q6Cc1zbmMD81npGD9ISu6ngLpicSEdyDP67eq0M5dCG3\nil9EZojIPhHJFJFFrTwuIvKM6/GdIjK+xeO+IrJNRP7ZUcFV52pqMvzuw12E9wrgF5cOtTqO6qaC\ne/hx3yVJbM4u5bPdOktXV2m3+EXEF1gMzARGAHNEZESL1WYCSa6PBcCzLR6/F+dE7cpDvLMln225\nZfx6ZjK9g/SEruo8syfEMDiyF49+speGxiar43gFd/b4JwKZxpgsY0w9sBKY1WKdWcAy47QRCBOR\ngQAiEg1cCbzYgblVJyqrqefRT/YyIb4PPxwfZXUc1c35+fqwaGYyWcXVrEzLszqOV3Cn+KOA5n8b\n+a5l7q7zF+BBQP8r9xBPfLqP8hMN/GHWKJ1AXXWJS5L7MTEhnKc/36+TtXSBTj25KyJXAUXGmC1u\nrLtARNJFJL24uLgzY6lT+C6/nBWbcrk5NZ7kgaFWx1FeQkR46IpkSqrqef6rg1bH6fbcKf4CIKbZ\n19GuZe6sMwW4RkSycR4iukhEXmvthxhjlhhjUowxKZGROgCYFZqaDL/9cBcRwT2479Ikq+MoLzMm\nJoyrxwzihXVZFJbXWh2nW3On+NOAJBFJEJEAYDawqsU6q4D5rqt7JgPlxpgjxphfG2OijTHxru2+\nMMbc1JG/gOo4b6bnsSOvjIeuSNZ5dJUlHrx8GE1N8NRnelNXZ2q3+I0xDuBuYA3OK3PeMsZkiMhC\nEVnoWm01kAVkAi8AP+ukvKqTHK+u50+f7GViQjizxg6yOo7yUjHhPZmfGsfbW/LZW1hhdZxuS+x4\n00RKSopJT0+3OoZX+fV73/FWeh6r75nGsAEhVsdRXqyspp7pj61lXGwfXr1totVxPIaIbDHGpLiz\nrt65q9ieV8bKtFxuPS9eS19ZLqxnAHdfNISv9hfzzQGdn7czaPF7uUbXkMv9Qnpwn96hq2xifmo8\nUWFBPLJ6j87P2wm0+L3cG5tz2ZlfzkNXjiC4h5/VcZQCnPPzPjhjGLuPVPDB9pYXEaqzpcXvxY5V\n1fH4mn2kJvbl6tEDrY6j1H+4evQgzonqzRM6P2+H0+L3Yo99so/qOgd/mDVS79BVtuPjI/zmimQO\nl9fy8rfZVsfpVrT4vdTW3OO8mZ7H7VMTSOqvJ3SVPaUO7svFw/vx97WZlOr8vB1Gi98LNTYZfvfB\nLgaEBvLzi/UOXWVvi2YOp7rewTP/OmB1lG5Di98Lvb4ph4zDFfz2qmQ9oatsL6l/CDdMiOW1jTlk\nl+j8vB1Bi9/LlLhO6E4dEsGV5+gJXeUZfnFpEgF+Pjy2Zq/VUboFLX4v8+jHeznR0Mjvr9ETuspz\n9AsJZMH0RFZ/V8iWHJ2f92xp8XuR9OxS3tmSzx3TEhnSL9jqOEqdljunJRIZ0oNHVu/R+XnPkha/\nl3A0NvG7DzMY1DuQn180xOo4Sp22Xj38+MUlQ9mSc5w1GYVWx/FoWvxewBjDwx/tYc+RCn571Qh6\nBugJXeWZrk+JZki/YP70yT7qHTqp35nS4vcCT366n1fWZ3PblARmjhpgdRylzpifrw8PXZnMoZJq\nHnxnh47jc4a0+Lu5xWsz+dvaTOZMjOV3VyXrCV3l8S4c1o9fXT6MD7Yf5k+f6FU+Z0Lf83djL397\niMfX7OPasYN4+FqdOF11Hz+7YDCF5bU8/3UW/UIDuX1qgtWRPIoWfzf1Zlou/+8fu7l8ZH+e+PEY\nfH209FX/ydsnAAANUUlEQVT3ISL8/pqRFFfW8T//3E1kSA+uGaMzx7lLD/V0Qx9uL2DRe98xfWgk\nz8wZh5+v/jWr7sfXR/jL7LFMTAjnl29tZ32mTtriLrcaQURmiMg+EckUkUWtPC4i8ozr8Z0iMt61\nPEZE1orIbhHJEJF7O/oXUP/p04xC7n9rBxPiw3n+pnPp4edrdSSlOk2gvy8vzEshIaIXC5ZvIeNw\nudWRPEK7xS8ivsBiYCYwApgjIiNarDYTSHJ9LACedS13AL80xowAJgN3tbKt6iDrDhRz9+vbOCeq\nN0tvmUBQgJa+6v569/Tn1dsmEhLoxy0vp5FXWmN1JNtzZ49/IpBpjMkyxtQDK4FZLdaZBSwzThuB\nMBEZaIw5YozZCmCMqQT2AFEdmF+5bD5Uyp3L0hncL5hXb52og68przKwdxDLbptIvaOJm5du1iGc\n2+FO8UcBec2+zuf75d3uOiISD4wDNp1uSHVqO/LKuO2VNAaFBbH89on07ulvdSSlulxS/xBevDmF\ngrIT3PZKGjX1Dqsj2VaXnPUTkWDgXeA+Y0xFG+ssEJF0EUkvLi7uiljdwp4jFcxfupk+vfx5/Y7J\nRAT3sDqSUpaZEB/OM3PGsTO/jLtf34ajUe/ubY07xV8AxDT7Otq1zK11RMQfZ+mvMMa819YPMcYs\nMcakGGNSIiMj3cnu9bKKq5j30iaC/H15/Y7JDOgdaHUkpSx3+cgB/GHWKL7YW8Rv3v9OB3RrhTvF\nnwYkiUiCiAQAs4FVLdZZBcx3Xd0zGSg3xhwR5x1DLwF7jDFPdWhyL5dXWsPcFzdhDLx2xyRiwnta\nHUkp27hpchz3XJzEW+n5PPXZfqvj2E67ZwCNMQ4RuRtYA/gCS40xGSKy0PX4c8Bq4AogE6gBbnVt\nPgWYB3wnIttdy35jjFndsb+Gdyksr2Xui5uoqW9k5YLJOsSyUq34xSVJFFXU8tcvMukXGsi8yXFW\nR7INty79cBX16hbLnmv2uQHuamW7bwC9ZbQDHauqY+6LGzlWVceKOyeTPDDU6khK2ZKI8PC1oyiu\nrOO/P9xFZHAAM0bprHOgd+56lPKaBua9tJmCshMsvWUCY2PCrI6klK35+frwtxvHMzYmjHtWbmfz\noVKrI9mCFr+HqKpzcMsrmzlQVMnz81KYlNjX6khKeYSgAF+W3jyB6D5B3PFqGvuPVlodyXJa/B6g\ntqGRO15NY2d+OX+dM57zh+pVT0qdjj69Anj11okE+vty89LNHC47YXUkS2nx21xVnYOFr21h06FS\nnrp+DDN0IhWlzkhMeE9euXUiVbUObnl5M+U1DVZHsowWv001Nhne2JzLBY+v5ct9xTzyg3OYNVZH\nu1DqbIwYFMrz888lu6SGO5elU9vQaHUkS2jx29A3B0q48pl1/Pq974jv24sP7prCnImxVsdSqls4\nb3AET90whrScUn7+xjYqar1vz19H8rKRzKJKHlm9ly/2FhETHsTf545n5qgBOnOWUh3sqtGDKKms\n4/f/2M0Fj3/JfZckMWdiLP5eMneF2PF25pSUFJOenm51jC5TWl3P05/v57VNufT09+Xui4Zw83nx\nBPrrsMpKdaZdBeU8/NFuNmaVkhjRi0Uzh3PpiP4eubMlIluMMSluravFb506RyPL1ufwzBcHqK5z\ncOOkWH5xyVD66kBrSnUZYwxf7C3ikdV7OFhczcSEcB66IpkxHnafjBa/zRljWJNRyB8/3kvOsRou\nGBbJb65IZmj/EKujKeW1HI1NrEzL48+f7edYdT2zxg7iV5cPI7qPZ4yDpcVvYzvzy3j4n3vYnF3K\n0P7BPHTlCL0uXykbqaxt4PmvsnhhXRYGuHVKPD+7YAi9g+w9z4UWvw0dKT/B42v28d7WAvr2CuD+\ny4ZyQ0qMToSulE0dLjvBk5/u571t+YQF+XPvxUnMnRxn2xPAWvw2Ul3n4Pmvs1jy9UGamuC2qQnc\ndeFgQgLtvfeglHLaVVDOI6v3sP7gMRIievFfM4Zz+Uj7nQDW4reBspp6PtlVyFOf7aeoso6rRg/k\nv2YM13HzlfJAxhjnjZSr93CgqIoJ8X146MoRthooUYvfAtV1DtKyS1l/8BjrD5aQcbgCY2BMTBj/\nfVUy58aFWx1RKXWWHI1Nrsld9lFSVc/VYwbx4OXDbLFDp8XfBeocjWzPLePbg8fYcLCEbbllOJoM\nAb4+jIsNY8qQCM4b3JfxsX3w8bHXW0Kl1NmpqnOw5KuDLFmXRVMT3DAhhqlJEYyP7UNkiDWXY2vx\nd4LGJkPG4XK+zXTu0adll1Lb0ISPwDlRvTnPVfQpceEEBeiNV0p5g8LyWp78dB8fbC+godHZpTHh\nQYyL6cP42DDGx/UheWBol5wQ1uLvAMYYDhRVsT6zhPUHj7Ex6xgVtQ4AhvYP5rzBzqKflNjX9pd5\nKaU6V21DIxmHy9maU8bW3ONszT3O0Yo6AAL9fRgdFca4uDDGx/bptHcFHV78IjIDeBrnnLsvGmMe\nbfG4uB6/Auecu7cYY7a6s21ruqr4a+odHK2o42hFbbOPOg6XnSAt+zglVc6/uNjwnpw3uC+pg/ty\n3uAIy97KKaU8gzGGw+W1bM1x/iewLbeMjMPlnfqu4HSKv91B2kTEF1gMXArkA2kissoYs7vZajOB\nJNfHJOBZYJKb23a4ekcTxVV1FJbXUnSy0CvrOFpey9HK2n+XfaVrD765IH9fBvQOZMqQvkwZHEHq\n4L62OHGjlPIcIkJUWBBRYUFcPWYQ8P13BZsOHWPVjsPA/70rGB/XhwcvH9bp5wXdGZ1zIpBpjMkC\nEJGVwCygeXnPApa5Jl3fKCJhIjIQiHdj2w7R1GS4ZvE3HCmr5Vh1/fce9/cV+oUE0j+0B0n9gpk6\nJIL+oc6vm/8Z3MPPdtfnKqU8X6C/L+fGhf/7Cr/W3hWs3VvEopnDOz2LO8UfBeQ1+zof5159e+tE\nublth/DxEZL6hXBOVBgD/qPQnZ/36RmgV9copWyjtXcFTU1dc87VNuPxi8gCYAFAbOyZTTry5xvG\ndmQkpZTqUl21c+rO2YQCIKbZ19GuZe6s4862ABhjlhhjUowxKZGROmiZUkp1FneKPw1IEpEEEQkA\nZgOrWqyzCpgvTpOBcmPMETe3VUop1YXaPdRjjHGIyN3AGpyXZC41xmSIyELX488Bq3FeypmJ83LO\nW0+1baf8JkoppdyiN3AppVQ3cDrX8dtzYGmllFKdRotfKaW8jBa/Ukp5GS1+pZTyMrY8uSsixUCO\nxTEigBKLM5wOT8sLmrmreFpmT8sL9sgcZ4xx6yYoWxa/HYhIurtnyO3A0/KCZu4qnpbZ0/KC52XW\nQz1KKeVltPiVUsrLaPG3bYnVAU6Tp+UFzdxVPC2zp+UFD8usx/iVUsrL6B6/Ukp5Ga8rfhFZKiJF\nIrKr2bJwEflMRA64/uzTxrYzRGSfiGSKyCIL8z4uIntFZKeIvC8iYW1smy0i34nIdhHpssGP2sj8\nexEpcGXZLiJXtLFtlz/Hp8j8ZrO82SKyvY1tu/x5FpEYEVkrIrtFJENE7nUtt/Nrua3Mtn09nyKz\nrV/P7TLGeNUHMB0YD+xqtuwxYJHr80XAn1rZzhc4CCQCAcAOYIRFeS8D/Fyf/6m1vK7HsoEImzzH\nvwceaGc7S57jtjK3ePxJ4L/t8jwDA4Hxrs9DgP3ACJu/ltvKbNvX8yky2/r13N6H1+3xG2O+Bkpb\nLJ4FvOr6/FXg2lY2/ffcw8aYeuDk/MGdqrW8xphPjTEnZ4rfiHOCG9to4zl2hyXPMZw6szgnYb4e\neKMrsrjDGHPEGLPV9XklsAfnVKd2fi23mtnOr+dTPM/usOz13B6vK/429DfOiWMACoH+razT1rzC\nVrsN+LiNxwzwuYhscU1tabWfu97OL23jEIRdn+NpwFFjzIE2Hrf0eRaReGAcsAkPeS23yNycbV/P\nrWT21NezFn9LxvkezSMudRKRhwAHsKKNVaYaY8YCM4G7RGR6l4X7vmdxvuUdCxzBeejEU8zh1Hv7\nlj3PIhIMvAvcZ4ypaP6YXV/LbWW28+u5lcye/HrW4nc5KiIDAVx/FrWyjtvzB3cFEbkFuAqY6/oH\n/j3GmALXn0XA+zjfelrCGHPUGNNojGkCXmgji62eYwAR8QN+CLzZ1jpWPc8i4o+zjFYYY95zLbb1\na7mNzLZ+PbeW2VNfzydp8TutAm52fX4z8GEr69hm/mARmQE8CFxjjKlpY51eIhJy8nOcJ9B2tbZu\nVzhZRi4/aCOLbZ7jZi4B9hpj8lt70Krn2XXe4SVgjzHmqWYP2fa13FZmO7+eT5HZU1/PTlafXe7q\nD5xv2Y8ADTiPud0O9AX+BRwAPgfCXesOAlY32/YKnGf1DwIPWZg3E+exw+2uj+da5sX5NnSH6yOj\nq/KeIvNy4DtgJ84X/0C7PMdtZXYtfwVY2GJdy59nYCrOwzg7m70OrrD5a7mtzLZ9PZ8is61fz+19\n6J27SinlZfRQj1JKeRktfqWU8jJa/Eop5WW0+JVSysto8SullJfR4ldKKS+jxa+UUl5Gi18ppbzM\n/wf5SXT+nzGTWQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "for i in size:\n", " k.append(toss(p,n))\n", "\n", "k.sort()\n", "\n", "for i in k:\n", " pk = binom.pmf(i,n,p)\n", " pmf.append(pk)\n", "\n", "plt.plot(k,pmf) " ] }, { "cell_type": "code", "execution_count": 38, "metadata": { "tags": [ "s1", "hid", "l1" ] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "continue\n" ] } ], "source": [ "ref_tmp_var = False\n", "\n", "try:\n", " head = k[:10]\n", " prob = pmf[:10]\n", " ks = [9, 9, 10, 10, 11, 11, 11, 11, 12, 12]\n", " pmfs = [round(x,4) for x in [0.00063412401653505015, 0.00063412401653505015, 0.0019974906520854141, 0.0019974906520854141, 0.005447701778414739, 0.005447701778414739, 0.005447701778414739, 0.005447701778414739, 0.012938291723735, 0.012938291723735]]\n", " if ks == head and pmfs == [round(a,4) for a in prob]:\n", " ref_assert_var = True\n", " ref_tmp_var = True\n", " else:\n", " ref_assert_var = False\n", " print('Please follow the instructions given and use the same variables provided in the instructions. ')\n", "except Exception:\n", " print('Please follow the instructions given and use the same variables provided in the instructions. ')\n", "\n", "assert ref_tmp_var" ] } ], "metadata": { "executed_sections": [], "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.1" } }, "nbformat": 4, "nbformat_minor": 2 }