Visualizing expected value

A visual interpretation of the expectation of a function of a random variable (discrete or continuous) as the area over the unit interval

· Jacob Hoover Vigly ·


Let \(X\) be a random variable taking values in (some subset of) \(\mathbb{R}\), write its law as \(P_X\), and its cumulative distribution function as \(F_X(x) = P_X\big((-\infty, x]\big) = P(X \le x)\).If we want to be really formal to define all the notation, \(X : \Omega \to \mathbb{R}\) is a measurable map on a probability space \((\Omega, \mathcal{F}, P)\), and its law is the pushforward \(P_X \defeq X_* P\), so that \(P_X(A) = P(X \in A)\) for every Borel set \(A\).

The expected value of a function \(g(X)\) is given by

\[ \E[g(X)] \defeq \int_{\mathbb{R}} g \dee{P_X} \]

Writing \(p_X\) for the probability density function (or probability mass function), this general definition simplifies to a sum or integral of \(g\) weighted by \(p_X\):In what follows, we’ll consider both discrete and continuous cases, and write \(p_X\) for either pdf or pmf. This same notation reflects the fact that either can be defined formally as \(p_X = \dee{P_X} / \dee{\mu}\), the density of the random variable’s law, \(P_X\), with respect to some reference measure \(\mu\) on \(\mathbb{R}\), which is the counting measure for the discrete case, or the Lebesgue measure for the continuous case, so either the integral or the sum is the same: \(\E[g(X)] = \int_{\mathbb{R}} g \, p_X \dee{\mu}\).

\[ \E[g(X)] = \begin{cases} \displaystyle\int_{-\infty}^{\infty} g(x)\,p_X(x) \dee{x} & \text{if } X \text{ is continuous,}\\[1.2ex] \displaystyle\sum_{x \in \supp{X}} g(x)\,p_X(x) & \text{if } X \text{ is discrete.} \end{cases} \]

We can equivalently write the expectation as an integral over the unit interval, using a substitution by \(F_X^{-1}(u) = \inf\{x : F_X(x) \ge u\}\), the inverse CDF (aka quantile function).If \(X\) is discrete, \(F_X\) is a staircase that jumps by \(p_X(x)\) at each outcome \(x\). The infimum simply picks the first \(x\) at which the staircase reaches \(u\), so \(F_X^{-1}(u) = x\) for every \(u\) in the interval \(\big(P(X < x),\, P(X \le x)\big]\). These intervals partition \([0, 1]\) into blocks, one per outcome, each as wide as its probability, in order of \(x\); \(F_X^{-1}\) is the staircase that takes the value \(x\) on the block for \(x\).

\[ \E[g(X)] = \int_0^1 g\big(F_X^{-1}(u)\big) \dee{u} \]

More on that substitution...

This is the substitution \(x = F_X^{-1}(u)\): where \(X\) has a density, \(\dee{u} = p_X(x) \dee{x}\), so the weight \(p_X(x) \dee{x}\) becomes the width \(\dee{u}\). In general, for \(U \sim \mathrm{Uniform}[0, 1]\),

\[ F_X^{-1}(u) \le x \iff u \le F_X(x) \\\Longrightarrow\quad P\big(F_X^{-1}(U) \le x\big) = P\big(U \le F_X(x)\big) = F_X(x), \]

so \(F_X^{-1}(U)\) has the law of \(X\) (this is the identity behind “inverse-transform sampling”), and \(\E[g(X)] = \E\big[g\big(F_X^{-1}(U)\big)\big]\).

I wrote up a note about this kind of density transformation in a previous post (see examples 3, 4, and 5).

The fact that we’ve expressed the expectation as an integral over the unit interval of the function \(g\), suggests a nice visual interpretation of the expectation as an area. We use the map \(x \mapsto u = F_X(x)\) to rescale the real line into \([0, 1]\) so that each outcome takes up as much room as its probability. Then evaluate \(g\) over the rescaled axis, and the area under \(g \circ F_X^{-1}\) is \(\E[g(X)]\).

Below, play around with visualizing the expectation as an area under a curve. Set the density or mass function \(p_X\), and the function \(g\) to integrate, and scrub the \(x\) or \(u\) axis to see the area build up.

\(g(x) =\)

Why not just visualize the expectation as area over the whole line?

Directly from the definition of expectation, it’s also visualizable as an area over the whole support of \(X\), of the function \(g(x)\,p_X(x)\). Its part up to \(x\), shaded below, is the area up to \(u = F_X(x)\) above.

But, compared with the unit interval visualization above, I think this simple visualization isn’t as useful for some of the following reasons: