Theorem 5.1 (CDF Method). If Y = g ( X ) Y = g(X) Y = g ( X ) and g g g is monotone, then
F Y ( y ) = P ( g ( X ) ≤ y ) = { F X ( g − 1 ( y ) ) if g is i n c r e a s i n g 1 − F X ( g − 1 ( y ) ) if g is d e c r e a s i n g F_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} F Y ( y ) = P ( g ( X ) ≤ y ) = { F X ( g − 1 ( y )) 1 − F X ( g − 1 ( y )) if g is in cr e a s in g if g is d ecr e a s in g
Theorem 5.2 (Change of Variables). If Y = g ( X ) Y = g(X) Y = g ( X ) where g g g is differentiable and strictly monotone, then
f Y ( y ) = f X ( g − 1 ( y ) ) ⋅ ∣ d d y g − 1 ( y ) ∣ f_Y(y) = f_X(g^{-1}(y)) \cdot \left|\frac{d}{dy} g^{-1}(y)\right| f Y ( y ) = f X ( g − 1 ( y )) ⋅ d y d g − 1 ( y )
Worked Example: Distribution of $X^2$ where $X \sim N(0, 1)$ Solution. Let Y = X 2 Y = X^2 Y = X 2 where X ∼ N ( 0 , 1 ) X \sim N(0, 1) X ∼ N ( 0 , 1 ) . For y ≥ 0 y \geq 0 y ≥ 0 :
F Y ( y ) = P ( X 2 ≤ y ) = P ( − y ≤ X ≤ y ) = Φ ( y ) − Φ ( − y ) = 2 Φ ( y ) − 1 F_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 F Y ( y ) = P ( X 2 ≤ y ) = P ( − y ≤ X ≤ y ) = Φ ( y ) − Φ ( − y ) = 2Φ ( y ) − 1
f Y ( y ) = d d y [ 2 Φ ( y ) − 1 ] = 2 ϕ ( y ) ⋅ 1 2 y = 1 2 π y e − y / 2 f_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} f Y ( y ) = d y d [ 2Φ ( y ) − 1 ] = 2 ϕ ( y ) ⋅ 2 y 1 = 2 π y 1 e − y /2
This is the PDF of the χ 2 ( 1 ) \chi^2(1) χ 2 ( 1 ) distribution. ■ \blacksquare ■
Theorem 5.3. If X X X and Y Y Y are independent continuous random variables, the PDF of Z = X + Y Z = X + Y Z = X + Y is
f Z ( z ) = ( f X ∗ f Y ) ( z ) = ∫ − ∞ ∞ f X ( x ) f Y ( z − x ) d x f_Z(z) = (f_X * f_Y)(z) = \int_{-\infty}^{\infty} f_X(x)\, f_Y(z - x)\, dx f Z ( z ) = ( f X ∗ f Y ) ( z ) = ∫ − ∞ ∞ f X ( x ) f Y ( z − x ) d x
Proof. F Z ( z ) = P ( X + Y ≤ z ) = ∬ x + y ≤ z f X , Y ( x , y ) d x d y = ∫ − ∞ ∞ f X ( x ) [ ∫ − ∞ z − x f Y ( y ) d y ] d x = ∫ − ∞ ∞ f X ( x ) F Y ( z − x ) d x F_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 F Z ( z ) = P ( X + Y ≤ z ) = ∬ x + y ≤ z f X , Y ( x , y ) d x d y = ∫ − ∞ ∞ f X ( x ) [ ∫ − ∞ z − x f Y ( y ) d y ] d x = ∫ − ∞ ∞ f X ( x ) F Y ( z − x ) d x .
Differentiating: f Z ( z ) = ∫ − ∞ ∞ f X ( x ) f Y ( z − x ) d x f_Z(z) = \int_{-\infty}^{\infty} f_X(x)\, f_Y(z - x)\, dx f Z ( z ) = ∫ − ∞ ∞ f X ( x ) f Y ( z − x ) d x . ■ \blacksquare ■
Corollary 5.4. The sum of independent normals is normal: if X ∼ N ( μ 1 , σ 1 2 ) X \sim N(\mu_1, \sigma_1^2) X ∼ N ( μ 1 , σ 1 2 ) and Y ∼ N ( μ 2 , σ 2 2 ) Y \sim N(\mu_2, \sigma_2^2) Y ∼ N ( μ 2 , σ 2 2 ) are independent, then X + Y ∼ N ( μ 1 + μ 2 , σ 1 2 + σ 2 2 ) X + Y \sim N(\mu_1 + \mu_2, \sigma_1^2 + \sigma_2^2) X + Y ∼ N ( μ 1 + μ 2 , σ 1 2 + σ 2 2 ) .
Proof. The convolution of two Gaussian PDFs is Gaussian. This follows from the MGF: M X + Y ( t ) = M X ( t ) M Y ( t ) = exp ( ( μ 1 + μ 2 ) t + ( σ 1 2 + σ 2 2 ) t 2 / 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) M X + Y ( t ) = M X ( t ) M Y ( t ) = exp (( μ 1 + μ 2 ) t + ( σ 1 2 + σ 2 2 ) t 2 /2 ) Which is the MGF of N ( μ 1 + μ 2 , σ 1 2 + σ 2 2 ) N(\mu_1 + \mu_2, \sigma_1^2 + \sigma_2^2) N ( μ 1 + μ 2 , σ 1 2 + σ 2 2 ) . ■ \blacksquare ■
Confusing PDF and CDF. PDF f ( x ) f(x) f ( x ) : probability density; CDF F ( x ) = P ( X ≤ x ) = ∫ − ∞ x f ( t ) d t F(x) = P(X \leq x) = \int_{-\infty}^x f(t)\, dt F ( x ) = P ( X ≤ x ) = ∫ − ∞ x f ( t ) d t . Fix: F " ( x ) = f ( x ) F"(x) = f(x) F " ( x ) = f ( x ) ; P ( a < X < b ) = F ( b ) − F ( a ) 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 n n n . Fix: X ˉ n → d N ( μ , σ 2 / n ) \bar{X}_n \xrightarrow{d} N(\mu, \sigma^2/n) X ˉ n d N ( μ , σ 2 / n ) as n → ∞ n \to \infty n → ∞ .Confusing type I and type II errors. Type I: rejecting H 0 H_0 H 0 when it is true (α \alpha α ). Type II: failing to reject H 0 H_0 H 0 when it is false (β \beta β ). Fix: Type I = false positive; Type II = false negative. Decreasing one increases the other.Problem. X ∼ N ( 100 , 15 2 ) X \sim N(100, 15^2) X ∼ N ( 100 , 1 5 2 ) . Find P ( X > 130 ) P(X > 130) P ( X > 130 ) .
Solution. Z = 130 − 100 15 = 2.0 Z = \frac{130 - 100}{15} = 2.0 Z = 15 130 − 100 = 2.0 . P ( X > 130 ) = P ( Z > 2 ) = 1 − Φ ( 2 ) ≈ 1 − 0.9772 = 0.0228 P(X > 130) = P(Z > 2) = 1 - \Phi(2) \approx 1 - 0.9772 = 0.0228 P ( X > 130 ) = P ( Z > 2 ) = 1 − Φ ( 2 ) ≈ 1 − 0.9772 = 0.0228 .
■ \blacksquare ■
Problem. Test H 0 : μ = 50 H_0: \mu = 50 H 0 : μ = 50 vs H 1 : μ > 50 H_1: \mu > 50 H 1 : μ > 50 given x ˉ = 53 \bar{x} = 53 x ˉ = 53 , s = 8 s = 8 s = 8 , n = 25 n = 25 n = 25 , α = 0.05 \alpha = 0.05 α = 0.05 .
Solution. t = 53 − 50 8 / 25 = 3 1.6 = 1.875 t = \frac{53 - 50}{8/\sqrt{25}} = \frac{3}{1.6} = 1.875 t = 8/ 25 53 − 50 = 1.6 3 = 1.875 . Critical value: t 0.05 , 24 = 1.711 t_{0.05, 24} = 1.711 t 0.05 , 24 = 1.711 . Since 1.875 > 1.711 1.875 > 1.711 1.875 > 1.711 , reject H 0 H_0 H 0 at the 5% level.
■ \blacksquare ■
Continuous distributions: PDF integrates to 1; CDF gives cumulative probability. Normal distribution: X ∼ N ( μ , σ 2 ) X \sim N(\mu, \sigma^2) X ∼ N ( μ , σ 2 ) ; standardise: Z = ( X − μ ) / σ Z = (X - \mu)/\sigma Z = ( X − μ ) / σ . Central limit theorem: sample mean is approximately normal for large n n n . Hypothesis testing: state H 0 H_0 H 0 and H 1 H_1 H 1 , choose significance level, compute test statistic, compare with critical value. Topic Site Link [Probability] A-Level View [Probability] IB View [Probability] DSE View [Probability] University View
Convolution satisfies several algebraic properties that simplify calculations:
Commutativity: f X ∗ f Y = f Y ∗ f X f_X * f_Y = f_Y * f_X f X ∗ f Y = f Y ∗ f X . The order of summation does not matter.
Associativity: ( f X ∗ f Y ) ∗ f Z = f X ∗ ( f Y ∗ f Z ) (f_X * f_Y) * f_Z = f_X * (f_Y * f_Z) ( f X ∗ f Y ) ∗ f Z = f X ∗ ( f Y ∗ f Z ) . Multiple sums can be evaluated in any order.
Scaling: If Z = a X Z = aX Z = a X where a > 0 a > 0 a > 0 , then f Z ( z ) = 1 a f X ( z / a ) f_Z(z) = \frac{1}{a}f_X(z/a) f Z ( z ) = a 1 f X ( z / a ) .
Location-scale family: If X X X has PDF f ( x ) f(x) f ( x ) , then Y = a X + b Y = aX + b Y = a X + b has PDF f Y ( y ) = 1 a f ( y − b a ) f_Y(y) = \frac{1}{a}f\!\left(\frac{y - b}{a}\right) f Y ( y ) = a 1 f ( a y − b ) .
The moment generating function (MGF) is M X ( t ) = E [ e t X ] M_X(t) = E[e^{tX}] M X ( t ) = E [ e tX ] .
Key property: If M X ( t ) = M Y ( t ) M_X(t) = M_Y(t) M X ( t ) = M Y ( t ) for all t t t in a neighbourhood of 0, then X X X and Y Y Y have the same distribution (MGFs uniquely determine distributions).
Linear transformations: If Y = a X + b Y = aX + b Y = a X + b , then M Y ( t ) = e b t M X ( a t ) M_Y(t) = e^{bt}M_X(at) M Y ( t ) = e b t M X ( a t ) .
Sums of independent variables: If X X X and Y Y Y are independent, M X + Y ( t ) = M X ( t ) ⋅ M Y ( t ) M_{X+Y}(t) = M_X(t)\cdot M_Y(t) M X + Y ( t ) = M X ( t ) ⋅ M Y ( t ) .
Worked Example: MGF of the Uniform Distribution Let X ∼ Uniform ( 0 , 1 ) X \sim \text{Uniform}(0, 1) X ∼ Uniform ( 0 , 1 ) . The MGF is:
M X ( t ) = E [ e t X ] = ∫ 0 1 e t x d x = e t − 1 t , t ≠ 0 M_X(t) = E[e^{tX}] = \int_0^1 e^{tx}\,dx = \frac{e^t - 1}{t}, \quad t \neq 0 M X ( t ) = E [ e tX ] = ∫ 0 1 e t x d x = t e t − 1 , t = 0
Differentiating: M X ′ ( 0 ) = E [ X ] = 1 / 2 M_X'(0) = E[X] = 1/2 M X ′ ( 0 ) = E [ X ] = 1/2 and M X ′ ′ ( 0 ) = E [ X 2 ] = 1 / 3 M_X''(0) = E[X^2] = 1/3 M X ′′ ( 0 ) = E [ X 2 ] = 1/3 .
Thus Var ( X ) = 1 / 3 − 1 / 4 = 1 / 12 \text{Var}(X) = 1/3 - 1/4 = 1/12 Var ( X ) = 1/3 − 1/4 = 1/12 , confirming the known result. ■ \blacksquare ■
Operation Resulting Distribution Key Formula Y = g ( X ) Y = g(X) Y = g ( X ) monotonef Y ( y ) = f X ( g − 1 ( y ) ) ∥ d g − 1 / d y ∥ f_Y(y) = f_X(g^{-1}(y))\|d g^{-1}/dy\| f Y ( y ) = f X ( g − 1 ( y )) ∥ d g − 1 / d y ∥ Change of variables Z = X + Y Z = X + Y Z = X + Y independentf Z = f X ∗ f Y f_Z = f_X * f_Y f Z = f X ∗ f Y Convolution integral Y = a X + b Y = aX + b Y = a X + b f Y ( y ) = 1 a f X ( y − b a ) f_Y(y) = \frac{1}{a}f_X\!\left(\frac{y-b}{a}\right) f Y ( y ) = a 1 f X ( a y − b ) Location-scale X ∼ N ( μ 1 , σ 1 2 ) X \sim N(\mu_1,\sigma_1^2) X ∼ N ( μ 1 , σ 1 2 ) , Y ∼ N ( μ 2 , σ 2 2 ) Y \sim N(\mu_2,\sigma_2^2) Y ∼ N ( μ 2 , σ 2 2 ) X + Y ∼ N ( μ 1 + μ 2 , σ 1 2 + σ 2 2 ) X+Y \sim N(\mu_1+\mu_2,\sigma_1^2+\sigma_2^2) X + Y ∼ N ( μ 1 + μ 2 , σ 1 2 + σ 2 2 ) Normal sum X ∼ Poisson ( λ 1 ) X \sim \text{Poisson}(\lambda_1) X ∼ Poisson ( λ 1 ) , Y ∼ Poisson ( λ 2 ) Y \sim \text{Poisson}(\lambda_2) Y ∼ Poisson ( λ 2 ) X + Y ∼ Poisson ( λ 1 + λ 2 ) X+Y \sim \text{Poisson}(\lambda_1+\lambda_2) X + Y ∼ Poisson ( λ 1 + λ 2 ) Poisson sum M a X + b ( t ) M_{aX+b}(t) M a X + b ( t ) e b t M X ( a t ) e^{bt}M_X(at) e b t M X ( a t ) MGF transformation