Skip to content

Transformations and Convolutions

5.1 Distribution of a Function of a Random Variable

Section titled “5.1 Distribution of a Function of a Random Variable”

Theorem 5.1 (CDF Method). If Y=g(X)Y = g(X) and gg is monotone, then

FY(y)=P(g(X)y)={FX(g1(y))ifg isincreasing1FX(g1(y))ifg isdecreasingF_Y(y) = P(g(X) \leq y) = \begin{cases} F_X(g^{-1}(y)) & \text{if} g \text{ is} increasing \\ 1 - F_X(g^{-1}(y)) & \text{if} g \text{ is} decreasing \end{cases}

Theorem 5.2 (Change of Variables). If Y=g(X)Y = g(X) where gg is differentiable and strictly monotone, then

fY(y)=fX(g1(y))ddyg1(y)f_Y(y) = f_X(g^{-1}(y)) \cdot \left|\frac{d}{dy} g^{-1}(y)\right|

Worked Example: Distribution of $X^2$ where $X \sim N(0, 1)$

Solution. Let Y=X2Y = X^2 where XN(0,1)X \sim N(0, 1). For y0y \geq 0:

FY(y)=P(X2y)=P(yXy)=Φ(y)Φ(y)=2Φ(y)1F_Y(y) = P(X^2 \leq y) = P(-\sqrt{y} \leq X \leq \sqrt{y}) = \Phi(\sqrt{y}) - \Phi(-\sqrt{y}) = 2\Phi(\sqrt{y}) - 1

fY(y)=ddy[2Φ(y)1]=2ϕ(y)12y=12πyey/2f_Y(y) = \frac{d}{dy}[2\Phi(\sqrt{y}) - 1] = 2\phi(\sqrt{y}) \cdot \frac{1}{2\sqrt{y}} = \frac{1}{\sqrt{2\pi y}}\, e^{-y/2}

This is the PDF of the χ2(1)\chi^2(1) distribution. \blacksquare

Theorem 5.3. If XX and YY are independent continuous random variables, the PDF of Z=X+YZ = X + Y is

fZ(z)=(fXfY)(z)=fX(x)fY(zx)dxf_Z(z) = (f_X * f_Y)(z) = \int_{-\infty}^{\infty} f_X(x)\, f_Y(z - x)\, dx

Proof. FZ(z)=P(X+Yz)=x+yzfX,Y(x,y)dxdy=fX(x)[zxfY(y)dy]dx=fX(x)FY(zx)dxF_Z(z) = P(X + Y \leq z) = \iint_{x+y \leq z} f_{X,Y}(x, y)\, dx\, dy = \int_{-\infty}^{\infty} f_X(x)\left[\int_{-\infty}^{z-x} f_Y(y)\, dy\right] dx = \int_{-\infty}^{\infty} f_X(x)\, F_Y(z - x)\, dx.

Differentiating: fZ(z)=fX(x)fY(zx)dxf_Z(z) = \int_{-\infty}^{\infty} f_X(x)\, f_Y(z - x)\, dx. \blacksquare

Corollary 5.4. The sum of independent normals is normal: if XN(μ1,σ12)X \sim N(\mu_1, \sigma_1^2) and YN(μ2,σ22)Y \sim N(\mu_2, \sigma_2^2) are independent, then X+YN(μ1+μ2,σ12+σ22)X + Y \sim N(\mu_1 + \mu_2, \sigma_1^2 + \sigma_2^2).

Proof. The convolution of two Gaussian PDFs is Gaussian. This follows from the MGF: MX+Y(t)=MX(t)MY(t)=exp((μ1+μ2)t+(σ12+σ22)t2/2)M_{X+Y}(t) = M_X(t)M_Y(t) = \exp((\mu_1 + \mu_2)t + (\sigma_1^2 + \sigma_2^2)t^2/2)Which is the MGF of N(μ1+μ2,σ12+σ22)N(\mu_1 + \mu_2, \sigma_1^2 + \sigma_2^2). \blacksquare

  • Confusing PDF and CDF. PDF f(x)f(x): probability density; CDF F(x)=P(Xx)=xf(t)dtF(x) = P(X \leq x) = \int_{-\infty}^x f(t)\, dt. Fix: F"(x)=f(x)F"(x) = f(x); P(a<X<b)=F(b)F(a)P(a < X < b) = F(b) - F(a).
  • Wrong central limit theorem application. The CLT applies to the sample mean, not individual observations, and requires sufficiently large nn. Fix: XˉndN(μ,σ2/n)\bar{X}_n \xrightarrow{d} N(\mu, \sigma^2/n) as nn \to \infty.
  • Confusing type I and type II errors. Type I: rejecting H0H_0 when it is true (α\alpha). Type II: failing to reject H0H_0 when it is false (β\beta). Fix: Type I = false positive; Type II = false negative. Decreasing one increases the other.

Problem. XN(100,152)X \sim N(100, 15^2). Find P(X>130)P(X > 130).

Solution. Z=13010015=2.0Z = \frac{130 - 100}{15} = 2.0. P(X>130)=P(Z>2)=1Φ(2)10.9772=0.0228P(X > 130) = P(Z > 2) = 1 - \Phi(2) \approx 1 - 0.9772 = 0.0228.

\blacksquare

Problem. Test H0:μ=50H_0: \mu = 50 vs H1:μ>50H_1: \mu > 50 given xˉ=53\bar{x} = 53, s=8s = 8, n=25n = 25, α=0.05\alpha = 0.05.

Solution. t=53508/25=31.6=1.875t = \frac{53 - 50}{8/\sqrt{25}} = \frac{3}{1.6} = 1.875. Critical value: t0.05,24=1.711t_{0.05, 24} = 1.711. Since 1.875>1.7111.875 > 1.711, reject H0H_0 at the 5% level.

\blacksquare

flowchart TD
A[5_Transformations And Convolutions] --> B[Key Definitions]
A --> C[Main Theorems]
A --> D[Applications]
B --> E[Core concepts]
C --> F[Proof techniques]
D --> G[Problem solving]
  • Continuous distributions: PDF integrates to 1; CDF gives cumulative probability.
  • Normal distribution: XN(μ,σ2)X \sim N(\mu, \sigma^2); standardise: Z=(Xμ)/σZ = (X - \mu)/\sigma.
  • Central limit theorem: sample mean is approximately normal for large nn.
  • Hypothesis testing: state H0H_0 and H1H_1, choose significance level, compute test statistic, compare with critical value.
TopicSiteLink
[Probability]A-LevelView
[Probability]IBView
[Probability]DSEView
[Probability]UniversityView

Convolution satisfies several algebraic properties that simplify calculations:

Commutativity: fXfY=fYfXf_X * f_Y = f_Y * f_X. The order of summation does not matter.

Associativity: (fXfY)fZ=fX(fYfZ)(f_X * f_Y) * f_Z = f_X * (f_Y * f_Z). Multiple sums can be evaluated in any order.

Scaling: If Z=aXZ = aX where a>0a > 0, then fZ(z)=1afX(z/a)f_Z(z) = \frac{1}{a}f_X(z/a).

Location-scale family: If XX has PDF f(x)f(x), then Y=aX+bY = aX + b has PDF fY(y)=1af ⁣(yba)f_Y(y) = \frac{1}{a}f\!\left(\frac{y - b}{a}\right).

5.4 Moment Generating Functions and Transformations

Section titled “5.4 Moment Generating Functions and Transformations”

The moment generating function (MGF) is MX(t)=E[etX]M_X(t) = E[e^{tX}].

Key property: If MX(t)=MY(t)M_X(t) = M_Y(t) for all tt in a neighbourhood of 0, then XX and YY have the same distribution (MGFs uniquely determine distributions).

Linear transformations: If Y=aX+bY = aX + b, then MY(t)=ebtMX(at)M_Y(t) = e^{bt}M_X(at).

Sums of independent variables: If XX and YY are independent, MX+Y(t)=MX(t)MY(t)M_{X+Y}(t) = M_X(t)\cdot M_Y(t).

Worked Example: MGF of the Uniform Distribution

Let XUniform(0,1)X \sim \text{Uniform}(0, 1). The MGF is:

MX(t)=E[etX]=01etxdx=et1t,t0M_X(t) = E[e^{tX}] = \int_0^1 e^{tx}\,dx = \frac{e^t - 1}{t}, \quad t \neq 0

Differentiating: MX(0)=E[X]=1/2M_X'(0) = E[X] = 1/2 and MX(0)=E[X2]=1/3M_X''(0) = E[X^2] = 1/3.

Thus Var(X)=1/31/4=1/12\text{Var}(X) = 1/3 - 1/4 = 1/12, confirming the known result. \blacksquare

When you transform a random variable, its distribution changes according to how the transformation stretches or compresses the number line. The change-of-variables formula accounts for the Jacobian — the local stretching factor. When you add two independent random variables, their distributions convolve: the density of the sum at a point is the integral of one density times a shifted version of the other, like sliding two waveforms past each other. Moment generating functions turn this convolution into multiplication, making sums of independent variables easy to handle. The central limit theorem is the grand payoff: no matter what distribution you start with, the sum of many independent copies approaches a normal distribution because convolutions smooth out irregularities.

OperationResulting DistributionKey Formula
Y=g(X)Y = g(X) monotonefY(y)=fX(g1(y))dg1/dyf_Y(y) = f_X(g^{-1}(y))\|d g^{-1}/dy\|Change of variables
Z=X+YZ = X + Y independentfZ=fXfYf_Z = f_X * f_YConvolution integral
Y=aX+bY = aX + bfY(y)=1afX ⁣(yba)f_Y(y) = \frac{1}{a}f_X\!\left(\frac{y-b}{a}\right)Location-scale
XN(μ1,σ12)X \sim N(\mu_1,\sigma_1^2), YN(μ2,σ22)Y \sim N(\mu_2,\sigma_2^2)X+YN(μ1+μ2,σ12+σ22)X+Y \sim N(\mu_1+\mu_2,\sigma_1^2+\sigma_2^2)Normal sum
XPoisson(λ1)X \sim \text{Poisson}(\lambda_1), YPoisson(λ2)Y \sim \text{Poisson}(\lambda_2)X+YPoisson(λ1+λ2)X+Y \sim \text{Poisson}(\lambda_1+\lambda_2)Poisson sum
MaX+b(t)M_{aX+b}(t)ebtMX(at)e^{bt}M_X(at)MGF transformation