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

  • 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

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