Theory of Statistics 1
MAR 2026

Bivariate & Trivariate Discrete Transformations

Direct Summation, Auxiliary Variables, and Convolutions

Bivariate Transformations

We have 2 discrete random variables (X1,X2)(X_1, X_2) and we know its distribution fX1,X2(x1,x2)f_{X_1, X_2}(x_1, x_2). We are interested in a random variable YY which is some function of X1,X2X_1, X_2: Y=g(X1,X2)Y = g(X_1, X_2)

We want to determine the probability distribution of YY.

Example 1

Consider the following bivariate probability distribution of X1X_1 and X2X_2:

X1\X2X_1 \backslash X_22-21122
1-11/61/61/61/61/121/12
001/121/121/121/1200
111/61/61/61/61/121/12

Find the probability distribution of Y=X1+X2Y = X_1 + X_2.

Bivariate Transformation: Y=X1+X2Y = X_1 + X_2

Hover over the joint table or the marginal table to trace the mappings.

Joint PMF fX1,X2(x1,x2)f_{X_1, X_2}(x_1, x_2)
x1\x2x_1 \backslash x_22-21122
1-1
0.167
Y=-3
0.167
Y=0
0.083
Y=1
00
0.083
Y=-2
0.083
Y=1
0.000
Y=2
11
0.167
Y=-1
0.167
Y=2
0.083
Y=3
Marginal PMF fY(y)f_Y(y)
yyfY(y)f_Y(y)
-3
0.167
-2
0.083
-1
0.167
0
0.167
1
0.167
2
0.167
3
0.083
Select any cell or marginal value to see the underlying summation theorem applied.
Method 1

The Direct Method (Not 1-1 Transformation)

Direct Summation

We map each (x1,x2)(x_1, x_2) pair to yy:

(x1,x2)(x_1, x_2)yyfX1,X2(x1,x2)f_{X_1, X_2}(x_1, x_2)
(1,2)(-1, -2)3-31/61/6
(1,1)(-1, 1)001/61/6
(1,2)(-1, 2)111/121/12
(0,2)(0, -2)2-21/121/12
(0,1)(0, 1)111/121/12
(0,2)(0, 2)2200
(1,2)(1, -2)1-11/61/6
(1,1)(1, 1)221/61/6
(1,2)(1, 2)331/121/12

By summing the probabilities for each unique yy, we get the distribution of YY.

Method 2

The Auxiliary Variable Method (1-1 Transformation)

Auxiliary Transformation

We introduce an auxiliary variable Y2=X2Y_2 = X_2, making the transformation (Y1=X1+X2,Y2=X2)(Y_1 = X_1 + X_2, Y_2 = X_2) one-to-one.

(x1,x2)(x_1, x_2)(y1,y2)(y_1, y_2)fY1,Y2(y1,y2)=fX1,X2(x1,x2)f_{Y_1, Y_2}(y_1, y_2) = f_{X_1, X_2}(x_1, x_2)
(1,2)(-1, -2)(3,2)(-3, -2)1/61/6
(1,1)(-1, 1)(0,1)(0, 1)1/61/6
(1,2)(-1, 2)(1,2)(1, 2)1/121/12
(0,2)(0, -2)(2,2)(-2, -2)1/121/12
(0,1)(0, 1)(1,1)(1, 1)1/121/12
(0,2)(0, 2)(2,2)(2, 2)00
(1,2)(1, -2)(1,2)(-1, -2)1/61/6
(1,1)(1, 1)(2,1)(2, 1)1/61/6
(1,2)(1, 2)(3,2)(3, 2)1/121/12

From fY1,Y2(y1,y2)f_{Y_1, Y_2}(y_1, y_2), we obtain fY1(y1)f_{Y_1}(y_1) by summing over y2y_2.

In general, let (X1,X2)(X_1, X_2) be a random vector with jpmf fX1,X2(x1,x2)f_{X_1, X_2}(x_1, x_2). S={(x1,x2):fX1,X2(x1,x2)>0}S = \{(x_1, x_2) : f_{X_1, X_2}(x_1, x_2) > 0\} is the support of (X1,X2)(X_1, X_2).

Let Y1=g1(X1,X2)Y_1 = g_1(X_1, X_2) and Y2=g2(X1,X2)Y_2 = g_2(X_1, X_2) define a one-to-one transformation that maps SS onto TT, where T={(y1,y2):fY1,Y2(y1,y2)>0}T = \{(y_1, y_2) : f_{Y_1, Y_2}(y_1, y_2) > 0\} is the support of (Y1,Y2)(Y_1, Y_2).

Then, the jpmf of (Y1,Y2)(Y_1, Y_2) is:

fY1,Y2(y1,y2)={fX1,X2(g11(y1,y2),g21(y1,y2)),for (y1,y2)T0,otherwisef_{Y_1, Y_2}(y_1, y_2) = \begin{cases} f_{X_1, X_2}(g_1^{-1}(y_1, y_2), g_2^{-1}(y_1, y_2)), & \text{for } (y_1, y_2) \in T \\ 0, & \text{otherwise} \end{cases}

Where x1=g11(y1,y2),x2=g21(y1,y2)x_1 = g_1^{-1}(y_1, y_2), x_2 = g_2^{-1}(y_1, y_2) is the single-valued inverse of y1=g1(x1,x2),y2=g2(x1,x2)y_1 = g_1(x_1, x_2), y_2 = g_2(x_1, x_2).

From fY1,Y2(y1,y2)f_{Y_1, Y_2}(y_1, y_2), we may obtain fY1(y1)f_{Y_1}(y_1) by summing over y2y_2 and fY2(y2)f_{Y_2}(y_2) by summing over y1y_1.

Example 2

Let X1X_1 and X2X_2 be 2 independent random variables having Poisson distributions with the parameters λ1\lambda_1 and λ2\lambda_2 respectively.

Find the probability distribution of the r.v. Y=X1+X2Y = X_1 + X_2.

Solution: Convolution of Poisson

By definition, YY takes values y{0,1,2,}y \in \{0, 1, 2, \dots\}.

Poisson Convolution Derivation
fY(y)\displaystyle f_Y(y)
=P(X1+X2=y)=x1=0yP(X1=x1,X2=yx1)\displaystyle = P(X_1 + X_2 = y) = \sum_{x_1=0}^y P(X_1 = x_1, X_2 = y - x_1)
=x1=0yfX1(x1)fX2(yx1)(since independent)\displaystyle = \sum_{x_1=0}^y f_{X_1}(x_1) f_{X_2}(y - x_1) \quad \text{(since independent)}
=x1=0yeλ1λ1x1x1!eλ2λ2yx1(yx1)!\displaystyle = \sum_{x_1=0}^y \frac{e^{-\lambda_1} \lambda_1^{x_1}}{x_1!} \frac{e^{-\lambda_2} \lambda_2^{y - x_1}}{(y - x_1)!}
=e(λ1+λ2)y!x1=0yy!x1!(yx1)!λ1x1λ2yx1(multiply & divide by y!)\displaystyle = \frac{e^{-(\lambda_1 + \lambda_2)}}{y!} \sum_{x_1=0}^y \frac{y!}{x_1!(y-x_1)!} \lambda_1^{x_1} \lambda_2^{y-x_1} \quad \text{(multiply \& divide by } y!)
=e(λ1+λ2)y!(λ1+λ2)y(via Binomial Theorem)\displaystyle = \frac{e^{-(\lambda_1 + \lambda_2)}}{y!} (\lambda_1 + \lambda_2)^y \quad \text{(via Binomial Theorem)}

This shows that YPoi(λ1+λ2)Y \sim Poi(\lambda_1 + \lambda_2).

Part III

Trivariate Case

Example

Let (X1,X2,X3)(X_1, X_2, X_3) have the joint discrete probability function given by:

(x1,x2,x3)(x_1, x_2, x_3)(0,0,0)(0,0,0)(0,0,1)(0,0,1)(0,1,1)(0,1,1)(1,0,1)(1,0,1)(1,1,0)(1,1,0)(1,1,1)(1,1,1)
f(x1,x2,x3)f(x_1, x_2, x_3)1/81/83/83/81/81/81/81/81/81/81/81/8

Find the joint probability function of Y1=X1+X2+X3Y_1 = X_1 + X_2 + X_3 & Y2=X3X2Y_2 = X_3 - X_2.

Solution

First, compute y1y_1 and y2y_2 for each possible state:

(x1,x2,x3)(x_1, x_2, x_3)(0,0,0)(0,0,0)(0,0,1)(0,0,1)(0,1,1)(0,1,1)(1,0,1)(1,0,1)(1,1,0)(1,1,0)(1,1,1)(1,1,1)
f(x1,x2,x3)f(x_1, x_2, x_3)1/81/83/83/81/81/81/81/81/81/81/81/8
y1=x1+x2+x3y_1 = x_1 + x_2 + x_3001122222233
y2=x3x2y_2 = x_3 - x_2001100111100

The joint pdf of (y1,y2)(y_1, y_2) is compiled by summing the probabilities for matching pairs of (y1,y2)(y_1, y_2):

y1\y2y_1 \backslash y_20011\sum
001/81/8001/81/8
11003/83/83/83/8
221/81/82/82/83/83/8
331/81/8001/81/8
\sum3/83/85/85/811