Common Probability Distributions in Finance › The Normal Distribution
Quantitative FinanceAll chapters 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.
A continuous probability distribution characterized by its symmetric, bell-shaped curve. It is fully specified by two parameters: the mean (center) and variance (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 has mean and variance , written , with probability density function:
Each component serves a specific purpose. Let's examine each component:
- : the value at which we evaluate the density
- : the mean (center of the distribution)
- : the standard deviation (controls the spread)
- : a normalization constant ensuring the density integrates to 1
- : one half of the squared standardized distance from to the mean
The heart of the normal distribution lies in the exponential term. The quantity measures the squared distance between the point and the center of the distribution . Dividing by scales this distance relative to the spread of the distribution. A value two units away from the mean is much more unusual when than when , 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 , and beyond five it is about . 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 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 arises from the Gaussian integral, one of the most famous results in calculus, while the in the denominator accounts for how wider distributions (larger ) must have lower peak heights to maintain the same total area under the curve.
Key parameters:
- Mean (): The center of the distribution and the expected value of . In finance, this represents the expected return.
- Variance (): Measures the spread of the distribution. The square root, , 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 and . Any normal variable with can be standardized:
For , this process turns a normal random variable into the canonical standard normal form. The operation 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 , gives the probability that a standard normal variable is less than . Before computers, statisticians relied on printed tables of 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.



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 is:
This formula translates a tail quantile into a loss estimate. For normally distributed returns with mean and standard deviation , we seek the threshold such that returns fall below it with probability , or . The standard normal quantile determines the lower-return threshold .
Each component of this formula plays a specific role:
- : the signed loss quantile at confidence level ; under a continuous model, losses exceed this threshold with probability
- : the expected return over the time horizon
- : the volatility (standard deviation of returns) over the time horizon
- : the -quantile of the standard normal distribution (e.g., for 95% confidence)
The formula works because under normality, the return quantile at probability is exactly . For example, the 95% VaR uses , and losses exceed the threshold with probability 5%. Taking the negative of this return quantile yields the signed loss quantile . It is positive when , 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 , but that floor is a reporting convention rather than the raw quantile.

Option pricing
The Black-Scholes model assumes that log-returns are normally distributed, as we will explore in a later chapter.
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)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 and outweighs the mean term, which grows with . Under that additive model, , so VaR need not increase with the horizon for every parameter choice. The square-root scaling follows from variance addition: the variance of independent daily increments is , and their sum has standard deviation .
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.