Common Probability Distributions in Finance › The Normal Distribution

Michael BrenndoerferOctober 21, 202514 min read
Quantitative FinanceAll chapters
Part 1 of 15Mathematical Foundations
Chapter 3Common Probability Distributions in Finance › The Normal Distribution
Reading section 2 of 8 · The Normal Distribution · 14 minSection 2 of 8 · 14 min

Book contents

Section 2 of 8 · 14 min read

Common Probability Distributions in Finance › The Normal DistributionOne section · Sections in this chapter2 / 8 ↓
Reading tools Definitions & inline links

Choose your expertise level to adjust how many terms are explained. Beginners see more tooltips, experts see fewer to maintain reading flow. Hover over underlined terms for instant definitions.

Article links

Make inline references clickable

The normal distribution, also called the Gaussian distribution, is a common baseline in financial modeling. One mathematical reason is the classical Central Limit Theorem: standardized sums of independent, identically distributed random variables with finite, strictly positive variance converge in distribution to a normal law. Thinking of a return as the aggregate effect of many small shocks can motivate a normal approximation, but market shocks need not be independent or identically distributed, so this is a modeling heuristic rather than a theorem about returns.

Normal Distribution

A continuous probability distribution characterized by its symmetric, bell-shaped curve. It is fully specified by two parameters: the mean μ\mu (center) and variance σ2\sigma^2 (spread). Also known as the Gaussian distribution.

Mathematical Properties

To understand what it means for a random variable to follow a normal distribution, we must first grasp the concept of a probability density function. Unlike discrete random variables that take on specific values with certain probabilities, a continuous random variable like one following a normal distribution can take any value along the real line. Probability over a region is the integral of the density over that region; density at one point is not itself a probability. A higher density at a point means that sufficiently small, equal-width neighborhoods around that point have greater probability than neighborhoods around a lower-density point.

A nondegenerate normal variable XX has mean μ\mu and variance σ2>0\sigma^2>0, written X∼N(μ,σ2)X \sim N(\mu, \sigma^2), with probability density function:

f(x)=1σ2πexp⁡(−(x−μ)22σ2),−∞<x<∞f(x) = \frac{1}{\sigma\sqrt{2\pi}} \exp\left(-\frac{(x-\mu)^2}{2\sigma^2}\right), \quad -\infty < x < \infty

Each component serves a specific purpose. Let's examine each component:

  • xx: the value at which we evaluate the density
  • μ\mu: the mean (center of the distribution)
  • σ\sigma: the standard deviation (controls the spread)
  • 1σ2π\frac{1}{\sigma\sqrt{2\pi}}: a normalization constant ensuring the density integrates to 1
  • (x−μ)2/(2σ2)(x-\mu)^2/(2\sigma^2): one half of the squared standardized distance from xx to the mean

The heart of the normal distribution lies in the exponential term. The quantity (x−μ)2(x-\mu)^2 measures the squared distance between the point xx and the center of the distribution μ\mu. Dividing by 2σ22\sigma^2 scales this distance relative to the spread of the distribution. A value two units away from the mean is much more unusual when σ=0.5\sigma = 0.5 than when σ=10\sigma = 10, because the standardized distance is larger. The negative sign in the exponent ensures that the function decreases as we move away from the mean, and the exponential function turns this squared distance into a probability density that decays smoothly and symmetrically.

The exponential of a negative quadratic creates the characteristic bell shape: values near the mean have high probability density, while values far from the mean have exponentially decreasing density. For a standard normal variable, the two-sided probability beyond four standard deviations is about 6.3×10−56.3 \times 10^{-5}, and beyond five it is about 5.7×10−75.7 \times 10^{-7}. Whether those benchmark probabilities adequately describe a financial return series is an empirical question tied to the asset, horizon, and sample. The quadratic term in the exponent is the key: because squaring grows faster than linear growth, the tails of the normal distribution thin out quickly.

The normalization constant 1σ2π\frac{1}{\sigma\sqrt{2\pi}} deserves special attention. This factor ensures that when we integrate the density function over all possible values from negative infinity to positive infinity, we obtain exactly 1. This is the total probability. The 2π\sqrt{2\pi} arises from the Gaussian integral, one of the most famous results in calculus, while the σ\sigma in the denominator accounts for how wider distributions (larger σ\sigma) must have lower peak heights to maintain the same total area under the curve.

Key parameters:

  • Mean (μ\mu): The center of the distribution and the expected value of XX. In finance, this represents the expected return.
  • Variance (σ2\sigma^2): Measures the spread of the distribution. The square root, σ\sigma, is the standard deviation, which in finance we call volatility.
  • Skewness: Zero for the normal distribution, meaning it is perfectly symmetric around the mean.
  • Kurtosis: Equal to 3 for the normal distribution (or excess kurtosis of 0). Kurtosis is the fourth standardized moment; by itself, it does not uniquely determine how probability is divided among the center, shoulders, and tails.

The symmetry of the normal distribution is both its elegance and its limitation. Symmetry is a modeling assumption: it assigns equal probabilities to equal-sized deviations above and below the mean. Some financial return series exhibit asymmetry. This possibility is one reason to treat the normal distribution as a baseline rather than a universal model of returns.

The standard normal distribution is the special case where μ=0\mu = 0 and σ2=1\sigma^2 = 1. Any normal variable with σ>0\sigma>0 can be standardized:

Z=X−μσ∼N(0,1)Z = \frac{X - \mu}{\sigma} \sim N(0, 1)

For σ>0\sigma>0, this process turns a normal random variable into the canonical standard normal form. The operation X−μσ\frac{X - \mu}{\sigma} first centers the distribution by subtracting the mean (so the new center is at zero), then scales it by dividing by the standard deviation (so the new spread equals one). This step preserves the essential shape of the distribution while placing it in a standard coordinate system.

This standardization is central because it allows us to use a single table of probabilities for any normal distribution. The cumulative distribution function of the standard normal, denoted Φ(z)\Phi(z), gives the probability that a standard normal variable is less than zz. Before computers, statisticians relied on printed tables of Φ(z)\Phi(z) values, and the standardization formula allowed them to answer questions about any normal distribution using this single reference. Even today, the standard normal is the universal reference point: when we speak of a "3-sigma event," we mean an outcome more than three standard deviations from the mean, regardless of what the actual mean and standard deviation are.

The 68-95-99.7 Rule

One of the most useful properties of the normal distribution is how probability mass is concentrated around the mean:

  • Approximately 68% of observations fall within one standard deviation of the mean
  • Approximately 95% fall within two standard deviations
  • Approximately 99.7% fall within three standard deviations

These percentages emerge directly from integrating the normal density function, but their approximate values (68, 95, and 99.7) are worth committing to memory because they provide instant intuition for assessing how unusual any particular observation is. If someone tells you that an event was a "two-sigma move," you immediately know that such events occur roughly 5% of the time, or about one trading day in twenty.

This rule provides quick mental estimates for risk. If daily returns are normal with a standard deviation of 1%, the probability of a deviation from their mean larger than 3 percentage points is about 0.27%. For independent daily draws, the expected waiting time is about 370 trading days, or roughly 1.5 trading years at 252 days per year. An absolute return larger than 3% has that same probability only when the mean is zero. These are model-based estimates; a different return distribution can produce very different frequencies.

Out[3]:
Visualization
Bell curve for the standard normal distribution with shaded standard-deviation probability regions.
The standard normal distribution with shaded standard-deviation bands. The central band contains about 68.3% of the probability mass, the two bands from one to two standard deviations add about 27.2% (95.4% cumulative), and the bands from two to three standard deviations add about 4.3% (99.7% cumulative).
Out[4]:
Visualization
Line plot of normal probability density functions with different means and fixed standard deviation.
Normal distributions with different means (σ fixed). Changing the mean shifts the distribution horizontally.
Line plot of normal probability density functions with fixed mean and different standard deviations.
Normal distributions with different standard deviations (μ fixed). Changing the standard deviation affects the spread and peak height.

Applications in Finance

The normal distribution appears throughout quantitative finance:

Portfolio theory

Portfolio theory often summarizes portfolios by expected return and variance. Normality is sufficient for those two moments to determine the full distribution of portfolio returns, but the mean-variance framework itself does not require normally distributed returns. When candidate portfolio returns are normal, expected utility for a fixed utility function can be evaluated from the mean and variance because those parameters determine the distribution, provided that expected utility exists and is finite.

Value at Risk (VaR)

Parametric normal VaR assumes normally distributed returns. Under that model, the signed VaR loss quantile at confidence level 0<α<10<\alpha<1 is:

VaRα=zασ−μ\text{VaR}_{\alpha} = z_{\alpha} \sigma - \mu

This formula translates a tail quantile into a loss estimate. For normally distributed returns with mean μ\mu and standard deviation σ>0\sigma>0, we seek the threshold such that returns fall below it with probability 1−α1-\alpha, or 100(1−α)%100(1-\alpha)\%. The standard normal quantile zαz_\alpha determines the lower-return threshold μ−zασ\mu-z_\alpha\sigma.

Each component of this formula plays a specific role:

  • VaRα\text{VaR}_\alpha: the signed loss quantile at confidence level α\alpha; under a continuous model, losses exceed this threshold with probability 1−α1-\alpha
  • μ\mu: the expected return over the time horizon
  • σ\sigma: the volatility (standard deviation of returns) over the time horizon
  • zαz_\alpha: the α\alpha-quantile of the standard normal distribution (e.g., z0.95≈1.645z_{0.95} \approx 1.645 for 95% confidence)

The formula works because under normality, the return quantile at probability 1−α1-\alpha is exactly μ−zασ\mu-z_\alpha\sigma. For example, the 95% VaR uses z0.95≈1.645z_{0.95}\approx1.645, and losses exceed the threshold with probability 5%. Taking the negative of this return quantile yields the signed loss quantile zασ−μz_\alpha\sigma-\mu. It is positive when zασ>μz_\alpha\sigma>\mu, zero at equality, and negative when the modeled lower-tail return quantile is still a gain. A reporting policy that disallows negative VaR can use max⁡(0,zασ−μ)\max(0,z_\alpha\sigma-\mu), but that floor is a reporting convention rather than the raw quantile.

Out[5]:
Visualization
Normal distribution with VaR threshold and shaded tail region showing Expected Shortfall.
Value at Risk (VaR) and Expected Shortfall (ES) visualized on a return distribution. With tail probability p=0.05, the plotted line is the lower-return threshold with probability p of a worse return; for these parameters, its negative is a positive 95% VaR loss threshold. ES is the average loss in that tail.

Option pricing

The Black-Scholes model assumes that log-returns are normally distributed, as we will explore in a later chapter.

In[6]:
Code
import numpy as np
from scipy import stats


def calculate_var_normal(mu, sigma, confidence_level=0.95, horizon=1):
    """
    Calculate Value at Risk assuming normal returns.

    Parameters:
        mu: Mean daily additive return increment (for example, a log return)
        sigma: Volatility (daily standard deviation)
        confidence_level: Confidence level (e.g., 0.95 for 95% VaR)
        horizon: Time horizon in days

    Returns:
        Signed VaR loss quantile; it is positive only when the modeled
        lower-tail return quantile is negative

    Raises:
        ValueError for invalid inputs; FloatingPointError if the loss
        calculation remains nonfinite after the scaled fallback

    Scaling avoids some intermediate overflow, not loss of relative accuracy
    when the mean and tail terms nearly cancel. The quantile coefficient is
    itself a double-precision approximation; extreme cancellation can even
    change the sign of a small residual. This is an instructional helper,
    not an arbitrary-range high-precision risk calculator.
    """
    if not (np.isfinite(mu) and np.isfinite(sigma) and sigma > 0):
        raise ValueError("mu must be finite and sigma must be finite and positive")
    if not (np.isfinite(confidence_level) and 0 < confidence_level < 1):
        raise ValueError("confidence_level must lie strictly between 0 and 1")
    if not (np.isfinite(horizon) and horizon > 0):
        raise ValueError("horizon must be finite and positive")

    # Use the loss quantile directly: 1-confidence can round to one.
    z = stats.norm.ppf(confidence_level)
    # Sum iid daily increments; this is not exact simple-return compounding.
    with np.errstate(over="ignore", invalid="ignore"):
        mu_scaled = mu * horizon
        sigma_scaled = sigma * np.sqrt(horizon)
        var = -(mu_scaled - z * sigma_scaled)
        if not np.isfinite(var):
            if z == 0:
                # At the median only the mean term matters, even if sigma_scaled overflows.
                var = -mu_scaled
            else:
                scale = max(abs(mu), sigma)
                var = (
                    z * (sigma / scale) * np.sqrt(horizon) - (mu / scale) * horizon
                ) * scale
    if not np.isfinite(var):
        raise FloatingPointError(
            "The normal VaR loss quantile overflows floating point"
        )
    return var


# Example: Daily return of 0.05% with 1.5% volatility
daily_return = 0.0005
daily_vol = 0.015

var_95_1day = calculate_var_normal(daily_return, daily_vol, 0.95, 1)
var_99_1day = calculate_var_normal(daily_return, daily_vol, 0.99, 1)
var_95_10day = calculate_var_normal(daily_return, daily_vol, 0.95, 10)
Out[7]:
Console
Portfolio with daily return 0.05% and volatility 1.50%
1-day 95% VaR: 2.42%
1-day 99% VaR: 3.44%
10-day 95% VaR: 7.30%

These VaR figures are loss thresholds, not maximum or expected losses. Under the assumed normal parameters, losses exceed the unrounded 1-day 95% threshold (about 2.42%) with probability 5%. The horizon calculation sums iid normal increments, such as log returns; it is not an exact model of compounded simple returns. The 10-day VaR is higher here because the volatility term grows with n\sqrt{n} and outweighs the mean term, which grows with nn. Under that additive model, VaRn=zασn−nμ\text{VaR}_n=z_\alpha\sigma\sqrt{n}-n\mu, so VaR need not increase with the horizon for every parameter choice. The square-root scaling follows from variance addition: the variance of nn independent daily increments is nσ2n\sigma^2, and their sum has standard deviation nσ\sqrt{n}\sigma.

The helper uses double-precision arithmetic. Its scaled fallback handles some overflowing intermediate terms, but it cannot promise relative accuracy when two large terms nearly cancel. Rounding in the normal quantile coefficient also matters in that regime: more careful subtraction alone cannot recover an accurate true-distribution residual. Use appropriate numerical error analysis when the answer depends on such a cancellation.