User:Boris Tsirelson/Sandbox1

Beliefs depend on the available information. This idea is formalized in probability theory by conditioning. Conditional probabilities, conditional expectations and conditional distributions are treated on three levels: discrete probabilities, probability density functions, and measure theory. Conditioning leads to a non-random result if the condition is completely specified; otherwise, if the condition is left random, the result of conditioning is also random.

This article concentrates on interrelations between various kinds of conditioning, shown mostly by examples. For systematic treatment (and corresponding literature) see more specialized articles mentioned below.

Conditioning on the discrete level

Example. A fair coin is tossed 10 times; the random variable X is the number of heads in these 10 tosses, and Y — the number of heads in the first 3 tosses. In spite of the fact that Y emerges before X it may happen that someone knows X but not Y.

Conditional probability

Given that X = 1, the conditional probability of the event Y = 0 is Template:Nowrap beginP ( Y = 0 | X = 1 ) = P ( Y = 0, X = 1 ) / P ( X = 1 ) = 0.7.Template:Nowrap end More generally,

$\mathbb {P} (Y=0|X=x)={\frac {\binom {7}{x}}{\binom {10}{x}}}={\frac {7!(10-x)!}{(7-x)!10!}}$

for x = 0, 1, 2, 3, 4, 5, 6, 7; otherwise (for x = 8, 9, 10), Template:Nowrap beginP ( Y = 0 | X = x ) = 0.Template:Nowrap end One may also treat the conditional probability as a random variable, — a function of the random variable X, namely,

$\mathbb {P} (Y=0|X)={\begin{cases}{\binom {7}{X}}/{\binom {10}{X}}&{\text{for }}X\leq 7,\\0&{\text{for }}X>7.\end{cases}}$

The expectation of this random variable is equal to the (unconditional) probability,

$\mathbb {E} (\mathbb {P} (Y=0|X))=\sum _{x}\mathbb {P} (Y=0|X=x)\mathbb {P} (X=x)=\mathbb {P} (Y=0),$

namely,

$\sum _{x=0}^{7}{\frac {\binom {7}{x}}{\binom {10}{x}}}\cdot {\frac {1}{2^{10}}}{\binom {10}{x}}={\frac {1}{8}},$

which is an instance of the law of total probability Template:Nowrap beginE ( P ( A | X ) ) = P ( A ).Template:Nowrap end

Thus, Template:Nowrap beginP ( Y = 0 | X = 1 )Template:Nowrap end may be treated as the value of the random variable Template:Nowrap beginP ( Y = 0 | X )Template:Nowrap end corresponding to X = 1. On the other hand, Template:Nowrap beginP ( Y = 0 | X = 1 )Template:Nowrap end is well-defined irrespective of other possible values of X.

Conditional expectation

Given that X = 1, the conditional expectation of the random variable Y is Template:Nowrap beginE ( Y | X = 1 ) = 0.3.Template:Nowrap end More generally,

$\mathbb {E} (Y|X=x)={\frac {3}{10}}x$

for x = 0, ..., 10. (In this example it appears to be a linear function, but in general it is nonlinear.) One may also treat the conditional expectation as a random variable, — a function of the random variable X, namely,

$\mathbb {E} (Y|X)={\frac {3}{10}}X.$

The expectation of this random variable is equal to the (unconditional) expectation of Y,

$\mathbb {E} (\mathbb {E} (Y|X))=\sum _{x}\mathbb {E} (Y|X=x)\mathbb {P} (X=x)=\mathbb {E} (Y),$

namely,

$\sum _{x=0}^{10}{\frac {3}{10}}x\cdot {\frac {1}{2^{10}}}{\binom {10}{x}}={\frac {3}{2}}\,,$   or simply   $\mathbb {E} {\Big (}{\frac {3}{10}}X{\Big )}={\frac {3}{10}}\mathbb {E} (X)={\frac {3}{10}}\cdot 5={\frac {3}{2}}\,,$

which is an instance of the law of total expectation Template:Nowrap beginE ( E ( Y | X ) ) = E ( Y ).Template:Nowrap end

The random variable Template:Nowrap beginE ( Y | X )Template:Nowrap end is the best predictor of Y given X. That is, it minimizes the mean square error Template:Nowrap beginE ( Y - f(X) )²Template:Nowrap end on the class of all random variables of the form f (X). This class of random variables remains intact if X is replaced, say, with 2X. Thus, Template:Nowrap beginE ( Y | 2X ) = E ( Y | X ).Template:Nowrap end It does not mean that Template:Nowrap beginE ( Y | 2X ) = 0.3 × 2X;Template:Nowrap end rather, Template:Nowrap beginE ( Y | 2X ) = 0.15 × 2X = 0.3 X.Template:Nowrap end In particular, Template:Nowrap beginE ( Y | 2X=2 ) = 0.3.Template:Nowrap end More generally, Template:Nowrap beginE ( Y | g(X) ) = E ( Y | X )Template:Nowrap end for every function g that is one-to-one on the set of all possible values of X. The values of X are irrelevant; what matters is the partition (denote it α_X)

$\Omega =\{X=x_{1}\}\uplus \{X=x_{2}\}\uplus \dots$

of the sample space Ω into disjoint sets $\{X=x_{n}\}.$ (Here $x_{1},x_{2},\dots$ are all possible values of X.) Given an arbitrary partition α of Ω, one may define the random variable Template:Nowrap beginE ( Y | α ).Template:Nowrap end Still, Template:Nowrap beginE ( E ( Y | α ) ) = E ( Y ).Template:Nowrap end

Conditional probability may be treated as a special case of conditional expectation. Namely, Template:Nowrap beginP ( A | X ) = E ( Y | X )Template:Nowrap end if Y is the indicator of A. Therefore the conditional probability also depends on the partition α_X generated by X rather than on X itself; Template:Nowrap beginP ( A | g(X) ) = P ( A | X ) = P ( A | α ),Template:Nowrap end Template:Nowrap beginα = α_X = α_g(X).Template:Nowrap end

On the other hand, conditioning on an event B is well-defined, provided that Template:Nowrap beginP ( B ) ≠ 0,Template:Nowrap end irrespective of any partition that may contain B as one of several parts.

Conditional distribution

Given X = x, the conditional distribution of Y is

$\mathbb {P} (Y=y|X=x)={\frac {{\binom {3}{y}}{\binom {7}{x-y}}}{\binom {10}{x}}}={\frac {{\binom {x}{y}}{\binom {10-x}{3-y}}}{\binom {10}{3}}}$

for Template:Nowrap begin0 ≤ y ≤ min ( 3, x ).Template:Nowrap end It is the hypergeometric distribution Template:Nowrap beginH ( x; 3, 7 ),Template:Nowrap end or equivalently, Template:Nowrap beginH ( 3; x, 10-x ).Template:Nowrap end The corresponding expectation 0.3 x, obtained from the general formula $n{\frac {R}{R+W}}$ for Template:Nowrap beginH ( n; R, W ),Template:Nowrap end is nothing but the conditional expectation Template:Nowrap beginE ( Y | X = x ) = 0.3 x.Template:Nowrap end

Treating Template:Nowrap beginH ( X; 3, 7 )Template:Nowrap end as a random distribution (a random vector in the four-dimensional space of all measures on {0,1,2,3}), one may take its expectation, getting the unconditional distribution of Y, — the binomial distribution Template:Nowrap beginBin ( 3, 0.5 ).Template:Nowrap end This fact amounts to the the equality

$\sum _{x=0}^{10}\mathbb {P} (Y=y|X=x)\mathbb {P} (X=x)=\mathbb {P} (Y=y)={\frac {1}{2^{3}}}{\binom {3}{y}}$

for y = 0,1,2,3; just the law of total probability.

Conditioning on the level of densities

Example. A point of the sphere x² + y² + z² = 1 is chosen at random according to the uniform distribution on the sphere . The random variables X, Y, Z are the coordinates of the random point. The joint density of X, Y, Z does not exist (since the sphere is of zero volume), but the joint density f_X,Y of X, Y exists,

$f_{X,Y}(x,y)={\begin{cases}{\frac {1}{2\pi {\sqrt {1-x^{2}-y^{2}}}}}&{\text{if }}x^{2}+y^{2}<1,\\0&{\text{otherwise}}.\end{cases}}$

(The density is non-constant because of a non-constant angle between the sphere and the plane.) The density of X may be calculated by integration,

$f_{X}(x)=\int _{-\infty }^{+\infty }f_{X,Y}(x,y)\,\mathrm {d} y=\int _{-{\sqrt {1-x^{2}}}}^{+{\sqrt {1-x^{2}}}}{\frac {\mathrm {d} y}{2\pi {\sqrt {1-x^{2}-y^{2}}}}}\,;$

surprisingly, the result does not depend on x in (-1,1),

$f_{X}(x)={\begin{cases}0.5&{\text{for }}-1<x<1,\\0&{\text{otherwise}},\end{cases}}$

which means that X is distributed uniformly on (-1,1). The same holds for Y and Z (and in fact, for Template:Nowrap beginaX + bY + cZTemplate:Nowrap end whenever Template:Nowrap begina² + b² + c² = 1).Template:Nowrap end

Conditional probability

Calculation

Given that X = 0.5, the conditional probability of the event Y ≤ 0.75 is the integral of the conditional density,

${\begin{aligned}&f_{Y|X=0.5}(y)={\frac {f_{X,Y}(0.5,y)}{f_{X}(0.5)}}={\begin{cases}{\frac {1}{\pi {\sqrt {0.75-y^{2}}}}}&{\text{for }}-{\sqrt {0.75}}<y<{\sqrt {0.75}},\\0&{\text{otherwise}}.\end{cases}}\\&\mathbb {P} (Y\leq 0.75|X=0.5)=\int _{-\infty }^{0.75}f_{Y|X=0.5}(y)\,\mathrm {d} y=\\&=\int _{-{\sqrt {0.75}}}^{0.75}{\frac {\mathrm {d} y}{\pi {\sqrt {0.75-y^{2}}}}}={\frac {1}{2}}+{\frac {1}{\pi }}\arcsin {\sqrt {0.75}}={\frac {5}{6}}\,.\end{aligned}}$

More generally,

$\mathbb {P} (Y\leq y|X=x)={\frac {1}{2}}+{\frac {1}{\pi }}\arcsin {\frac {y}{\sqrt {1-x^{2}}}}$

for all x and y such that -1 < x < 1 (otherwise the denominator f_X(x) vanishes) and $\textstyle -{\sqrt {1-x^{2}}}<y<{\sqrt {1-x^{2}}}$ (otherwise the conditional probability degenerates to 0 or 1). One may also treat the conditional probability as a random variable, — a function of the random variable X, namely,

$\mathbb {P} (Y\leq y|X)={\begin{cases}0&{\text{for }}X^{2}\geq 1-y^{2}{\text{ and }}y<0,\\{\frac {1}{2}}+{\frac {1}{\pi }}\arcsin {\frac {y}{\sqrt {1-X^{2}}}}&{\text{for }}X^{2}<1-y^{2},\\1&{\text{for }}X^{2}\geq 1-y^{2}{\text{ and }}y>0.\end{cases}}$

The expectation of this random variable is equal to the (unconditional) probability,

$\mathbb {E} (\mathbb {P} (Y\leq y|X))=\int _{-\infty }^{+\infty }\mathbb {P} (Y\leq y|X=x)f_{X}(x)\,\mathrm {d} x=\mathbb {P} (Y\leq y),$

which is an instance of the law of total probability Template:Nowrap beginE ( P ( A | X ) ) = P ( A ).Template:Nowrap end

Interpretation

The conditional probability Template:Nowrap beginP ( Y ≤ 0.75 | X = 0.5 )Template:Nowrap end cannot be interpreted as Template:Nowrap beginP ( Y ≤ 0.75, X = 0.5 ) / P ( X = 0.5 ),Template:Nowrap end since the latter gives 0/0. Accordingly, Template:Nowrap beginP ( Y ≤ 0.75 | X = 0.5 )Template:Nowrap end cannot be interpreted via empirical frequencies, since the exact value X = 0.5 has no chance to appear at random, not even once during an infinite sequence of independent trials.

The conditional probability can be interpreted as a limit,

${\begin{aligned}&\mathbb {P} (Y\leq 0.75|X=0.5)=\lim _{\varepsilon \to 0+}\mathbb {P} (Y\leq 0.75|0.5-\varepsilon <X<0.5+\varepsilon )=\\&=\lim _{\varepsilon \to 0+}{\frac {\mathbb {P} (Y\leq 0.75,0.5-\varepsilon <X<0.5+\varepsilon )}{\mathbb {P} (0.5-\varepsilon <X<0.5+\varepsilon )}}=\\&=\lim _{\varepsilon \to 0+}{\frac {\int _{0.5-\varepsilon }^{0.5+\varepsilon }\mathrm {d} x\int _{-\infty }^{0.75}\mathrm {d} y\,f_{X,Y}(x,y)}{\int _{0.5-\varepsilon }^{0.5+\varepsilon }\mathrm {d} x\,f_{X}(x)}}\,.\end{aligned}}$

Conditional expectation

The conditional expectation Template:Nowrap beginE ( Y | X = 0.5 )Template:Nowrap end is of little interest; it vanishes just by symmetry. It is more interesting to calculate Template:Nowrap beginE ( |Z| | X = 0.5 )Template:Nowrap end treating |Z| as a function of X, Y:

${\begin{aligned}&|Z|=h(X,Y)={\sqrt {1-X^{2}-Y^{2}}}\,;\\&\mathrm {E} (|Z||X=0.5)=\int _{-\infty }^{+\infty }h(0.5,y)f_{Y|X=0.5}(y)\,\mathrm {d} y=\\&=\int _{-{\sqrt {0.75}}}^{+{\sqrt {0.75}}}{\sqrt {0.75-y^{2}}}\cdot {\frac {\mathrm {d} y}{\pi {\sqrt {0.75-y^{2}}}}}={\frac {2}{\pi }}{\sqrt {0.75}}\,.\end{aligned}}$

More generally,

$\mathbb {E} (|Z||X=x)={\frac {2}{\pi }}{\sqrt {1-x^{2}}}$

for -1 < x < 1. One may also treat the conditional expectation as a random variable, — a function of the random variable X, namely,

$\mathbb {E} (|Z||X)={\frac {2}{\pi }}{\sqrt {1-X^{2}}}\,.$

The expectation of this random variable is equal to the (unconditional) expectation of |Z|,

$\mathbb {E} (\mathbb {E} (|Z||X))=\int _{-\infty }^{+\infty }\mathbb {E} (|Z||X=x)f_{X}(x)\,\mathrm {d} x=\mathbb {E} (|Z|)\,,$

namely,

$\int _{-1}^{+1}{\frac {2}{\pi }}{\sqrt {1-x^{2}}}\cdot {\frac {\mathrm {d} x}{2}}={\frac {1}{2}}\,,$

which is an instance of the law of total expectation Template:Nowrap beginE ( E ( Y | X ) ) = E ( Y ).Template:Nowrap end

The random variable Template:Nowrap beginE ( |Z| | X )Template:Nowrap end is the best predictor of |Z| given X. That is, it minimizes the mean square error Template:Nowrap beginE ( |Z| - f(X) )²Template:Nowrap end on the class of all random variables of the form f (X). Similarly to the discrete case, Template:Nowrap beginE ( |Z| | g(X) ) = E ( |Z| | X )Template:Nowrap end for every measurable function g that is one-to-one on (-1,1).

Conditional distribution

Given X = x, the conditional distribution of Y, given by the density f_Y|X=x(y), is the (rescaled) arcsin distribution; its cumulative distribution function is

$F_{Y|X=x}(y)=\mathbb {P} (Y\leq y|X=x)={\frac {1}{2}}+{\frac {1}{\pi }}\arcsin {\frac {y}{\sqrt {1-x^{2}}}}$

for all x and y such that x² + y² < 1.The corresponding expectation of h(x,Y) is nothing but the conditional expectation Template:Nowrap beginE ( h(X,Y) | X=x ).Template:Nowrap end The mixture of these conditional distributions, taken for all x (according to the distribution of X) is the unconditional distribution of Y. This fact amounts to the equalities

${\begin{aligned}&\int _{-\infty }^{+\infty }f_{Y|X=x}(y)f_{X}(x)\,\mathrm {d} x=f_{Y}(y)\,,\\&\int _{-\infty }^{+\infty }F_{Y|X=x}(y)f_{X}(x)\,\mathrm {d} x=F_{Y}(y)\,,\end{aligned}}$

the latter being the instance of the law of total probability mentioned above.

What conditioning is not

On the discrete level conditioning is possible only if the condition is of nonzero probability (one cannot divide by zero). On the level of densities, conditioning on X = x is possible even though Template:Nowrap beginP ( X = x ) = 0.Template:Nowrap end This success may create the illusion that conditioning is always possible. Regretfully, it is not, for several reasons presented below.

Geometric intuition: caution

The result Template:Nowrap beginP ( Y ≤ 0.75 | X = 0.5 ) = 5/6,Template:Nowrap end mentioned above, is geometrically evident in the following sense. The points (x,y,z) of the sphere x² + y² + z² = 1, satisfying the condition x = 0.5, are a circle y² + z² = 0.75 of radius ${\sqrt {0.75}}$ on the plane x = 0.5. The inequality y ≤ 0.75 holds on an arc. The length of the arc is 5/6 of the length of the circle, which is why the conditional probability is equal to 5/6.

This successful geometric explanation may create the illusion that the following question is trivial.

A point of a given sphere is chosen at random (uniformly). Given that the point lies on a given plane, what is its conditional distribution?

It may seem evident that the conditional distribution must be uniform on the given circle (the intersection of the given sphere and the given plane). Sometimes it really is, but in general it is not. Especially, Z is distributed uniformly on (-1,+1) and independent of the ratio Y/X, thus, Template:Nowrap beginP ( Z ≤ 0.5 | Y/X ) = 0.75.Template:Nowrap end On the other hand, the inequality z ≤ 0.5 holds on an arc of the circle Template:Nowrap beginx² + y² + z² = 1,Template:Nowrap end Template:Nowrap beginy = cxTemplate:Nowrap end (for any given c). The length of the arc is 2/3 of the length of the circle. However, the conditional probability is 3/4, not 2/3. This is a manifestation of the classical Borel paradox .

‘	Appeals to symmetry can be misleading if not formalized as invariance arguments.	’
—Pollard

Another example. A random rotation of the three-dimensional space is a rotation by a random angle around a random axis. Geometric intuition suggests that the angle is independent of the axis and distributed uniformly. However, the latter is wrong; small values of the angle are less probable.

The limiting procedure

Given an event B of zero probability, the formula $\textstyle \mathbb {P} (A|B)=\mathbb {P} (A\cap B)/\mathbb {P} (B)$ is useless, however, one can try $\textstyle \mathbb {P} (A|B)=\lim _{n\to \infty }\mathbb {P} (A\cap B_{n})/\mathbb {P} (B_{n})$ for an appropriate sequence of events B_n of nonzero probability such that B_n ↓ B (that is, $\textstyle B_{1}\supset B_{2}\supset \dots$ and $\textstyle B_{1}\cap B_{2}\cap \dots =B$). One example is given above. Two more examples are Brownian bridge and Brownian excursion.

In the latter two examples the law of total probability is irrelevant, since only a single event (the condition) is given. In contrast, in the example above the law of total probability applies, since the event X = 0.5 is included into a family of events X = x where x runs over (-1,1), and these events are a partition of the probability space.

In order to avoid paradoxes (such as the Borel's paradox), the following important distinction should be taken into account. If a given event is of nonzero probability then conditioning on it is well-defined (irrespective of any other events), as was noted above. In contrast, if the given event is of zero probability then conditioning on it is ill-defined unless some additional input is provided. Wrong choice of this additional input leads to wrong conditional probabilities (expectations, distributions). In this sense, "the concept of a conditional probability with regard to an isolated hypothesis whose probability equals 0 is inadmissible." (Kolmogorov; quoted in ).

The additional input may be (a) a symmetry (invariance group); (b) a sequence of events B_n such that B_n ↓ B, P ( B_n ) > 0; (c) a partition containing the given event. Measure-theoretic conditioning (below) investigates Case (c), discloses its relation to (b) in general and to (a) when applicable.

Some events of zero probability are beyond the reach of conditioning. An example: let X_n be independent random variables distributed uniformly on (0,1), and B the event Template:Nowrap begin"X_n → 0Template:Nowrap end as Template:Nowrap begin$\textstyle n\to \infty$";Template:Nowrap end what about Template:Nowrap beginP ( X_n < 0.5 | B ) ?Template:Nowrap end Does it tend to 1, or not? Another example: let X be a random variable distributed uniformly on (0,1), and B the event "X is a rational number"; what about Template:Nowrap beginP ( X = 1/n | B ) ?Template:Nowrap end The only answer is that, once again,

‘	the concept of a conditional probability with regard to an isolated hypothesis whose probability equals 0 is inadmissible.	’
—Kolmogorov, quoted in

Conditioning on the level of measure theory

Example. Let Y be a random variable distributed uniformly on (0,1), and X = f(Y) where f is a given function. Two cases are treated below: f = f₁ and f = f₂, where f₁ is the continuous piecewise-linear function

$f_{1}(y)={\begin{cases}3y&{\text{for }}0\leq y\leq 1/3,\\1.5(1-y)&{\text{for }}1/3\leq y\leq 2/3,\\0.5&{\text{for }}2/3\leq y\leq 1,\end{cases}}$

and f₂ is the Weierstrass function.

Geometric intuition: caution

Given X = 0.75, two values of Y are possible, 0.25 and 0.5. It may seem evident that both values are of conditional probability 0.5 just because one point is congruent to another point. However, this is an illusion; see below.

Conditional probability

The conditional probability Template:Nowrap beginP ( Y ≤ 1/3 | X )Template:Nowrap end may be defined as the best predictor of the indicator

$I={\begin{cases}1&{\text{if }}Y\leq 1/3,\\0&{\text{otherwise}},\end{cases}}$

given X. That is, it minimizes the mean square error Template:Nowrap beginE ( I - g(X) )²Template:Nowrap end on the class of all random variables of the form g (X).

In the case f = f₁ the corresponding function g = g₁ may be calculated explicitly,

$g_{1}(x)={\begin{cases}1&{\text{for }}0<x<0.5,\\0&{\text{for }}x=0.5,\\1/3&{\text{for }}0.5<x<1.\end{cases}}$

Alternatively, the limiting procedure may be used,

$g_{1}(x)=\lim _{\varepsilon \to 0+}\mathbb {P} (Y\leq 1/3|x-\varepsilon \leq X\leq x+\varepsilon )\,,$

giving the same result.

Thus, Template:Nowrap beginP ( Y ≤ 1/3 | X ) = g₁ (X).Template:Nowrap end The expectation of this random variable is equal to the (unconditional) probability, Template:Nowrap beginE ( P ( Y ≤ 1/3 | X ) ) = P ( Y ≤ 1/3 ),Template:Nowrap end namely,

$1\cdot \mathbb {P} (X<0.5)+0\cdot \mathbb {P} (X=0.5)+{\frac {1}{3}}\cdot \mathbb {P} (X>0.5)=1\cdot {\frac {1}{6}}+0\cdot {\frac {1}{3}}+{\frac {1}{3}}\cdot {\Big (}{\frac {1}{6}}+{\frac {1}{3}}{\Big )}={\frac {1}{3}}\,,$

which is an instance of the law of total probability Template:Nowrap beginE ( P ( A | X ) ) = P ( A ).Template:Nowrap end

In the case f = f₂ the corresponding function g = g₂ probably cannot be calculated explicitly. Nevertheless it exists, and can be computed numerically. Indeed, the space L₂ (Ω) of all square integrable random variables is a Hilbert space; the indicator I is a vector of this space; and random variables of the form g (X) are a (closed, linear) subspace. The orthogonal projection of this vector to this subspace is well-defined. It can be computed numerically, using finite-dimensional approximations to the infinite-dimensional Hilbert space.

Once again, the expectation of the random variable Template:Nowrap beginP ( Y ≤ 1/3 | X ) = g₂ (X)Template:Nowrap end is equal to the (unconditional) probability, Template:Nowrap beginE ( P ( Y ≤ 1/3 | X ) ) = P ( Y ≤ 1/3 ),Template:Nowrap end namely,

$\int _{0}^{1}g_{2}(f_{2}(y))\,\mathrm {d} y={\frac {1}{3}}\,.$

However, the Hilbert space approach treats g₂ as an equivalence class of functions rather than an individual function. Measurability of g₂ is ensured, but continuity (or even Riemann integrability) is not. The value g₂ (0.5) is determined uniquely, since the point 0.5 is an atom of the distribution of X. Other values x are not atoms, thus, corresponding values g₂ (x) are not determined uniquely. Once again, "the concept of a conditional probability with regard to an isolated hypothesis whose probability equals 0 is inadmissible." (Kolmogorov; quoted in ).

Alternatively, the same function g (be it g₁ or g₂) may be defined as the Radon-Nikodym derivative

$g={\frac {\mathrm {d} \nu }{\mathrm {d} \mu }}\,,$

where measures μ, ν are defined by

${\begin{aligned}\mu (B)&=\mathbb {P} (X\in B)\,,\\\nu (B)&=\mathbb {P} (X\in B,\,Y\leq 1/3)\end{aligned}}$

for all Borel sets $B\subset \mathbb {R} .$ That is, μ is the (unconditional) distribution of X, while ν is one third of its conditional distribution,

$\nu (B)=\mathbb {P} (X\in B|Y\leq 1/3)\mathbb {P} (Y\leq 1/3)={\frac {1}{3}}\mathbb {P} (X\in B|Y\leq 1/3)\,.$

Both approaches (via the Hilbert space, and via the Radon-Nikodym derivative) treat g as an equivalence class of functions; two functions g and g′ are treated as equivalent, if g (X) = g′ (X) almost surely. Accordingly, the conditional probability Template:Nowrap beginP ( Y ≤ 1/3 | X )Template:Nowrap end is treated as an equivalence class of random variables; as usual, two random variables are treated as equivalent if they are equal almost surely.

Conditional expectation

The conditional expectation Template:Nowrap beginE ( Y | X )Template:Nowrap end may be defined as the best predictor of Y given X. That is, it minimizes the mean square error Template:Nowrap beginE ( Y - h(X) )²Template:Nowrap end on the class of all random variables of the form h(X).

In the case f = f₁ the corresponding function h = h₁ may be calculated explicitly,

$h_{1}(x)={\begin{cases}x/3&{\text{for }}0<x<0.5,\\5/6&{\text{for }}x=0.5,\\(2-x)/3&{\text{for }}0.5<x<1,\end{cases}}$

Alternatively, the limiting procedure may be used,

$h_{1}(x)=\lim _{\varepsilon \to 0+}\mathbb {E} (Y|x-\varepsilon \leq X\leq x+\varepsilon )\,,$

giving the same result.

Thus, Template:Nowrap beginE ( Y | X ) = h₁ (X).Template:Nowrap end The expectation of this random variable is equal to the (unconditional) expectation, Template:Nowrap beginE ( E ( Y | X ) ) = E ( Y ),Template:Nowrap end namely,

${\begin{aligned}&\int _{0}^{1}h_{1}(f_{1}(y))\,\mathrm {d} y=\int _{0}^{1/6}{\frac {3y}{3}}\,\mathrm {d} y+\\&\quad +\int _{1/6}^{1/3}{\frac {2-3y}{3}}\,\mathrm {d} y+\int _{1/3}^{2/3}{\frac {2-1.5(1-y)}{3}}\,\mathrm {d} y+\int _{2/3}^{1}{\frac {5}{6}}\,\mathrm {d} y={\frac {1}{2}}\,,\end{aligned}}$

which is an instance of the law of total expectation Template:Nowrap beginE ( E ( Y | X ) ) = E ( Y ).Template:Nowrap end

In the case f = f₂ the corresponding function h = h₂ probably cannot be calculated explicitly. Nevertheless it exists, and can be computed numerically in the same way as g₂ above, — as the orthogonal projection in the Hilbert space. The law of total expectation holds, since the projection cannot change the scalar product by the constant 1 belonging to the subspace.

Alternatively, the same function h (be it h₁ or h₂) may be defined as the Radon-Nikodym derivative

$h={\frac {\mathrm {d} \nu }{\mathrm {d} \mu }}\,,$

where measures μ, ν are defined by

${\begin{aligned}\mu (B)&=\mathbb {P} (X\in B)\,,\\\nu (B)&=\mathbb {E} (Y,\,X\in B)\end{aligned}}$

for all Borel sets $B\subset \mathbb {R} .$ Here Template:Nowrap beginE ( Y; A )Template:Nowrap end is the restricted expectation, not to be confused with the conditional expectation Template:Nowrap beginE ( Y | A ) = E (Y; A ) / P ( A ).Template:Nowrap end

Conditional distribution

In the case f = f₁ the conditional cumulative distribution function may be calculated explicitly, similarly to g₁. The limiting procedure gives

${\begin{aligned}&F_{Y|X=0.75}(y)=\mathbb {P} (Y\leq y|X=0.75)=\\&=\lim _{\varepsilon \to 0+}\mathbb {P} (Y\leq y|0.75-\varepsilon \leq X\leq 0.75+\varepsilon )=\\&={\begin{cases}0&{\text{for }}-\infty <y<1/4,\\1/6&{\text{for }}y=1/4,\\1/3&{\text{for }}1/4<y<1/2,\\2/3&{\text{for }}y=1/2,\\1&{\text{for }}1/2<y<\infty ,\end{cases}}\end{aligned}}$

which cannot be correct, since a cumulative distribution function must be right-continuous!

This paradoxical result is explained by measure theory as follows. For a given y the corresponding Template:Nowrap beginF_Y|X=x(y) = P ( Y ≤ y | X = x )Template:Nowrap end is well-defined (via the Hilbert space or the Radon-Nikodym derivative) as an equivalence class of functions (of x). Treated as a function of y for a given x it is ill-defined unless some additional input is provided. Namely, a function (of x) must be chosen within every (or at least almost every) equivalence class. Wrong choice leads to wrong conditional cumulative distribution functions.

A right choice can be made as follows. First, Template:Nowrap beginF_Y|X=x(y) = P ( Y ≤ y | X = x )Template:Nowrap end is considered for rational numbers y only. (Any other dense countable set may be used equally well.) Thus, only a countable set of equivalence classes is used; all choices of functions within these classes are mutually equivalent, and the corresponding function of rational y is well-defined (for almost every x). Second, the function is extended from rational numbers to real numbers by right continuity.

In general the conditional distribution is defined for almost all x (according to the distribution of X), but sometimes the result is continuous in x, in which case individual values are acceptable. In the considered example this is the case; the correct result for x = 0.75,

${\begin{aligned}&F_{Y|X=0.75}(y)=\mathbb {P} (Y\leq y|X=0.75)=\\&={\begin{cases}0&{\text{for }}-\infty <y<1/4,\\1/3&{\text{for }}1/4\leq y<1/2,\\1&{\text{for }}1/2\leq y<\infty \end{cases}}\end{aligned}}$

shows that the conditional distribution of Y given X = 0.75 consists of two atoms, at 0.25 and 0.5, of probabilities 1/3 and 2/3 respectively.

Similarly, the conditional distribution may be calculated for all x in (0, 0.5) or (0.5, 1).

The value x = 0.5 is an atom of the distribution of X, thus, the corresponding conditional distribution is well-defined and may be calculated by elementary means (the denominator does not vanish); the conditional distribution of Y given X = 0.5 is uniform on (2/3, 1). Measure theory leads to the same result.

The mixture of all conditional distributions is the (unconditional) distribution of Y.

The conditional expectation Template:Nowrap beginE ( Y | X = x )Template:Nowrap end is nothing but the expectation with respect to the conditional distribution.

In the case f = f₂ the corresponding Template:Nowrap beginF_Y|X=x(y) = P ( Y ≤ y | X = x )Template:Nowrap end probably cannot be calculated explicitly. For a given y it is well-defined (via the Hilbert space or the Radon-Nikodym derivative) as an equivalence class of functions (of x). The right choice of functions within these equivalence classes may be made as above; it leads to correct conditional cumulative distribution functions, thus, conditional distributions. In general, conditional distributions need not be atomic or absolutely continuous (nor mixtures of both types). Probably, in the considered example they are singular (like the Cantor distribution).

Once again, the mixture of all conditional distributions is the (unconditional) distribution, and the conditional expectation is the expectation with respect to the conditional distribution.

See also

• Conditional probability
• Conditional expectation
• Conditional probability distribution
• Borel's paradox
• Regular conditional probability
• Disintegration theorem

References

• Durrett, Richard (1996), Probability: theory and examples (Second ed.)
• Pollard, David (2002), A user's guide to measure theoretic probability, Cambridge University Press