Cumulative Distribution Functions

Suppose you know a light bulb's lifetime is random. You want to answer a single practical question. What is the chance it survives at least 500 hours? To answer this, and infinitely many questions like it, you do not need to reconstruct the entire probabilistic mechanism from scratch each time. You need one function that packages every such answer at once. That function is the cumulative distribution function, or CDF. It is the single most versatile object in the theory of Random Variables. The CDF matters because it is the only description of a random variable's distribution. It works uniformly across discrete, continuous, and every other type of variable. This includes hybrid variables that are neither. Probability mass functions only make sense for variables taking countably many values. Densities (introduced later for Continuous Random Variables) require an integral to even exist. The CDF, by contrast, is defined for every real-valued random variable. There are no exceptions and no side conditions. It converts the abstract measure-theoretic content of a probability space into something concrete. It becomes a monotone real function that can be graphed, differentiated, integrated, and manipulated with ordinary calculus. The big picture is this: a random variable XX induces a probability measure on the real line. The CDF is a bookkeeping device. It records the mass this measure assigns to every half-line (−∞,x](-\infty, x]. Remarkably, this one family of numbers — indexed by x∈Rx \in \mathbb{R} — determines the entire distribution of XX. Two random variables with the same CDF are distributed identically. This holds for every practical and theoretical purpose. It holds even if they live on completely different probability spaces. This article gives the precise definition of the CDF. It proves the four properties that characterize CDFs among all functions. This means one can manufacture new distributions by writing down any function with these properties. The article shows exactly how the CDF specializes to the discrete case. There, it becomes a step function built from the PMF. It also foreshadows the continuous case, where it becomes an integral of a density. We close with worked numerical examples. We include a counterexample showing why right-continuity cannot be dropped. Finally, we provide a map of where CDFs resurface throughout probability and statistics.

Key takeaways
  • The CDF of a random variable XX is FX(x)=P(X≤x)F_X(x) = P(X \le x), defined for all x∈Rx \in \mathbb{R}, regardless of whether XX is discrete, continuous, or neither.
  • FXF_X is non-decreasing, right-continuous, with limits 00 at −∞-\infty and 11 at +∞+\infty — and any function with these four properties is the CDF of some random variable.
  • For a discrete variable, FXF_X is a right-continuous step function that jumps by pX(x)p_X(x) at each atom xx; the PMF is recovered from the jump sizes.
  • The CDF determines probabilities of all intervals: P(a<X≤b)=FX(b)−FX(a)P(a < X \le b) = F_X(b) - F_X(a), and with care, of open and closed intervals and single points too.
  • Two random variables have the same distribution if and only if they have the same CDF — this is the precise sense in which the CDF is a complete description.
  • The CDF is the bridge to continuous random variables: when FXF_X is absolutely continuous, its derivative is the probability density function.
  • Key idea: the CDF FX(x)=P(X≤x)F_X(x)=P(X\le x) is a single monotone, right-continuous function that fully encodes a random variable's distribution, discrete or continuous alike.

Where this lives today

The CDF is not a purely theoretical bookkeeping device — it is the working interface between probability theory and nearly every quantitative discipline that touches uncertainty.

Statistics. The empirical CDF, built from a sample by replacing P(X≤x)P(X\le x) with the sample proportion of observations ≤x\le x, is the basic nonparametric estimator of an unknown distribution, and underlies goodness-of-fit tests such as Kolmogorov–Smirnov and Cramér–von Mises. Percentiles, medians, and quantiles used throughout applied statistics are literally defined via the (generalized) inverse of the CDF.

Simulation and computer science. The inverse transform method — generate UU uniform on (0,1)(0,1) and output FX−1(U)F_X^{-1}(U) — is exactly the existence construction used in the converse direction of the characterization theorem above, and it is how most pseudo-random variate generators (for exponential, Weibull, or custom empirical distributions) are implemented in software.

Reliability engineering and survival analysis. The complement 1−FX(t)=P(X>t)1-F_X(t) = P(X>t), the survival function, describes the probability a component or patient survives past time tt; hazard rates and reliability curves are built directly from the CDF's structure.

Finance and risk management. Value-at-risk, a standard regulatory risk metric, is precisely a quantile — an inverse-CDF evaluation — of a portfolio's loss distribution.

Physics and information theory. Cumulative distributions of energy levels, particle counts, or noise amplitudes are used identically to how they are used here; entropy and information-theoretic quantities are often computed by first passing through the CDF or its density.

Ask the tutor: How does the derivative of the CDF give rise to the probability density function for continuous random variables?

DefinitionCumulative Distribution Function

Let XX be a random variable on a probability space (Ω,F,P)(\Omega, \mathcal{F}, P), meaning X:Ω→RX : \Omega \to \mathbb{R} is F\mathcal{F}-measurable. The cumulative distribution function (CDF) of XX is the function FX:R→[0,1]F_X : \mathbb{R} \to [0,1] defined by

FX(x)  =  P(X≤x)  =  P({ω∈Ω:X(ω)≤x}),x∈R.F_X(x) \;=\; P(X \le x) \;=\; P(\{\omega \in \Omega : X(\omega) \le x\}), \qquad x \in \mathbb{R}.

This is well-defined for every real xx because {X≤x}∈F\{X \le x\} \in \mathcal{F} by measurability of XX, so it has a well-defined probability.

Intuition

Think of the CDF as a running total. Imagine sweeping a vertical line from left to right across the real number line, and at each position xx recording all the probability mass that lies at or to the left of the line. As the line moves rightward, it can only accumulate more mass, never lose any — this is why FXF_X is non-decreasing. Far to the left, the line has swept up nothing, so FX(−∞)=0F_X(-\infty) = 0; far to the right, it has swept up everything, so FX(+∞)=1F_X(+\infty) = 1. For a discrete variable, mass sits at isolated points, so the running total is flat between these points and jumps abruptly at each one — a staircase. For a continuous variable, mass is spread out with no isolated atoms, so the running total climbs smoothly with no jumps at all. The CDF framework treats both pictures as instances of the same object.

Figure
The CDF as an accumulating sweep: a discrete staircase (left) vs. a continuous ramp (right)
x (discrete X) 1 0 1 2 4 jumps = pmf values x (continuous X) 1 0 smooth, no jumps

The four defining properties

TheoremCharacterization of CDFs

Let XX be a random variable with CDF FXF_X. Then FXF_X has the following four properties: 1. (Monotonicity) FXF_X is non-decreasing: if a≤ba \le b then FX(a)≤FX(b)F_X(a) \le F_X(b). 2. (Limits) lim⁡x→−∞FX(x)=0\displaystyle\lim_{x \to -\infty} F_X(x) = 0 and lim⁡x→+∞FX(x)=1\displaystyle\lim_{x \to +\infty} F_X(x) = 1. 3. (Right-continuity) For every x∈Rx \in \mathbb{R}, lim⁡h→0+FX(x+h)=FX(x)\displaystyle\lim_{h \to 0^+} F_X(x+h) = F_X(x). 4. (Left limits exist) For every x∈Rx \in \mathbb{R}, lim⁡h→0+FX(x−h)\displaystyle\lim_{h \to 0^+} F_X(x-h) exists, and equals P(X<x)P(X < x); moreover

FX(x)−lim⁡h→0+FX(x−h)=P(X=x).F_X(x) - \lim_{h\to 0^+} F_X(x-h) = P(X = x).

Conversely, if F:R→[0,1]F : \mathbb{R} \to [0,1] satisfies properties (1), (2), and (3), then there exists a probability space and a random variable XX on it with FX=FF_X = F.

Proof

We prove the forward direction (properties of any CDF) using only the axioms of probability (countable additivity, monotonicity of measures, continuity of measure along monotone sequences of events), and then sketch the converse (existence) via the standard quantile transform construction.

Property 1 (Monotonicity). Let a≤ba \le b. Then {X≤a}⊆{X≤b}\{X \le a\} \subseteq \{X \le b\} as subsets of Ω\Omega: if X(ω)≤aX(\omega) \le a and a≤ba \le b, then X(ω)≤bX(\omega) \le b. By monotonicity of probability measures (if A⊆BA \subseteq B then P(A)≤P(B)P(A) \le P(B), which follows from countable additivity applied to B=A∪(B∖A)B = A \cup (B\setminus A), a disjoint union), we get

FX(a)=P(X≤a)≤P(X≤b)=FX(b).F_X(a) = P(X\le a) \le P(X \le b) = F_X(b).

Property 2 (Limits at infinity). Let An={X≤n}A_n = \{X \le n\} for n=1,2,3,…n = 1, 2, 3, \dots ranging over the integers, or more generally consider any sequence xn→+∞x_n \to +\infty and let An={X≤xn}A_n = \{X \le x_n\} with xnx_n increasing. The events AnA_n are increasing (An⊆An+1A_n \subseteq A_{n+1}), and ⋃n=1∞An=Ω\bigcup_{n=1}^\infty A_n = \Omega because every real number X(ω)X(\omega) is eventually exceeded by xnx_n as n→∞n \to \infty. By continuity of probability measures from below (a consequence of countable additivity: if An↑AA_n \uparrow A then P(An)→P(A)P(A_n) \to P(A)),

FX(xn)=P(An)→P(Ω)=1.F_X(x_n) = P(A_n) \to P(\Omega) = 1.

Since this holds along every sequence xn→∞x_n \to \infty, and FXF_X is monotone (Property 1), the full limit lim⁡x→∞FX(x)=1\lim_{x\to\infty} F_X(x) = 1 follows: monotone functions have limits at +∞+\infty equal to sup⁡xFX(x)\sup_x F_X(x), and since some sequence achieves this sup, the sup is 11. The argument for −∞-\infty is symmetric: let xn↓−∞x_n \downarrow -\infty, so An={X≤xn}A_n = \{X \le x_n\} is decreasing with ⋂nAn=∅\bigcap_n A_n = \varnothing (no real number is ≤\le every xnx_n as xn→−∞x_n \to -\infty). By continuity of measure from above, FX(xn)=P(An)→P(∅)=0F_X(x_n) = P(A_n) \to P(\varnothing) = 0.

Property 3 (Right-continuity). Fix x∈Rx \in \mathbb{R} and let hn↓0h_n \downarrow 0 be any sequence of positive reals decreasing to 00. Let An={X≤x+hn}A_n = \{X \le x + h_n\}. Since hnh_n is decreasing, AnA_n is a decreasing sequence of events, and

⋂n=1∞An={X≤x},\bigcap_{n=1}^\infty A_n = \{X \le x\},

because ω\omega satisfies X(ω)≤x+hnX(\omega) \le x + h_n for all nn if and only if X(ω)≤xX(\omega) \le x (taking n→∞n\to\infty so hn→0h_n \to 0; if X(ω)≤x+hnX(\omega) \le x+h_n for every nn then X(ω)≤inf⁡n(x+hn)=xX(\omega) \le \inf_n (x+h_n) = x, and conversely X(ω)≤xX(\omega)\le x trivially gives X(ω)≤x+hnX(\omega) \le x+h_n for all nn). By continuity of measure from above,

FX(x+hn)=P(An)→P(X≤x)=FX(x).F_X(x+h_n) = P(A_n) \to P(X \le x) = F_X(x).

Since this holds for every sequence hn↓0h_n \downarrow 0, we conclude lim⁡h→0+FX(x+h)=FX(x)\lim_{h\to 0^+} F_X(x+h) = F_X(x).

Property 4 (Left limits and jumps). Fix xx and let hn↓0h_n \downarrow 0. Let Bn={X≤x−hn}B_n = \{X \le x - h_n\}, an increasing sequence of events with

⋃n=1∞Bn={X<x},\bigcup_{n=1}^\infty B_n = \{X < x\},

since ω\omega satisfies X(ω)≤x−hnX(\omega) \le x - h_n for some nn iff X(ω)<xX(\omega) < x (if X(ω)<xX(\omega) < x, choose nn large enough that hn<x−X(ω)h_n < x - X(\omega); conversely X(ω)≤x−hn<xX(\omega) \le x - h_n < x for any nn). By continuity of measure from below,

FX(x−hn)=P(Bn)→P(X<x).F_X(x - h_n) = P(B_n) \to P(X < x).

This limit exists (it equals P(X<x)P(X<x)) and does not depend on the sequence hnh_n, establishing the first claim. For the jump size,

FX(x)−lim⁡h→0+FX(x−h)=P(X≤x)−P(X<x)=P(X=x),F_X(x) - \lim_{h\to 0^+}F_X(x-h) = P(X\le x) - P(X < x) = P(X = x),

using finite additivity: {X≤x}={X<x}⊔{X=x}\{X \le x\} = \{X < x\} \sqcup \{X = x\} is a disjoint union, so P(X≤x)=P(X<x)+P(X=x)P(X\le x) = P(X<x) + P(X=x).

Converse (existence, quantile transform). Suppose F:R→[0,1]F:\mathbb{R}\to[0,1] satisfies (1)–(3). Define the generalized inverse (quantile function)

F−1(u)=inf⁡{x∈R:F(x)≥u},u∈(0,1).F^{-1}(u) = \inf\{x \in \mathbb{R} : F(x) \ge u\}, \qquad u \in (0,1).

Take Ω=(0,1)\Omega = (0,1), F\mathcal{F} the Borel σ\sigma-algebra on (0,1)(0,1), and PP = Lebesgue measure (i.e. let UU be a random variable uniformly distributed on (0,1)(0,1) under PP). Define X=F−1(U)X = F^{-1}(U). We claim FX=FF_X = F. The key fact is: for all u∈(0,1)u \in (0,1) and x∈Rx \in \mathbb{R},

F−1(u)≤x  ⟺  u≤F(x).F^{-1}(u) \le x \iff u \le F(x).

(⇐\Leftarrow): if u≤F(x)u \le F(x), then xx belongs to the set {y:F(y)≥u}\{y : F(y) \ge u\}, so by definition of infimum, F−1(u)≤xF^{-1}(u) \le x. (⇒\Rightarrow): if F−1(u)≤xF^{-1}(u) \le x, then since FF is right-continuous and non-decreasing, F(x)≥F(F−1(u))F(x) \ge F(F^{-1}(u)); and F(F−1(u))≥uF(F^{-1}(u)) \ge u because F−1(u)F^{-1}(u) is the infimum of the set {y:F(y)≥u}\{y: F(y)\ge u\}, and by right-continuity of FF this infimum is itself attained in the set (a decreasing sequence yn↓F−1(u)y_n \downarrow F^{-1}(u) with F(yn)≥uF(y_n)\ge u gives, by right-continuity, F(F−1(u))=lim⁡F(yn)≥uF(F^{-1}(u)) = \lim F(y_n) \ge u). Hence F(x)≥uF(x) \ge u. Using this equivalence,

FX(x)=P(X≤x)=P(F−1(U)≤x)=P(U≤F(x))=F(x),F_X(x) = P(X \le x) = P(F^{-1}(U) \le x) = P(U \le F(x)) = F(x),

where the last equality uses that UU is uniform on (0,1)(0,1), so P(U≤t)=tP(U \le t) = t for t∈[0,1]t \in [0,1], and F(x)∈[0,1]F(x) \in [0,1]. This proves FX=FF_X = F, completing the converse. ■\blacksquare

∎
CorollaryInterval probabilities from the CDF

For any random variable XX with CDF FXF_X and any reals a<ba < b,

P(a<X≤b)=FX(b)−FX(a).P(a < X \le b) = F_X(b) - F_X(a).

More generally, using Property 4 above,

P(a≤X≤b)=FX(b)−lim⁡h→0+FX(a−h),P(X=x)=FX(x)−lim⁡h→0+FX(x−h).P(a \le X \le b) = F_X(b) - \lim_{h\to0^+}F_X(a-h), \qquad P(X=x) = F_X(x)-\lim_{h\to 0^+}F_X(x-h).
Proof

Write {X≤b}={X≤a}⊔{a<X≤b}\{X \le b\} = \{X \le a\} \sqcup \{a < X \le b\} as a disjoint union (every ω\omega with X(ω)≤bX(\omega)\le b either has X(ω)≤aX(\omega)\le a or a<X(ω)≤ba < X(\omega) \le b, exclusively). By finite additivity,

FX(b)=P(X≤b)=P(X≤a)+P(a<X≤b)=FX(a)+P(a<X≤b),F_X(b) = P(X\le b) = P(X \le a) + P(a < X \le b) = F_X(a) + P(a<X\le b),

and rearranging gives the first formula. The remaining formulas follow identically from Property 4 of the theorem, using the corresponding decompositions {X≤b}={X<a}⊔{a≤X≤b}\{X\le b\} = \{X<a\}\sqcup\{a\le X\le b\} and {X≤x}={X<x}⊔{X=x}\{X\le x\}=\{X<x\}\sqcup\{X=x\}. ■\blacksquare

∎

Worked examples

Example 1 — a discrete variable with unequal, asymmetric masses. Let XX be the number of heads-runs-of-length-2-or-more observed in a specific finite experiment, and suppose (from a prior calculation using the PMF of XX) that XX takes values in {−1,2,5,7}\{-1, 2, 5, 7\} with

P(X=−1)=0.15,P(X=2)=0.35,P(X=5)=0.10,P(X=7)=0.40.P(X=-1) = 0.15,\quad P(X=2) = 0.35, \quad P(X=5) = 0.10, \quad P(X=7) = 0.40.

(Check: 0.15+0.35+0.10+0.40=1.000.15+0.35+0.10+0.40 = 1.00, so this is a valid PMF.) The CDF is built by accumulating these masses left to right:

FX(x)={0x<−10.15−1≤x<20.502≤x<50.605≤x<71.00x≥7.F_X(x) = \begin{cases} 0 & x < -1\\ 0.15 & -1 \le x < 2\\ 0.50 & 2 \le x < 5\\ 0.60 & 5 \le x < 7\\ 1.00 & x \ge 7. \end{cases}

Check right-continuity at each atom: at x=2x=2, FX(2)=0.50F_X(2) = 0.50 and lim⁡h→0+FX(2+h)=0.50\lim_{h\to0^+}F_X(2+h)=0.50 — matches, confirming the half-open intervals [2,5)[2,5) etc. are correctly placed with the jump landing at the atom. Check the left-limit rule: lim⁡h→0+FX(2−h)=0.15\lim_{h\to0^+}F_X(2-h) = 0.15, so the jump at x=2x=2 is 0.50−0.15=0.35=P(X=2)0.50 - 0.15 = 0.35 = P(X=2), exactly the PMF value, as guaranteed by Property 4. Now use the corollary to compute P(1≤X≤6)P(1 \le X \le 6) without re-deriving anything from the PMF directly:

P(1≤X≤6)=FX(6)−lim⁡h→0+FX(1−h)=0.60−0.15=0.45.P(1 \le X \le 6) = F_X(6) - \lim_{h\to0^+}F_X(1-h) = 0.60 - 0.15 = 0.45.

Sanity check against the PMF directly: the values of XX in [1,6][1,6] are 22 and 55, with total probability 0.35+0.10=0.450.35+0.10=0.45. Matches.

Example 2 — mixed (hybrid) random variable, demonstrating the CDF's universality. Consider a warranty payout YY (in units of $1000) structured as follows: with probability 0.30.3 the claim is denied and Y=0Y=0 exactly; otherwise (probability 0.70.7), YY is a continuous amount uniformly distributed on (0,4)(0, 4), but this branch's probability is itself scaled by 0.70.7. Formally,

P(Y=0)=0.3,P(Y∈(a,b])=0.7⋅b−a4 for 0≤a<b≤4.P(Y = 0) = 0.3, \qquad P(Y \in (a,b]) = 0.7 \cdot \frac{b-a}{4} \text{ for } 0 \le a < b \le 4.

This YY is neither purely discrete (it has a continuous part on (0,4)(0,4)) nor purely continuous (it has an atom at 00) — exactly the kind of object where PMFs and densities alone are inadequate, but the CDF still works without modification. For x<0x<0: FY(x)=0F_Y(x) = 0. For x=0x = 0: FY(0)=P(Y≤0)=P(Y=0)=0.3F_Y(0) = P(Y\le 0) = P(Y=0) = 0.3 (there is no probability of Y<0Y<0). For 0<x≤40 < x \le 4: FY(x)=P(Y=0)+P(0<Y≤x)=0.3+0.7⋅x4=0.3+0.175xF_Y(x) = P(Y=0) + P(0<Y\le x) = 0.3 + 0.7\cdot\frac{x}{4} = 0.3 + 0.175x. For x>4x > 4: FY(x)=0.3+0.7⋅1=1.0F_Y(x) = 0.3 + 0.7\cdot 1 = 1.0. Check continuity/limits: FY(4)=0.3+0.175(4)=0.3+0.7=1.0F_Y(4) = 0.3+0.175(4) = 0.3+0.7=1.0 ✓, matching the value for x>4x>4, so FYF_Y is continuous at x=4x=4 (no atom there). At x=0x=0, FY(0)=0.3F_Y(0)=0.3, while lim⁡h→0+FY(h)=0.3+0.175(0)=0.3\lim_{h\to0^+}F_Y(h) = 0.3 + 0.175(0)=0.3 as h→0+h\to0^+ — actually let's verify right-continuity precisely: lim⁡h→0+FY(0+h)=lim⁡h→0+(0.3+0.175h)=0.3=FY(0)\lim_{h\to 0^+} F_Y(0+h) = \lim_{h\to0^+}\left(0.3+0.175h\right) = 0.3 = F_Y(0). Right-continuous, as required. But the left limit at 00 is lim⁡h→0+FY(0−h)=lim⁡h→0+0=0≠FY(0)=0.3\lim_{h\to0^+}F_Y(0-h) = \lim_{h\to0^+}0 = 0 \ne F_Y(0) = 0.3, so there is a jump of size 0.30.3 at x=0x=0, correctly recording the atom P(Y=0)=0.3P(Y=0)=0.3. Now compute P(1.5<Y≤3.5)P(1.5 < Y \le 3.5) using only the CDF:

FY(3.5)−FY(1.5)=(0.3+0.175×3.5)−(0.3+0.175×1.5)=0.6125−0.5625=0.05.F_Y(3.5) - F_Y(1.5) = (0.3+0.175\times3.5) - (0.3+0.175\times1.5) = 0.6125 - 0.5625 = 0.05.

Direct check: this event lies entirely in the continuous branch, with length 22 out of range 44, scaled by 0.70.7: 0.7×24=0.350.7 \times \frac{2}{4} = 0.35. Wait — these disagree (0.050.05 vs 0.350.35); let's recheck the arithmetic. 0.175×3.5=0.61250.175\times3.5 = 0.6125, and 0.175×1.5=0.26250.175\times 1.5 = 0.2625, not 0.56250.5625. Recomputing: FY(1.5)=0.3+0.2625=0.5625F_Y(1.5) = 0.3+0.2625 = 0.5625 is correct as written; the error is in the direct check. Recompute directly: 3.5−1.54=24=0.5\frac{3.5-1.5}{4}=\frac{2}{4}=0.5, times 0.70.7 gives 0.350.35. But via CDF we got FY(3.5)−FY(1.5)=0.6125−0.5625=0.05F_Y(3.5)-F_Y(1.5) = 0.6125-0.5625=0.05. These still disagree, so recheck FY(3.5)F_Y(3.5): 0.175×3.5=0.61250.175\times 3.5 = 0.6125, so FY(3.5)=0.3+0.6125=0.9125F_Y(3.5)=0.3+0.6125=0.9125, not 0.61250.6125 — the earlier line dropped the +0.3+0.3. Correcting:

FY(3.5)=0.3+0.175(3.5)=0.3+0.6125=0.9125,FY(1.5)=0.3+0.175(1.5)=0.3+0.2625=0.5625,F_Y(3.5) = 0.3+0.175(3.5) = 0.3+0.6125 = 0.9125,\qquad F_Y(1.5)=0.3+0.175(1.5)=0.3+0.2625=0.5625,
FY(3.5)−FY(1.5)=0.9125−0.5625=0.35,F_Y(3.5)-F_Y(1.5) = 0.9125-0.5625 = 0.35,

which now matches the direct check of 0.350.35 exactly. This illustrates precisely why one always re-verifies CDF arithmetic against the underlying probability structure — a dropped constant is easy to miss, and the corollary formula is only as good as the correctly-evaluated FYF_Y.

Counterexamplewhy right-continuity cannot be dropped

Drop Property 3 (right-continuity) from the characterization theorem, keeping monotonicity and the correct limits at ±∞\pm\infty, and consider the function

G(x)={0x<01/20≤x<11x≥1modified toG~(x)={0x<01/20<x≤11x>1,G~(0)=0, G~(1)=1/2.G(x) = \begin{cases} 0 & x < 0\\ 1/2 & 0 \le x < 1 \\ 1 & x \ge 1 \end{cases} \quad\text{modified to}\quad \tilde G(x) = \begin{cases} 0 & x < 0\\ 1/2 & 0 < x \le 1\\ 1 & x > 1, \end{cases} \qquad \tilde G(0) = 0, \ \tilde G(1) = 1/2.

Here G~\tilde G is non-decreasing and satisfies G~(−∞)=0\tilde G(-\infty)=0, G~(+∞)=1\tilde G(+\infty)=1, but at x=0x=0: lim⁡h→0+G~(h)=1/2≠G~(0)=0\lim_{h\to0^+}\tilde G(h) = 1/2 \ne \tilde G(0) = 0. So G~\tilde G is left-continuous at 00, not right-continuous — it fails Property 3. Suppose, for contradiction, G~\tilde G were the CDF of some random variable XX, i.e. G~(x)=P(X≤x)\tilde G(x) = P(X\le x) for all xx. Then P(X≤0)=G~(0)=0P(X \le 0) = \tilde G(0) = 0, so X>0X > 0 almost surely. But also P(X≤h)=1/2P(X \le h) = 1/2 for every h∈(0,1]h \in (0,1], in particular for arbitrarily small h>0h>0, so

P(0<X≤h)=G~(h)−G~(0)=1/2−0=1/2for all h∈(0,1].P(0 < X \le h) = \tilde G(h) - \tilde G(0) = 1/2 - 0 = 1/2 \quad \text{for all } h \in (0,1].

Let hn=1/n↓0h_n = 1/n \downarrow 0. The events Cn={0<X≤hn}C_n = \{0 < X \le h_n\} are decreasing with ⋂nCn={X≤0}∩{X>0}=∅\bigcap_n C_n = \{X \le 0\} \cap \{X>0\}=\varnothing (no positive real number is ≤\le every 1/n1/n). By continuity of probability measures from above (a genuine axiom-derived property that any actual probability P(X≤x)P(X\le x), viewed as xx varies, must obey — this is exactly what Property 3 encodes), we would need

P(Cn)→P(∅)=0,P(C_n) \to P(\varnothing) = 0,

but P(Cn)=1/2P(C_n) = 1/2 for every nn, a constant sequence that cannot converge to 00. This is a contradiction. Hence no random variable can have G~\tilde G as its CDF: right-continuity is not a cosmetic convention but a mathematical consequence of the axioms of probability (specifically, continuity of measures along decreasing sequences of events with empty intersection), and dropping it from the hypothesis list breaks the existence half of the characterization theorem. This is exactly why the definition FX(x)=P(X≤x)F_X(x) = P(X\le x) uses a closed half-line (−∞,x](-\infty, x] rather than an open one (−∞,x)(-\infty, x) — the closed convention is what forces right-continuity to hold automatically.

Scope, edge cases & common traps

SituationWhat to watch for
Discrete XXFXF_X is a right-continuous step function; jump at xx equals pX(x)p_X(x); flat (constant) between atoms.
Continuous XX (foreshadowing)FXF_X is continuous everywhere (no jumps, since single points have probability 00); if additionally FXF_X is differentiable, fX=FX′f_X = F_X' is the density.
Mixed/hybrid XXFXF_X has jumps at atoms and smooth increase elsewhere; the CDF handles this with no special casing, unlike PMF/PDF, which cannot describe such XX alone.
Computing P(X=x)P(X=x) from FXF_XMust use the left-limit gap, FX(x)−FX(x−)F_X(x) - F_X(x^-), not just read off FX(x)F_X(x), unless XX is known to be discrete at that point.
Open vs. closed interval endpointsP(X<b)≠FX(b)P(X<b) \ne F_X(b) in general; P(X<b)=FX(b−)P(X<b) = F_X(b^-). Confusing these is the single most common computational error.
Converse directionAny function satisfying (1)-(3) is a CDF of some (possibly newly constructed) random variable — a genuine existence theorem, not just a checklist.

The mathematician behind it

Origins · 1 minAndrey Nikolaevich Kolmogorov and how this was discovered
AN
Andrey Nikolaevich Kolmogorov
1903–1987 · Soviet Union

Kolmogorov gave probability theory its modern axiomatic foundation in his 1933 monograph Grundbegriffe der Wahrscheinlichkeitsrechnung, building probability on measure theory. The cumulative distribution function, though its basic idea traces to earlier work of Chebyshev, Markov, and the theory of distribution of mass, became rigorous and universal only once random variables were defined as measurable functions on a probability space, exactly the framework in which the CDF's four defining properties (monotonicity, limits, right-continuity, controlled jumps) can be proved from axioms rather than merely observed empirically.

History
- 1901–1909 — Early work on distribution functions of sums of independent variables by Lyapunov and Markov, in the context of the Central Limit Theorem, uses proto-CDF ideas without a fully axiomatic footing.
- 1920s — Measure theory (Lebesgue, Borel, Carathéodory) matures, providing the tools needed to treat arbitrary distributions, not just those with densities.
- 1933 — Kolmogorov's Grundbegriffe establishes the modern axiomatic framework: probability spaces, measurable random variables, and FX(x)=P(X≤x)F_X(x)=P(X\le x) as the canonical description of a distribution.
- 1930s–1940s — The quantile transform / inverse-CDF construction becomes a standard tool for simulating and constructing random variables with prescribed distributions.
- Present — The CDF underlies modern statistical inference (empirical CDFs, Kolmogorov–Smirnov tests) and simulation (inverse transform sampling), directly descending from this foundational formalism.
Why is the CDF defined with ≤\le rather than <<?

The choice of ≤\le (a closed half-line (−∞,x](-\infty,x]) is what forces right-continuity to hold automatically as a theorem, rather than requiring it as a separate assumption. If one instead defined F~(x)=P(X<x)\tilde F(x) = P(X<x), the resulting function would be left-continuous, not right-continuous, by the mirror-image argument in the counterexample above. Both conventions carry the same information (since P(X<x)=FX(x−)P(X<x) = F_X(x^-)), but the ≤\le convention is the universal standard because it makes the characterization theorem's third property a clean consequence of continuity of measure from above, rather than an awkward continuity-from-below argument that behaves differently at atoms.

Does the CDF always determine the PMF or density?

It determines the PMF whenever XX is discrete: the PMF value at xx is exactly the jump size FX(x)−lim⁡h→0+FX(x−h)F_X(x) - \lim_{h\to0^+}F_X(x-h), as proved in Property 4. It determines the density only when FXF_X is absolutely continuous, in which case the density is (almost everywhere) the derivative FX′(x)F_X'(x); this is developed fully when Continuous Random Variables are introduced. In general (for arbitrary, possibly pathological distributions), the CDF always exists and always determines the distribution, but a PMF or density need not exist — this is precisely why the CDF, not the PMF or density, is the universally applicable object.

If two random variables have the same CDF, are they equal?

Not necessarily as functions on a sample space — they can live on entirely different probability spaces — but they are equal in distribution, meaning P(X∈B)=P(Y∈B)P(X\in B) = P(Y\in B) for every Borel set B⊆RB \subseteq \mathbb{R}. This is because the CDF determines the full induced probability measure on R\mathbb{R} (a standard result from measure theory: two finite measures on R\mathbb{R} that agree on all sets (−∞,x](-\infty,x] agree on the entire Borel σ\sigma-algebra, since these half-lines generate it and form a π\pi-system). So "same CDF" is exactly the mathematical meaning of "same distribution."

How does the CDF relate to the PMF for a discrete random variable?

For discrete XX with PMF pXp_X, the CDF is obtained by accumulating mass: FX(x)=∑k≤xpX(k)F_X(x) = \sum_{k\le x} p_X(k), summing over all atoms kk of XX that are ≤x\le x. This is always a right-continuous step function, constant between consecutive atoms and jumping exactly at each atom by the corresponding PMF value. Conversely, given the CDF of a known-discrete XX, the PMF is recovered via pX(x)=FX(x)−lim⁡h→0+FX(x−h)p_X(x) = F_X(x) - \lim_{h\to0^+}F_X(x-h), the jump-size formula from Property 4 — so the two descriptions carry identical information for discrete variables, but only the CDF generalizes beyond them.

Can a CDF have infinitely many jumps?

Yes. A discrete random variable taking countably infinitely many values (e.g. a Poisson or geometric random variable) has a CDF with a countable infinity of jumps, one at each atom, with jump sizes summing to at most 11 (they sum to exactly 11 if XX is purely discrete). Monotonicity plus boundedness in [0,1][0,1] forces the set of jump discontinuities of any CDF to be at most countable — this is a general fact about monotone bounded functions, not special to probability, but it explains why a CDF can never have uncountably many atoms of positive probability.

Is every non-decreasing, right-continuous function with limits 00 and 11 automatically a valid CDF?

Yes — this is exactly the converse half of the characterization theorem proved above, using the quantile transform construction X=F−1(U)X = F^{-1}(U) with UU uniform on (0,1)(0,1). This converse is what makes the CDF characterization so powerful in practice: to define a new probability distribution, one only needs to write down a function satisfying properties (1)–(3), with no need to separately verify measurability, countable additivity, or any other axiom — those are automatically inherited from the construction.

Why can't we just always use the PMF or density instead of the CDF?

Because PMFs require countably many atoms (undefined for genuinely continuous variables), and densities require absolute continuity of the distribution (undefined whenever there are atoms, as in the mixed Example 2 above). The CDF has no such restriction: FX(x)=P(X≤x)F_X(x)=P(X\le x) is well-defined for literally every random variable, discrete, continuous, mixed, or pathological (e.g. singular continuous distributions like the Cantor distribution, which have a continuous CDF yet no density). This universality is why the CDF, not the PMF or PDF, is the definition used to state general theorems about "the distribution of XX."

Before this

Where this leads