binomial-example.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuMAAAEGCAYAAADPHJsIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAWc0lEQVR4nO3dfZSeZX0n8O/PEDddoEjRg1RqVRYFgxmUNxWjlFPARY8VUBbN7hHpYYvtonXdrbLStVvrqlst1R6qC4WCPTGAgOjxZX05C0KVVRJKAoKASqrBEBGhEDUIeO0f8xCGMJlMJjPPNcl8Pufk5Jnrfvs9c5175vtcc933Xa21AAAAw/ek3gUAAMBcJYwDAEAnwjgAAHQijAMAQCfCOAAAdLLTFpa71QoAAGy7Gq/RyDgAAHQijAMAQCfCOAAAdCKMAwBAJ1u6gBMAgO3EQw89lDVr1mTDhg29S5mzFixYkL333jvz58+f1PrV2oQ3THE3FQCA7cQdd9yRXXfdNXvssUeqxr15BzOotZZ77rknDzzwQJ797GdvutjdVAAAdmQbNmwQxDuqquyxxx5b9ZcJYRwAYAciiPe1td9/YRwAADoRxgEAdlDr7t8w4b/LP//l7PHUp+Xyz395i+uuu39yUy/uuuuunHTSSdlnn31y0EEH5dhjj81tt9027rovfelLkySrV6/OJz/5yY3ty5cvz1vf+tYpvecjjjgiy5cvH3fZT37yk8yfPz8f//jHp7TvmSCMA3TyL0/Ztn8A2+Lr11ydU9+0JOdeuDSHL375pNbfktZajjvuuBxxxBH53ve+lxUrVuT9739/1q1b97j1Hn744STJN77xjSRPDOMHH3xwPvrRj27N25mUT33qU3nxi1+cZcuWbXadRx55ZNqPOxFhHABgjplKED/1TUu2uN6VV16Z+fPn57TTTtvYNjIyksWLF+eqq67K4sWL85rXvCbPf/7zkyS77LJLkuRd73pXrrnmmhx44IE566yzctVVV+XVr351kmT9+vV585vfnBe84AVZtGhRLrvssiTJW97ylhx88MFZuHBh3vOe90zqfS9btiwf/vCHc+edd2bNmjUb23fZZZe84x3vyMjISK699tqsWLEir3jFK3LQQQflmGOOydq1a5Mk5557bg455JCMjIzkhBNOyM9//vNJHXciwjgAwBwy1SB+7oVLt7juTTfdlIMOOmizy6+//vp85CMfecK0lQ984ANZvHhxbrjhhrz97W9/3LL3vve92W233XLjjTdm1apVOfLII5Mk73vf+7J8+fKsWrUqX/va17Jq1aoJa/vhD3+YtWvX5tBDD82JJ56Yiy++eOOyn/3sZznssMOycuXKHHbYYTn99NNz6aWXZsWKFTnllFPy7ne/O0ly/PHH57rrrsvKlSuz//7757zzztvi92RLPPQHAGCO2JYgPpn1t+TQQw8d7/7bE/rqV7+aiy66aOPXu+++e5LkkksuyTnnnJOHH344a9euzc0335xFixZtdj8XX3xxTjzxxCTJSSedlFNOOSXveMc7kiTz5s3LCSeckCS59dZbc9NNN+Woo45KMjptZa+99koy+mHjzDPPzH333Zf169fnmGOO2ar3Mh5hHABgDhhGEF+4cGEuvfTSzS7feeedJ13vRO6444586EMfynXXXZfdd989J5988hbv7b1s2bLcddddWbp0dIT/Rz/6UW6//fbsu+++WbBgQebNm5dkdN77woULc+211z5hHyeffHKuuOKKjIyM5IILLshVV121ze/FNBUAgB3csEbEjzzyyDz44IM555xzNratWrUq11xzzYTb7brrrnnggQfGXXbUUUfl7LPP3vj1vffem/vvvz8777xzdtttt6xbty5f/OIXJ9z/bbfdlvXr1+fOO+/M6tWrs3r16pxxxhnjXsj5vOc9L3fffffGMP7QQw/l29/+dpLkgQceyF577ZWHHnpoY6jfVsI4AMAOas9fX5A9f31Bjn/V0bnnJ3fn+FcdvbFton+bW39Lqiqf/vSn89WvfjX77LNPFi5cmDPOOCNPf/rTJ9xu0aJFmTdvXkZGRnLWWWc9btmZZ56Ze++9NwcccEBGRkZy5ZVXZmRkJC984Quz33775Y1vfGMOP/zwCfe/bNmyHHfccY9rO+GEE8YN409+8pNz6aWX5p3vfGdGRkZy4IEHbrzry3vf+94cdthhOfzww7Pffvtt8fsxGdVam2j5hAsBmLptvT3hbvdNTx3AjuOWW27J/vvv37uMOW8z/TDuozmNjAMAQCfCOAAAdCKMAwDsQLYwBZkZtrXff7c2BNgKP/jtqW/7zH+evjrGNe5sxEnyuxt2CAsWLMg999yTPfbYI1Xb8kOBqWit5Z577smCBVu+2PVRwjgAwA5i7733zpo1a3L33Xf3LmXOWrBgQfbee+9Jry+MAwDsIObPn7/VT7ikL3PGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADrZqXcBADPpxkXbtv0LVk1PHdud2sbt27RUAbDDMzIOAACdCOMAANCJMA4AAJ0I4wAA0IkwDgAAnQjjAADQiTAOAACdCOMAANCJMA4AAJ0I4wAA0IkwDgAAnQjjAADQiTAOAACdCOMAANCJMA4AAJ0I4wAA0IkwDgAAnQjjAADQiTAOAACdCOMAANCJMA4AAJ0I4wAA0IkwDgAAnQjjAADQiTAOAACdCOMAANCJMA4AAJ0I4wAA0IkwDgAAnQjjAADQiTAOAACdCOMAANCJMA4AAJ0I4wAA0IkwDgAAnQjjAADQiTAOAACdCOMAANCJMA4AAJ0I4wAA0MlOvQsA2NQ/Lp76ti+7ZvrqYBrVNmzbpq0KgFnHyDgAAHQijAMAQCfCOAAAdCKMAwBAJ8I4AAB0IowDAEAnwjgAAHQijAMAQCfCOAAAdCKMAwBAJ8I4AAB0IowDAEAnwjgAAHQijAMAQCfCOAAAdCKMAwBAJ8I4AAB0IowDAEAnwjgAAHQijAMAQCfCOAAAdCKMAwBAJ8I4AAB0IowDAEAnwjgAAHQijAMAQCfCOAAAdCKMAwBAJ8I4AAB0IowDAEAnwjgAAHQijAMAQCfCOAAAdCKMAwBAJ8I4AAB0IowDAEAnwjgAAHQijAMAQCfCOAAAdCKMAwBAJzv1LgDY/n3+1du2/as+Nz11MEfUNm7fpqUKgGlhZBwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOhHGAQCgE2EcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOtmpdwFAH59cMvVt37h0+uqA7mobt2/TUgUwRxkZBwCAToRxAADoRBgHAIBOhHEAAOhEGAcAgE6EcQAA6EQYBwCAToRxAADoRBgHAIBOhHEAAOhEGAcAgE6EcQAA6EQYBwCAToRxAADoRBgHAIBOhHEAAOhEGAcAgE6EcQAA6EQYBwCAToRxAADoRBgHAIBOhHEAAOhkp94FAJPzsT/ctu3f8rfTUweweet23bDx9dcfvjqnbliScxcszeE7vXyL23794atz/C+OnsnygFnIyDgATLOpBPFTNywZQmXAbCOMA8A0mmoQP3fB0iFUB8w2wjgATJNtCeKTWR/Y8QjjADBNBHFgawnjADBNBHFgawnjADBNBHFgawnjADAkgjiwKWEcAIZAEAfGI4wDwAwTxIHNEcYBYAYJ4sBEhHEAmCGCOLAlwjgAzICpPAAImHuEcQCYZlN9Eicw9+zUuwDYkf3puzYkSVb/4OpcfsWSHP/apXnWM7f8i3n1D67OP3zy6JkuD5gBUw3i5y5Y+oRl63bdsE373/OBBZMvHOjCyDjMsKkE8cuvMEIG26NtCeIeGARzkzAOM2iqQfz41z5xhAyY/WZTEDcHHbYPwjjMoKkG8cmsD8w+symIm4MO2wdhHGaQIA5zy2wK4uPNQQdmH2EcZpAgDoxlTjmwKWEcOhLEYe6YSlAWxGHH59aGMMbb/udjtxFb8/2r84VlS3LsG5Zm7+ds+Rfbmu9fncv+bvK3IxTEYe6YalCeziA+9jaJUxlxP/4XbrcKM8HIOIxjKkH8C8smf7GUIA5zx7aMWM+GqSwuBoWZJYzDJqYaxI99w+QulhLEYe6Y6akjLgaF7Z9pKmz3TvmbJz6h7lFrb786V56/JL9zytLste8Tf/Gcf/rjn063LUF8MusL4jB3TGkqSCY/FaT3xaCbPh10a/fv6aAwShhn6F534egP8LtvuTr/72+X5MV/uDRP23/LP7gfXf/B+++e1HG2FMTHI4gD02Gqwff4TO7n22y7GHRbP3hsGuy3dv+CPduzaq1tfmHV/0ny1Gk+5lOT/GSa98nW0w+zg37oTx/MDvqhP30wO+iH/maqD37SWnvlpo0ThvGZUFXLW2sHD/WgPIF+mB30Q3/6YHbQD/3pg9lBP/Q37D5wAScAAHQijAMAQCc9wvg5HY7JE+mH2UE/9KcPZgf90J8+mB30Q39D7YOhzxkHAABGmaYCAACdCOMAANDJUMN4Vb2yqm6tqu9W1buGeWweU1Wrq+rGqrqhqpb3rmcuqKrzq+rHVXXTmLbfqKqvVNXtg/9371njXLCZfvizqrpzcD7cUFXH9qxxR1dVv1VVV1bVzVX17ap626Dd+TBEE/SD82FIqmpBVX2rqlYO+uB/DNqfXVXfHGSli6vqyb1r3ZFN0A8XVNUdY86FA2eshmHNGa+qeUluS3JUkjVJrkvyhtbazUMpgI2qanWSg1trHiowJFX18iTrk3yitXbAoO1/Jflpa+0Dgw+nu7fW3tmzzh3dZvrhz5Ksb619qGdtc0VV7ZVkr9ba9VW1a5IVSV6b5OQ4H4Zmgn44Mc6HoaiqSrJza219Vc1P8o9J3pbkPye5vLV2UVV9PMnK1trHeta6I5ugH05L8rnW2qUzXcMwR8YPTfLd1tr3W2u/THJRkt8b4vGhm9ba1Ul+uknz7yW5cPD6woz+ImQGbaYfGKLW2trW2vWD1w8kuSXJM+J8GKoJ+oEhaaPWD76cP/jXkhyZ5NEA6FyYYRP0w9AMM4w/I8kPx3y9Jk78XlqSL1fViqr6j72LmcP2bK2tHby+K8mePYuZ4/5TVa0aTGMxPWJIqupZSV6Y5JtxPnSzST8kzoehqap5VXVDkh8n+UqS7yW5r7X28GAVWWkINu2H1tqj58L7BufCWVX1r2bq+C7gnJte1lp7UZJ/m+SPBn+6p6M2Ol/MfUb7+FiSfZIcmGRtkg/3LWduqKpdklyW5I9ba/ePXeZ8GJ5x+sH5MESttUdaawcm2TujMwj261zSnLRpP1TVAUnOyGh/HJLkN5LM2LS5YYbxO5P81piv9x60MWSttTsH//84yacz+gOA4Vs3mLf56PzNH3euZ05qra0b/CD+VZJz43yYcYN5mZclWdpau3zQ7HwYsvH6wfnQR2vtviRXJnlJkqdU1U6DRbLSEI3ph1cOpnK11tqDSf4+M3guDDOMX5dk38FVwk9OclKSzw7x+CSpqp0HF+ukqnZOcnSSmybeihny2SRvGrx+U5LPdKxlzno0AA4cF+fDjBpcLHVekltaa381ZpHzYYg21w/Oh+GpqqdV1VMGr38toze4uCWjYfB1g9WcCzNsM/3wnTGDA5XRefszdi4M9Qmcg1sk/XWSeUnOb629b2gHJ0lSVc/J6Gh4kuyU5JP6YeZV1bIkRyR5apJ1Sd6T5IoklyR5ZpJ/TnJia83FhTNoM/1wREb/JN+SrE7yB2PmLjPNquplSa5JcmOSXw2a/1tG5ys7H4Zkgn54Q5wPQ1FVizJ6gea8jA6OXtJa+/PB7+mLMjo14p+S/PvB6CwzYIJ++L9JnpakktyQ5LQxF3pObw3DDOMAAMBjXMAJAACdCOMAANCJMA4AAJ0I4wAA0IkwDgAAnQjjwJxVVa2qPjzm6/9SVX82Tfu+oKpet+U1t/k4r6+qW6rqyk3aj6iqz83QMbfq9l5Vdengdm2bW35RVe277ZUBbH+EcWAuezDJ8VX11N6FjDXm6XuT8ftJTm2t/c5M1bMtqmphknmtte9PsNrHkvzJkEoCmFWEcWAuezjJOUnevumCTUe2Hx0NHow4f62qPlNV36+qD1TVkqr6VlXdWFX7jNnN71bV8qq6rapePdh+XlX9ZVVdV1WrquoPxuz3mqr6bJKbx6nnDYP931RVHxy0/fckL0tyXlX95Tjvb5fBqPR3qmrp4ElyqaqDBu9hRVV9acyT5k4d1LWyqi6rqn89aH92VV07OP5fjKlpr6q6uqpuGNS1eJwalmTwBMGqes1g3Ruq6taqumOwzjWD79XWfAgB2CEI48Bcd3aSJVW121ZsM5LktCT7J/kPSZ7bWjs0yd8lOX3Mes9KcmiSVyX5eFUtyOhI9r+01g5JckiSU6vq2YP1X5Tkba215449WFX9ZpIPJjkyo09HPKSqXtta+/Mky5Msaa3913HqfGGSP07y/CTPSXJ4Vc1P8jdJXtdaOyjJ+UkefQrv5a21Q1prIxl9LPfvD9o/kuRjrbUXJBn7NMY3JvlSa+3AwffkhnFqODzJiiRprX22tXbgYP2VST40aP9Vku8O9gEwpxiFAOa01tr9VfWJJG9N8otJbnbdo48Ir6rvJfnyoP3GJGOni1wyCJq3V9X3k+yX5Ogki8aMuu+WZN8kv0zyrdbaHXmiQ5Jc1Vq7e3DMpUlenuSKLdT5rdbamsE2N2T0w8F9SQ5I8pXBQPm8PBawDxiMfD8lyS5JvjRoPzzJCYPX/5DRDwZJcl2S8wcB/4rW2nhhfK8kd49tqKo/SfKL1trZY5p/nOQ3MwjuAHOFMA6Q/HWS65P8/Zi2hzP462FVPSnJk8cse3DM61+N+fpXefzP1bbJcVqSSnJ6a+1LYxdU1RFJfja18jdrbJ2PDGqrJN9urb1knPUvSPLa1trKqjo5yRFjlm36XtJau7qqXp7Rkf8LquqvWmuf2GS1XyRZ8OgXVfW7SV6f0Q8TYy3I5D8MAewwTFMB5rzW2k+TXJLHpmUkyeokBw1evybJ/Cns+vVV9aTBPPLnJLk1o6PNbxmMJqeqnltVO29hP99K8oqqempVzUvyhiRfm0I9GdTwtKp6yeD48wcXWSbJrknWDmpbMmabryc5afB6Y3tV/XaSda21czM6RedF4xzvliT/Zsz6Zyd5fWtt0+D93CQ3TfE9AWy3hHGAUR9OMvauKudmNACvTPKSTG3U+gcZDdJfTHJaa21DRkPrzUmur6qbkvzvbOGvlIMpMe9KcmVG51qvaK19Zgr1pLX2yySvS/LBwXu7IclLB4v/NMk3Mxq+vzNms7cl+aOqujHJM8a0H5FkZVX9U5J/l9G55Zv6fB4bYT85yR5JrhhcxPmFJKmqPTM6beWuqbwngO1ZtfaEvzwCwLSoql/L6IeIw1trj2xmnbcnub+1dt5QiwOYBYyMAzBjBtNR3pPHj6hv6r4kFw6nIoDZxcg4AAB0YmQcAAA6EcYBAKATYRwAADoRxgEAoBNhHAAAOvn/6KzdLbdQ5g4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 936x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "from scipy.stats import binom, nbinom\n",
    "from dsbay.distribution import DistributionMass\n",
    "\n",
    "N = 34\n",
    "z = 11\n",
    "theta = 0.5\n",
    "x = [_z for _z in range(1,N+1)]\n",
    "y = [binom.pmf(_z, N, theta) for _z in range(1,N+1)]\n",
    "\n",
    "d = DistributionMass(x,y)\n",
    "ax = d.plot()\n",
    "_ = ax.set_xlabel(\"Number of heads (z)\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}