Normal Approximation to the Binomial Distribution and Inverse Exponential Functions

Transitioning from Discrete to Continuous: Normal Approximation to the Binomial

  • Conceptual Overview

    • In Chapter 5, the binomial distribution is defined as a discrete distribution with parameters nn (finite number of trials) and π\pi (constant probability of success).

    • Calculating cumulative probabilities for discrete distributions requires the summation of individual probabilities. For example, to find P(X50)P(X \leq 50) when n=100n = 100 and π=0.3\pi = 0.3, one must calculate:         P(X=0)+P(X=1)+P(X=2)++P(X=50)P(X = 0) + P(X = 1) + P(X = 2) + \dots + P(X = 50)

    • This is described as a "very long and tedious calculation" when performed manually.

  • The Symmetrical Nature of Large Trials

    • Plotting a binomial distribution with a large nn reveals that the discrete bars begin to resemble a bell-shaped curve.

    • Because the distribution becomes symmetrical and bell-shaped as nn increases, the normal distribution (a continuous distribution) can be used to approximate these probabilities, significantly shortening the calculation process.

Rules and Requirements for Using the Approximation

  • Minimum Threshold Checks

    • The normal approximation is not always valid. To use it, the binomial distribution must pass two tests based on a "rule of thumb" to ensure accuracy. If these values are below 5, the approximation breakdown and loses accuracy.

    • The criteria are:

      1. n×π>5n \times \pi > 5

      2. n×(1π)>5n \times (1 - \pi) > 5

  • Assessment Importance

    • In exams offering calculation marks (such as second-semester assessments), explicitly writing down these checks and verifying they hold can be worth two marks.

Formulating the Parameters for Approximation

  • Variable notation

    • To distinguish the binomial random variable from its normal replacement, the approximated variable is denoted as xax_a (x approximate).

  • Equating the Mean and Standard Deviation

    • For the approximation to be valid, the normal distribution must share the same center and spread as the original binomial distribution.

    • Mean (μ\mu): The expected value of the binomial distribution is calculated as:         μ=n×π\mu = n \times \pi

    • Standard Deviation (σ\sigma): Calculated as:         σ=n×π×(1π)\sigma = \sqrt{n \times \pi \times (1 - \pi)}

  • Standardization to Z-Values

    • To find the probability using the standard normal distribution table or a calculator, the approximated variable must be standardized:         Z=xa(n×π)n×π×(1π)Z = \frac{x_a - (n \times \pi)}{\sqrt{n \times \pi \times (1 - \pi)}}

    • The resulting ZZ variable is normal with a mean of 0 and a standard deviation of 1: ZN(0,1)Z \sim N(0, 1).

The Continuity Correction

  • Theory of Continuity Correction

    • Discrete variables represent points, while continuous variables represent the area under a curve. Moving from one to the other is not a direct "one-to-one" transition.

    • To calculate the area under the curve effectively, we "pad" the discrete value by adding or subtracting 0.50.5. This creates a small slice of area representing the discrete point.

  • Correction Scenarios

    • Exact Value: P(X=x)P(X = x) becomes P(x0.5<xa<x+0.5)P(x - 0.5 < x_a < x + 0.5).

    • Less Than or Equal To: P(Xx)P(X \leq x) becomes P(xa<x+0.5)P(x_a < x + 0.5).

    • Greater Than or Equal To: P(Xx)P(X \geq x) becomes P(xa>x0.5)P(x_a > x - 0.5).

Practical Application: Airline Overbooking Example

  • Scenario Details

    • A flight accommodates 300 passengers300 \text{ passengers}.

    • The airline accepts n=350 bookingsn = 350 \text{ bookings}, hoping no more than 300 people show up.

    • The probability that a passenger shows up for their booking is π=0.8\pi = 0.8.

    • Goal: Find the probability that more than 300 passengers show up: P(X>300)P(X > 300).

  • Step 1: Verification Checks

    • n×π=350×0.8=280n \times \pi = 350 \times 0.8 = 280 (Greater than 5)

    • n×(1π)=350×0.2=70n \times (1 - \pi) = 350 \times 0.2 = 70 (Greater than 5)

    • Conclusion: Normal approximation is valid.

  • Step 2: Calculating Parameters

    • μ=280\mu = 280

    • σ=350×0.8×0.2=56\sigma = \sqrt{350 \times 0.8 \times 0.2} = \sqrt{56}

  • Step 3: Probability Transformation and Continuity Correction

    • Original binomial target: P(X>300)P(X > 300).

    • Rewrite as cumulative: 1P(X300)1 - P(X \leq 300).

    • Apply continuity correction for the normal variable: 1P(xa<300.5)1 - P(x_a < 300.5).

  • Step 4: Standardization and Result

    • Standardize to Z: Z=300.528056=20.5562.74Z = \frac{300.5 - 280}{\sqrt{56}} = \frac{20.5}{\sqrt{56}} \approx 2.74

    • Using an HP calculator (Shift ZPZ \rightarrow P): P(Z<2.74)=0.9969P(Z < 2.74) = 0.9969

    • Final Probability: 10.9969=0.00311 - 0.9969 = 0.0031.

  • Accuracy Comparison

    • Excel exact binomial calculation (1BINOM.DIST(300,350,0.8,TRUE)1 - \text{BINOM.DIST}(300, 350, 0.8, \text{TRUE})): 0.002240.00224

    • The approximation result (0.00310.0031) is close to the exact value (0.002240.00224). Accuracy typically increases with higher nn.

Inverse Exponential Distribution: Finding Percentiles

  • Standard Probability Problem

    • In the standard exponential distribution, we know the time tt and solve for the probability α\alpha:         α=P(X<t)=1eλ×t\alpha = P(X < t) = 1 - e^{-\lambda \times t}

  • The Inverse Problem

    • One may be given a specific probability (e.g., the 80th80^{th} percentile, where α=0.8\alpha = 0.8) and asked to find the unknown point pp on the curve associated with that probability.

  • Mathematical Derivation for p

    1. Start with the identity: α=1eλ×p\alpha = 1 - e^{-\lambda \times p}

    2. Rearrange: eλ×p=1αe^{-\lambda \times p} = 1 - \alpha

    3. Take the natural logarithm (ln) of both sides to cancel ee:         λ×p=ln(1α)-\lambda \times p = \ln(1 - \alpha)

    4. Solve for the point pp:         p=1λ×ln(1α)p = -\frac{1}{\lambda} \times \ln(1 - \alpha)

  • Example Calculation

    • Given λ=2\lambda = 2 and α=0.8\alpha = 0.8:         p=12×ln(10.8)p = -\frac{1}{2} \times \ln(1 - 0.8)         p=0.5×ln(0.2)0.8047p = -0.5 \times \ln(0.2) \approx 0.8047

Questions & Discussion

  • Question: Where does the value of 5 come from in the verification checks?

  • Response: It is a standard rule of thumb developed through observation. When values are lower than 5, the model's accuracy deteriorates significantly.

  • Question: On the HP calculator, how do you perform the natural log (ln\ln)?

  • Response: On the HP calculator, the process is often the reverse of the Casio. You must input the number first and then press the function button. (The instructor noted they are still working on mastering the HP calculator interface for this specific function).

  • Question: In the Excel function BINOM.DIST, what does "True" mean?

  • Response: "True" is used for cumulative probabilities. If set to "False," it would calculate the probability of the variable being exactly equal to a specific value.