Let (Ω,F,λ)(\Omega, \mathcal{F}, \lambda) be a probability space, and X,Y:Ω→RX, Y : \Omega \rightarrow \mathbb{R} be independent random variables.

Define the distributions of XX and YY by

μ(A)≐Pr⁡[X∈A]≐λ(X−1(A)),∀A∈B,ν(B)≐Pr⁡[Y∈B]≐λ(Y−1(B)),∀B∈B. \begin{aligned} \mu(A) \doteq \Pr[ X \in A ] \doteq \lambda \Big(X^{-1}(A) \Big), \quad \forall A \in \mathcal{B}, \\ \nu(B) \doteq \Pr[ Y \in B ] \doteq \lambda \Big(Y^{-1}(B) \Big), \quad \forall B \in \mathcal{B}. \end{aligned}

Hence, both μ\mu and ν\nu are probability measures on (R,B)(\mathbb{R}, \mathcal{B}). Denote μ×ν\mu \times \nu the product measure of μ\mu and ν\nu on (R2,B2)(\mathbb{R}^2, \mathcal{B}^2).

By independence of XX and YY, ∀A,B∈B\forall A, B \in \mathcal{B},

Pr⁡[X∈A,Y∈B]=Pr⁡[X∈A]⋅Pr⁡[Y∈B]=μ(A)⋅ν(B)=(μ×ν)(A×B). \begin{aligned} \Pr[ X \in A, Y \in B ] &= \Pr[ X \in A ] \cdot \Pr[ Y \in B ] \\ &= \mu(A) \cdot \nu(B) \\ &= (\mu \times \nu) (A \times B). \end{aligned}

Since S≐{A×B:A,B∈B}\mathcal{S} \doteq \{ A \times B: A, B \in \mathcal{B} \} generates B2\mathcal{B}^2, i.e., σ(S)=B2\sigma(\mathcal{S}) = \mathcal{B}^2, applying π\pi-λ\lambda theorem, we get

Pr⁡[(X,Y)∈M]=(μ×ν)(M),∀M∈B2. \Pr[ (X, Y) \in M ] = (\mu \times \nu) (M), \quad \forall M \in \mathcal{B}^2.

It follows that μ×ν\mu \times \nu is the distribution of (X,Y)(X, Y).

We list some interesting corollaries based on Fubini's theorem.

Corollary One. For each M∈B2M \in \mathcal{B}^2, it holds that

Pr⁡[(X,Y)∈M]=∫R2Pr⁡[(x,Y)∈M] μ(dx). \Pr[ (X, Y) \in M ] = \int_{\mathbb{R}^2} \Pr[ (x, Y) \in M ] \, \mu( dx ).

Similarly, for each A∈B,M∈B2A \in \mathcal{B}, M \in \mathcal{B}^2,

Pr⁡[X∈A,(X,Y)∈M]=∫APr⁡[(x,Y)∈M] μ(dx). \Pr[ X \in A, (X, Y) \in M ] = \int_{A} \Pr[ (x, Y) \in M ] \, \mu( dx ).

Proof. Via Fubini's theorem,

Pr⁡[(X,Y)∈M]=∫R2I[(x,y)∈M]  (μ×ν)(d(x,y)).=∫R∫RI[(x,y)∈M]  ν(dy)  μ(dx)=∫Rν({y:(x,y)∈M})  μ(dx)=∫RPr⁡[(x,Y)∈M]  μ(dx), \begin{aligned} \Pr[ (X, Y) \in M ] &= \int_{\mathbb{R}^2} \mathbb{I}[ (x, y) \in M ] \, \, (\mu \times \nu)( d (x, y) ). \\ &= \int_{\mathbb{R}} \int_{\mathbb{R}} \mathbb{I}[ (x, y) \in M ] \, \, \nu( dy ) \, \, \mu( dx ) \\ &= \int_{\mathbb{R}} \nu( \{ y : (x, y) \in M \} ) \, \, \mu( dx ) \\ &= \int_{\mathbb{R}} \Pr[ (x, Y) \in M ] \, \, \mu( dx ), \end{aligned}

where Pr⁡[(x,Y)∈M]\Pr[ (x, Y) \in M ] is a shorthand for

Y∈Mx≐{y:(x,y)∈M}. Y \in M_x \doteq \{ y : (x, y) \in M \}.

Substituting the function I[(x,y)∈M]\mathbb{I}[ (x, y) \in M ] with

I[x∈A]⋅I[(x,y)∈M]. \mathbb{I}[ x \in A ] \cdot \mathbb{I}[ (x, y) \in M ].

proves the second claim.

□\square

Another corollary is a direct application of the Fubini's theorem.

Corollary Two. Assume that f:R2→Rf: \mathbb{R}^2 \rightarrow \mathbb{R} is a function measurable B2/B\mathcal{B}^2 / \mathcal{B}. Define

g(x)≐E[f(x,Y)]=∫Rf(x,y)  ν(dy),∀x∈R. g(x) \doteq \mathbb{E}[ f(x, Y) ] = \int_{\mathbb{R}} f(x, y) \,\, \nu (dy), \quad \forall x \in \mathbb{R}.

Then we have:

E[g(X)]=E[f(X,Y)], \mathbb{E}[ g(X) ] = \mathbb{E} [ f(X, Y) ],

i.e.,

∫R(∫Rf(x,y)  ν(dy))  μ(dx)=∫R2f(x,y)  (μ×ν)(d(x,y)) \int_{\mathbb{R}} \left( \int_{\mathbb{R}} f(x, y) \,\, \nu (dy) \right) \, \, \mu (dx) = \int_{\mathbb{R}^2} f(x, y) \, \, (\mu \times \nu)(d (x,y))

Remarks:

  1. If XX and YY are not independent, then the above claim may not hold, since the product measure μ×ν\mu \times \nu may not be the joint distribution of (X,Y)(X, Y).

  2. In the above, we operate on the product space (R×R,B×B,μ×ν)=(R2,B2,μ×ν)(\mathbb{R} \times \mathbb{R}, \mathcal{B} \times \mathcal{B}, \mu \times \nu) = (\mathbb{R}^2, \mathcal{B}^2, \mu \times \nu), instead of (Ω×Ω,F×F,λ×λ)=(Ω2,F2,λ2)(\Omega \times \Omega, \mathcal{F} \times \mathcal{F}, \lambda \times \lambda) = (\Omega^2, \mathcal{F}^2, \lambda^2).

  3. There are various interpretations of E[f(X,Y)]\mathbb{E} [ f(X, Y) ]:

  • f(X,Y)f(X, Y) can be viewed as a function from Ω\Omega to R\mathbb{R}, given by ω→f(X(ω),Y(ω))\omega \rightarrow f( X(\omega), Y(\omega) ). It is easy to see that it is measurable F/B\mathcal{F} / \mathcal{B}, and therefore is a random variable. Then

    E[f(X,Y)]=∫Ωf(X(ω),Y(ω))  λ(dω). \mathbb{E} [ f(X, Y) ] = \int_\Omega f( X(\omega), Y(\omega) ) \, \, \lambda (d \omega).
  • Alternative, we can view f(x,y)f(x, y) as a measurable function from the probability space (R2,B2,μ×ν)(\mathbb{R}^2, \mathcal{B}^2, \mu \times \nu), to R\mathbb{R}. It is a random variable, and

    E[f(X,Y)]=∫R2f(x,y)  (μ×ν)(d(x,y)). \mathbb{E} [ f(X, Y) ] = \int_{\mathbb{R}^2} f(x, y) \, \, (\mu \times \nu)(d (x,y)).
  • It can be proven that the two integrals are equal:

    • If f(x,y)=I[(x,y)∈A×B]f(x, y) = \mathbb{I}[ (x, y) \in A \times B], for some A,B∈BA, B \in \mathcal{B}, then by independence of XX and YY, both integrals equal μ(A)⋅ν(B)\mu(A) \cdot \nu(B).

    • Since S≐{A×B:A,B∈B}\mathcal{S} \doteq \{ A \times B: A, B \in \mathcal{B} \} generates B2\mathcal{B}^2, i.e., σ(S)=B2\sigma(\mathcal{S}) = \mathcal{B}^2, applying π\pi-λ\lambda theorem, if f(x,y)=I[(x,y)∈M]f(x, y) = \mathbb{I}[ (x, y) \in M] for some M∈B2M \in \mathcal{B}^2, the two integral are equal.

    • By linearity of integration, if f(x,y)f(x, y) is B2\mathcal{B}^2 simple function, the two integral are equal.

    • By monotone convergence theorem, if f(x,y)f(x, y) is a non-negative measurable function, the two integral are equal.

    • By decomposing ff into non-negative and negative parts, if f(x,y)f(x, y) is an integrable function, the two integral are equal.

Reference

[1] P. Billingsley, Probability and measure, 3rd ed. New York: Wiley, 1995.