{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Distributions" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "Probability distributions reflect the probabilities of occurence for the possible outcomes of a function / data source. \n", "
\n", "\n", "
\n", "Probability distributions on \n", "wikipedia.\n", "If you want a more general refresher on probability / distributions, check out this\n", "article.\n", "
" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Probability Distributions\n", "\n", "Typically, given a data source, we want to think about and check what kind of probability distribution our data sample appears to follow. More specifically, we are trying to infer the probability distribution that the data generator follows, asking the question: what function could it be replaced by?\n", "\n", "Checking the distribution of data is important, as we typically want to apply statistical tests to our data, and many statistical tests come with underlying assumptions about the distributions of the data they are applied to. Ensuring we apply appropriate statistical methodology requires thinking about and checking the distribution of our data. \n", "\n", "Informally, we can start by visualizing our data, and seeing what 'shape' it takes, and which distribution it appears to follow. More formally, we can statistically test whether a sample of data follows a particular distribution.\n", "\n", "Here we will start be visualizing some of the most common distributions. `Scipy` (scipy.stats) has code for working with, and generating different distributions. We will generate synthetic data from different underlying distributions, and do a quick survey of how they look, plotting histograms of the generated data. \n", "\n", "You can use this notebook to explore different parameters to get a feel for these distributions. For further exploration, explore plotting the probability density functions of each distribution. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Uniform Distribution" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "A uniform distribution is a distribution in which each possible value is equally probable. \n", "
\n", "\n", "
\n", "Uniform distribution on\n", "wikipedia.\n", "
" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "from scipy.stats import uniform" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "data = uniform.rvs(size=10000)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAPQElEQVR4nO3df4ylVX3H8fdHVrT+XGAHQ3eXLsa1lZg0koliTax1jRVsWP6ABlPrSjbdxKK1YlrX9g8aTRrsL1pSg90W6tJYhVJTNpbWEH7EtukSF7HIjxq2aGEKlbHA9gexSv32j3vQYZnduXufmTsO5/1KJvd5znOee87ZGT733PPc+5CqQpLUh+esdgckSdNj6EtSRwx9SeqIoS9JHTH0Jakj61a7A0ezYcOG2rJly2p3Q5LWlNtvv/2bVTWz2LEf6NDfsmULBw4cWO1uSNKakuRfj3TM5R1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SerID/Q3cqWlbNn916vW9tcvffuqtS1NytCXJrRaLzi+2GgIl3ckqSOGviR1xOUdSWPrbUnr2XjNaMmZfpKrkjyS5K4FZScmuTHJfe3xhFaeJJcnOZjkziRnLDhnR6t/X5IdKzIaSdJRjTPT/yTwh8DVC8p2AzdV1aVJdrf9DwFnAVvbz+uAK4DXJTkRuASYBQq4Pcm+qnpsuQaiZ+esRNLyWnKmX1VfAB49rHg7sLdt7wXOXVB+dY3sB9YnOQX4aeDGqnq0Bf2NwNuWYwCSpPFNuqb/sqp6GKCqHk5ycivfCDy4oN5cKztS+TMk2QXsAjj11FMn7N7qWs0ZtyQdzXJ/eieLlNVRyp9ZWLWnqmaranZmZtH/xaMkaUKTzvS/keSUNss/BXiklc8BmxfU2wQ81MrfdFj5rRO2LakzvntePpOG/j5gB3Bpe7x+Qfl7k3yG0YXcQ+2F4fPAbz71KR/grcCHJ+/2ePxDkaSnWzL0k3ya0Sx9Q5I5Rp/CuRS4NslO4AHg/Fb9BuBs4CDwBHAhQFU9muSjwBdbvY9U1eEXhyVJK2zJ0K+qdxzh0LZF6hZw0RGe5yrgqmPqnaRn8B2shvA2DJLUEUNfkjpi6EtSRwx9SeqIoS9JHfHWyloWfqJEWhuc6UtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjoyKPSTfCDJ3UnuSvLpJM9PclqS25Lcl+SaJMe3us9r+wfb8S3LMQBJ0vgmDv0kG4FfAmar6tXAccAFwMeAy6pqK/AYsLOdshN4rKpeAVzW6kmSpmjo8s464IeSrANeADwMvBm4rh3fC5zbtre3fdrxbUkysH1J0jGYOPSr6t+A3wEeYBT2h4Dbgcer6slWbQ7Y2LY3Ag+2c59s9U+atH1J0rEbsrxzAqPZ+2nADwMvBM5apGo9dcpRji183l1JDiQ5MD8/P2n3JEmLGLK88xbga1U1X1XfAT4L/ASwvi33AGwCHmrbc8BmgHb8pcCjhz9pVe2pqtmqmp2ZmRnQPUnS4YaE/gPAmUle0NbmtwH3ALcA57U6O4Dr2/a+tk87fnNVPWOmL0laOUPW9G9jdEH2S8BX2nPtAT4EXJzkIKM1+yvbKVcCJ7Xyi4HdA/otSZrAuqWrHFlVXQJccljx/cBrF6n7LeD8Ie1JkobxG7mS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHVkUOgnWZ/kuiT/nOTeJK9PcmKSG5Pc1x5PaHWT5PIkB5PcmeSM5RmCJGlcQ2f6fwD8bVX9GPDjwL3AbuCmqtoK3NT2Ac4CtrafXcAVA9uWJB2jiUM/yUuANwJXAlTVt6vqcWA7sLdV2wuc27a3A1fXyH5gfZJTJu65JOmYDZnpvxyYB/40yR1J/iTJC4GXVdXDAO3x5FZ/I/DggvPnWtnTJNmV5ECSA/Pz8wO6J0k63JDQXwecAVxRVa8B/ofvL+UsJouU1TMKqvZU1WxVzc7MzAzoniTpcENCfw6Yq6rb2v51jF4EvvHUsk17fGRB/c0Lzt8EPDSgfUnSMZo49Kvq34EHk/xoK9oG3APsA3a0sh3A9W17H/Cu9imeM4FDTy0DSZKmY93A898HfCrJ8cD9wIWMXkiuTbITeAA4v9W9ATgbOAg80epKkqZoUOhX1ZeB2UUObVukbgEXDWlPkjSM38iVpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjg0M/yXFJ7kjyubZ/WpLbktyX5Jokx7fy57X9g+34lqFtS5KOzXLM9N8P3Ltg/2PAZVW1FXgM2NnKdwKPVdUrgMtaPUnSFA0K/SSbgLcDf9L2A7wZuK5V2Quc27a3t33a8W2tviRpSobO9H8f+FXgu23/JODxqnqy7c8BG9v2RuBBgHb8UKv/NEl2JTmQ5MD8/PzA7kmSFpo49JP8DPBIVd2+sHiRqjXGse8XVO2pqtmqmp2ZmZm0e5KkRawbcO4bgHOSnA08H3gJo5n/+iTr2mx+E/BQqz8HbAbmkqwDXgo8OqB9SdIxmnimX1UfrqpNVbUFuAC4uap+DrgFOK9V2wFc37b3tX3a8Zur6hkzfUnSylmJz+l/CLg4yUFGa/ZXtvIrgZNa+cXA7hVoW5J0FEOWd76nqm4Fbm3b9wOvXaTOt4Dzl6M9SdJk/EauJHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdmTj0k2xOckuSe5PcneT9rfzEJDcmua89ntDKk+TyJAeT3JnkjOUahCRpPENm+k8CH6yqVwFnAhclOR3YDdxUVVuBm9o+wFnA1vazC7hiQNuSpAlMHPpV9XBVfalt/xdwL7AR2A7sbdX2Aue27e3A1TWyH1if5JSJey5JOmbLsqafZAvwGuA24GVV9TCMXhiAk1u1jcCDC06ba2WHP9euJAeSHJifn1+O7kmSmsGhn+RFwF8Cv1xV/3m0qouU1TMKqvZU1WxVzc7MzAztniRpgUGhn+S5jAL/U1X12Vb8jaeWbdrjI618Dti84PRNwEND2pckHZshn94JcCVwb1X93oJD+4AdbXsHcP2C8ne1T/GcCRx6ahlIkjQd6wac+wbg54GvJPlyK/s14FLg2iQ7gQeA89uxG4CzgYPAE8CFA9qWJE1g4tCvqr9n8XV6gG2L1C/goknbkyQN5zdyJakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqyNRDP8nbknw1ycEku6fdviT1bKqhn+Q44OPAWcDpwDuSnD7NPkhSz6Y9038tcLCq7q+qbwOfAbZPuQ+S1K11U25vI/Dggv054HULKyTZBexqu/+d5KsTtrUB+OaE565VjrkPjrkD+digMf/IkQ5MO/SzSFk9badqD7BncEPJgaqaHfo8a4lj7oNj7sNKjXnayztzwOYF+5uAh6bcB0nq1rRD/4vA1iSnJTkeuADYN+U+SFK3prq8U1VPJnkv8HngOOCqqrp7hZobvES0BjnmPjjmPqzImFNVS9eSJD0r+I1cSeqIoS9JHVnzob/UbR2SPC/JNe34bUm2TL+Xy2uMMV+c5J4kdya5KckRP7O7Vox7+44k5yWpJGv+433jjDnJz7bf9d1J/nzafVxuY/xtn5rkliR3tL/vs1ejn8slyVVJHkly1xGOJ8nl7d/jziRnDG60qtbsD6OLwf8CvBw4Hvgn4PTD6vwi8Im2fQFwzWr3ewpj/ingBW37PT2MudV7MfAFYD8wu9r9nsLveStwB3BC2z95tfs9hTHvAd7Ttk8Hvr7a/R445jcCZwB3HeH42cDfMPqO05nAbUPbXOsz/XFu67Ad2Nu2rwO2JVnsS2JrxZJjrqpbquqJtruf0fch1rJxb9/xUeC3gG9Ns3MrZJwx/wLw8ap6DKCqHplyH5fbOGMu4CVt+6Ws8e/5VNUXgEePUmU7cHWN7AfWJzllSJtrPfQXu63DxiPVqaongUPASVPp3coYZ8wL7WQ0U1jLlhxzktcAm6vqc9Ps2Aoa5/f8SuCVSf4hyf4kb5ta71bGOGP+DeCdSeaAG4D3Tadrq+ZY/3tf0rRvw7Dclrytw5h11pKxx5PkncAs8JMr2qOVd9QxJ3kOcBnw7ml1aArG+T2vY7TE8yZG7+b+Lsmrq+rxFe7bShlnzO8APllVv5vk9cCftTF/d+W7tyqWPb/W+kx/nNs6fK9OknWM3hIe7e3UD7qxbmWR5C3ArwPnVNX/TqlvK2WpMb8YeDVwa5KvM1r73LfGL+aO+7d9fVV9p6q+BnyV0YvAWjXOmHcC1wJU1T8Cz2d0M7Znq2W/dc1aD/1xbuuwD9jRts8Dbq52hWSNWnLMbanjjxgF/lpf54UlxlxVh6pqQ1VtqaotjK5jnFNVB1anu8tinL/tv2J00Z4kGxgt99w/1V4ur3HG/ACwDSDJqxiF/vxUezld+4B3tU/xnAkcqqqHhzzhml7eqSPc1iHJR4ADVbUPuJLRW8CDjGb4F6xej4cbc8y/DbwI+It2zfqBqjpn1To90JhjflYZc8yfB96a5B7g/4Bfqar/WL1eDzPmmD8I/HGSDzBa5nj3Wp7EJfk0o+W5De06xSXAcwGq6hOMrlucDRwEngAuHNzmGv73kiQdo7W+vCNJOgaGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SerI/wM+mvwjd40ptgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.hist(data);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Normal Distribution" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "The Normal (also Gaussian, or 'Bell Curve') distribution, is a distribution defined by it's mean and standard deviation.\n", "
\n", "\n", "
\n", "Normal distribution on \n", "wikipedia.\n", "
" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "from scipy.stats import norm" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "data = norm.rvs(size=10000)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAR80lEQVR4nO3df4xl5X3f8fen4B+xE2exGVtkd92hzcoNsZwYjTCppcryungBi6VVkEBJvHJWWkXFrVOnspdaCmosS2u5ConVlGprtsYqxUGOI1aGFG8wFqoUCIuNMXjtMCKUnUC8Yy0mSVHibvLtH/fZ+mb2zs6POzt3huf9kq7mnO95ztzvwO5nnn3uueemqpAk9eEfTLoBSdL6MfQlqSOGviR1xNCXpI4Y+pLUkfMn3cDZXHjhhTU9PT3pNiRpU3n00Ue/V1VTo45t6NCfnp7m6NGjk25DkjaVJP97sWMu70hSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkc29DtypXNtev89qz73mQNXr2En0vpwpi9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkeWDP0kh5KcSPLEiGP/LkklubDtJ8mnk8wmeTzJpUNj9yR5qj32rO2PIUlajuXM9D8L7FpYTLId+OfAs0PlK4Ed7bEPuLWNfT1wM/AO4DLg5iQXjNO4JGnllgz9qnoQODni0C3AR4Aaqu0GPlcDDwFbklwEvBc4UlUnq+oF4AgjfpFIks6tVa3pJ7kG+LOq+saCQ1uB40P7c622WH3U996X5GiSo/Pz86tpT5K0iBWHfpLXAB8Dfn3U4RG1Okv9zGLVwaqaqaqZqamplbYnSTqL1dxw7R8DFwPfSAKwDfhakssYzOC3D43dBjzX6u9aUP/qKp5b2jC8WZs2oxXP9Kvqm1X1xqqarqppBoF+aVX9OXAYeH+7iudy4MWqeh64D7giyQXtBdwrWk2StI6Wc8nmncAfAW9JMpdk71mG3ws8DcwC/xX4VwBVdRL4OPBIe/xGq0mS1tGSyztVdcMSx6eHtgu4cZFxh4BDK+xPkrSGfEeuJHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6shq3pErbRjjvCtW6pEzfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSPee0eagHHuGfTMgavXsBP1ZjkfjH4oyYkkTwzVPpXk20keT/L7SbYMHbspyWyS7yR571B9V6vNJtm/9j+KJGkpy1ne+Sywa0HtCPDWqnob8CfATQBJLgGuB366nfOfk5yX5Dzgd4ArgUuAG9pYSdI6WjL0q+pB4OSC2per6lTbfQjY1rZ3A5+vqr+pqj8FZoHL2mO2qp6uqh8An29jJUnraC1eyP1l4A/a9lbg+NCxuVZbrH6GJPuSHE1ydH5+fg3akySdNlboJ/kYcAq443RpxLA6S/3MYtXBqpqpqpmpqalx2pMkLbDqq3eS7AHeB+ysqtMBPgdsHxq2DXiubS9WlyStk1XN9JPsAj4KXFNVLw0dOgxcn+RVSS4GdgB/DDwC7EhycZJXMnix9/B4rUuSVmrJmX6SO4F3ARcmmQNuZnC1zquAI0kAHqqqX6mqJ5PcBXyLwbLPjVX1t+37fBC4DzgPOFRVT56Dn0eSdBZLhn5V3TCifNtZxn8C+MSI+r3AvSvqTpK0prwNgyR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI6v+jFxprUzvv2fSLUjdcKYvSR0x9CWpI0uGfpJDSU4keWKo9vokR5I81b5e0OpJ8ukks0keT3Lp0Dl72vinkuw5Nz+OJOlsljPT/yywa0FtP3B/Ve0A7m/7AFcCO9pjH3ArDH5JADcD7wAuA24+/YtCkrR+lgz9qnoQOLmgvBu4vW3fDlw7VP9cDTwEbElyEfBe4EhVnayqF4AjnPmLRJJ0jq12Tf9NVfU8QPv6xlbfChwfGjfXaovVz5BkX5KjSY7Oz8+vsj1J0ihr/UJuRtTqLPUzi1UHq2qmqmampqbWtDlJ6t1qQ/+7bdmG9vVEq88B24fGbQOeO0tdkrSOVhv6h4HTV+DsAe4eqr+/XcVzOfBiW/65D7giyQXtBdwrWk2StI6WfEdukjuBdwEXJpljcBXOAeCuJHuBZ4Hr2vB7gauAWeAl4AMAVXUyyceBR9q436iqhS8OS5LOsSVDv6puWOTQzhFjC7hxke9zCDi0ou4kSWvKd+RKUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRPy5R2mTG/XjJZw5cvUadaDNypi9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHVkrNBP8m+TPJnkiSR3Jnl1kouTPJzkqSS/m+SVbeyr2v5sOz69Fj+AJGn5Vh36SbYC/waYqaq3AucB1wOfBG6pqh3AC8Dedspe4IWq+kngljZOkrSOxl3eOR/4kSTnA68BngfeDXyhHb8duLZt7277tOM7k2TM55ckrcCqQ7+q/gz4j8CzDML+ReBR4PtVdaoNmwO2tu2twPF27qk2/g0Lv2+SfUmOJjk6Pz+/2vYkSSOMs7xzAYPZ+8XATwCvBa4cMbROn3KWYz8sVB2sqpmqmpmamlpte5KkEcZZ3nkP8KdVNV9V/xf4IvBPgS1tuQdgG/Bc254DtgO04z8OnBzj+SVJKzRO6D8LXJ7kNW1tfifwLeAB4OfbmD3A3W37cNunHf9KVZ0x05cknTvjrOk/zOAF2a8B32zf6yDwUeDDSWYZrNnf1k65DXhDq38Y2D9G35KkVRjrM3Kr6mbg5gXlp4HLRoz9a+C6cZ5PkjQe35ErSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdGesdudJp0/vvmXQLkpbBmb4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHVkrNBPsiXJF5J8O8mxJD+X5PVJjiR5qn29oI1Nkk8nmU3yeJJL1+ZHkCQt17gz/d8G/mdV/RPgZ4BjwH7g/qraAdzf9gGuBHa0xz7g1jGfW5K0QqsO/SSvA/4ZcBtAVf2gqr4P7AZub8NuB65t27uBz9XAQ8CWJBetunNJ0oqNM9P/R8A88N+SfD3JZ5K8FnhTVT0P0L6+sY3fChwfOn+u1f6eJPuSHE1ydH5+foz2JEkLjRP65wOXArdW1duB/8MPl3JGyYhanVGoOlhVM1U1MzU1NUZ7kqSFxgn9OWCuqh5u+19g8Evgu6eXbdrXE0Pjtw+dvw14boznlySt0KpDv6r+HDie5C2ttBP4FnAY2NNqe4C72/Zh4P3tKp7LgRdPLwNJktbHuB+i8q+BO5K8Enga+ACDXyR3JdkLPAtc18beC1wFzAIvtbGS1tk4H3jzzIGr17ATTcJYoV9VjwEzIw7tHDG2gBvHeT5J0nh8R64kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUkbFDP8l5Sb6e5Ett/+IkDyd5Ksnvtg9NJ8mr2v5sOz497nNLklZmLWb6HwKODe1/ErilqnYALwB7W30v8EJV/SRwSxsnSVpHY4V+km3A1cBn2n6AdwNfaENuB65t27vbPu34zjZekrROzh/z/N8CPgL8WNt/A/D9qjrV9ueArW17K3AcoKpOJXmxjf/e8DdMsg/YB/DmN795zPa0EtP775l0C5LOsVXP9JO8DzhRVY8Ol0cMrWUc+2Gh6mBVzVTVzNTU1GrbkySNMM5M/53ANUmuAl4NvI7BzH9LkvPbbH8b8FwbPwdsB+aSnA/8OHByjOeXJK3Qqmf6VXVTVW2rqmngeuArVfULwAPAz7dhe4C72/bhtk87/pWqOmOmL0k6d87FdfofBT6cZJbBmv1trX4b8IZW/zCw/xw8tyTpLMZ9IReAqvoq8NW2/TRw2Ygxfw1ctxbPJ0laHd+RK0kdMfQlqSOGviR1ZE3W9CX1YZw38D1z4Oo17ESr5Uxfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjqw69JNsT/JAkmNJnkzyoVZ/fZIjSZ5qXy9o9ST5dJLZJI8nuXStfghJ0vKMM9M/BfxaVf0UcDlwY5JLgP3A/VW1A7i/7QNcCexoj33ArWM8tyRpFVYd+lX1fFV9rW3/JXAM2ArsBm5vw24Hrm3bu4HP1cBDwJYkF626c0nSiq3Jmn6SaeDtwMPAm6rqeRj8YgDe2IZtBY4PnTbXagu/174kR5McnZ+fX4v2JEnN2J+Rm+RHgd8DfrWq/iLJokNH1OqMQtVB4CDAzMzMGcd1duN8hqmkl7+xQj/JKxgE/h1V9cVW/m6Si6rq+bZ8c6LV54DtQ6dvA54b5/klbR5+qPrGMM7VOwFuA45V1W8OHToM7Gnbe4C7h+rvb1fxXA68eHoZSJK0PsaZ6b8T+CXgm0kea7V/DxwA7kqyF3gWuK4duxe4CpgFXgI+MMZzS5JWYdWhX1X/i9Hr9AA7R4wv4MbVPp8kaXy+I1eSOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1ZOzbMGhteRsFSeeSM31J6oihL0kdMfQlqSOGviR1xNCXpI549Y6kDc978a8dZ/qS1BFDX5I6YuhLUkcMfUnqiC/kngPeSkHSRmXoS3pZG3cS9nK7+mfdQz/JLuC3gfOAz1TVgfXuYTmcrUt6OVrXNf0k5wG/A1wJXALckOSS9exBknq23jP9y4DZqnoaIMnngd3At87FkzlblzSuSeXIuVpWWu/Q3wocH9qfA94xPCDJPmBf2/2rJN9Zp95W40Lge5NuYhk2Q5+boUfYHH1uhh5hc/Q5sR7zyRUNX9jnP1xs4HqHfkbU6u/tVB0EDq5PO+NJcrSqZibdx1I2Q5+boUfYHH1uhh5hc/S5GXqElfW53tfpzwHbh/a3Ac+tcw+S1K31Dv1HgB1JLk7ySuB64PA69yBJ3VrX5Z2qOpXkg8B9DC7ZPFRVT65nD2tsUyxDsTn63Aw9wuboczP0CJujz83QI6ygz1TV0qMkSS8L3ntHkjpi6EtSRwz9MSX5eJLHkzyW5MtJfmLSPY2S5FNJvt16/f0kWybd00JJrkvyZJK/S7KhLpNLsivJd5LMJtk/6X5GSXIoyYkkT0y6l8Uk2Z7kgSTH2v/rD026p1GSvDrJHyf5RuvzP0y6p8UkOS/J15N8aTnjDf3xfaqq3lZVPwt8Cfj1STe0iCPAW6vqbcCfADdNuJ9RngD+JfDgpBsZtoluH/JZYNekm1jCKeDXquqngMuBGzfof8u/Ad5dVT8D/CywK8nlE+5pMR8Cji13sKE/pqr6i6Hd17LgzWYbRVV9uapOtd2HGLxHYkOpqmNVtRHfgf3/bx9SVT8ATt8+ZEOpqgeBk5Pu42yq6vmq+lrb/ksGYbV1sl2dqQb+qu2+oj023N/tJNuAq4HPLPccQ38NJPlEkuPAL7BxZ/rDfhn4g0k3sYmMun3IhguqzSbJNPB24OHJdjJaWzZ5DDgBHKmqjdjnbwEfAf5uuScY+suQ5A+TPDHisRugqj5WVduBO4APbtQ+25iPMfgn9h0btccNaMnbh2hlkvwo8HvAry741/KGUVV/25ZttwGXJXnrpHsaluR9wImqenQl5/khKstQVe9Z5tD/AdwD3HwO21nUUn0m2QO8D9hZE3qDxgr+W24k3j5kDSV5BYPAv6OqvjjpfpZSVd9P8lUGr5dspBfJ3wlck+Qq4NXA65L896r6xbOd5Ex/TEl2DO1eA3x7Ur2cTfvwmo8C11TVS5PuZ5Px9iFrJEmA24BjVfWbk+5nMUmmTl/hluRHgPewwf5uV9VNVbWtqqYZ/Jn8ylKBD4b+WjjQliceB65g8Er6RvSfgB8DjrTLS//LpBtaKMm/SDIH/BxwT5L7Jt0TDG4fwmDZ7j4GLzzetRFvH5LkTuCPgLckmUuyd9I9jfBO4JeAd7c/h4+1mepGcxHwQPt7/QiDNf1lXRK50XkbBknqiDN9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I68v8AHOHXVcCwS80AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.hist(data, bins=20);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Bernouilli Distribution" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "The Bernouilli Distribution is a binary distribution - it takes only two values (0 or 1), with some probably 'p'. \n", "
\n", "\n", "
\n", "Bernouilli distribution on wikipedia\n", "wikipedia.\n", "
" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "from scipy.stats import bernoulli" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "data = bernoulli.rvs(0.5, size=1000)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAOrUlEQVR4nO3ca4xdV3nG8f9DTKAXikM8iSLbrUEYiQgJiEapK6QWMEKJqXA+JFVQaVxk1YKmFRWV2rR86PUDqVRSRUK0VoNwUIGktDQWTS+pk4gWNYFJE3IhpTFpmowcxUOTuEURlJS3H84yndgznu2Zc2Yyy/+fdHTWXnvN7Hf5jJ/Zs84+O1WFJKkvL1nrAiRJ42e4S1KHDHdJ6pDhLkkdMtwlqUMb1roAgE2bNtW2bdvWugxJWlfuueeeb1bV1EL7XhThvm3bNmZmZta6DElaV5L8x2L7XJaRpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOvSg+oSpJa2nbNX+9Zsd+7CPvmsj39cxdkjpkuEtShwx3SerQoHBP8liSB5Lcl2Sm9b0qyW1JHmnP57T+JLk+yeEk9ye5aJITkCSd7HTO3N9WVW+qqum2fQ1wqKq2A4faNsClwPb22Ad8fFzFSpKGWcmyzG7gQGsfAC6b139jjdwFbExywQqOI0k6TUPDvYC/T3JPkn2t7/yqehKgPZ/X+jcDT8z72tnW9wJJ9iWZSTIzNze3vOolSQsaep37W6rqSJLzgNuS/OspxmaBvjqpo2o/sB9genr6pP2SpOUbdOZeVUfa81Hg88DFwFPHl1va89E2fBbYOu/LtwBHxlWwJGlpS4Z7kh9K8orjbeCdwIPAQWBPG7YHuKW1DwJXtatmdgDHji/fSJJWx5BlmfOBzyc5Pv7TVfW3Sb4C3JxkL/A4cEUbfyuwCzgMPAe8b+xVS5JOaclwr6pHgTcu0P+fwM4F+gu4eizVSZKWxU+oSlKH1v1dIXu8m5skrZRn7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHRoc7knOSnJvki+07VcnuTvJI0luSnJ2639Z2z7c9m+bTOmSpMWczpn7B4GH521fC1xXVduBZ4C9rX8v8ExVvRa4ro2TJK2iQeGeZAvwLuBP23aAtwOfa0MOAJe19u62Tdu/s42XJK2SoWfufwT8GvC9tn0u8GxVPd+2Z4HNrb0ZeAKg7T/Wxr9Akn1JZpLMzM3NLbN8SdJClgz3JD8NHK2qe+Z3LzC0Buz7/46q/VU1XVXTU1NTg4qVJA2zYcCYtwDvTrILeDnwI4zO5Dcm2dDOzrcAR9r4WWArMJtkA/BK4OmxVy5JWtSSZ+5V9RtVtaWqtgFXArdX1c8CdwCXt2F7gFta+2Dbpu2/vapOOnOXJE3OSq5z/3XgQ0kOM1pTv6H13wCc2/o/BFyzshIlSadryLLM91XVncCdrf0ocPECY74NXDGG2iRJy+QnVCWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6tCS4Z7k5Um+nOSrSR5K8jut/9VJ7k7ySJKbkpzd+l/Wtg+3/dsmOwVJ0omGnLl/B3h7Vb0ReBNwSZIdwLXAdVW1HXgG2NvG7wWeqarXAte1cZKkVbRkuNfIt9rmS9ujgLcDn2v9B4DLWnt326bt35kkY6tYkrSkQWvuSc5Kch9wFLgN+AbwbFU934bMAptbezPwBEDbfww4d4HvuS/JTJKZubm5lc1CkvQCg8K9qv63qt4EbAEuBl6/0LD2vNBZep3UUbW/qqaranpqampovZKkAU7rapmqeha4E9gBbEyyoe3aAhxp7VlgK0Db/0rg6XEUK0kaZsjVMlNJNrb2DwDvAB4G7gAub8P2ALe09sG2Tdt/e1WddOYuSZqcDUsP4QLgQJKzGP0yuLmqvpDka8Bnk/w+cC9wQxt/A/CpJIcZnbFfOYG6JUmnsGS4V9X9wJsX6H+U0fr7if3fBq4YS3WSpGXxE6qS1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHVoyXBPsjXJHUkeTvJQkg+2/lcluS3JI+35nNafJNcnOZzk/iQXTXoSkqQXGnLm/jzwq1X1emAHcHWSC4FrgENVtR041LYBLgW2t8c+4ONjr1qSdEpLhntVPVlV/9La/w08DGwGdgMH2rADwGWtvRu4sUbuAjYmuWDslUuSFnVaa+5JtgFvBu4Gzq+qJ2H0CwA4rw3bDDwx78tmW9+J32tfkpkkM3Nzc6dfuSRpUYPDPckPA38B/EpV/dephi7QVyd1VO2vqumqmp6amhpahiRpgEHhnuSljIL9z6rqL1v3U8eXW9rz0dY/C2yd9+VbgCPjKVeSNMSQq2UC3AA8XFUfnbfrILCntfcAt8zrv6pdNbMDOHZ8+UaStDo2DBjzFuDngAeS3Nf6fhP4CHBzkr3A48AVbd+twC7gMPAc8L6xVixJWtKS4V5V/8TC6+gAOxcYX8DVK6xLkrQCfkJVkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOLRnuST6R5GiSB+f1vSrJbUkeac/ntP4kuT7J4ST3J7loksVLkhY25Mz9k8AlJ/RdAxyqqu3AobYNcCmwvT32AR8fT5mSpNOxZLhX1ReBp0/o3g0caO0DwGXz+m+skbuAjUkuGFexkqRhlrvmfn5VPQnQns9r/ZuBJ+aNm219J0myL8lMkpm5ubllliFJWsi431DNAn210MCq2l9V01U1PTU1NeYyJOnMttxwf+r4ckt7Ptr6Z4Gt88ZtAY4svzxJ0nIsN9wPAntaew9wy7z+q9pVMzuAY8eXbyRJq2fDUgOSfAZ4K7ApySzwW8BHgJuT7AUeB65ow28FdgGHgeeA902gZknSEpYM96p6zyK7di4wtoCrV1qUJGll/ISqJHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdmki4J7kkydeTHE5yzSSOIUla3NjDPclZwMeAS4ELgfckuXDcx5EkLW4SZ+4XA4er6tGq+h/gs8DuCRxHkrSIDRP4npuBJ+ZtzwI/fuKgJPuAfW3zW0m+vszjbQK+ucyvXZFcuxZHBdZwzmvIOZ8Zzrg559oVzfnHFtsxiXDPAn11UkfVfmD/ig+WzFTV9Eq/z3rinM8MzvnMMKk5T2JZZhbYOm97C3BkAseRJC1iEuH+FWB7klcnORu4Ejg4geNIkhYx9mWZqno+yS8BfwecBXyiqh4a93HmWfHSzjrknM8MzvnMMJE5p+qk5XBJ0jrnJ1QlqUOGuyR1aN2E+1K3NEjysiQ3tf13J9m2+lWO14A5fyjJ15Lcn+RQkkWveV0vht66IsnlSSrJur9sbsick/xMe60fSvLp1a5x3Ab8bP9okjuS3Nt+vnetRZ3jkuQTSY4meXCR/Ulyffv3uD/JRSs+aFW96B+M3pj9BvAa4Gzgq8CFJ4z5ReCPW/tK4Ka1rnsV5vw24Adb+wNnwpzbuFcAXwTuAqbXuu5VeJ23A/cC57Tt89a67lWY837gA619IfDYWte9wjn/JHAR8OAi+3cBf8Poc0I7gLtXesz1cuY+5JYGu4EDrf05YGeShT5QtV4sOeequqOqnmubdzH6TMF6NvTWFb8H/AHw7dUsbkKGzPkXgI9V1TMAVXV0lWsctyFzLuBHWvuVrPPPylTVF4GnTzFkN3BjjdwFbExywUqOuV7CfaFbGmxebExVPQ8cA85dleomY8ic59vL6Df/erbknJO8GdhaVV9YzcImaMjr/DrgdUm+lOSuJJesWnWTMWTOvw28N8kscCvwy6tT2po53f/vS5rE7QcmYcgtDQbd9mAdGTyfJO8FpoGfmmhFk3fKOSd5CXAd8POrVdAqGPI6b2C0NPNWRn+d/WOSN1TVsxOubVKGzPk9wCer6g+T/ATwqTbn702+vDUx9vxaL2fuQ25p8P0xSTYw+lPuVH8GvdgNuo1DkncAHwbeXVXfWaXaJmWpOb8CeANwZ5LHGK1NHlznb6oO/dm+paq+W1X/DnydUdivV0PmvBe4GaCq/hl4OaObivVq7LdtWS/hPuSWBgeBPa19OXB7tXcq1qkl59yWKP6EUbCv93VYWGLOVXWsqjZV1baq2sbofYZ3V9XM2pQ7FkN+tv+K0ZvnJNnEaJnm0VWtcryGzPlxYCdAktczCve5Va1ydR0ErmpXzewAjlXVkyv6jmv9LvJpvNu8C/g3Ru+yf7j1/S6j/9wwevH/HDgMfBl4zVrXvApz/gfgKeC+9ji41jVPes4njL2TdX61zMDXOcBHga8BDwBXrnXNqzDnC4EvMbqS5j7gnWtd8wrn+xngSeC7jM7S9wLvB94/7zX+WPv3eGAcP9fefkCSOrRelmUkSafBcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkd+j+LtQ7VBICcpAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.hist(data);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Gamma Distribution" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "The Gamma Distribution is continuous probably distribution defined by two parameters.\n", "
\n", "\n", "
\n", "Gamma distribution on\n", "wikipedia.\n", "
\n", "\n", "Given different parameters, gamma distributions can look quite different. Explore different parameters. \n", "\n", "The exponential distribution is technically a special case of the Gamma Distribution, but is also implemented separately in scipy as 'expon'." ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "from scipy.stats import gamma" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "data = gamma.rvs(a=1, size=100000)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAD4CAYAAAAO9oqkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAPCUlEQVR4nO3df6zd9V3H8edLuul+hiKFYNt40TRzSBxjDVRJzByuFFhWTCSB6LiZJDULU2aWuKJ/1GzOdFGnI05M3TpKRJCwEZoBKzd1yWIyZi8b4ce62YZVuGulF8sYSiKib/84n5scy7m9557be889u89HcnLOeZ/P93vep2nvq5/P93u+N1WFJGll+7FhNyBJGj7DQJJkGEiSDANJEoaBJAlYNewGBnX22WfX2NjYsNuQpJHy6KOPPl9Va06uj2wYjI2NMTk5Oew2JGmkJPnXXnWXiSRJhoEkyTCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CSxAh/A3khxrY/0LN+ZOfVS9yJJC0PzgwkSYaBJMkwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEkYBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEYSJIwDCRJGAaSJAwDSRJ9hEGS9Um+muRgkqeS3NzqZyWZSHKo3a9u9SS5NcnhJI8nubhrX+Nt/KEk4131dyV5om1za5IsxoeVJPXWz8zgVeCjVfV2YBNwU5ILgO3A/qraAOxvzwGuBDa02zbgNuiEB7ADuBS4BNgxEyBtzLau7bYs/KNJkvo1ZxhU1bGq+mZ7/BJwEFgLbAX2tGF7gGva463AHdXxCHBmkvOAK4CJqjpRVS8AE8CW9tpbq+rrVVXAHV37kiQtgXkdM0gyBrwT+AZwblUdg05gAOe0YWuBZ7s2m2q1U9WnetR7vf+2JJNJJqenp+fTuiTpFPoOgyRvBr4IfKSqfniqoT1qNUD9tcWqXVW1sao2rlmzZq6WJUl96isMkryOThDcWVVfauXn2hIP7f54q08B67s2XwccnaO+rkddkrRE+jmbKMDngYNV9emul/YCM2cEjQP3d9VvaGcVbQJebMtI+4DNSVa3A8ebgX3ttZeSbGrvdUPXviRJS2BVH2MuAz4APJHksVb7A2AncE+SG4FngGvbaw8CVwGHgZeBDwJU1YkknwAOtHEfr6oT7fGHgNuBNwAPtZskaYnMGQZV9U/0XtcHuLzH+AJummVfu4HdPeqTwIVz9SJJWhx+A1mSZBhIkgwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEmiv99nsGKMbX+gZ/3IzquXuBNJWlrODCRJhoEkyTCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEn0EQZJdic5nuTJrtofJfl+ksfa7aqu125JcjjJd5Nc0VXf0mqHk2zvqp+f5BtJDiX5hySvP50fUJI0t35mBrcDW3rU/6KqLmq3BwGSXABcB/x82+avk5yR5Azgs8CVwAXA9W0swKfavjYALwA3LuQDSZLmb84wqKqvASf63N9W4O6q+q+q+h5wGLik3Q5X1dNV9QpwN7A1SYD3APe27fcA18zzM0iSFmghxww+nOTxtoy0utXWAs92jZlqtdnqPwn8oKpePaneU5JtSSaTTE5PTy+gdUlSt0HD4DbgZ4GLgGPAn7d6eoytAeo9VdWuqtpYVRvXrFkzv44lSbNaNchGVfXczOMkfwt8uT2dAtZ3DV0HHG2Pe9WfB85MsqrNDrrHS5KWyEAzgyTndT39NWDmTKO9wHVJfjzJ+cAG4J+BA8CGdubQ6+kcZN5bVQV8Ffj1tv04cP8gPUmSBjfnzCDJXcC7gbOTTAE7gHcnuYjOks4R4LcBquqpJPcA3wZeBW6qqv9p+/kwsA84A9hdVU+1t/gYcHeSPwa+BXz+tH06SVJf5gyDqrq+R3nWH9hV9Ungkz3qDwIP9qg/TedsI0nSkPgNZEmSYSBJMgwkSRgGkiQMA0kShoEkiQG/gbzSjG1/oGf9yM6rl7gTSVoczgwkSYaBJMkwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEkYBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgSQJWDbuBUTa2/YGe9SM7r17iTiRpYZwZSJIMA0lSH2GQZHeS40me7KqdlWQiyaF2v7rVk+TWJIeTPJ7k4q5txtv4Q0nGu+rvSvJE2+bWJDndH1KSdGr9zAxuB7acVNsO7K+qDcD+9hzgSmBDu20DboNOeAA7gEuBS4AdMwHSxmzr2u7k95IkLbI5w6CqvgacOKm8FdjTHu8Brumq31EdjwBnJjkPuAKYqKoTVfUCMAFsaa+9taq+XlUF3NG1L0nSEhn0mMG5VXUMoN2f0+prgWe7xk212qnqUz3qPSXZlmQyyeT09PSArUuSTna6DyD3Wu+vAeo9VdWuqtpYVRvXrFkzYIuSpJMNGgbPtSUe2v3xVp8C1neNWwccnaO+rkddkrSEBg2DvcDMGUHjwP1d9RvaWUWbgBfbMtI+YHOS1e3A8WZgX3vtpSSb2llEN3TtS5K0ROb8BnKSu4B3A2cnmaJzVtBO4J4kNwLPANe24Q8CVwGHgZeBDwJU1YkknwAOtHEfr6qZg9IfonPG0huAh9pNkrSE5gyDqrp+lpcu7zG2gJtm2c9uYHeP+iRw4Vx9SJIWj99AliQZBpIkw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEmijwvVaf7Gtj/Qs35k59VL3Ikk9ceZgSTJMJAkGQaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAk4S+3WVL+0htJy5UzA0mSYSBJMgwkSRgGkiQMA0kShoEkiQWGQZIjSZ5I8liSyVY7K8lEkkPtfnWrJ8mtSQ4neTzJxV37GW/jDyUZX9hHkiTN1+mYGfxKVV1UVRvb8+3A/qraAOxvzwGuBDa02zbgNuiEB7ADuBS4BNgxEyCSpKWxGMtEW4E97fEe4Jqu+h3V8QhwZpLzgCuAiao6UVUvABPAlkXoS5I0i4WGQQEPJ3k0ybZWO7eqjgG0+3NafS3wbNe2U602W/01kmxLMplkcnp6eoGtS5JmLPRyFJdV1dEk5wATSb5zirHpUatT1F9brNoF7ALYuHFjzzGjaLbLVICXqpC0NBY0M6iqo+3+OHAfnTX/59ryD+3+eBs+Bazv2nwdcPQUdUnSEhk4DJK8KclbZh4Dm4Engb3AzBlB48D97fFe4IZ2VtEm4MW2jLQP2JxkdTtwvLnVJElLZCHLROcC9yWZ2c/fV9VXkhwA7klyI/AMcG0b/yBwFXAYeBn4IEBVnUjyCeBAG/fxqjqxgL4kSfM0cBhU1dPAO3rU/x24vEe9gJtm2dduYPegvUiSFsZvIEuSDANJkmEgScIwkCRhGEiSMAwkSSz8chRaZLNdqsLLVEg6nZwZSJIMA0mSYSBJwjCQJGEYSJIwDCRJeGrpyPKUU0mnkzMDSZJhIEkyDCRJGAaSJAwDSRKGgSQJTy39keMpp5IG4cxAkmQYSJIMA0kSHjNYMTyWIOlUnBlIkgwDSZLLRCuey0eSwJmBJAnDQJKEy0SahctH0srizECS5MxA8+OMQfrRZBjotDAkpNHmMpEkyZmBFpczBmk0LJswSLIF+AxwBvC5qto55Ja0iGYLidkYHtLiWhZhkOQM4LPAe4Ep4ECSvVX17eF2puVivuEBBog0H8siDIBLgMNV9TRAkruBrYBhoIENEiCngyGkUbRcwmAt8GzX8yng0pMHJdkGbGtP/yPJdwd8v7OB5wfcdphGtW9YQb3nU4vYyfysmD/zZWQU+v7pXsXlEgbpUavXFKp2AbsW/GbJZFVtXOh+ltqo9g32Pgyj2jeMbu+j2jcsn1NLp4D1Xc/XAUeH1IskrTjLJQwOABuSnJ/k9cB1wN4h9yRJK8ayWCaqqleTfBjYR+fU0t1V9dQivuWCl5qGZFT7BnsfhlHtG0a391Htm1S9ZmlekrTCLJdlIknSEBkGkqSVFQZJtiT5bpLDSbYPu59+JVmf5KtJDiZ5KsnNw+5pPpKckeRbSb487F7mI8mZSe5N8p32Z/+Lw+6pX0l+r/1deTLJXUl+Ytg99ZJkd5LjSZ7sqp2VZCLJoXa/epg9zmaW3v+0/X15PMl9Sc4cZo/zsWLCoOuSF1cCFwDXJ7lguF317VXgo1X1dmATcNMI9Q5wM3Bw2E0M4DPAV6rq54B3MCKfIcla4HeBjVV1IZ2TMq4bblezuh3YclJtO7C/qjYA+9vz5eh2Xtv7BHBhVf0C8C/ALUvd1KBWTBjQdcmLqnoFmLnkxbJXVceq6pvt8Ut0fiitHW5X/UmyDrga+Nywe5mPJG8Ffhn4PEBVvVJVPxhuV/OyCnhDklXAG1mm39upqq8BJ04qbwX2tMd7gGuWtKk+9eq9qh6uqlfb00fofGdqJKykMOh1yYuR+IHaLckY8E7gG8PtpG9/Cfw+8L/DbmSefgaYBr7Qlrg+l+RNw26qH1X1feDPgGeAY8CLVfXwcLual3Or6hh0/iMEnDPkfgb1W8BDw26iXyspDPq65MVyluTNwBeBj1TVD4fdz1ySvA84XlWPDruXAawCLgZuq6p3Av/J8l2u+H/aGvtW4Hzgp4A3JfnN4Xa1siT5QzrLu3cOu5d+raQwGOlLXiR5HZ0guLOqvjTsfvp0GfD+JEfoLMu9J8nfDbelvk0BU1U1MwO7l044jIJfBb5XVdNV9d/Al4BfGnJP8/FckvMA2v3xIfczL0nGgfcBv1Ej9EWulRQGI3vJiyShs3Z9sKo+Pex++lVVt1TVuqoao/Pn/Y9VNRL/Q62qfwOeTfK2Vrqc0bmk+jPApiRvbH93LmdEDn43e4Hx9ngcuH+IvcxL+yVdHwPeX1UvD7uf+VgxYdAO6sxc8uIgcM8iX/LidLoM+ACd/1k/1m5XDbupFeB3gDuTPA5cBPzJkPvpS5vN3At8E3iCzr/zZXmZhCR3AV8H3pZkKsmNwE7gvUkO0fmFV8vytx7O0vtfAW8BJtq/078ZapPz4OUoJEkrZ2YgSZqdYSBJMgwkSYaBJAnDQJKEYSBJwjCQJAH/B44YIEWUWvlYAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.hist(data, 50);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Beta Distribution " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "The Beta Distribution is a distribution defined on the interval [0, 1], defined by two shape parameters. \n", "
\n", "\n", "
\n", "Beta distribution on \n", "wikipedia.\n", "
" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "from scipy.stats import beta" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "data = beta.rvs(1, 1, size=1000)" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD5CAYAAAA+0W6bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAOH0lEQVR4nO3db6xkdX3H8fdHQG0rLdC9kA2wXW2wYWNSIDeEhsSiqEGaCCa0gUTcppuuWmk0+oToA0ntA9oUSJoQ7BoI20ZBilo2FttSCqEYoV1gCwsbCtKtXdmwSxCENFWBbx/MWbNe7t05e+fPnd/e9yuZzJkzZ2a+v3vu/fDbM99zSFUhSWrPm1a6AEnS8hjgktQoA1ySGmWAS1KjDHBJapQBLkmNOnrYBkneCtwHvKXb/vaq+kKStwO3AicADwOXV9VPDvVea9asqfXr149ctCStJg899NDzVTW3cP3QAAd+DLy3ql5Jcgxwf5JvA58BrquqW5N8CdgE3HCoN1q/fj3bt29fRvmStHol+e/F1g89hFIDr3QPj+luBbwXuL1bvxW4eAx1SpJ66nUMPMlRSXYA+4C7gO8BL1bVq90me4CTJ1OiJGkxvQK8ql6rqjOAU4CzgdMX22yx1ybZnGR7ku379+9ffqWSpJ9zWF0oVfUicC9wDnBckgPH0E8Bnl3iNVuqar6q5ufm3nAMXpK0TEMDPMlckuO65V8A3gfsAu4BLuk22wjcMakiJUlv1KcLZS2wNclRDAL/tqr6VpIngFuT/CnwCHDjBOuUJC0wNMCr6lHgzEXWP8PgeLgkaQV4JqYkNcoAl6RG9TkGPhPWX/n3i67fffXvTLkSSZoNzsAlqVEGuCQ1ygCXpEYZ4JLUKANckhplgEtSo5ppI5RmzaRbW4/k1tkjeWzT5AxckhplgEtSowxwSWqUAS5JjTLAJalRBrgkNco2QmmVspVvuFn/GTkDl6RGGeCS1CgDXJIaZYBLUqMMcElqlAEuSY2yjVDSzJj1tr1Z4wxckhplgEtSowxwSWrU0ABPcmqSe5LsSvJ4kk91669K8oMkO7rbhZMvV5J0QJ8vMV8FPltVDyc5FngoyV3dc9dV1V9MrjxJ0lKGBnhV7QX2dssvJ9kFnDzpwiRJh3ZYbYRJ1gNnAg8C5wJXJPkosJ3BLP2Hi7xmM7AZYN26dSOWq2lYqpULbOeSZknvLzGTvA34OvDpqvoRcAPw68AZDGbo1yz2uqraUlXzVTU/Nzc3hpIlSdAzwJMcwyC8v1JV3wCoqueq6rWqeh34MnD25MqUJC3UpwslwI3Arqq69qD1aw/a7MPAzvGXJ0laSp9j4OcClwOPJdnRrfsccFmSM4ACdgMfm0iFkqRF9elCuR/IIk/dOf5yJEl9eSamJDXKqxFOgFdUW93c/7Nr0vtm2i24zsAlqVEGuCQ1ygCXpEYZ4JLUKANckhplgEtSowxwSWqUfeCaKfZQaxwO1Y99JHEGLkmNMsAlqVEGuCQ1ygCXpEYZ4JLUKANckhplG6HEkd2+uFpa6lYjZ+CS1CgDXJIaZYBLUqMMcElqlAEuSY0ywCWpUc23ER5u+9e0/6/RfT77SGhV08qzXXD1cQYuSY0ywCWpUQa4JDVqaIAnOTXJPUl2JXk8yae69SckuSvJU9398ZMvV5J0QJ8Z+KvAZ6vqdOAc4JNJNgBXAndX1WnA3d1jSdKUDA3wqtpbVQ93yy8Du4CTgYuArd1mW4GLJ1WkJOmNDquNMMl64EzgQeCkqtoLg5BPcuISr9kMbAZYt27dKLVOnG1+kvqYlZbN3l9iJnkb8HXg01X1o76vq6otVTVfVfNzc3PLqVGStIheAZ7kGAbh/ZWq+ka3+rkka7vn1wL7JlOiJGkxfbpQAtwI7Kqqaw96ahuwsVveCNwx/vIkSUvpcwz8XOBy4LEkO7p1nwOuBm5Lsgn4PvC7kylRkrSYoQFeVfcDWeLp88dbjiSpL8/ElKRGNX81wiPZuNoax9keebjtU623YM5Ku9gsWM7PovX9P+ucgUtSowxwSWqUAS5JjTLAJalRBrgkNcoAl6RG2UY4AlvMhpu1VshZ3GezWFMrVvvPzhm4JDXKAJekRhngktQoA1ySGmWAS1KjDHBJatQR20bYUnvR4dZ6JLfUjcssjm0Wa1LbnIFLUqMMcElqlAEuSY0ywCWpUQa4JDXKAJekRhngktSoI7YPfJzs351d7pvVwf28OGfgktQoA1ySGjU0wJPclGRfkp0HrbsqyQ+S7OhuF062TEnSQn1m4DcDFyyy/rqqOqO73TnesiRJwwwN8Kq6D3hhCrVIkg7DKMfAr0jyaHeI5fixVSRJ6mW5bYQ3AF8Eqru/BviDxTZMshnYDLBu3bplfpykA1pqqWup1hYtawZeVc9V1WtV9TrwZeDsQ2y7parmq2p+bm5uuXVKkhZYVoAnWXvQww8DO5faVpI0GUMPoSS5BTgPWJNkD/AF4LwkZzA4hLIb+NgEa5QkLWJogFfVZYusvnECtUiSDoNnYkpSowxwSWqUVyPUirC9TBqdM3BJapQBLkmNMsAlqVEGuCQ1ygCXpEYZ4JLUKNsIjyC25knL0+rfjjNwSWqUAS5JjTLAJalRBrgkNcoAl6RGGeCS1CgDXJIaZYBLUqMMcElqlAEuSY0ywCWpUQa4JDXKAJekRnk1Qkk/p9Ur861GzsAlqVEGuCQ1ygCXpEYNDfAkNyXZl2TnQetOSHJXkqe6++MnW6YkaaE+M/CbgQsWrLsSuLuqTgPu7h5LkqZoaIBX1X3ACwtWXwRs7Za3AhePuS5J0hDLPQZ+UlXtBejuT1xqwySbk2xPsn3//v3L/DhJ0kIT/xKzqrZU1XxVzc/NzU364yRp1VhugD+XZC1Ad79vfCVJkvpYboBvAzZ2yxuBO8ZTjiSprz5thLcA3wV+I8meJJuAq4H3J3kKeH/3WJI0RUOvhVJVly3x1PljrkWSdBg8E1OSGmWAS1KjDHBJapQBLkmNMsAlqVEGuCQ1ygCXpEYZ4JLUKANckhplgEtSowxwSWqUAS5JjTLAJalRBrgkNcoAl6RGGeCS1CgDXJIaZYBLUqMMcElqlAEuSY0ywCWpUQa4JDXKAJekRhngktQoA1ySGmWAS1KjDHBJatTRo7w4yW7gZeA14NWqmh9HUZKk4UYK8M57qur5MbyPJOkweAhFkho1aoAX8E9JHkqyebENkmxOsj3J9v3794/4cZKkA0YN8HOr6izgg8Ank7x74QZVtaWq5qtqfm5ubsSPkyQdMFKAV9Wz3f0+4JvA2eMoSpI03LIDPMkvJTn2wDLwAWDnuAqTJB3aKF0oJwHfTHLgfb5aVf8wlqokSUMtO8Cr6hngN8dYiyTpMNhGKEmNMsAlqVEGuCQ1ygCXpEYZ4JLUKANckhplgEtSowxwSWqUAS5JjTLAJalRBrgkNcoAl6RGGeCS1CgDXJIaZYBLUqMMcElqlAEuSY0ywCWpUQa4JDXKAJekRhngktQoA1ySGmWAS1KjDHBJapQBLkmNMsAlqVEjBXiSC5I8meTpJFeOqyhJ0nDLDvAkRwHXAx8ENgCXJdkwrsIkSYc2ygz8bODpqnqmqn4C3ApcNJ6yJEnDjBLgJwP/c9DjPd06SdIUHD3Ca7PIunrDRslmYHP38JUkTy7z89YAzy/zta1yzKuDY14F8mcjjfnXFls5SoDvAU496PEpwLMLN6qqLcCWET4HgCTbq2p+1PdpiWNeHRzz6jCJMY9yCOXfgdOSvD3Jm4FLgW3jKUuSNMyyZ+BV9WqSK4B/BI4Cbqqqx8dWmSTpkEY5hEJV3QncOaZahhn5MEyDHPPq4JhXh7GPOVVv+N5RktQAT6WXpEbNXIAPOz0/yVuSfK17/sEk66df5Xj1GPNnkjyR5NEkdydZtKWoJX0vw5DkkiSVpPmOhT5jTvJ73b5+PMlXp13juPX43V6X5J4kj3S/3xeuRJ3jkuSmJPuS7Fzi+ST5y+7n8WiSs0b6wKqamRuDL0O/B7wDeDPwH8CGBdv8EfClbvlS4GsrXfcUxvwe4Be75U+shjF32x0L3Ac8AMyvdN1T2M+nAY8Ax3ePT1zpuqcw5i3AJ7rlDcDula57xDG/GzgL2LnE8xcC32ZwHs05wIOjfN6szcD7nJ5/EbC1W74dOD/JYicVtWLomKvqnqr63+7hAwx67lvW9zIMXwT+HPi/aRY3IX3G/IfA9VX1Q4Cq2jflGsetz5gL+OVu+VdY5FySllTVfcALh9jkIuCva+AB4Lgka5f7ebMW4H1Oz//ZNlX1KvAS8KtTqW4yDveSBJsY/Be8ZUPHnORM4NSq+tY0C5ugPvv5ncA7k3wnyQNJLphadZPRZ8xXAR9JsodBR9sfT6e0FTPWS5CM1EY4AX1Oz+91Cn9Deo8nyUeAeeC3J1rR5B1yzEneBFwH/P60CpqCPvv5aAaHUc5j8K+sf03yrqp6ccK1TUqfMV8G3FxV1yT5LeBvujG/PvnyVsRY82vWZuB9Ts//2TZJjmbwz65D/ZNl1vW6JEGS9wGfBz5UVT+eUm2TMmzMxwLvAu5NspvBscJtjX+R2fd3+46q+mlV/RfwJINAb1WfMW8CbgOoqu8Cb2VwnZQjVa+/975mLcD7nJ6/DdjYLV8C/Et13w40auiYu8MJf8UgvFs/LgpDxlxVL1XVmqpaX1XrGRz3/1BVbV+Zcseiz+/23zH4wpokaxgcUnlmqlWOV58xfx84HyDJ6QwCfP9Uq5yubcBHu26Uc4CXqmrvst9tpb+1XeJb2v9k8O3157t1f8LgDxgGO/hvgaeBfwPesdI1T2HM/ww8B+zobttWuuZJj3nBtvfSeBdKz/0c4FrgCeAx4NKVrnkKY94AfIdBh8oO4AMrXfOI470F2Av8lMFsexPwceDjB+3j67ufx2Oj/l57JqYkNWrWDqFIknoywCWpUQa4JDXKAJekRhngktQoA1ySGmWAS1KjDHBJatT/A1q13sKbcOdlAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.hist(data, 50);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Poisson Distribution" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "The Poisson Distribution that models events in fixed intervals of time, given a known average rate (and independent occurences).\n", "
\n", "\n", "
\n", "Poisson distribution on\n", "wikipedia.\n", "
" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "from scipy.stats import poisson" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "data = poisson.rvs(mu=1, size=100000)" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAD4CAYAAAAO9oqkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAATgUlEQVR4nO3df4xd5X3n8fendmhIutQmDIja1pqqVhqKFH6MwF2kqht3jQ1VzB9FMtrWI2RpVoh0k1Wl1uk/VqGRiFQ1raUUyQoudjcb6qWJsBIT13IS9R8gHgKFgIM8JRTP2sXTjiG0qGFJv/vHfby9su947ozHc/3j/ZKuzjnf85xznyMQnznPee4hVYUk6dL2U4PugCRp8AwDSZJhIEkyDCRJGAaSJGDxoDswV1dddVWtXLly0N2QpAvGc889949VNdRr3wUbBitXrmRsbGzQ3ZCkC0aSv59un8NEkiTDQJJkGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEniAv4F8tlYueUbA/ne1x++ayDfK0kz8c5AkmQYSJIMA0kShoEkCcNAkkQfYZDko0le6Pr8KMlnklyZZH+Sw225tLVPkm1JxpO8mOTmrnONtPaHk4x01W9J8lI7ZluSnJvLlST1MmMYVNWrVXVjVd0I3AK8C3wN2AIcqKpVwIG2DbAeWNU+o8AjAEmuBLYCtwG3AltPBkhrM9p13Lp5uTpJUl9mO0y0Bvi7qvp7YAOws9V3Ane39Q3Arup4BliS5FrgDmB/VU1V1QlgP7Cu7buiqp6uqgJ2dZ1LkrQAZhsGG4GvtPVrquoYQFte3erLgCNdx0y02pnqEz3qkqQF0ncYJLkM+CTwv2dq2qNWc6j36sNokrEkY5OTkzN0Q5LUr9ncGawHvldVb7btN9sQD215vNUngBVdxy0Hjs5QX96jfpqq2l5Vw1U1PDQ0NIuuS5LOZDZhcC//PkQEsAc4OSNoBHiyq76pzSpaDbzdhpH2AWuTLG0PjtcC+9q+d5KsbrOINnWdS5K0APp6UV2SDwH/BfhvXeWHgd1JNgNvAPe0+l7gTmCczsyj+wCqairJQ8DB1u7Bqppq6/cDjwGXA0+1jyRpgfQVBlX1LvCRU2r/RGd20altC3hgmvPsAHb0qI8BN/TTF0nS/PMXyJIkw0CSZBhIkjAMJEkYBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSfQZBkmWJHkiyQ+SHEryy0muTLI/yeG2XNraJsm2JONJXkxyc9d5Rlr7w0lGuuq3JHmpHbMtSeb/UiVJ0+n3zuBPgW9W1S8CHwcOAVuAA1W1CjjQtgHWA6vaZxR4BCDJlcBW4DbgVmDryQBpbUa7jlt3dpclSZqNGcMgyRXArwCPAlTVe1X1FrAB2Nma7QTubusbgF3V8QywJMm1wB3A/qqaqqoTwH5gXdt3RVU9XVUF7Oo6lyRpAfRzZ/DzwCTw50meT/KlJB8GrqmqYwBteXVrvww40nX8RKudqT7Ro36aJKNJxpKMTU5O9tF1SVI/+gmDxcDNwCNVdRPwL/z7kFAvvcb7aw7104tV26tquKqGh4aGztxrSVLf+gmDCWCiqp5t20/QCYc32xAPbXm8q/2KruOXA0dnqC/vUZckLZAZw6Cq/gE4kuSjrbQGeAXYA5ycETQCPNnW9wCb2qyi1cDbbRhpH7A2ydL24HgtsK/teyfJ6jaLaFPXuSRJC2Bxn+1+G/hyksuA14D76ATJ7iSbgTeAe1rbvcCdwDjwbmtLVU0leQg42No9WFVTbf1+4DHgcuCp9pEkLZC+wqCqXgCGe+xa06NtAQ9Mc54dwI4e9THghn76Ikmaf/4CWZJkGEiSDANJEoaBJAnDQJKEYSBJwjCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kSfYZBkteTvJTkhSRjrXZlkv1JDrfl0lZPkm1JxpO8mOTmrvOMtPaHk4x01W9p5x9vx2a+L1SSNL3Z3Bn856q6saqG2/YW4EBVrQIOtG2A9cCq9hkFHoFOeABbgduAW4GtJwOktRntOm7dnK9IkjRrZzNMtAHY2dZ3And31XdVxzPAkiTXAncA+6tqqqpOAPuBdW3fFVX1dFUVsKvrXJKkBdBvGBTw10meSzLaatdU1TGAtry61ZcBR7qOnWi1M9UnetRPk2Q0yViSscnJyT67LkmayeI+291eVUeTXA3sT/KDM7TtNd5fc6ifXqzaDmwHGB4e7tlGkjR7fd0ZVNXRtjwOfI3OmP+bbYiHtjzemk8AK7oOXw4cnaG+vEddkrRAZgyDJB9O8h9OrgNrge8De4CTM4JGgCfb+h5gU5tVtBp4uw0j7QPWJlnaHhyvBfa1fe8kWd1mEW3qOpckaQH0M0x0DfC1NttzMfC/quqbSQ4Cu5NsBt4A7mnt9wJ3AuPAu8B9AFU1leQh4GBr92BVTbX1+4HHgMuBp9pHkrRAZgyDqnoN+HiP+j8Ba3rUC3hgmnPtAHb0qI8BN/TRX0nSOeAvkCVJhoEkyTCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkMYswSLIoyfNJvt62r0vybJLDSf4yyWWt/tNte7ztX9l1js+2+qtJ7uiqr2u18SRb5u/yJEn9mM2dwaeBQ13bnwe+UFWrgBPA5lbfDJyoql8AvtDakeR6YCPwS8A64M9awCwCvgisB64H7m1tJUkLpK8wSLIcuAv4UtsO8AngidZkJ3B3W9/Qtmn717T2G4DHq+rHVfVDYBy4tX3Gq+q1qnoPeLy1lSQtkH7vDP4E+F3g39r2R4C3qur9tj0BLGvry4AjAG3/2639/6+fcsx09dMkGU0ylmRscnKyz65LkmYyYxgk+XXgeFU9113u0bRm2Dfb+unFqu1VNVxVw0NDQ2fotSRpNhb30eZ24JNJ7gQ+CFxB505hSZLF7a//5cDR1n4CWAFMJFkM/Cww1VU/qfuY6eqSpAUw451BVX22qpZX1Uo6D4C/VVX/Ffg28But2QjwZFvf07Zp+79VVdXqG9tso+uAVcB3gYPAqjY76bL2HXvm5eokSX3p585gOr8HPJ7kD4HngUdb/VHgL5KM07kj2AhQVS8n2Q28ArwPPFBVPwFI8ilgH7AI2FFVL59FvyRJszSrMKiq7wDfaeuv0ZkJdGqbfwXumeb4zwGf61HfC+ydTV8kSfPHXyBLkgwDSZJhIEnCMJAkYRhIkjAMJEkYBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgSaKPMEjywSTfTfK3SV5O8getfl2SZ5McTvKXSS5r9Z9u2+Nt/8quc3221V9NckdXfV2rjSfZMv+XKUk6k37uDH4MfKKqPg7cCKxLshr4PPCFqloFnAA2t/abgRNV9QvAF1o7klwPbAR+CVgH/FmSRUkWAV8E1gPXA/e2tpKkBTJjGFTHP7fND7RPAZ8Anmj1ncDdbX1D26btX5Mkrf54Vf24qn4IjAO3ts94Vb1WVe8Bj7e2kqQF0tczg/YX/AvAcWA/8HfAW1X1fmsyASxr68uAIwBt/9vAR7rrpxwzXb1XP0aTjCUZm5yc7KfrkqQ+9BUGVfWTqroRWE7nL/mP9WrWlplm32zrvfqxvaqGq2p4aGho5o5Lkvoyq9lEVfUW8B1gNbAkyeK2azlwtK1PACsA2v6fBaa666ccM11dkrRA+plNNJRkSVu/HPg14BDwbeA3WrMR4Mm2vqdt0/Z/q6qq1Te22UbXAauA7wIHgVVtdtJldB4y75mPi5Mk9WfxzE24FtjZZv38FLC7qr6e5BXg8SR/CDwPPNraPwr8RZJxOncEGwGq6uUku4FXgPeBB6rqJwBJPgXsAxYBO6rq5Xm7QknSjGYMg6p6EbipR/01Os8PTq3/K3DPNOf6HPC5HvW9wN4++itJOgf8BbIkyTCQJBkGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEkYBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEYSJIwDCRJGAaSJPoIgyQrknw7yaEkLyf5dKtfmWR/ksNtubTVk2RbkvEkLya5uetcI6394SQjXfVbkrzUjtmWJOfiYiVJvfVzZ/A+8DtV9TFgNfBAkuuBLcCBqloFHGjbAOuBVe0zCjwCnfAAtgK3AbcCW08GSGsz2nXcurO/NElSv2YMg6o6VlXfa+vvAIeAZcAGYGdrthO4u61vAHZVxzPAkiTXAncA+6tqqqpOAPuBdW3fFVX1dFUVsKvrXJKkBTCrZwZJVgI3Ac8C11TVMegEBnB1a7YMONJ12ESrnak+0aPe6/tHk4wlGZucnJxN1yVJZ9B3GCT5GeCvgM9U1Y/O1LRHreZQP71Ytb2qhqtqeGhoaKYuS5L61FcYJPkAnSD4clV9tZXfbEM8tOXxVp8AVnQdvhw4OkN9eY+6JGmB9DObKMCjwKGq+uOuXXuAkzOCRoAnu+qb2qyi1cDbbRhpH7A2ydL24HgtsK/teyfJ6vZdm7rOJUlaAIv7aHM78FvAS0leaLXfBx4GdifZDLwB3NP27QXuBMaBd4H7AKpqKslDwMHW7sGqmmrr9wOPAZcDT7WPJGmBpDOB58IzPDxcY2Njczp25ZZvzHNvzn+vP3zXoLsgacCSPFdVw732+QtkSZJhIEkyDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEn0EQZJdiQ5nuT7XbUrk+xPcrgtl7Z6kmxLMp7kxSQ3dx0z0tofTjLSVb8lyUvtmG1JMt8XKUk6s37uDB4D1p1S2wIcqKpVwIG2DbAeWNU+o8Aj0AkPYCtwG3ArsPVkgLQ2o13HnfpdkqRzbMYwqKq/AaZOKW8Adrb1ncDdXfVd1fEMsCTJtcAdwP6qmqqqE8B+YF3bd0VVPV1VBezqOpckaYHM9ZnBNVV1DKAtr271ZcCRrnYTrXam+kSPuiRpAc33A+Re4/01h3rvkyejScaSjE1OTs6xi5KkU801DN5sQzy05fFWnwBWdLVbDhydob68R72nqtpeVcNVNTw0NDTHrkuSTjXXMNgDnJwRNAI82VXf1GYVrQbebsNI+4C1SZa2B8drgX1t3ztJVrdZRJu6ziVJWiCLZ2qQ5CvArwJXJZmgMyvoYWB3ks3AG8A9rfle4E5gHHgXuA+gqqaSPAQcbO0erKqTD6XvpzNj6XLgqfaRJC2gGcOgqu6dZteaHm0LeGCa8+wAdvSojwE3zNQPSdK54y+QJUmGgSTJMJAkYRhIkjAMJEkYBpIk+phaqovDyi3fGMj3vv7wXQP5Xkmz452BJMkwkCQZBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEYSJIwDCRJGAaSJHw3kc6xQb0TCXwvkjQb3hlIks6fMEiyLsmrScaTbBl0fyTpUnJehEGSRcAXgfXA9cC9Sa4fbK8k6dJxvjwzuBUYr6rXAJI8DmwAXhlor3RB8//hIPXvfAmDZcCRru0J4LZTGyUZBUbb5j8neXWO33cV8I9zPPZ857UNWD4/50MviOubI6/t/PAfp9txvoRBetTqtELVdmD7WX9ZMlZVw2d7nvOR13bhupivz2s7/50Xzwzo3Ams6NpeDhwdUF8k6ZJzvoTBQWBVkuuSXAZsBPYMuE+SdMk4L4aJqur9JJ8C9gGLgB1V9fI5/MqzHmo6j3ltF66L+fq8tvNcqk4bmpckXWLOl2EiSdIAGQaSpEsrDC7mV14kWZHk20kOJXk5yacH3af5lmRRkueTfH3QfZlPSZYkeSLJD9o/v18edJ/mU5L/0f6d/H6SryT54KD7NFdJdiQ5nuT7XbUrk+xPcrgtlw6yj3N1yYTBJfDKi/eB36mqjwGrgQcususD+DRwaNCdOAf+FPhmVf0i8HEuomtMsgz478BwVd1AZ4LIxsH26qw8Bqw7pbYFOFBVq4ADbfuCc8mEAV2vvKiq94CTr7y4KFTVsar6Xlt/h85/UJYNtlfzJ8ly4C7gS4Puy3xKcgXwK8CjAFX1XlW9NdhezbvFwOVJFgMf4gL+DVFV/Q0wdUp5A7Czre8E7l7QTs2TSykMer3y4qL5j2W3JCuBm4BnB9uTefUnwO8C/zbojsyznwcmgT9vQ2BfSvLhQXdqvlTV/wH+CHgDOAa8XVV/PdhezbtrquoYdP4oA64ecH/m5FIKg75eeXGhS/IzwF8Bn6mqHw26P/Mhya8Dx6vquUH35RxYDNwMPFJVNwH/wgU6zNBLGz/fAFwH/Bzw4SS/OdheqZdLKQwu+ldeJPkAnSD4clV9ddD9mUe3A59M8jqd4b1PJPmfg+3SvJkAJqrq5F3cE3TC4WLxa8APq2qyqv4v8FXgPw24T/PtzSTXArTl8QH3Z04upTC4qF95kSR0xp0PVdUfD7o/86mqPltVy6tqJZ1/bt+qqovir8uq+gfgSJKPttIaLq5Xt78BrE7yofbv6BouogfkzR5gpK2PAE8OsC9zdl68jmIhDOCVFwvtduC3gJeSvNBqv19VewfYJ/Xnt4Evtz9SXgPuG3B/5k1VPZvkCeB7dGa8Pc8F/PqGJF8BfhW4KskEsBV4GNidZDOd8LtncD2cO19HIUm6pIaJJEnTMAwkSYaBJMkwkCRhGEiSMAwkSRgGkiTg/wHPuwjKXGtJvgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.hist(data);" ] } ], "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.4" } }, "nbformat": 4, "nbformat_minor": 2 }