Discrete_Vs_Continuous_Prob_Distributions.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First lets define some terms for clarity\n",
    "\n",
    "**The sample space Ω**\n",
    "The sample space is the set of all possible outcomes of the experiment,\n",
    "usually denoted by Ω. For example, two successive coin tosses have\n",
    "a sample space of {hh, tt, ht, th}, where “h” denotes “heads” and “t”\n",
    "denotes “tails”.\n",
    "\n",
    "**The event space A**\n",
    "The event space is the space of potential results of the experiment. A\n",
    "subset A of the sample space Ω is in the event space A if at the end\n",
    "of the experiment we can observe whether a particular outcome ω ∈ Ω\n",
    "is in A. The event space A is obtained by considering the collection of\n",
    "subsets of Ω, and for discrete probability distributions (Section 6.2.1)\n",
    "A is often the power set of Ω.\n",
    "\n",
    "**The probability P**\n",
    "With each event A ∈ A, we associate a number P (A) that measures the\n",
    "probability or degree of belief that the event will occur. P (A) is called\n",
    "the probability of A.\n",
    "\n",
    "The probability of a single event must lie in the interval [0, 1], and the\n",
    "total probability over all outcomes in the sample space Ω must be 1, i.e.,\n",
    "P (Ω) = 1. Given a probability space (Ω, A, P ), we want to use it to model\n",
    "some real-world phenomenon. In machine learning, we often avoid explic-\n",
    "itly referring to the probability space, but instead refer to probabilities on\n",
    "quantities of interest, which **we denote by T as the target space** and refer to elements of of T as states.\n",
    "\n",
    "The term `probability` relates is to an `event` and `probability distribution` relates is to a `random variable`.\n",
    "\n",
    "It is a convention that the term `probability mass function` refers to the `probability distribution` of a `discrete random variable` and the term `probability density function` refers to the probability function of a `continuous random variable`.\n",
    "\n",
    "### Understanding Probability Density\n",
    "\n",
    "First a quick reference on PMF, PDF and CDF\n",
    "\n",
    "![](https://imgur.com/jQ5MuCL.png)\n",
    "\n",
    "\n",
    "In order to understand the heart of modern probability, we need to extend the concept of integration from basic calculus.\n",
    "\n",
    "To begin, let us consider the following piecewise function\n",
    "\n",
    "![](https://imgur.com/igRkfyD.png)\n",
    "\n",
    "\n",
    "Applying the fundamental Riemann integration of Calculus we get\n",
    "\n",
    "![](https://imgur.com/Fp9lkG9.png)\n",
    "\n",
    "\n",
    "which has the usual interpretation as the area of the two rectangles that make up f (x).\n",
    "\n",
    "The question is given f (x) = 1, what is the set of x values for which this is true? For our example, this is true whenever x ∈ (0, 1]. So now we have a correspondence between the values of the function (namely, 1 and 2) and the sets of x values for which this is true, namely, {(0, 1]} and {(1, 2]}, respectively. To compute the integral, we simply take the function values\n",
    "(i.e., 1,2) and some way of measuring the size of the corresponding interval.\n",
    "\n",
    "Since areas can be defined by definite integrals, we can also define the probability of an event occuring within an interval [a, b] by the definite integral\n",
    "\n",
    "![](https://imgur.com/FDABhQX.png)\n",
    "\n",
    "\n",
    "where f(x) is called the probability density function (pdf).\n",
    "\n",
    "A function f(x) is called a **probability density function** if\n",
    "\n",
    "1.  f(x)≥0 for all x\n",
    "2.  The area under the graph of f(x) over all the real line is exactly 1\n",
    "3.  The probability that x is in the interval [a, b] is\n",
    "\n",
    "![](https://imgur.com/HBl73is.png)\n",
    "\n",
    "\n",
    "i.e. the area under the graph of f(x) from a to b.\n",
    "\n",
    "In the problem above, the probability density function f(x) is called a **uniform (flat) probability density function (pdf)**.\n",
    "\n",
    "> So fundamentally, what does a probability density at point 𝒙 mean?\n",
    "\n",
    "> Probability density function’s value at some specific point does not give you probability; it is a measure of how dense the distribution is around that value. It means how much probability is concentrated per unit length (**d**𝒙) near 𝒙, or how dense the probability is near 𝒙.\n",
    "\n",
    "**For discrete random variables, we look up the value of a PMF at a single point to find its probability P(𝐗=𝒙) \n",
    "For continuous random variables, we take an integral of a PDF over a certain interval** to find its probability that **X** will fall in that interval.\n",
    "\n",
    "### Discrete Random Variable\n",
    "\n",
    "#### First what is a Random Variable\n",
    "\n",
    "Given a random experiment with sample space S,a **random variable** X is a set function that assigns one and only one real number to each element s that belongs in the sample space S.\n",
    "\n",
    "The set of all possible values of the random variable X, denoted x, I am calling here as the **support**, or **space**, of X.\n",
    "\n",
    "Note that the capital letters at the end of the alphabet, such as W,X,Y, and Z typically represent the definition of the random variable. The corresponding lowercase letters, such as w,x,y, and z, represent the random variable’s possible values.\n",
    "\n",
    "#### And now what is a Discrete Random Variable\n",
    "\n",
    "By a discrete random variable, it is meant a function (or a mapping), say X, from a sample space Ω, into the set of real numbers. Symbolically, if ω ∈Ω, then X (ω ) = x, where x is a real number.\n",
    "\n",
    "A random variable X is a **discrete random variable** if:\n",
    "\n",
    "*   there are a finite number of possible outcomes of X, or\n",
    "*   there are a countably infinite number of possible outcomes of X.\n",
    "\n",
    "**A countably infinite number of possible outcomes means that there is a one-to-one correspondence between the outcomes and the set of integers.**\n",
    "\n",
    "No such one-to-one correspondence exists for an uncountably infinite number of possible outcomes.\n",
    "\n",
    "**For a value x of the set of possible outcomes of the random variable X , i.e., x ∈ T , p(x) denotes the probability that random variable X has the outcome x.**\n",
    "\n",
    "**For discrete random variables, this is written as P (X = x), which is known as the probability mass function. The pmf is often referred to as the distribution”. For continuous variables, p(x) is called the probability density function (often referred to as a density).**\n",
    "\n",
    "When we say probability distribution it may pertain to a discrete random variable or a continuous random variable, depending on the context.\n",
    "\n",
    "When the random variable is discrete, probability distribution means, how the total probability is distributed over various possible values of the random variable. Consider the experiment of tossing two unbiased coins simultaneously. Then, sample space S associated with this experiment is:\n",
    "\n",
    "S = {HH,HT,TH,TT}\n",
    "\n",
    "**If we define a random variable X as: the number of heads on this sample space S,** then we will have\n",
    "\n",
    "X(HH)=2,\n",
    "\n",
    "X(HT)=X(TH)=1,\n",
    "\n",
    "X(TT)=0\n",
    "\n",
    "X(HH)=2,\n",
    "\n",
    "X(HT)=X(TH)=1,\n",
    "\n",
    "X(TT)=0.\n",
    "\n",
    "The probability distribution of XX is then given by\n",
    "\n",
    "![](https://imgur.com/SYnwiaR.png)\n",
    "\n",
    "\n",
    "For a discrete random variable, we consider events of the type {X=x} and compute probabilities of such events to describe the distribution of the random variable.\n",
    "\n",
    "> The Probability Mass Function of a Discrete Random Variable expresses the probability of the variable being equal to each specific value in the range of all potential discrete values defi ned.The sum of these probabilities over all possible values equals 100%.\n",
    "\n",
    "In mathematical form, the probability that a discrete random variable X takes on a particular value x, that is, P(X=x), is frequently denoted f(x). The function f(x) is typically called the **probability mass function**\n",
    "\n",
    "Let X be a discrete random variable with possible values denoted x1, x2, xi, x1, x2, xi,…. The probability mass function of X, denoted p\n",
    "\n",
    "![](https://imgur.com/Ur0SOp6.png)\n",
    "\n",
    "\n",
    "The same above in more general mathematical form, the **probability mass function**, P(X=x)=f(x), of a discrete random variable X is a function that satisfies the following properties:\n",
    "\n",
    "![](https://imgur.com/QvyzSZ3.png)\n",
    "\n",
    "\n",
    "First item basically says that, for every element x in the support S, all of the probabilities must be positive. Note that if x does not belong in the support S, then f(x)=0. The second item basically says that if you add up the probabilities for all of the possible x values in the support S, then the sum must equal 1. And, the third item says to determine the probability associated with the event A, you just sum up the probabilities of the x values in A.\n",
    "\n",
    "Since f(x) is a function, it can be presented:\n",
    "\n",
    "*   in tabular form\n",
    "*   in graphical form\n",
    "*   as a formula\n",
    "\n",
    "#### Some more daily life examples of discrete random variables\n",
    "\n",
    "If a random variable can take only a finite number of discrete values, then it is\n",
    "discrete.\n",
    "\n",
    "A fair die is a small cube with a natural number from 1 to 6 engraved on each side equally spaced without repetition. The fairness means that a die is made so that its weight is equally spread and, thus, all six faces are equally likely to face when rolled. So, if rolled, the set of numbers { 1,2,3,4,5,6} is the sample space of this experiment.\n",
    "\n",
    "\n",
    "### Continuous Random Variables\n",
    "\n",
    "A **continuous random variable** differs from a discrete random variable in that it takes on an uncountably infinite number of possible outcomes.\n",
    "\n",
    "While for a **discrete random variable** X that takes on a finite or countably infinite number of possible values, we determined P(X=x) for all of the possible values of X, and called it the probability mass function (“p.m.f.”). For **continuous random variables**, the probability that X takes on any particular value x is 0. That is, finding P(X=x) for a continuous random variable X is not going to work. Instead, we’ll need to find the probability that X falls in some interval (a,b), that is, we’ll need to find P(a<X<b). We’ll do that using a probability density function (“p.d.f.”).\n",
    "\n",
    "> The Probability Density Function of a Continuous Random Variable expresses\n",
    "> the rate of change in the probability distribution over the range of potential continuous values defined, and expresses the relative likelihood of getting one value in comparison with another.\n",
    "\n",
    "A nondiscrete random variable X is said to be absolutely continuous, or simply continuous, if its distribution function may be represented as\n",
    "\n",
    "![](https://imgur.com/CRGhMLE.png)\n",
    "\n",
    "\n",
    "It follows from the above that if X is a continuous random variable, then the probability that X takes on any one particular value is zero, whereas the interval probability that X lies between two different values, say, a and b,\n",
    "is given by\n",
    "\n",
    "![](https://imgur.com/P4MMU8a.png)\n",
    "\n",
    "\n",
    "A function $f(x)$ that satisfies the above requirements is called a probability function or probability distribution for a continuous random variable, but it is more often called a probability density function or simply density function. Any function f (x) satisfying Properties 1 and 2 above will automatically be a density function, and required probabilities can then be obtained from the more general form below\n",
    "\n",
    "---\n",
    "\n",
    "### Probability Density Function\n",
    "\n",
    "A function f : RD → R is called a probability density function (pdf ) if\n",
    "1. ∀x ∈ RD : f (x) > 0\n",
    "2. Its integral exists and\n",
    "\n",
    "![](https://imgur.com/b2ES6UE.png)\n",
    "\n",
    "\n",
    "As you can see, the definition for the p.d.f. of a continuous random variable differs from the definition for the p.m.f. of a discrete random variable by simply changing the summations that appeared in the discrete case to integrals in the continuous case.\n",
    "\n",
    "Now at the start of this article we discussed how density histogram (representing frequency) is defined so that the area of each rectangle equals the relative frequency of the corresponding class, and the area of the entire histogram equals 1. That suggests then that finding the probability that a **continuous random variable X** falls in some interval of values involves finding the area under the curve f(x) sandwiched by the endpoints of the interval.\n",
    "\n",
    "So from a large sample space of Pizza, the probability that a randomly selected Pizza weighs between 0.20 and 0.30 pounds is then this area: (which is what the definite Integral formulae above calculates )\n",
    "\n",
    "![](https://imgur.com/IG9J22d.png)\n",
    "\n",
    "\n",
    "**Some examples of well known discrete probability distributions include:**\n",
    "\n",
    "*   Poisson distribution.\n",
    "*   Bernoulli and binomial distributions.\n",
    "*   Multinoulli and multinomial distributions.\n",
    "*   Discrete uniform distribution.\n",
    "*   The Geometric Distribution\n",
    "*   The Negative-Binomial Distribution\n",
    "*   The Hypergeometric Distribution\n",
    "\n",
    "Some examples of common domains with well-known discrete probability distributions include:\n",
    "\n",
    "*   The probabilities of dice rolls form a discrete uniform distribution.\n",
    "*   The probabilities of coin flips form a Bernoulli distribution.\n",
    "*   The probabilities car colors form a multinomial distribution.\n",
    "\n",
    "---\n",
    "\n",
    "### A quick summary\n",
    "\n",
    "![](https://cdn-images-1.medium.com/max/800/1*51lPbnyMX69R6neu-v8aVA.png)\n",
    "\n",
    "\n",
    "Now lets see an simple actual exmaple of Discrete Probability Distribution. Quickly revisit the definition\n",
    "\n",
    "The **probability distribution** of a discrete random variable X is a list of each possible value of X together with the probability that X takes that value in one trial of the experiment.\n",
    "\n",
    "I start with a simple experiment, tossing a fair coin 10 times, and measured how many successes/heads I observe. I can use the number of successes (heads) observed in many ways to understand the basics of probability. For example, I could simply count how many times we see 0 heads, 1 head, 2 heads with our fair coin toss, and so on. Or here, I am just denoting the outcome with ‘H’ or ‘T’ for each experiment.\n",
    "\n",
    "Now a quick and simple Math example of PDF\n",
    "\n",
    "Let X be a continuous random variable whose probability density function is:\n",
    "\n",
    "![](https://imgur.com/gTZBX0F.png)\n",
    "\n",
    "\n",
    "which is clearly not a probability! In the continuous case, f(x) is instead the height of the curve at X=x, so that the total area under the curve is 1. In the continuous case, it is areas under the curve that define the probabilities.\n",
    "\n",
    "**What is P(X=1/2)?**\n",
    "\n",
    "It is a straightforward integration to see that the probability is 0:\n",
    "\n",
    "![](https://imgur.com/Fv6Oefb.png)\n",
    "\n",
    "\n",
    "In general, if X is continuous, the probability that X takes on any specific value x is 0. That is, when X is continuous, P(X=x)=0 for all x in the support.\n",
    "\n",
    "An implication of the fact that **P(X=x)=0 for all x when X is continuous** is that you can be less precise about the endpoints of intervals when finding probabilities of continuous random variables. That is:\n",
    "\n",
    "![](https://imgur.com/s5beiPR.png)\n",
    "\n",
    "\n",
    "for any constants a and b.\n",
    "\n",
    "**Further explanation of the above principle**\n",
    "\n",
    "The probability of observing any single value of the continuous random variable is 0 since the number of possible outcomes of a continuous random variable is uncountable and infinite. That is, for a continuous random variable, we must calculate a probability over an interval rather than at a particular point. This is why the probability for a continuous random variable can be interpreted as an area under the curve on an interval. In other words, we cannot describe the probability distribution of a continuous random variable by giving probability of single values of the random variable as we did for a discrete random variable. This property can also be seen from the fact that\n",
    "\n",
    "![](https://imgur.com/VxmYRGV.png)\n",
    "\n",
    "\n",
    "for any real **c**\n",
    "\n",
    "### Why do I need to Integrate over the PDF to get the Probability\n",
    "\n",
    "In the case of of continuous random variable, we should not ask for the probability that **X** is exactly a single number (since that probability is zero). Instead, we need to think about the probability that **x** is close to a single number.\n",
    "\n",
    "We capture the notion of being close to a number with a probability density function which is normally denoted by **P(x)**. If the probability density around a point **x** is large, that means the random variable **X** is likely to be close to **x**. If, on the other hand, **P(x)=0** in some interval, then **X** won’t be in that interval.\n",
    "\n",
    "So building on the Integration concept of Calculus\n",
    "\n",
    "If the probability of **X** being exactly at point 𝒙 is zero, how about an extremely small interval around the point **𝒙**? Say, **[𝒙, 𝒙+d𝒙]?**\n",
    "\n",
    "Let’s assume **d𝒙** is infinitesimally small with a value of 0.00000000001.\n",
    "\n",
    "Then the probability that **X** will fall in **[𝒙, 𝒙+d𝒙]** is the **Area** under the curve **f(𝒙)** sandwiched by **[𝒙, 𝒙+d𝒙]**.\n",
    "\n",
    "**The Area Under a Curve — Integral Calculus Basics**\n",
    "\n",
    "The area under a curve between two points can be found by doing a definite integral between the two points. To find the area under the curve y = f(x) between x = a and x = b, integrate y = f(x) between the limits of a and b.\n",
    "\n",
    "To translate the probability density P(x) into a probability, imagine that Ix is some small interval around the point x. Then, assuming P is continuous, the probability that X is in that interval will depend both on the density P(x) and the length of the interval\n",
    "\n",
    "P ( X ∈ Ix) ≈ P (x) × Length of Ix\n",
    "\n",
    "We don’t have a true equality here, because the density P may vary over the interval Ix. But, the approximation becomes better and better as the interval Ix shrinks around the point x, as P will be come closer and closer to a constant inside that small interval. The probability P ( X ∈ Ix ) approaches zero as Ix shrinks down to an infinitesemally small value to the point x (consistent with our above result for single numbers), but the information about X is contained in the rate that this probability goes to zero as Ix shrinks.\n",
    "\n",
    "So, to determine the probability that X is in any subset A of the real numbers, we simply add up the values of P(x) in the subset. By “add up,” we mean integrate the function P(x) over the set A.\n",
    "\n",
    "### Cumulative Distribution Function\n",
    "\n",
    "> **The Cumulative Distribution Function of a Discrete Random Variable** expresses the theoretical or observed probability of that variable being less than or equal to any given value. **It equates to the sum of the probabilities** of achieving that value and each successive lower value.\n",
    "\n",
    "#### Example of the Cumulative Distribution Function for Rolling a Single Die\n",
    "\n",
    "![](https://imgur.com/gmOzMVa.png)\n",
    "\n",
    "\n",
    "And now the same for Continuous Random Variable\n",
    "\n",
    "> **The Cumulative Distribution Function of a Continuous Random Variable**\n",
    "> expresses the theoretical or observed probability of that variable being less than or equal to any given value. **It equates to the area under the Probability Density Function curve** to the left of the value in question.\n",
    "\n",
    "Now implementing some very basic PDF with Python and Scipy\n",
    "\n",
    "Another Jupyter Notebook to understand how PDF is different from Probability."
   ]
  }
 ],
 "metadata": {
  "language_info": {
   "name": "python"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}