This topic explores fundamental concepts that shape our understanding of the world.
The exponential distribution models the time between events in a Poisson process, while the theory Of continuous random variables extends probability to quantities that can take any value in an Interval.
Board Paper Notes AQA Paper 2 Continuous RVs; limited exponential coverage Edexcel S3, S4 Exponential distribution in S4; continuous RVs in S3 OCR (A) Paper 2 Continuous RVs and exponential CIE (9231) S2 Both continuous RVs and exponential covered
Both are memoryless. The Poisson process links all three distributions: Poisson counts events, Exponential measures inter-arrival times, and geometric counts trials until the first event. Definition. A probability density function (PDF) f ( x ) f(x) f ( x ) of a continuous random variable X X X Is a non-negative function satisfying:
f ( x ) ≥ 0 f o r a l l x , ∫ − ∞ ∞ f ( x ) d x = 1 f(x) \geq 0 \quad \mathrm{for all } x, \qquad \int_{-\infty}^{\infty}f(x)\,dx = 1 f ( x ) ≥ 0 forall x , ∫ − ∞ ∞ f ( x ) d x = 1
Probabilities are found by integration:
P ( a ≤ X ≤ b ) = ∫ a b f ( x ) d x P(a \leq X \leq b) = \int_a^b f(x)\,dx P ( a ≤ X ≤ b ) = ∫ a b f ( x ) d x
$P(a \leq X \leq b) = P(a < X < b)$ — the inequalities at individual points do not matter. ### 1.2 Cumulative distribution function
Definition. The cumulative distribution function (CDF) is
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
Properties:
F ( − ∞ ) = 0 F(-\infty) = 0 F ( − ∞ ) = 0 , F ( ∞ ) = 1 F(\infty) = 1 F ( ∞ ) = 1 F F F is non-decreasingf ( x ) = F " ( x ) f(x) = F"(x) f ( x ) = F " ( x ) where F F F is differentiableP ( a < X ≤ b ) = F ( b ) − F ( a ) P(a < X \leq b) = F(b) - F(a) P ( a < X ≤ b ) = F ( b ) − F ( a ) Definition. The expected value of a continuous random variable X X X is
E ( X ) = ∫ − ∞ ∞ x f ( x ) d x \boxed{E(X) = \int_{-\infty}^{\infty}x\,f(x)\,dx} E ( X ) = ∫ − ∞ ∞ x f ( x ) d x
For a function g ( X ) g(X) g ( X ) :
E ( g ( X ) ) = ∫ − ∞ ∞ g ( x ) f ( x ) d x E(g(X)) = \int_{-\infty}^{\infty}g(x)\,f(x)\,dx E ( g ( X )) = ∫ − ∞ ∞ g ( x ) f ( x ) d x
Definition.
V a r ( X ) = E ( X 2 ) − [ E ( X ) ] 2 = ∫ − ∞ ∞ x 2 f ( x ) d x − ( ∫ − ∞ ∞ x f ( x ) d x ) 2 \boxed{\mathrm{Var}(X) = E(X^2) - [E(X)]^2 = \int_{-\infty}^{\infty}x^2\,f(x)\,dx - \left(\int_{-\infty}^{\infty}x\,f(x)\,dx\right)^2} Var ( X ) = E ( X 2 ) − [ E ( X ) ] 2 = ∫ − ∞ ∞ x 2 f ( x ) d x − ( ∫ − ∞ ∞ x f ( x ) d x ) 2
The linear properties carry over from the discrete case:
E ( a X + b ) = a E ( X ) + b , V a r ( a X + b ) = a 2 V a r ( X ) E(aX + b) = aE(X) + b, \qquad \mathrm{Var}(aX + b) = a^2\,\mathrm{Var}(X) E ( a X + b ) = a E ( X ) + b , Var ( a X + b ) = a 2 Var ( X )
Definition. The median m m m satisfies F ( m ) = 0.5 F(m) = 0.5 F ( m ) = 0.5 I.e., ∫ − ∞ m f ( x ) d x = 0.5 \int_{-\infty}^{m}f(x)\,dx = 0.5 ∫ − ∞ m f ( x ) d x = 0.5 .
Definition. The mode is the value of x x x at which f ( x ) f(x) f ( x ) is maximised.
Definition. The lower quartile Q 1 Q_1 Q 1 satisfies F ( Q 1 ) = 0.25 F(Q_1) = 0.25 F ( Q 1 ) = 0.25 and the upper quartile Q 3 Q_3 Q 3 satisfies F ( Q 3 ) = 0.75 F(Q_3) = 0.75 F ( Q 3 ) = 0.75 .
The interquartile range is I Q R = Q 3 − Q 1 \mathrm{IQR} = Q_3 - Q_1 IQR = Q 3 − Q 1 .
Definition. A continuous random variable X X X follows an exponential distribution with rate Parameter λ \lambda λ (where λ > 0 \lambda > 0 λ > 0 ), written X ∼ E x p ( λ ) X \sim \mathrm{Exp}(\lambda) X ∼ Exp ( λ ) If
f ( x ) = λ e − λ x , x ≥ 0 \boxed{f(x) = \lambda e^{-\lambda x}, \quad x \geq 0} f ( x ) = λ e − λ x , x ≥ 0
And f ( x ) = 0 f(x) = 0 f ( x ) = 0 for x < 0 x < 0 x < 0 .
F ( x ) = P ( X ≤ x ) = ∫ 0 x λ e − λ t d t = [ − e − λ t ] 0 x = 1 − e − λ x F(x) = P(X \leq x) = \int_0^x \lambda e^{-\lambda t}\,dt = \left[-e^{-\lambda t}\right]_0^x = 1 - e^{-\lambda x} F ( x ) = P ( X ≤ x ) = ∫ 0 x λ e − λ t d t = [ − e − λ t ] 0 x = 1 − e − λ x
F ( x ) = 1 − e − λ x , x ≥ 0 \boxed{F(x) = 1 - e^{-\lambda x}, \quad x \geq 0} F ( x ) = 1 − e − λ x , x ≥ 0
E ( X ) = ∫ 0 ∞ x ⋅ λ e − λ x d x \begin{aligned} E(X) &= \int_0^{\infty}x\cdot\lambda e^{-\lambda x}\,dx \end{aligned} E ( X ) = ∫ 0 ∞ x ⋅ λ e − λ x d x Using integration by parts with u = x u = x u = x , d v = λ e − λ x d x dv = \lambda e^{-\lambda x}\,dx d v = λ e − λ x d x :
d u = d x du = dx d u = d x , v = − e − λ x v = -e^{-\lambda x} v = − e − λ x .
E ( X ) = [ − x e − λ x ] 0 ∞ + ∫ 0 ∞ e − λ x d x = lim x → ∞ ( − x e − λ x ) + 0 + [ − 1 λ e − λ x ] 0 ∞ = 0 + 1 λ = 1 λ ■ \begin{aligned} E(X) &= \left[-xe^{-\lambda x}\right]_0^{\infty} + \int_0^{\infty}e^{-\lambda x}\,dx \\ &= \lim_{x\to\infty}(-xe^{-\lambda x}) + 0 + \left[-\frac{1}{\lambda}e^{-\lambda x}\right]_0^{\infty} \\ &= 0 + \frac{1}{\lambda} = \frac{1}{\lambda} \quad \blacksquare \end{aligned} E ( X ) = [ − x e − λ x ] 0 ∞ + ∫ 0 ∞ e − λ x d x = x → ∞ lim ( − x e − λ x ) + 0 + [ − λ 1 e − λ x ] 0 ∞ = 0 + λ 1 = λ 1 ■ Note: lim x → ∞ x e − λ x = 0 \lim_{x\to\infty}xe^{-\lambda x} = 0 lim x → ∞ x e − λ x = 0 by L’Hôpital’s rule (exponential decay dominates).
First compute E ( X 2 ) E(X^2) E ( X 2 ) :
E ( X 2 ) = ∫ 0 ∞ x 2 ⋅ λ e − λ x d x E(X^2) = \int_0^{\infty}x^2\cdot\lambda e^{-\lambda x}\,dx E ( X 2 ) = ∫ 0 ∞ x 2 ⋅ λ e − λ x d x
Integration by parts twice with u = x 2 u = x^2 u = x 2 , d v = λ e − λ x d x dv = \lambda e^{-\lambda x}\,dx d v = λ e − λ x d x :
d u = 2 x d x du = 2x\,dx d u = 2 x d x , v = − e − λ x v = -e^{-\lambda x} v = − e − λ x .
E ( X 2 ) = [ − x 2 e − λ x ] 0 ∞ + ∫ 0 ∞ 2 x e − λ x d x = 0 + 2 ⋅ 1 λ ⋅ E ( X ) = 2 λ ⋅ 1 λ = 2 λ 2 \begin{aligned} E(X^2) &= \left[-x^2 e^{-\lambda x}\right]_0^{\infty} + \int_0^{\infty}2x\,e^{-\lambda x}\,dx \\ &= 0 + 2\cdot\frac{1}{\lambda}\cdot E(X) = \frac{2}{\lambda}\cdot\frac{1}{\lambda} = \frac{2}{\lambda^2} \end{aligned} E ( X 2 ) = [ − x 2 e − λ x ] 0 ∞ + ∫ 0 ∞ 2 x e − λ x d x = 0 + 2 ⋅ λ 1 ⋅ E ( X ) = λ 2 ⋅ λ 1 = λ 2 2 V a r ( X ) = E ( X 2 ) − [ E ( X ) ] 2 = 2 λ 2 − 1 λ 2 = 1 λ 2 ■ \mathrm{Var}(X) = E(X^2) - [E(X)]^2 = \frac{2}{\lambda^2} - \frac{1}{\lambda^2} = \frac{1}{\lambda^2} \quad \blacksquare Var ( X ) = E ( X 2 ) − [ E ( X ) ] 2 = λ 2 2 − λ 2 1 = λ 2 1 ■
E ( X ) = 1 λ , V a r ( X ) = 1 λ 2 , σ = 1 λ \boxed{E(X) = \frac{1}{\lambda}, \qquad \mathrm{Var}(X) = \frac{1}{\lambda^2}, \qquad \sigma = \frac{1}{\lambda}} E ( X ) = λ 1 , Var ( X ) = λ 2 1 , σ = λ 1
Theorem. The exponential distribution is the only continuous memoryless distribution:
P ( X > s + t ∣ X > s ) = P ( X > t ) P(X > s + t \mid X > s) = P(X > t) P ( X > s + t ∣ X > s ) = P ( X > t )
P ( X > s + t ∣ X > s ) = P ( X > s + t ) P ( X > s ) = e − λ ( s + t ) e − λ s = e − λ t = P ( X > t ) ■ \begin{aligned} P(X > s + t \mid X > s) &= \frac{P(X > s+t)}{P(X > s)} = \frac{e^{-\lambda(s+t)}}{e^{-\lambda s}} \\ &= e^{-\lambda t} = P(X > t) \quad \blacksquare \end{aligned} P ( X > s + t ∣ X > s ) = P ( X > s ) P ( X > s + t ) = e − λ s e − λ ( s + t ) = e − λ t = P ( X > t ) ■ This uses P ( X > x ) = 1 − F ( x ) = e − λ x P(X > x) = 1 - F(x) = e^{-\lambda x} P ( X > x ) = 1 − F ( x ) = e − λ x .
Exponentially distributed lifetime has been working for $s$ hours, the remaining lifetime has the Same distribution as a brand new component. This means exponential lifetimes imply no "wear out" Effect — which is why it is more appropriate for electronic components than mechanical ones. ### 2.6 Link to Poisson processes
A Poisson process with rate λ \lambda λ satisfies:
The number of events in an interval of length t t t follows P o ( λ t ) \mathrm{Po}(\lambda t) Po ( λ t ) The time between consecutive events follows E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) The inter-arrival times are independent and identically distributed Proof sketch that inter-arrival times are exponential. Let T T T be the time until the first Event. P ( T > t ) = P ( n o e v e n t s i n [ 0 , t ] ) = P ( N ( t ) = 0 ) = e − λ t ( λ t ) 0 0 ! = e − λ t P(T > t) = P(\mathrm{no events in }[0,t]) = P(N(t) = 0) = \dfrac{e^{-\lambda t}(\lambda t)^0}{0!} = e^{-\lambda t} P ( T > t ) = P ( noeventsin [ 0 , t ]) = P ( N ( t ) = 0 ) = 0 ! e − λ t ( λ t ) 0 = e − λ t .
So P ( T ≤ t ) = 1 − e − λ t P(T \leq t) = 1 - e^{-\lambda t} P ( T ≤ t ) = 1 − e − λ t Which is the CDF of E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) . ■ \blacksquare ■
For the exponential distribution:
F ( x ) = 1 − e − λ x = p ⟹ x = − 1 λ ln ( 1 − p ) F(x) = 1 - e^{-\lambda x} = p \implies x = -\frac{1}{\lambda}\ln(1-p) F ( x ) = 1 − e − λ x = p ⟹ x = − λ 1 ln ( 1 − p )
The median is x 0.5 = − ln ( 0.5 ) λ = ln 2 λ x_{0.5} = -\dfrac{\ln(0.5)}{\lambda} = \dfrac{\ln 2}{\lambda} x 0.5 = − λ ln ( 0.5 ) = λ ln 2 .
Example. X ∼ E x p ( 0.5 ) X \sim \mathrm{Exp}(0.5) X ∼ Exp ( 0.5 ) . Find P ( X > 3 ) P(X > 3) P ( X > 3 ) , P ( 1 < X < 4 ) P(1 < X < 4) P ( 1 < X < 4 ) And the median.
P ( X > 3 ) = e − 0.5 × 3 = e − 1.5 ≈ 0.2231 P(X > 3) = e^{-0.5 \times 3} = e^{-1.5} \approx 0.2231 P ( X > 3 ) = e − 0.5 × 3 = e − 1.5 ≈ 0.2231 .
P ( 1 < X < 4 ) = F ( 4 ) − F ( 1 ) = ( 1 − e − 2 ) − ( 1 − e − 0.5 ) = e − 0.5 − e − 2 ≈ 0.6065 − 0.1353 = 0.4712 P(1 < X < 4) = F(4) - F(1) = (1-e^{-2}) - (1-e^{-0.5}) = e^{-0.5} - e^{-2} \approx 0.6065 - 0.1353 = 0.4712 P ( 1 < X < 4 ) = F ( 4 ) − F ( 1 ) = ( 1 − e − 2 ) − ( 1 − e − 0.5 ) = e − 0.5 − e − 2 ≈ 0.6065 − 0.1353 = 0.4712 .
Median = ln 2 0.5 = 2 ln 2 ≈ 1.386 = \dfrac{\ln 2}{0.5} = 2\ln 2 \approx 1.386 = 0.5 ln 2 = 2 ln 2 ≈ 1.386 .
Example. X X X has PDF f ( x ) = 3 8 x 2 f(x) = \dfrac{3}{8}x^2 f ( x ) = 8 3 x 2 for 0 ≤ x ≤ 2 0 \leq x \leq 2 0 ≤ x ≤ 2 and f ( x ) = 0 f(x) = 0 f ( x ) = 0 otherwise.
Verify: ∫ 0 2 3 8 x 2 d x = 3 8 ⋅ 8 3 = 1 \int_0^2 \dfrac{3}{8}x^2\,dx = \dfrac{3}{8}\cdot\dfrac{8}{3} = 1 ∫ 0 2 8 3 x 2 d x = 8 3 ⋅ 3 8 = 1 . ✓ \checkmark ✓
E ( X ) = ∫ 0 2 x ⋅ 3 8 x 2 d x = 3 8 ∫ 0 2 x 3 d x = 3 8 ⋅ 16 4 = 3 2 E(X) = \int_0^2 x\cdot\dfrac{3}{8}x^2\,dx = \dfrac{3}{8}\int_0^2 x^3\,dx = \dfrac{3}{8}\cdot\dfrac{16}{4} = \dfrac{3}{2} E ( X ) = ∫ 0 2 x ⋅ 8 3 x 2 d x = 8 3 ∫ 0 2 x 3 d x = 8 3 ⋅ 4 16 = 2 3 .
E ( X 2 ) = 3 8 ∫ 0 2 x 4 d x = 3 8 ⋅ 32 5 = 12 5 E(X^2) = \dfrac{3}{8}\int_0^2 x^4\,dx = \dfrac{3}{8}\cdot\dfrac{32}{5} = \dfrac{12}{5} E ( X 2 ) = 8 3 ∫ 0 2 x 4 d x = 8 3 ⋅ 5 32 = 5 12 .
V a r ( X ) = 12 5 − ( 3 2 ) 2 = 12 5 − 9 4 = 3 20 \mathrm{Var}(X) = \dfrac{12}{5} - \left(\dfrac{3}{2}\right)^2 = \dfrac{12}{5} - \dfrac{9}{4} = \dfrac{3}{20} Var ( X ) = 5 12 − ( 2 3 ) 2 = 5 12 − 4 9 = 20 3 .
Median: ∫ 0 m 3 8 x 2 d x = 0.5 ⟹ m 3 8 = 0.5 ⟹ m 3 = 4 ⟹ m = 4 3 ≈ 1.587 \int_0^m \dfrac{3}{8}x^2\,dx = 0.5 \implies \dfrac{m^3}{8} = 0.5 \implies m^3 = 4 \implies m = \sqrt[3]{4} \approx 1.587 ∫ 0 m 8 3 x 2 d x = 0.5 ⟹ 8 m 3 = 0.5 ⟹ m 3 = 4 ⟹ m = 3 4 ≈ 1.587 .
Example. f ( x ) = 2 x f(x) = 2x f ( x ) = 2 x for 0 ≤ x ≤ 1 0 \leq x \leq 1 0 ≤ x ≤ 1 . Find F ( x ) F(x) F ( x ) .
For 0 ≤ x ≤ 1 0 \leq x \leq 1 0 ≤ x ≤ 1 : F ( x ) = ∫ 0 x 2 t d t = x 2 F(x) = \int_0^x 2t\,dt = x^2 F ( x ) = ∫ 0 x 2 t d t = x 2 .
For x < 0 x < 0 x < 0 : F ( x ) = 0 F(x) = 0 F ( x ) = 0 . For x > 1 x > 1 x > 1 : F ( x ) = 1 F(x) = 1 F ( x ) = 1 .
Example. The lifetime of a component is modelled by X ∼ E x p ( λ ) X \sim \mathrm{Exp}(\lambda) X ∼ Exp ( λ ) . A sample of 10 components gives a mean lifetime of 420 hours. Test at the 5% level whether λ = 0.005 \lambda = 0.005 λ = 0.005 Against H 1 : λ ≠ 0.005 H_1: \lambda \neq 0.005 H 1 : λ = 0.005 .
Under H 0 H_0 H 0 : E ( X ) = 1 / λ = 200 E(X) = 1/\lambda = 200 E ( X ) = 1/ λ = 200 hours. Since n n n is large, use the approximate normal Distribution of X ˉ \bar{X} X ˉ :
X ˉ ∼ N ( 1 λ , 1 n λ 2 ) = N ( 200 , 4000 ) \bar{X} \sim N\!\left(\frac{1}{\lambda}, \frac{1}{n\lambda^2}\right) = N(200, 4000) X ˉ ∼ N ( λ 1 , n λ 2 1 ) = N ( 200 , 4000 ) Approximately.
z = 420 − 200 4000 = 220 63.25 = 3.48 z = \dfrac{420 - 200}{\sqrt{4000}} = \dfrac{220}{63.25} = 3.48 z = 4000 420 − 200 = 63.25 220 = 3.48 .
∣ z ∣ = 3.48 > 1.96 |z| = 3.48 > 1.96 ∣ z ∣ = 3.48 > 1.96 So reject H 0 H_0 H 0 .
Problem 1 $X \sim \mathrm{Exp}(2)$. Find $P(X > 1)$, $P(0.5 < X < 2)$And the 90th percentile.Solution 1 $P(X > 1) = e^{-2(1)} = e^{-2} \approx 0.1353$.P ( 0.5 < X < 2 ) = ( 1 − e − 4 ) − ( 1 − e − 1 ) = e − 1 − e − 4 ≈ 0.3679 − 0.0183 = 0.3496 P(0.5 < X < 2) = (1-e^{-4}) - (1-e^{-1}) = e^{-1} - e^{-4} \approx 0.3679 - 0.0183 = 0.3496 P ( 0.5 < X < 2 ) = ( 1 − e − 4 ) − ( 1 − e − 1 ) = e − 1 − e − 4 ≈ 0.3679 − 0.0183 = 0.3496 .
90th percentile: F ( x ) = 0.9 ⟹ 1 − e − 2 x = 0.9 ⟹ x = − ln ( 0.1 ) 2 ≈ 1.151 F(x) = 0.9 \implies 1 - e^{-2x} = 0.9 \implies x = -\dfrac{\ln(0.1)}{2} \approx 1.151 F ( x ) = 0.9 ⟹ 1 − e − 2 x = 0.9 ⟹ x = − 2 ln ( 0.1 ) ≈ 1.151 .
If you get this wrong, revise: Percentiles — Section 2.7.
Problem 2 A continuous random variable $X$ has PDF $f(x) = \dfrac{3x^2}{8}$ for $0 \leq x \leq 2$. Find $E(X)$, $\mathrm{Var}(X)$And the median.Solution 2 $E(X) = \int_0^2 x\cdot\dfrac{3x^2}{8}\,dx = \dfrac{3}{8}\cdot\left[\dfrac{x^4}{4}\right]_0^2 = \dfrac{3}{8}\cdot 4 = 1.5$.E ( X 2 ) = ∫ 0 2 x 2 ⋅ 3 x 2 8 d x = 3 8 ⋅ [ x 5 5 ] 0 2 = 3 8 ⋅ 32 5 = 2.4 E(X^2) = \int_0^2 x^2\cdot\dfrac{3x^2}{8}\,dx = \dfrac{3}{8}\cdot\left[\dfrac{x^5}{5}\right]_0^2 = \dfrac{3}{8}\cdot\dfrac{32}{5} = 2.4 E ( X 2 ) = ∫ 0 2 x 2 ⋅ 8 3 x 2 d x = 8 3 ⋅ [ 5 x 5 ] 0 2 = 8 3 ⋅ 5 32 = 2.4 .
V a r ( X ) = 2.4 − 1.5 2 = 2.4 − 2.25 = 0.15 \mathrm{Var}(X) = 2.4 - 1.5^2 = 2.4 - 2.25 = 0.15 Var ( X ) = 2.4 − 1. 5 2 = 2.4 − 2.25 = 0.15 .
Median: ∫ 0 m 3 x 2 8 d x = 0.5 ⟹ m 3 8 = 0.5 ⟹ m = 4 3 ≈ 1.587 \int_0^m \dfrac{3x^2}{8}\,dx = 0.5 \implies \dfrac{m^3}{8} = 0.5 \implies m = \sqrt[3]{4} \approx 1.587 ∫ 0 m 8 3 x 2 d x = 0.5 ⟹ 8 m 3 = 0.5 ⟹ m = 3 4 ≈ 1.587 .
If you get this wrong, revise: Median, mode, and quartiles — Section 1.5.
Problem 3 Prove the memoryless property of the exponential distribution.Solution 3 $P(X > s+t \mid X > s) = \dfrac{P(X > s+t)}{P(X > s)} = \dfrac{e^{-\lambda(s+t)}}{e^{-\lambda s}} = e^{-\lambda t} = P(X > t)$. $\blacksquare$This uses the survival function P ( X > x ) = e − λ x P(X > x) = e^{-\lambda x} P ( X > x ) = e − λ x .
If you get this wrong, revise: The memoryless property — Section 2.5.
Problem 4 Calls arrive at a switchboard as a Poisson process with rate $\lambda = 4$ per hour. Find the probability that the time between two consecutive calls exceeds 30 minutes.Solution 4 The inter-arrival time $T \sim \mathrm{Exp}(4)$ (rate in hours).P ( T > 0.5 ) = e − 4 × 0.5 = e − 2 ≈ 0.1353 P(T > 0.5) = e^{-4 \times 0.5} = e^{-2} \approx 0.1353 P ( T > 0.5 ) = e − 4 × 0.5 = e − 2 ≈ 0.1353 .
If you get this wrong, revise: Link to Poisson processes — Section 2.6.
Problem 5 $X$ has PDF $f(x) = \dfrac{1}{2}x$ for $0 \leq x \leq 2$. Find the CDF, $E(X)$And $\mathrm{Var}(X)$.Solution 5 CDF: $F(x) = \int_0^x \dfrac{t}{2}\,dt = \dfrac{x^2}{4}$ for $0 \leq x \leq 2$. $F(x) = 0$ for $x < 0$, $F(x) = 1$ for $x > 2$.E ( X ) = ∫ 0 2 x ⋅ x 2 d x = 1 2 [ x 3 3 ] 0 2 = 4 3 E(X) = \int_0^2 x\cdot\dfrac{x}{2}\,dx = \dfrac{1}{2}\left[\dfrac{x^3}{3}\right]_0^2 = \dfrac{4}{3} E ( X ) = ∫ 0 2 x ⋅ 2 x d x = 2 1 [ 3 x 3 ] 0 2 = 3 4 .
E ( X 2 ) = 1 2 [ x 4 4 ] 0 2 = 1 2 ⋅ 4 = 2 E(X^2) = \dfrac{1}{2}\left[\dfrac{x^4}{4}\right]_0^2 = \dfrac{1}{2}\cdot 4 = 2 E ( X 2 ) = 2 1 [ 4 x 4 ] 0 2 = 2 1 ⋅ 4 = 2 .
V a r ( X ) = 2 − ( 4 3 ) 2 = 2 − 16 9 = 2 9 \mathrm{Var}(X) = 2 - \left(\dfrac{4}{3}\right)^2 = 2 - \dfrac{16}{9} = \dfrac{2}{9} Var ( X ) = 2 − ( 3 4 ) 2 = 2 − 9 16 = 9 2 .
If you get this wrong, revise: Expected value — Section 1.3.
Problem 6 The lifetime of a light bulb follows $X \sim \mathrm{Exp}(0.01)$ (in hours). Given that the bulb has been working for 500 hours, find the probability it lasts at least another 200 hours.Solution 6 By the memoryless property: $P(X > 500+200 \mid X > 500) = P(X > 200) = e^{-0.01 \times 200} = e^{-2} \approx 0.1353$.If you get this wrong, revise: The memoryless property — Section 2.5.
Problem 7 A continuous random variable $X$ has CDF $F(x) = \dfrac{x^3}{27}$ for $0 \leq x \leq 3$. Find the PDF, $E(X)$And the upper quartile.Solution 7 PDF: $f(x) = F'(x) = \dfrac{x^2}{9}$ for $0 \leq x \leq 3$.E ( X ) = ∫ 0 3 x ⋅ x 2 9 d x = 1 9 [ x 4 4 ] 0 3 = 1 9 ⋅ 81 4 = 9 4 = 2.25 E(X) = \int_0^3 x\cdot\dfrac{x^2}{9}\,dx = \dfrac{1}{9}\left[\dfrac{x^4}{4}\right]_0^3 = \dfrac{1}{9}\cdot\dfrac{81}{4} = \dfrac{9}{4} = 2.25 E ( X ) = ∫ 0 3 x ⋅ 9 x 2 d x = 9 1 [ 4 x 4 ] 0 3 = 9 1 ⋅ 4 81 = 4 9 = 2.25 .
Upper quartile: F ( Q 3 ) = 0.75 ⟹ Q 3 3 27 = 0.75 ⟹ Q 3 3 = 20.25 ⟹ Q 3 ≈ 2.725 F(Q_3) = 0.75 \implies \dfrac{Q_3^3}{27} = 0.75 \implies Q_3^3 = 20.25 \implies Q_3 \approx 2.725 F ( Q 3 ) = 0.75 ⟹ 27 Q 3 3 = 0.75 ⟹ Q 3 3 = 20.25 ⟹ Q 3 ≈ 2.725 .
If you get this wrong, revise: Cumulative distribution function — Section 1.2.
Problem 8 Prove that $E(X) = 1/\lambda$ for $X \sim \mathrm{Exp}(\lambda)$Using integration by parts.Solution 8 $E(X) = \int_0^{\infty}x\lambda e^{-\lambda x}\,dx$.Let u = x$$dv = \lambda e^{-\lambda x}\,dx So du = dx$$v = -e^{-\lambda x} .
E ( X ) = [ − x e − λ x ] 0 ∞ + ∫ 0 ∞ e − λ x d x = 0 + [ − 1 λ e − λ x ] 0 ∞ = 1 λ E(X) = \left[-xe^{-\lambda x}\right]_0^{\infty} + \int_0^{\infty}e^{-\lambda x}\,dx = 0 + \left[-\dfrac{1}{\lambda}e^{-\lambda x}\right]_0^{\infty} = \dfrac{1}{\lambda} E ( X ) = [ − x e − λ x ] 0 ∞ + ∫ 0 ∞ e − λ x d x = 0 + [ − λ 1 e − λ x ] 0 ∞ = λ 1 . ■ \blacksquare ■
If you get this wrong, revise: Proof that E ( X ) = 1 λ E(X) = \frac{1}{\lambda} E ( X ) = λ 1 — Section 2.3.
Problem 9 Buses arrive at a stop as a Poisson process with rate 6 per hour. Find the probability that a passenger waits between 5 and 15 minutes for a bus.Solution 9 Waiting time $T \sim \mathrm{Exp}(6)$ (rate per hour).P ( 1 / 12 < T < 1 / 4 ) = F ( 1 / 4 ) − F ( 1 / 12 ) = ( 1 − e − 1.5 ) − ( 1 − e − 0.5 ) = e − 0.5 − e − 1.5 P(1/12 < T < 1/4) = F(1/4) - F(1/12) = (1-e^{-1.5}) - (1-e^{-0.5}) = e^{-0.5} - e^{-1.5} P ( 1/12 < T < 1/4 ) = F ( 1/4 ) − F ( 1/12 ) = ( 1 − e − 1.5 ) − ( 1 − e − 0.5 ) = e − 0.5 − e − 1.5
≈ 0.6065 − 0.2231 = 0.3834 \approx 0.6065 - 0.2231 = 0.3834 ≈ 0.6065 − 0.2231 = 0.3834 .
If you get this wrong, revise: Link to Poisson processes — Section 2.6.
Problem 10 $X$ has PDF $f(x) = 4x^3$ for $0 \leq x \leq 1$. Find $P(X > 0.5)$$E(X)$$\mathrm{Var}(X)$And the mode.Solution 10 $P(X > 0.5) = \int_{0.5}^1 4x^3\,dx = \left[x^4\right]_{0.5}^1 = 1 - 0.0625 = 0.9375$.E ( X ) = ∫ 0 1 4 x 4 d x = [ 4 x 5 5 ] 0 1 = 4 5 E(X) = \int_0^1 4x^4\,dx = \left[\dfrac{4x^5}{5}\right]_0^1 = \dfrac{4}{5} E ( X ) = ∫ 0 1 4 x 4 d x = [ 5 4 x 5 ] 0 1 = 5 4 .
E ( X 2 ) = ∫ 0 1 4 x 5 d x = [ 4 x 6 6 ] 0 1 = 2 3 E(X^2) = \int_0^1 4x^5\,dx = \left[\dfrac{4x^6}{6}\right]_0^1 = \dfrac{2}{3} E ( X 2 ) = ∫ 0 1 4 x 5 d x = [ 6 4 x 6 ] 0 1 = 3 2 .
V a r ( X ) = 2 3 − ( 4 5 ) 2 = 2 3 − 16 25 = 50 − 48 75 = 2 75 \mathrm{Var}(X) = \dfrac{2}{3} - \left(\dfrac{4}{5}\right)^2 = \dfrac{2}{3} - \dfrac{16}{25} = \dfrac{50 - 48}{75} = \dfrac{2}{75} Var ( X ) = 3 2 − ( 5 4 ) 2 = 3 2 − 25 16 = 75 50 − 48 = 75 2 .
Mode: f ( x ) = 4 x 3 f(x) = 4x^3 f ( x ) = 4 x 3 is increasing on [ 0 , 1 ] [0,1] [ 0 , 1 ] So the mode is at x = 1 x = 1 x = 1 .
If you get this wrong, revise: Median, mode, and quartiles — Section 1.5.
Theorem. If X ∼ E x p ( λ ) X \sim \mathrm{Exp}(\lambda) X ∼ Exp ( λ ) Then for all s , t > 0 s, t > 0 s , t > 0 :
P ( X > s + t ∣ X > s ) = P ( X > t ) P(X > s + t \mid X > s) = P(X > t) P ( X > s + t ∣ X > s ) = P ( X > t )
By definition of conditional probability:
P ( X > s + t ∣ X > s ) = P ( X > s + t ∩ X > s ) P ( X > s ) = P ( X > s + t ) P ( X > s ) P(X > s + t \mid X > s) = \frac{P(X > s + t \,\cap\, X > s)}{P(X > s)} = \frac{P(X > s + t)}{P(X > s)} P ( X > s + t ∣ X > s ) = P ( X > s ) P ( X > s + t ∩ X > s ) = P ( X > s ) P ( X > s + t )
Since X > s + t X > s + t X > s + t implies X > s X > s X > s .
Using the survival function S ( x ) = P ( X > x ) = e − λ x S(x) = P(X > x) = e^{-\lambda x} S ( x ) = P ( X > x ) = e − λ x :
P ( X > s + t ) P ( X > s ) = e − λ ( s + t ) e − λ s = e − λ t = P ( X > t ) ■ \frac{P(X > s + t)}{P(X > s)} = \frac{e^{-\lambda(s+t)}}{e^{-\lambda s}} = e^{-\lambda t} = P(X > t) \quad \blacksquare P ( X > s ) P ( X > s + t ) = e − λ s e − λ ( s + t ) = e − λ t = P ( X > t ) ■
Theorem. If a continuous random variable X X X on ( 0 , ∞ ) (0, \infty) ( 0 , ∞ ) satisfies P ( X > s + t ∣ X > s ) = P ( X > t ) P(X > s+t \mid X > s) = P(X > t) P ( X > s + t ∣ X > s ) = P ( X > t ) for all s , t > 0 s, t > 0 s , t > 0 Then X ∼ E x p ( λ ) X \sim \mathrm{Exp}(\lambda) X ∼ Exp ( λ ) for some λ > 0 \lambda > 0 λ > 0 .
Let G ( t ) = P ( X > t ) G(t) = P(X > t) G ( t ) = P ( X > t ) . The memoryless condition gives:
G ( s + t ) = G ( s ) G ( t ) G(s + t) = G(s)G(t) G ( s + t ) = G ( s ) G ( t )
This is Cauchy’s functional equation. Since G G G is non-increasing and 0 ≤ G ≤ 1 0 \leq G \leq 1 0 ≤ G ≤ 1 The only Solutions are:
G ( t ) = e − λ t G(t) = e^{-\lambda t} G ( t ) = e − λ t
For some λ ≥ 0 \lambda \geq 0 λ ≥ 0 . Since G G G is non-trivial (not identically 1), λ > 0 \lambda > 0 λ > 0 . Therefore:
P ( X ≤ t ) = 1 − e − λ t P(X \leq t) = 1 - e^{-\lambda t} P ( X ≤ t ) = 1 − e − λ t
Which is the CDF of E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) . ■ \blacksquare ■
The memoryless property means:
If a light bulb has been on for 100 hours, the probability it lasts another 50 hours is the same as a new bulb lasting 50 hours. If you have waited 20 minutes for a bus, your expected additional wait time is the same as if you had just arrived. This property makes exponential models appropriate for random failure mechanisms (electronic components) but inappropriate for wear-out mechanisms (mechanical parts). Theorem. In a Poisson process with rate λ \lambda λ The time between consecutive events follows E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) .
Let T T T be the time from an arbitrary starting point until the first event.
P ( T > t ) = P ( n o e v e n t s i n [ 0 , t ] ) P(T > t) = P(\mathrm{no events in }[0,t]) P ( T > t ) = P ( noeventsin [ 0 , t ])
Since the number of events in [ 0 , t ] [0,t] [ 0 , t ] follows P o ( λ t ) \mathrm{Po}(\lambda t) Po ( λ t ) :
P ( N ( t ) = 0 ) = e − λ t ( λ t ) 0 0 ! = e − λ t P(N(t) = 0) = \frac{e^{-\lambda t}(\lambda t)^0}{0!} = e^{-\lambda t} P ( N ( t ) = 0 ) = 0 ! e − λ t ( λ t ) 0 = e − λ t
Therefore P ( T ≤ t ) = 1 − e − λ t P(T \leq t) = 1 - e^{-\lambda t} P ( T ≤ t ) = 1 − e − λ t Which is the CDF of E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) . ■ \blacksquare ■
If T 1 , T 2 , … , T n T_1, T_2, \ldots, T_n T 1 , T 2 , … , T n are n n n independent inter-arrival times, each ∼ E x p ( λ ) \sim \mathrm{Exp}(\lambda) ∼ Exp ( λ ) Then the total time until the n n n -th event is:
S n = T 1 + T 2 + ⋯ + T n ∼ G a m m a ( n , λ ) S_n = T_1 + T_2 + \cdots + T_n \sim \mathrm{Gamma}(n, \lambda) S n = T 1 + T 2 + ⋯ + T n ∼ Gamma ( n , λ )
This connects the exponential to the gamma distribution.
Example. Calls arrive at a call centre as a Poisson process at rate 5 per hour.
(a) Find the probability that the time between two consecutive calls exceeds 20 minutes.
T ∼ E x p ( 5 ) T \sim \mathrm{Exp}(5) T ∼ Exp ( 5 ) (rate per hour). P ( T > 1 / 3 ) = e − 5 / 3 ≈ 0.1889 P(T > 1/3) = e^{-5/3} \approx 0.1889 P ( T > 1/3 ) = e − 5/3 ≈ 0.1889 .
(b) Find the probability that at least 3 calls arrive in the next 30 minutes.
N ( 0.5 ) ∼ P o ( 2.5 ) N(0.5) \sim \mathrm{Po}(2.5) N ( 0.5 ) ∼ Po ( 2.5 ) . P ( N ≥ 3 ) = 1 − P ( N ≤ 2 ) = 1 − e − 2.5 ( 1 + 2.5 + 2.5 2 / 2 ) = 1 − 0.0821 × 9.125 ≈ 1 − 0.749 = 0.251 P(N \geq 3) = 1 - P(N \leq 2) = 1 - e^{-2.5}(1 + 2.5 + 2.5^2/2) = 1 - 0.0821 \times 9.125 \approx 1 - 0.749 = 0.251 P ( N ≥ 3 ) = 1 − P ( N ≤ 2 ) = 1 − e − 2.5 ( 1 + 2.5 + 2. 5 2 /2 ) = 1 − 0.0821 × 9.125 ≈ 1 − 0.749 = 0.251 .
(c) Find the median inter-arrival time.
Median = ln 2 λ = ln 2 5 ≈ 0.139 h o u r s ≈ 8.3 m i n u t e s = \dfrac{\ln 2}{\lambda} = \dfrac{\ln 2}{5} \approx 0.139\,\mathrm{hours} \approx 8.3\,\mathrm{minutes} = λ ln 2 = 5 ln 2 ≈ 0.139 hours ≈ 8.3 minutes .
Definition. X ∼ U ( a , b ) X \sim U(a, b) X ∼ U ( a , b ) if:
f ( x ) = 1 b − a , a ≤ x ≤ b f(x) = \frac{1}{b - a}, \quad a \leq x \leq b f ( x ) = b − a 1 , a ≤ x ≤ b
And f ( x ) = 0 f(x) = 0 f ( x ) = 0 otherwise.
E ( X ) = ∫ a b x ⋅ 1 b − a d x = 1 b − a [ x 2 2 ] a b = b 2 − a 2 2 ( b − a ) = ( b − a ) ( b + a ) 2 ( b − a ) = a + b 2 ■ E(X) = \int_a^b x \cdot \frac{1}{b-a}\,dx = \frac{1}{b-a}\left[\frac{x^2}{2}\right]_a^b = \frac{b^2 - a^2}{2(b-a)} = \frac{(b-a)(b+a)}{2(b-a)} = \frac{a+b}{2} \quad \blacksquare E ( X ) = ∫ a b x ⋅ b − a 1 d x = b − a 1 [ 2 x 2 ] a b = 2 ( b − a ) b 2 − a 2 = 2 ( b − a ) ( b − a ) ( b + a ) = 2 a + b ■
E ( X 2 ) = ∫ a b x 2 ⋅ 1 b − a d x = 1 b − a [ x 3 3 ] a b = b 3 − a 3 3 ( b − a ) = a 2 + a b + b 2 3 E(X^2) = \int_a^b x^2 \cdot \frac{1}{b-a}\,dx = \frac{1}{b-a}\left[\frac{x^3}{3}\right]_a^b = \frac{b^3 - a^3}{3(b-a)} = \frac{a^2 + ab + b^2}{3} E ( X 2 ) = ∫ a b x 2 ⋅ b − a 1 d x = b − a 1 [ 3 x 3 ] a b = 3 ( b − a ) b 3 − a 3 = 3 a 2 + ab + b 2
V a r ( X ) = a 2 + a b + b 2 3 − ( a + b ) 2 4 = 4 ( a 2 + a b + b 2 ) − 3 ( a + b ) 2 12 \mathrm{Var}(X) = \frac{a^2 + ab + b^2}{3} - \frac{(a+b)^2}{4} = \frac{4(a^2 + ab + b^2) - 3(a+b)^2}{12} Var ( X ) = 3 a 2 + ab + b 2 − 4 ( a + b ) 2 = 12 4 ( a 2 + ab + b 2 ) − 3 ( a + b ) 2
= 4 a 2 + 4 a b + 4 b 2 − 3 a 2 − 6 a b − 3 b 2 12 = a 2 − 2 a b + b 2 12 = \frac{4a^2 + 4ab + 4b^2 - 3a^2 - 6ab - 3b^2}{12} = \frac{a^2 - 2ab + b^2}{12} = 12 4 a 2 + 4 ab + 4 b 2 − 3 a 2 − 6 ab − 3 b 2 = 12 a 2 − 2 ab + b 2
V a r ( X ) = ( b − a ) 2 12 ■ \boxed{\mathrm{Var}(X) = \frac{(b-a)^2}{12}} \quad \blacksquare Var ( X ) = 12 ( b − a ) 2 ■
F ( x ) = { 0 x < a x − a b − a a ≤ x ≤ b 1 x > b F(x) = \begin{cases} 0 & x \lt a \\ \dfrac{x - a}{b - a} & a \leq x \leq b \\ 1 & x > b \end{cases} F ( x ) = ⎩ ⎨ ⎧ 0 b − a x − a 1 x < a a ≤ x ≤ b x > b
A mixture distribution arises when a random variable is selected from one of several Sub-populations. If X X X comes from distribution 1 with probability p p p and from distribution 2 with Probability 1 − p 1-p 1 − p :
f ( x ) = p f 1 ( x ) + ( 1 − p ) f 2 ( x ) f(x) = p\,f_1(x) + (1-p)\,f_2(x) f ( x ) = p f 1 ( x ) + ( 1 − p ) f 2 ( x )
E ( X ) = p E ( X 1 ) + ( 1 − p ) E ( X 2 ) E(X) = p\,E(X_1) + (1-p)\,E(X_2) E ( X ) = p E ( X 1 ) + ( 1 − p ) E ( X 2 )
V a r ( X ) = p V a r ( X 1 ) + ( 1 − p ) V a r ( X 2 ) + p ( 1 − p ) [ E ( X 1 ) − E ( X 2 ) ] 2 \mathrm{Var}(X) = p\,\mathrm{Var}(X_1) + (1-p)\,\mathrm{Var}(X_2) + p(1-p)[E(X_1) - E(X_2)]^2 Var ( X ) = p Var ( X 1 ) + ( 1 − p ) Var ( X 2 ) + p ( 1 − p ) [ E ( X 1 ) − E ( X 2 ) ] 2
The variance formula includes an extra term from the difference in means — this is the law of total Variance.
Example. A machine produces components. With probability 0.7 it is correctly calibrated, Producing components with lifetime ∼ E x p ( 0.01 ) \sim \mathrm{Exp}(0.01) ∼ Exp ( 0.01 ) . With probability 0.3 it is faulty, Producing components with lifetime ∼ E x p ( 0.002 ) \sim \mathrm{Exp}(0.002) ∼ Exp ( 0.002 ) . Find the overall PDF, the expected Lifetime, and P ( X > 100 ) P(X > 100) P ( X > 100 ) .
f ( x ) = 0.7 ( 0.01 e − 0.01 x ) + 0.3 ( 0.002 e − 0.002 x ) = 0.007 e − 0.01 x + 0.0006 e − 0.002 x f(x) = 0.7(0.01\,e^{-0.01x}) + 0.3(0.002\,e^{-0.002x}) = 0.007\,e^{-0.01x} + 0.0006\,e^{-0.002x} f ( x ) = 0.7 ( 0.01 e − 0.01 x ) + 0.3 ( 0.002 e − 0.002 x ) = 0.007 e − 0.01 x + 0.0006 e − 0.002 x
E ( X ) = 0.7 × 100 + 0.3 × 500 = 70 + 150 = 220 h o u r s E(X) = 0.7 \times 100 + 0.3 \times 500 = 70 + 150 = 220\,\mathrm{hours} E ( X ) = 0.7 × 100 + 0.3 × 500 = 70 + 150 = 220 hours .
P ( X > 100 ) = 0.7 e − 0.01 × 100 + 0.3 e − 0.002 × 100 = 0.7 e − 1 + 0.3 e − 0.2 ≈ 0.7 ( 0.3679 ) + 0.3 ( 0.8187 ) P(X > 100) = 0.7\,e^{-0.01 \times 100} + 0.3\,e^{-0.002 \times 100} = 0.7e^{-1} + 0.3e^{-0.2} \approx 0.7(0.3679) + 0.3(0.8187) P ( X > 100 ) = 0.7 e − 0.01 × 100 + 0.3 e − 0.002 × 100 = 0.7 e − 1 + 0.3 e − 0.2 ≈ 0.7 ( 0.3679 ) + 0.3 ( 0.8187 )
≈ 0.2575 + 0.2456 = 0.5031 \approx 0.2575 + 0.2456 = 0.5031 ≈ 0.2575 + 0.2456 = 0.5031 .
The PDF f ( x ) f(x) f ( x ) gives the density of probability at x x x . It is not a probability itself — f ( x ) f(x) f ( x ) Can be greater than 1. The CDF F ( x ) F(x) F ( x ) gives the accumulated probability up to x x x And always Satisfies 0 ≤ F ( x ) ≤ 1 0 \leq F(x) \leq 1 0 ≤ F ( x ) ≤ 1 .
Common error: writing P ( X = a ) = f ( a ) P(X = a) = f(a) P ( X = a ) = f ( a ) for a continuous RV. This is wrong — P ( X = a ) = 0 P(X = a) = 0 P ( X = a ) = 0 always. Probabilities are areas under the PDF, not values of the PDF.
For a discrete RV, P ( X = a ) > 0 P(X = a) > 0 P ( X = a ) > 0 for specific values, and probabilities sum to 1.
For a continuous RV, P ( X = a ) = 0 P(X = a) = 0 P ( X = a ) = 0 for any single value, and probabilities integrate to 1.
Never use summation for a continuous RV or integration for a discrete RV (unless using the CDF Formalism).
X ∼ E x p ( λ ) X \sim \mathrm{Exp}(\lambda) X ∼ Exp ( λ ) has mean 1 / λ 1/\lambda 1/ λ . A common error is to confuse λ \lambda λ (the Rate) with the mean. If the mean lifetime is 200 hours, then λ = 1 / 200 = 0.005 \lambda = 1/200 = 0.005 λ = 1/200 = 0.005 Not λ = 200 \lambda = 200 λ = 200 .
Check: a larger λ \lambda λ means shorter lifetimes on average (events happen more frequently).
If events occur at rate λ = 5 \lambda = 5 λ = 5 per hour, then the inter-arrival time is E x p ( 5 ) \mathrm{Exp}(5) Exp ( 5 ) Measured in hours . If you want the probability of waiting more than 20 minutes, convert to Hours: t = 1 / 3 t = 1/3 t = 1/3 hours. Using t = 20 t = 20 t = 20 directly would give P ( T > 20 ) = e − 100 P(T > 20) = e^{-100} P ( T > 20 ) = e − 100 Which is Essentially zero and wrong.
Q1. $X \sim \mathrm{Exp}(\lambda)$. Find the value of $\lambda$ such that $P(X > 2) = 0.3$And hence find $E(X)$ and the 80th percentile. P ( X > 2 ) = e − 2 λ = 0.3 ⟹ − 2 λ = ln ( 0.3 ) ⟹ λ = − ln ( 0.3 ) 2 = 1.204 2 = 0.602 P(X > 2) = e^{-2\lambda} = 0.3 \implies -2\lambda = \ln(0.3) \implies \lambda = \dfrac{-\ln(0.3)}{2} = \dfrac{1.204}{2} = 0.602 P ( X > 2 ) = e − 2 λ = 0.3 ⟹ − 2 λ = ln ( 0.3 ) ⟹ λ = 2 − ln ( 0.3 ) = 2 1.204 = 0.602 .
E ( X ) = 1 / 0.602 ≈ 1.661 E(X) = 1/0.602 \approx 1.661 E ( X ) = 1/0.602 ≈ 1.661 .
80th percentile: F ( x ) = 0.8 ⟹ 1 − e − 0.602 x = 0.8 ⟹ e − 0.602 x = 0.2 ⟹ x = − ln ( 0.2 ) 0.602 = 1.609 0.602 ≈ 2.673 F(x) = 0.8 \implies 1 - e^{-0.602x} = 0.8 \implies e^{-0.602x} = 0.2 \implies x = \dfrac{-\ln(0.2)}{0.602} = \dfrac{1.609}{0.602} \approx 2.673 F ( x ) = 0.8 ⟹ 1 − e − 0.602 x = 0.8 ⟹ e − 0.602 x = 0.2 ⟹ x = 0.602 − ln ( 0.2 ) = 0.602 1.609 ≈ 2.673 .
Q2. A continuous random variable $X$ has PDF $f(x) = kx(4-x)$ for $0 \leq x \leq 4$. Find $k$, $E(X)$, $\mathrm{Var}(X)$And the median. ∫ 0 4 k x ( 4 − x ) d x = k ∫ 0 4 ( 4 x − x 2 ) d x = k [ 2 x 2 − x 3 3 ] 0 4 = k ( 32 − 64 / 3 ) = k ( 32 / 3 ) = 1 ⟹ k = 3 / 32 \int_0^4 kx(4-x)\,dx = k\int_0^4 (4x - x^2)\,dx = k\left[2x^2 - \dfrac{x^3}{3}\right]_0^4 = k(32 - 64/3) = k(32/3) = 1 \implies k = 3/32 ∫ 0 4 k x ( 4 − x ) d x = k ∫ 0 4 ( 4 x − x 2 ) d x = k [ 2 x 2 − 3 x 3 ] 0 4 = k ( 32 − 64/3 ) = k ( 32/3 ) = 1 ⟹ k = 3/32 .
E ( X ) = 3 32 ∫ 0 4 ( 4 x 2 − x 3 ) d x = 3 32 [ 4 x 3 3 − x 4 4 ] 0 4 = 3 32 ( 256 3 − 64 ) = 3 32 ( 256 − 192 3 ) = 64 32 = 2 E(X) = \dfrac{3}{32}\int_0^4 (4x^2 - x^3)\,dx = \dfrac{3}{32}\left[\dfrac{4x^3}{3} - \dfrac{x^4}{4}\right]_0^4 = \dfrac{3}{32}\left(\dfrac{256}{3} - 64\right) = \dfrac{3}{32}\left(\dfrac{256 - 192}{3}\right) = \dfrac{64}{32} = 2 E ( X ) = 32 3 ∫ 0 4 ( 4 x 2 − x 3 ) d x = 32 3 [ 3 4 x 3 − 4 x 4 ] 0 4 = 32 3 ( 3 256 − 64 ) = 32 3 ( 3 256 − 192 ) = 32 64 = 2 .
E ( X 2 ) = 3 32 ∫ 0 4 ( 4 x 3 − x 4 ) d x = 3 32 [ x 4 − x 5 5 ] 0 4 = 3 32 ( 256 − 204.8 ) = 3 × 51.2 32 = 4.8 E(X^2) = \dfrac{3}{32}\int_0^4 (4x^3 - x^4)\,dx = \dfrac{3}{32}\left[x^4 - \dfrac{x^5}{5}\right]_0^4 = \dfrac{3}{32}(256 - 204.8) = \dfrac{3 \times 51.2}{32} = 4.8 E ( X 2 ) = 32 3 ∫ 0 4 ( 4 x 3 − x 4 ) d x = 32 3 [ x 4 − 5 x 5 ] 0 4 = 32 3 ( 256 − 204.8 ) = 32 3 × 51.2 = 4.8 .
V a r ( X ) = 4.8 − 4 = 0.8 \mathrm{Var}(X) = 4.8 - 4 = 0.8 Var ( X ) = 4.8 − 4 = 0.8 .
Median: 3 32 ∫ 0 m ( 4 x − x 2 ) d x = 0.5 ⟹ 3 32 ( 2 m 2 − m 3 3 ) = 0.5 \dfrac{3}{32}\int_0^m (4x - x^2)\,dx = 0.5 \implies \dfrac{3}{32}\left(2m^2 - \dfrac{m^3}{3}\right) = 0.5 32 3 ∫ 0 m ( 4 x − x 2 ) d x = 0.5 ⟹ 32 3 ( 2 m 2 − 3 m 3 ) = 0.5 .
2 m 2 − m 3 3 = 16 3 ⟹ 6 m 2 − m 3 = 16 2m^2 - \dfrac{m^3}{3} = \dfrac{16}{3} \implies 6m^2 - m^3 = 16 2 m 2 − 3 m 3 = 3 16 ⟹ 6 m 2 − m 3 = 16 .
By inspection: m = 2 m = 2 m = 2 gives 24 − 8 = 16 24 - 8 = 16 24 − 8 = 16 . So the median is 2 2 2 (equal to the mean, reflecting the Symmetry of the PDF about x = 2 x = 2 x = 2 ).
Q3. Emails arrive at a server as a Poisson process with rate 12 per hour. Find the probability that the time between two consecutive emails is between 2 and 5 minutes, and find the probability that the third email arrives within 10 minutes. Inter-arrival time T ∼ E x p ( 12 ) T \sim \mathrm{Exp}(12) T ∼ Exp ( 12 ) per hour.
P ( 2 / 60 < T < 5 / 60 ) = e − 12 ( 2 / 60 ) − e − 12 ( 5 / 60 ) = e − 0.4 − e − 1 ≈ 0.6703 − 0.3679 = 0.3024 P(2/60 \lt T \lt 5/60) = e^{-12(2/60)} - e^{-12(5/60)} = e^{-0.4} - e^{-1} \approx 0.6703 - 0.3679 = 0.3024 P ( 2/60 < T < 5/60 ) = e − 12 ( 2/60 ) − e − 12 ( 5/60 ) = e − 0.4 − e − 1 ≈ 0.6703 − 0.3679 = 0.3024 .
For the third email: the total time S 3 = T 1 + T 2 + T 3 ∼ G a m m a ( 3 , 12 ) S_3 = T_1 + T_2 + T_3 \sim \mathrm{Gamma}(3, 12) S 3 = T 1 + T 2 + T 3 ∼ Gamma ( 3 , 12 ) .
Alternatively, use the Poisson count: P ( a t l e a s t 3 i n 10 m i n ) = P ( N ( 1 / 6 ) ≥ 3 ) P(\mathrm{at least 3 in 10 min}) = P(N(1/6) \geq 3) P ( atleast3in10min ) = P ( N ( 1/6 ) ≥ 3 ) where N ( 1 / 6 ) ∼ P o ( 2 ) N(1/6) \sim \mathrm{Po}(2) N ( 1/6 ) ∼ Po ( 2 ) .
= 1 − P ( N ≤ 2 ) = 1 − e − 2 ( 1 + 2 + 2 ) = 1 − 5 e − 2 ≈ 1 − 0.6767 = 0.3233 = 1 - P(N \leq 2) = 1 - e^{-2}(1 + 2 + 2) = 1 - 5e^{-2} \approx 1 - 0.6767 = 0.3233 = 1 − P ( N ≤ 2 ) = 1 − e − 2 ( 1 + 2 + 2 ) = 1 − 5 e − 2 ≈ 1 − 0.6767 = 0.3233 .
Q4. $X \sim U(0, a)$. Given that $E(X) = 3$ and $\mathrm{Var}(X) = 3$Find $a$ and the 90th percentile. E ( X ) = a / 2 = 3 ⟹ a = 6 E(X) = a/2 = 3 \implies a = 6 E ( X ) = a /2 = 3 ⟹ a = 6 .
Check: V a r ( X ) = ( b − a ) 2 / 12 = 36 / 12 = 3 \mathrm{Var}(X) = (b-a)^2/12 = 36/12 = 3 Var ( X ) = ( b − a ) 2 /12 = 36/12 = 3 . ✓ \checkmark ✓
90th percentile: F ( x ) = ( x − 0 ) / 6 = 0.9 ⟹ x = 5.4 F(x) = (x - 0)/6 = 0.9 \implies x = 5.4 F ( x ) = ( x − 0 ) /6 = 0.9 ⟹ x = 5.4 .
Q5. The lifetime of Component A follows $\mathrm{Exp}(0.02)$ and Component B follows $\mathrm{Exp}(0.05)$ (in hours). They are connected in series, so the system fails when either component fails. Assuming independence, find the PDF of the system lifetime, the expected system lifetime, and $P(\mathrm{system lasts} > 50\,\mathrm{hours})$. System lifetime T = min ( X A , X B ) T = \min(X_A, X_B) T = min ( X A , X B ) .
P ( T > t ) = P ( X A > t ) P ( X B > t ) = e − 0.02 t ⋅ e − 0.05 t = e − 0.07 t P(T > t) = P(X_A > t)P(X_B > t) = e^{-0.02t} \cdot e^{-0.05t} = e^{-0.07t} P ( T > t ) = P ( X A > t ) P ( X B > t ) = e − 0.02 t ⋅ e − 0.05 t = e − 0.07 t .
So T ∼ E x p ( 0.07 ) T \sim \mathrm{Exp}(0.07) T ∼ Exp ( 0.07 ) .
E ( T ) = 1 / 0.07 ≈ 14.3 h o u r s E(T) = 1/0.07 \approx 14.3\,\mathrm{hours} E ( T ) = 1/0.07 ≈ 14.3 hours .
P ( T > 50 ) = e − 0.07 × 50 = e − 3.5 ≈ 0.0302 P(T > 50) = e^{-0.07 \times 50} = e^{-3.5} \approx 0.0302 P ( T > 50 ) = e − 0.07 × 50 = e − 3.5 ≈ 0.0302 .
Note: the minimum of independent exponential RVs is itself exponential, with rate equal to the sum Of the individual rates.
Q6. A random variable $X$ has PDF $f(x) = \dfrac{2x}{9}$ for $0 \leq x \leq 3$. Find the CDF, $E(X)$, $\mathrm{Var}(X)$The median, and $P(1 \lt X \lt 2)$. CDF: F ( x ) = ∫ 0 x 2 t 9 d t = x 2 9 F(x) = \int_0^x \dfrac{2t}{9}\,dt = \dfrac{x^2}{9} F ( x ) = ∫ 0 x 9 2 t d t = 9 x 2 for 0 ≤ x ≤ 3 0 \leq x \leq 3 0 ≤ x ≤ 3 . F ( x ) = 0 F(x) = 0 F ( x ) = 0 for x < 0 x \lt 0 x < 0 , F ( x ) = 1 F(x) = 1 F ( x ) = 1 for x > 3 x > 3 x > 3 .
E ( X ) = ∫ 0 3 x ⋅ 2 x 9 d x = 2 9 [ x 3 3 ] 0 3 = 2 9 × 9 = 2 E(X) = \int_0^3 x \cdot \dfrac{2x}{9}\,dx = \dfrac{2}{9}\left[\dfrac{x^3}{3}\right]_0^3 = \dfrac{2}{9} \times 9 = 2 E ( X ) = ∫ 0 3 x ⋅ 9 2 x d x = 9 2 [ 3 x 3 ] 0 3 = 9 2 × 9 = 2 .
E ( X 2 ) = 2 9 [ x 4 4 ] 0 3 = 2 9 × 81 4 = 9 2 = 4.5 E(X^2) = \dfrac{2}{9}\left[\dfrac{x^4}{4}\right]_0^3 = \dfrac{2}{9} \times \dfrac{81}{4} = \dfrac{9}{2} = 4.5 E ( X 2 ) = 9 2 [ 4 x 4 ] 0 3 = 9 2 × 4 81 = 2 9 = 4.5 .
V a r ( X ) = 4.5 − 4 = 0.5 \mathrm{Var}(X) = 4.5 - 4 = 0.5 Var ( X ) = 4.5 − 4 = 0.5 .
Median: m 2 9 = 0.5 ⟹ m 2 = 4.5 ⟹ m = 4.5 ≈ 2.12 \dfrac{m^2}{9} = 0.5 \implies m^2 = 4.5 \implies m = \sqrt{4.5} \approx 2.12 9 m 2 = 0.5 ⟹ m 2 = 4.5 ⟹ m = 4.5 ≈ 2.12 .
P ( 1 < X < 2 ) = F ( 2 ) − F ( 1 ) = 4 9 − 1 9 = 1 3 P(1 \lt X \lt 2) = F(2) - F(1) = \dfrac{4}{9} - \dfrac{1}{9} = \dfrac{1}{3} P ( 1 < X < 2 ) = F ( 2 ) − F ( 1 ) = 9 4 − 9 1 = 3 1 .
Problem. The lifetime T T T of a component follows an exponential distribution with mean 200 Hours. Given that the component has survived 150 hours, find the probability it survives a further 100 hours.
Solution. By the memoryless property of the exponential distribution:
P ( T > 150 + 100 ∣ T > 150 ) = P ( T > 100 ) P(T > 150+100 \mid T > 150) = P(T > 100) P ( T > 150 + 100 ∣ T > 150 ) = P ( T > 100 )
λ = 1 200 = 0.005 \lambda = \dfrac{1}{200} = 0.005 λ = 200 1 = 0.005 .
P ( T > 100 ) = e − 0.005 × 100 = e − 0.5 = 0.607 P(T > 100) = e^{-0.005 \times 100} = e^{-0.5} = \boxed{0.607} P ( T > 100 ) = e − 0.005 × 100 = e − 0.5 = 0.607 (3 s.f.).
Problem. X ∼ U ( 0 , 10 ) X \sim \mathrm{U}(0, 10) X ∼ U ( 0 , 10 ) . Find P ( X > 6 ∣ X > 3 ) P(X > 6 \mid X > 3) P ( X > 6 ∣ X > 3 ) .
Solution. P ( X > 6 ∣ X > 3 ) = P ( X > 6 ) P ( X > 3 ) = 0.4 0.7 = 4 7 ≈ 0.571 P(X > 6 \mid X > 3) = \dfrac{P(X > 6)}{P(X > 3)} = \dfrac{0.4}{0.7} = \dfrac{4}{7} \approx \boxed{0.571} P ( X > 6 ∣ X > 3 ) = P ( X > 3 ) P ( X > 6 ) = 0.7 0.4 = 7 4 ≈ 0.571 .
Alternatively: conditional on X > 3 X > 3 X > 3 The distribution is U ( 3 , 10 ) \mathrm{U}(3, 10) U ( 3 , 10 ) So P ( X > 6 ∣ X > 3 ) = 10 − 6 10 − 3 = 4 7 P(X > 6 \mid X > 3) = \dfrac{10-6}{10-3} = \dfrac{4}{7} P ( X > 6 ∣ X > 3 ) = 10 − 3 10 − 6 = 7 4 .
Problem. X X X and Y Y Y are independent with X ∼ E x p ( λ 1 ) X \sim \mathrm{Exp}(\lambda_1) X ∼ Exp ( λ 1 ) and Y ∼ E x p ( λ 2 ) Y \sim \mathrm{Exp}(\lambda_2) Y ∼ Exp ( λ 2 ) . Find P ( X < Y ) P(X < Y) P ( X < Y ) .
Solution. Using the joint density and integration:
P ( X < Y ) = ∫ 0 ∞ ∫ x ∞ λ 1 e − λ 1 x ⋅ λ 2 e − λ 2 y d y d x P(X < Y) = \int_0^{\infty} \int_x^{\infty} \lambda_1 e^{-\lambda_1 x} \cdot \lambda_2 e^{-\lambda_2 y}\,dy\,dx P ( X < Y ) = ∫ 0 ∞ ∫ x ∞ λ 1 e − λ 1 x ⋅ λ 2 e − λ 2 y d y d x
= ∫ 0 ∞ λ 1 e − λ 1 x ⋅ e − λ 2 x d x = λ 1 ∫ 0 ∞ e − ( λ 1 + λ 2 ) x d x = \int_0^{\infty} \lambda_1 e^{-\lambda_1 x} \cdot e^{-\lambda_2 x}\,dx = \lambda_1 \int_0^{\infty} e^{-(\lambda_1+\lambda_2)x}\,dx = ∫ 0 ∞ λ 1 e − λ 1 x ⋅ e − λ 2 x d x = λ 1 ∫ 0 ∞ e − ( λ 1 + λ 2 ) x d x
= λ 1 λ 1 + λ 2 = \frac{\lambda_1}{\lambda_1 + \lambda_2} = λ 1 + λ 2 λ 1
For example, if λ 1 = λ 2 \lambda_1 = \lambda_2 λ 1 = λ 2 : P ( X < Y ) = 1 2 P(X < Y) = \dfrac{1}{2} P ( X < Y ) = 2 1 (by symmetry).
Problem. A continuous random variable X X X has PDF f ( x ) = k x ( 4 − x ) f(x) = kx(4-x) f ( x ) = k x ( 4 − x ) for 0 ≤ x ≤ 4 0 \leq x \leq 4 0 ≤ x ≤ 4 . Find k k k The CDF, and P ( 1 < X < 3 ) P(1 < X < 3) P ( 1 < X < 3 ) .
Solution. ∫ 0 4 k x ( 4 − x ) d x = k [ 2 x 2 − x 3 3 ] 0 4 = k ( 32 − 64 3 ) = 32 k 3 = 1 ⟹ k = 3 32 \displaystyle\int_0^4 kx(4-x)\,dx = k\!\left[2x^2 - \frac{x^3}{3}\right]_0^4 = k\!\left(32 - \frac{64}{3}\right) = \frac{32k}{3} = 1 \implies k = \frac{3}{32} ∫ 0 4 k x ( 4 − x ) d x = k [ 2 x 2 − 3 x 3 ] 0 4 = k ( 32 − 3 64 ) = 3 32 k = 1 ⟹ k = 32 3 .
CDF: F ( x ) = 3 32 ( 2 x 2 − x 3 3 ) = 3 x 2 16 − x 3 32 F(x) = \dfrac{3}{32}\!\left(2x^2 - \dfrac{x^3}{3}\right) = \dfrac{3x^2}{16} - \dfrac{x^3}{32} F ( x ) = 32 3 ( 2 x 2 − 3 x 3 ) = 16 3 x 2 − 32 x 3 For 0 ≤ x ≤ 4 0 \leq x \leq 4 0 ≤ x ≤ 4 .
P ( 1 < X < 3 ) = F ( 3 ) − F ( 1 ) = ( 27 16 − 27 32 ) − ( 3 16 − 1 32 ) = 27 32 − 5 32 = 11 16 P(1 < X < 3) = F(3) - F(1) = \left(\dfrac{27}{16} - \dfrac{27}{32}\right) - \left(\dfrac{3}{16} - \dfrac{1}{32}\right) = \dfrac{27}{32} - \dfrac{5}{32} = \boxed{\dfrac{11}{16}} P ( 1 < X < 3 ) = F ( 3 ) − F ( 1 ) = ( 16 27 − 32 27 ) − ( 16 3 − 32 1 ) = 32 27 − 32 5 = 16 11 .
Problem. X 1 , X 2 , … , X 50 X_1, X_2, \ldots, X_{50} X 1 , X 2 , … , X 50 are independent, each following E x p ( 0.1 ) \mathrm{Exp}(0.1) Exp ( 0.1 ) . Approximate P ( X ‾ > 12 ) P(\overline{X} > 12) P ( X > 12 ) .
Solution. E ( X i ) = 10 E(X_i) = 10 E ( X i ) = 10 , V a r ( X i ) = 100 \mathrm{Var}(X_i) = 100 Var ( X i ) = 100 . E ( X ‾ ) = 10 E(\overline{X}) = 10 E ( X ) = 10 V a r ( X ‾ ) = 100 50 = 2 \mathrm{Var}(\overline{X}) = \dfrac{100}{50} = 2 Var ( X ) = 50 100 = 2 .
By the CLT, X ‾ ≈ N ( 10 , 2 ) \overline{X} \approx N(10, 2) X ≈ N ( 10 , 2 ) approximately.
P ( X ‾ > 12 ) = P ( Z > 12 − 10 2 ) = P ( Z > 1.414 ) = 1 − 0.9214 = 0.0786 P(\overline{X} > 12) = P\!\left(Z > \frac{12-10}{\sqrt{2}}\right) = P(Z > 1.414) = 1 - 0.9214 = \boxed{0.0786} P ( X > 12 ) = P ( Z > 2 12 − 10 ) = P ( Z > 1.414 ) = 1 − 0.9214 = 0.0786
Problem. X X X has PDF f X ( x ) = 2 x f_X(x) = 2x f X ( x ) = 2 x for 0 < x < 1 0 < x < 1 0 < x < 1 . Find the PDF of Y = X 2 Y = X^2 Y = X 2 .
Solution. For 0 < y < 1 0 < y < 1 0 < y < 1 : F Y ( y ) = P ( Y ≤ y ) = P ( X 2 ≤ y ) = P ( X ≤ y ) = ( y ) 2 = y F_Y(y) = P(Y \leq y) = P(X^2 \leq y) = P(X \leq \sqrt{y}) = (\sqrt{y})^2 = y F Y ( y ) = P ( Y ≤ y ) = P ( X 2 ≤ y ) = P ( X ≤ y ) = ( y ) 2 = y .
f Y ( y ) = d d y F Y ( y ) = 1 for 0 < y < 1 f_Y(y) = \frac{d}{dy}F_Y(y) = \boxed{1} \quad \text{for } 0 < y < 1 f Y ( y ) = d y d F Y ( y ) = 1 for 0 < y < 1
So Y ∼ U ( 0 , 1 ) Y \sim \mathrm{U}(0,1) Y ∼ U ( 0 , 1 ) .
Problem. X X X has the triangular distribution with PDF f ( x ) = { 2 x 0 ≤ x ≤ 1 2 ( 2 − x ) 1 < x ≤ 2 f(x) = \begin{cases}2x & 0 \leq x \leq 1 \\ 2(2-x) & 1 < x \leq 2\end{cases} f ( x ) = { 2 x 2 ( 2 − x ) 0 ≤ x ≤ 1 1 < x ≤ 2 . Find the mode, Median, mean, and variance.
Solution. Mode: The PDF peaks at x = 1 x = 1 x = 1 So mode = 1 = \boxed{1} = 1 .
Median: For m ≤ 1 m \leq 1 m ≤ 1 : ∫ 0 m 2 x d x = m 2 \displaystyle\int_0^m 2x\,dx = m^2 ∫ 0 m 2 x d x = m 2 . Set m 2 = 0.5 ⟹ m = 1 2 ≈ 0.707 m^2 = 0.5 \implies m = \dfrac{1}{\sqrt{2}} \approx 0.707 m 2 = 0.5 ⟹ m = 2 1 ≈ 0.707 .
Mean: E ( X ) = ∫ 0 1 2 x 2 d x + ∫ 1 2 2 x ( 2 − x ) d x = 2 3 + [ 2 x 2 − 2 x 3 3 ] 1 2 = 2 3 + 4 3 = 2 E(X) = \displaystyle\int_0^1 2x^2\,dx + \int_1^2 2x(2-x)\,dx = \dfrac{2}{3} + \left[2x^2 - \dfrac{2x^3}{3}\right]_1^2 = \dfrac{2}{3} + \dfrac{4}{3} = \boxed{2} E ( X ) = ∫ 0 1 2 x 2 d x + ∫ 1 2 2 x ( 2 − x ) d x = 3 2 + [ 2 x 2 − 3 2 x 3 ] 1 2 = 3 2 + 3 4 = 2 .
Wait, let me recalculate: ∫ 1 2 2 x ( 2 − x ) d x = ∫ 1 2 ( 4 x − 2 x 2 ) d x = [ 2 x 2 − 2 x 3 3 ] 1 2 = ( 8 − 16 3 ) − ( 2 − 2 3 ) = 8 3 − 4 3 = 4 3 \displaystyle\int_1^2 2x(2-x)\,dx = \int_1^2 (4x - 2x^2)\,dx = \left[2x^2 - \dfrac{2x^3}{3}\right]_1^2 = (8-\dfrac{16}{3}) - (2-\dfrac{2}{3}) = \dfrac{8}{3} - \dfrac{4}{3} = \dfrac{4}{3} ∫ 1 2 2 x ( 2 − x ) d x = ∫ 1 2 ( 4 x − 2 x 2 ) d x = [ 2 x 2 − 3 2 x 3 ] 1 2 = ( 8 − 3 16 ) − ( 2 − 3 2 ) = 3 8 − 3 4 = 3 4 .
E ( X ) = 2 3 + 4 3 = 2 E(X) = \dfrac{2}{3} + \dfrac{4}{3} = \boxed{2} E ( X ) = 3 2 + 3 4 = 2 . This is the midpoint of [ 0 , 2 ] [0,2] [ 0 , 2 ] As expected for a Symmetric triangular distribution.
E ( X 2 ) = ∫ 0 1 2 x 3 d x + ∫ 1 2 ( 4 x 2 − 2 x 3 ) d x = 1 2 + [ 4 x 3 3 − x 4 2 ] 1 2 = 1 2 + 32 3 − 8 − 4 3 + 1 2 = 8 3 E(X^2) = \displaystyle\int_0^1 2x^3\,dx + \int_1^2 (4x^2 - 2x^3)\,dx = \dfrac{1}{2} + \left[\dfrac{4x^3}{3} - \dfrac{x^4}{2}\right]_1^2 = \dfrac{1}{2} + \dfrac{32}{3} - 8 - \dfrac{4}{3} + \dfrac{1}{2} = \dfrac{8}{3} E ( X 2 ) = ∫ 0 1 2 x 3 d x + ∫ 1 2 ( 4 x 2 − 2 x 3 ) d x = 2 1 + [ 3 4 x 3 − 2 x 4 ] 1 2 = 2 1 + 3 32 − 8 − 3 4 + 2 1 = 3 8 .
V a r ( X ) = 8 3 − 4 = − 4 3 \mathrm{Var}(X) = \dfrac{8}{3} - 4 = \boxed{-\dfrac{4}{3}} Var ( X ) = 3 8 − 4 = − 3 4 ? This is impossible. Let me recheck E ( X 2 ) E(X^2) E ( X 2 ) .
E ( X 2 ) = ∫ 0 1 2 x 3 d x + ∫ 1 2 x 2 ⋅ 2 ( 2 − x ) d x = 1 2 + 28 3 − 4 3 ⋅ . . . E(X^2) = \displaystyle\int_0^1 2x^3\,dx + \int_1^2 x^2 \cdot 2(2-x)\,dx = \dfrac{1}{2} + \dfrac{28}{3} - \dfrac{4}{3} \cdot ... E ( X 2 ) = ∫ 0 1 2 x 3 d x + ∫ 1 2 x 2 ⋅ 2 ( 2 − x ) d x = 2 1 + 3 28 − 3 4 ⋅ ...
Actually E ( X 2 ) = 1 2 + ∫ 1 2 ( 4 x 2 − 2 x 3 ) d x = 1 2 + [ 4 3 x 3 − x 4 2 ] 1 2 = 1 2 + ( 32 3 − 8 ) − ( 4 3 − 1 2 ) = 1 2 + 8 3 − 4 3 + 1 2 = 1 + 4 3 = 7 3 E(X^2) = \dfrac{1}{2} + \int_1^2 (4x^2 - 2x^3)\,dx = \dfrac{1}{2} + \left[\dfrac{4}{3}x^3 - \dfrac{x^4}{2}\right]_1^2 = \dfrac{1}{2} + (\dfrac{32}{3}-8) - (\dfrac{4}{3}-\dfrac{1}{2}) = \dfrac{1}{2}+\dfrac{8}{3}-\dfrac{4}{3}+\dfrac{1}{2} = 1+\dfrac{4}{3} = \dfrac{7}{3} E ( X 2 ) = 2 1 + ∫ 1 2 ( 4 x 2 − 2 x 3 ) d x = 2 1 + [ 3 4 x 3 − 2 x 4 ] 1 2 = 2 1 + ( 3 32 − 8 ) − ( 3 4 − 2 1 ) = 2 1 + 3 8 − 3 4 + 2 1 = 1 + 3 4 = 3 7 .
V a r ( X ) = 7 3 − 4 = − 5 3 \mathrm{Var}(X) = \dfrac{7}{3} - 4 = -\dfrac{5}{3} Var ( X ) = 3 7 − 4 = − 3 5 . This still cannot be right. The issue is E ( X ) = 1 E(X) = 1 E ( X ) = 1 (not 2) since the distribution is on [ 0 , 2 ] [0,2] [ 0 , 2 ] with peak at 1.
Let me redo: E ( X ) = 2 3 + 4 3 = 2 E(X) = \dfrac{2}{3} + \dfrac{4}{3} = 2 E ( X ) = 3 2 + 3 4 = 2 . But the distribution is symmetric about x = 1 x=1 x = 1 So E ( X ) E(X) E ( X ) should be 1 1 1 .
Rechecking the second integral: ∫ 1 2 2 ( 2 − x ) x d x \int_1^2 2(2-x)x\,dx ∫ 1 2 2 ( 2 − x ) x d x . At x = 1 x=1 x = 1 : 2 ( 1 ) ( 1 ) = 2 2(1)(1) = 2 2 ( 1 ) ( 1 ) = 2 . At x = 2 x=2 x = 2 : 0 0 0 . This Integral should give 2 / 3 2/3 2/3 by symmetry.
∫ 1 2 ( 4 x − 2 x 2 ) d x = [ 2 x 2 − 2 x 3 3 ] 1 2 = ( 8 − 16 3 ) − ( 2 − 2 3 ) = 8 3 − 4 3 = 4 3 \int_1^2 (4x-2x^2)\,dx = [2x^2-\frac{2x^3}{3}]_1^2 = (8-\frac{16}{3})-(2-\frac{2}{3}) = \frac{8}{3}-\frac{4}{3} = \frac{4}{3} ∫ 1 2 ( 4 x − 2 x 2 ) d x = [ 2 x 2 − 3 2 x 3 ] 1 2 = ( 8 − 3 16 ) − ( 2 − 3 2 ) = 3 8 − 3 4 = 3 4 .
Total: 2 3 + 4 3 = 2 \frac{2}{3}+\frac{4}{3} = 2 3 2 + 3 4 = 2 . But the range is [ 0 , 2 ] [0,2] [ 0 , 2 ] and the function is symmetric about x = 1 x=1 x = 1 . The mean of a symmetric distribution on [ 0 , 2 ] [0,2] [ 0 , 2 ] about x = 1 x=1 x = 1 is 1 1 1 . There must be a Normalization error. Let me verify: ∫ 0 1 2 x d x = 1 \int_0^1 2x\,dx = 1 ∫ 0 1 2 x d x = 1 and ∫ 1 2 2 ( 2 − x ) d x = [ 4 x − x 2 ] 1 2 = ( 8 − 4 ) − ( 4 − 1 ) = 1 \int_1^2 2(2-x)\,dx = [4x-x^2]_1^2 = (8-4)-(4-1) = 1 ∫ 1 2 2 ( 2 − x ) d x = [ 4 x − x 2 ] 1 2 = ( 8 − 4 ) − ( 4 − 1 ) = 1 . Total area = 2 ≠ 1 = 2 \neq 1 = 2 = 1 .
The PDF should be f ( x ) = x f(x) = x f ( x ) = x for 0 ≤ x ≤ 1 0 \leq x \leq 1 0 ≤ x ≤ 1 and f ( x ) = 2 − x f(x) = 2-x f ( x ) = 2 − x for 1 < x ≤ 2 1 < x \leq 2 1 < x ≤ 2 . Then E ( X ) = ∫ 0 1 x 2 d x + ∫ 1 2 x ( 2 − x ) d x = 1 3 + 2 3 = 1 E(X) = \int_0^1 x^2\,dx + \int_1^2 x(2-x)\,dx = \frac{1}{3}+\frac{2}{3} = 1 E ( X ) = ∫ 0 1 x 2 d x + ∫ 1 2 x ( 2 − x ) d x = 3 1 + 3 2 = 1 . ✓
Pitfall Correct Approach Confusing the rate λ \lambda λ with the mean 1 λ \dfrac{1}{\lambda} λ 1 for exponential distributions E ( X ) = 1 λ E(X) = \dfrac{1}{\lambda} E ( X ) = λ 1 ; the rate parameter is λ \lambda λ Forgetting that the total area under a PDF must equal 1 Always verify: ∫ − ∞ ∞ f ( x ) d x = 1 \displaystyle\int_{-\infty}^{\infty} f(x)\,dx = 1 ∫ − ∞ ∞ f ( x ) d x = 1 Applying the exponential memoryless property to other distributions Only the exponential distribution has this property Using P ( a < X < b ) = f ( b ) − f ( a ) P(a < X < b) = f(b) - f(a) P ( a < X < b ) = f ( b ) − f ( a ) This is for CDFs, not PDFs. Use ∫ a b f ( x ) d x \displaystyle\int_a^b f(x)\,dx ∫ a b f ( x ) d x
Calls arrive at a call centre at a rate of 4 per hour. Find the probability that the time between Two consecutive calls exceeds 45 minutes.
Solution Time between calls T ∼ E x p ( 4 ) T \sim \mathrm{Exp}(4) T ∼ Exp ( 4 ) (rate = 4 = 4 = 4 per hour).
P ( T > 0.75 ) = e − 4 × 0.75 = e − 3 ≈ 0.0498 P(T > 0.75) = e^{-4 \times 0.75} = e^{-3} \approx \boxed{0.0498} P ( T > 0.75 ) = e − 4 × 0.75 = e − 3 ≈ 0.0498 .
X X X is a continuous random variable with PDF f ( x ) = 3 4 ( 2 x − x 2 ) f(x) = \dfrac{3}{4}(2x - x^2) f ( x ) = 4 3 ( 2 x − x 2 ) for 0 ≤ x ≤ 2 0 \leq x \leq 2 0 ≤ x ≤ 2 . Find E ( X ) E(X) E ( X ) , V a r ( X ) \mathrm{Var}(X) Var ( X ) And the median.
Solution E ( X ) = 3 4 ∫ 0 2 ( 2 x 2 − x 3 ) d x = 3 4 [ 2 x 3 3 − x 4 4 ] 0 2 = 3 4 ( 16 3 − 4 ) = 3 4 ⋅ 4 3 = 1 E(X) = \dfrac{3}{4}\displaystyle\int_0^2 (2x^2-x^3)\,dx = \dfrac{3}{4}\!\left[\dfrac{2x^3}{3}-\dfrac{x^4}{4}\right]_0^2 = \dfrac{3}{4}\!\left(\dfrac{16}{3}-4\right) = \dfrac{3}{4}\cdot\dfrac{4}{3} = 1 E ( X ) = 4 3 ∫ 0 2 ( 2 x 2 − x 3 ) d x = 4 3 [ 3 2 x 3 − 4 x 4 ] 0 2 = 4 3 ( 3 16 − 4 ) = 4 3 ⋅ 3 4 = 1 .
E ( X 2 ) = 3 4 ∫ 0 2 ( 2 x 3 − x 4 ) d x = 3 4 [ x 4 2 − x 5 5 ] 0 2 = 3 4 ( 8 − 32 5 ) = 3 4 ⋅ 8 5 = 6 5 E(X^2) = \dfrac{3}{4}\displaystyle\int_0^2 (2x^3-x^4)\,dx = \dfrac{3}{4}\!\left[\dfrac{x^4}{2}-\dfrac{x^5}{5}\right]_0^2 = \dfrac{3}{4}\!\left(8-\dfrac{32}{5}\right) = \dfrac{3}{4}\cdot\dfrac{8}{5} = \dfrac{6}{5} E ( X 2 ) = 4 3 ∫ 0 2 ( 2 x 3 − x 4 ) d x = 4 3 [ 2 x 4 − 5 x 5 ] 0 2 = 4 3 ( 8 − 5 32 ) = 4 3 ⋅ 5 8 = 5 6 .
V a r ( X ) = 6 5 − 1 = 1 5 \mathrm{Var}(X) = \dfrac{6}{5}-1 = \dfrac{1}{5} Var ( X ) = 5 6 − 1 = 5 1 .
Median m m m : 3 4 ( m 2 − m 3 3 ) = 1 2 \dfrac{3}{4}\!\left(m^2-\dfrac{m^3}{3}\right) = \dfrac{1}{2} 4 3 ( m 2 − 3 m 3 ) = 2 1 . By inspection or Numerical methods: m ≈ 0.908 m \approx 0.908 m ≈ 0.908 .
Prove that for X ∼ E x p ( λ ) X \sim \mathrm{Exp}(\lambda) X ∼ Exp ( λ ) The memoryless property holds: P ( X > s + t ∣ X > s ) = P ( X > t ) P(X > s+t \mid X > s) = P(X > t) P ( X > s + t ∣ X > s ) = P ( X > t ) .
Solution P ( X > s + t ∣ X > s ) = P ( X > s + t ) P ( X > s ) = e − λ ( s + t ) e − λ s = e − λ t = P ( X > t ) P(X > s+t \mid X > s) = \frac{P(X > s+t)}{P(X > s)} = \frac{e^{-\lambda(s+t)}}{e^{-\lambda s}} = e^{-\lambda t} = P(X > t) P ( X > s + t ∣ X > s ) = P ( X > s ) P ( X > s + t ) = e − λ s e − λ ( s + t ) = e − λ t = P ( X > t )
■ \blacksquare ■
The exponential distribution models inter-arrival times in a Poisson process. See Poisson and Geometric Distributions .
Finding CDFs, means, and variances of continuous random variables requires integration. See Further Calculus .
The Central Limit Theorem connects the exponential and uniform distributions to the normal Distribution. See Chi-Squared Tests .
Distribution PDF E ( X ) E(X) E ( X ) V a r ( X ) \mathrm{Var}(X) Var ( X ) E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) λ e − λ x \lambda e^{-\lambda x} λ e − λ x , x ≥ 0 x \geq 0 x ≥ 0 1 λ \dfrac{1}{\lambda} λ 1 1 λ 2 \dfrac{1}{\lambda^2} λ 2 1 U ( a , b ) \mathrm{U}(a,b) U ( a , b ) 1 b − a \dfrac{1}{b-a} b − a 1 , a ≤ x ≤ b a \leq x \leq b a ≤ x ≤ b a + b 2 \dfrac{a+b}{2} 2 a + b ( b − a ) 2 12 \dfrac{(b-a)^2}{12} 12 ( b − a ) 2
Property Exponential Uniform Memoryless Yes No CDF 1 − e − λ x 1 - e^{-\lambda x} 1 − e − λ x x − a b − a \dfrac{x-a}{b-a} b − a x − a Median ln 2 λ \dfrac{\ln 2}{\lambda} λ ln 2 a + b 2 \dfrac{a+b}{2} 2 a + b
The lifetime of a light bulb follows an exponential distribution with mean 500 hours. Find: (a) the Probability it lasts more than 600 hours; (b) the probability it lasts between 400 and 600 hours; (c) the median lifetime.
Solution λ = 1 500 = 0.002 \lambda = \dfrac{1}{500} = 0.002 λ = 500 1 = 0.002 .
(a) P ( X > 600 ) = e − 1.2 ≈ 0.301 P(X > 600) = e^{-1.2} \approx \boxed{0.301} P ( X > 600 ) = e − 1.2 ≈ 0.301 .
(b) P ( 400 < X < 600 ) = e − 0.8 − e − 1.2 ≈ 0.449 − 0.301 = 0.148 P(400 < X < 600) = e^{-0.8} - e^{-1.2} \approx 0.449 - 0.301 = \boxed{0.148} P ( 400 < X < 600 ) = e − 0.8 − e − 1.2 ≈ 0.449 − 0.301 = 0.148 .
(c) Median m m m : e − 0.002 m = 0.5 ⟹ m = ln 2 0.002 = 346.6 hours e^{-0.002m} = 0.5 \implies m = \dfrac{\ln 2}{0.002} = \boxed{346.6\,\text{hours}} e − 0.002 m = 0.5 ⟹ m = 0.002 ln 2 = 346.6 hours .
Prove that for X ∼ U ( a , b ) X \sim \mathrm{U}(a,b) X ∼ U ( a , b ) , V a r ( X ) = ( b − a ) 2 12 \mathrm{Var}(X) = \dfrac{(b-a)^2}{12} Var ( X ) = 12 ( b − a ) 2 .
Solution E ( X ) = a + b 2 E(X) = \dfrac{a+b}{2} E ( X ) = 2 a + b .
E ( X 2 ) = 1 b − a ∫ a b x 2 d x = b 3 − a 3 3 ( b − a ) = a 2 + a b + b 2 3 E(X^2) = \dfrac{1}{b-a}\displaystyle\int_a^b x^2\,dx = \dfrac{b^3-a^3}{3(b-a)} = \dfrac{a^2+ab+b^2}{3} E ( X 2 ) = b − a 1 ∫ a b x 2 d x = 3 ( b − a ) b 3 − a 3 = 3 a 2 + ab + b 2 .
V a r ( X ) = a 2 + a b + b 2 3 − ( a + b ) 2 4 = 4 a 2 + 4 a b + 4 b 2 − 3 a 2 − 6 a b − 3 b 2 12 = a 2 − 2 a b + b 2 12 = ( b − a ) 2 12 \mathrm{Var}(X) = \dfrac{a^2+ab+b^2}{3} - \dfrac{(a+b)^2}{4} = \dfrac{4a^2+4ab+4b^2-3a^2-6ab-3b^2}{12} = \dfrac{a^2-2ab+b^2}{12} = \boxed{\dfrac{(b-a)^2}{12}} Var ( X ) = 3 a 2 + ab + b 2 − 4 ( a + b ) 2 = 12 4 a 2 + 4 ab + 4 b 2 − 3 a 2 − 6 ab − 3 b 2 = 12 a 2 − 2 ab + b 2 = 12 ( b − a ) 2 . ■ \blacksquare ■
The MGF of a random variable X X X is M X ( t ) = E ( e t X ) M_X(t) = E(e^{tX}) M X ( t ) = E ( e tX ) .
Properties:
M X ( 0 ) = 1 M_X(0) = 1 M X ( 0 ) = 1 M X ′ ( 0 ) = E ( X ) M_X'(0) = E(X) M X ′ ( 0 ) = E ( X ) M X ′ ′ ( 0 ) = E ( X 2 ) M_X''(0) = E(X^2) M X ′′ ( 0 ) = E ( X 2 ) 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)M_Y(t) M X + Y ( t ) = M X ( t ) M Y ( t ) MGFs:
E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) : M ( t ) = λ λ − t M(t) = \dfrac{\lambda}{\lambda-t} M ( t ) = λ − t λ for t < λ t < \lambda t < λ U ( a , b ) \mathrm{U}(a,b) U ( a , b ) : M ( t ) = e b t − e a t ( b − a ) t M(t) = \dfrac{e^{bt}-e^{at}}{(b-a)t} M ( t ) = ( b − a ) t e b t − e a t For any continuous random variable with PDF f ( x ) f(x) f ( x ) :
P ( a < X ≤ b ) = F ( b ) − F ( a ) P(a < X \leq b) = F(b) - F(a) P ( a < X ≤ b ) = F ( b ) − F ( a ) where F ( x ) = ∫ − ∞ x f ( t ) d t F(x) = \displaystyle\int_{-\infty}^x f(t)\,dt F ( x ) = ∫ − ∞ x f ( t ) d t .
For i.i.d. Random variables X 1 , … , X n X_1, \ldots, X_n X 1 , … , X n The order statistics are X ( 1 ) ≤ X ( 2 ) ≤ ⋯ ≤ X ( n ) X_{(1)} \leq X_{(2)} \leq \cdots \leq X_{(n)} X ( 1 ) ≤ X ( 2 ) ≤ ⋯ ≤ X ( n ) .
For X ∼ U ( 0 , 1 ) X \sim \mathrm{U}(0,1) X ∼ U ( 0 , 1 ) : X ( k ) ∼ B e t a ( k , n − k + 1 ) X_{(k)} \sim \mathrm{Beta}(k, n-k+1) X ( k ) ∼ Beta ( k , n − k + 1 ) .
Find the MGF of X ∼ E x p ( λ ) X \sim \mathrm{Exp}(\lambda) X ∼ Exp ( λ ) and use it to find E ( X ) E(X) E ( X ) and V a r ( X ) \mathrm{Var}(X) Var ( X ) .
Solution M ( t ) = ∫ 0 ∞ e t x λ e − λ x d x = λ ∫ 0 ∞ e − ( λ − t ) x d x = λ λ − t M(t) = \displaystyle\int_0^{\infty} e^{tx}\lambda e^{-\lambda x}\,dx = \lambda\displaystyle\int_0^{\infty} e^{-(\lambda-t)x}\,dx = \frac{\lambda}{\lambda-t} M ( t ) = ∫ 0 ∞ e t x λ e − λ x d x = λ ∫ 0 ∞ e − ( λ − t ) x d x = λ − t λ For t < λ t < \lambda t < λ .
M ′ ( t ) = λ ( λ − t ) 2 M'(t) = \dfrac{\lambda}{(\lambda-t)^2} M ′ ( t ) = ( λ − t ) 2 λ . M ′ ( 0 ) = 1 λ = E ( X ) M'(0) = \dfrac{1}{\lambda} = E(X) M ′ ( 0 ) = λ 1 = E ( X ) . ✓
M ′ ′ ( t ) = 2 λ ( λ − t ) 3 M''(t) = \dfrac{2\lambda}{(\lambda-t)^3} M ′′ ( t ) = ( λ − t ) 3 2 λ . M ′ ′ ( 0 ) = 2 λ 2 = E ( X 2 ) M''(0) = \dfrac{2}{\lambda^2} = E(X^2) M ′′ ( 0 ) = λ 2 2 = E ( X 2 ) .
V a r ( X ) = 2 λ 2 − 1 λ 2 = 1 λ 2 \mathrm{Var}(X) = \dfrac{2}{\lambda^2} - \dfrac{1}{\lambda^2} = \boxed{\dfrac{1}{\lambda^2}} Var ( X ) = λ 2 2 − λ 2 1 = λ 2 1 . ✓
Prove that if X ∼ U ( 0 , 1 ) X \sim \mathrm{U}(0,1) X ∼ U ( 0 , 1 ) Then Y = − 1 λ ln X Y = -\dfrac{1}{\lambda}\ln X Y = − λ 1 ln X follows E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) .
Solution F Y ( y ) = P ( Y ≤ y ) = P ( − 1 λ ln X ≤ y ) = P ( ln X ≥ − λ y ) = P ( X ≥ e − λ y ) F_Y(y) = P(Y \leq y) = P\!\left(-\dfrac{1}{\lambda}\ln X \leq y\right) = P(\ln X \geq -\lambda y) = P(X \geq e^{-\lambda y}) F Y ( y ) = P ( Y ≤ y ) = P ( − λ 1 ln X ≤ y ) = P ( ln X ≥ − λ y ) = P ( X ≥ e − λ y ) .
= 1 − e − λ y = 1 - e^{-\lambda y} = 1 − e − λ y for y ≥ 0 y \geq 0 y ≥ 0 .
This is the CDF of E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) . ■ \blacksquare ■