All revision notes topics

Probability generating functionsEdexcel A-Level Further Maths: Revision notes

Section 1

Definition and basic properties

For a discrete random variable XX taking non-negative integer values, the probability generating function (PGF) is GX(t)=E(tX)=∑rP(X=r) tr.\mathrm{G}_X(t)=\mathrm{E}(t^X)=\sum_{r}\mathrm{P}(X=r)\,t^r. The coefficient of trt^r is P(X=r)\mathrm{P}(X=r), so you can read off probabilities from a PGF written as a polynomial or power series. Since probabilities sum to 1, GX(1)=1\mathrm{G}_X(1)=1 always, which is a useful check. Example: if P(X=0)=0.2\mathrm{P}(X=0)=0.2, P(X=1)=0.5\mathrm{P}(X=1)=0.5 and P(X=2)=0.3\mathrm{P}(X=2)=0.3, then GX(t)=0.2+0.5t+0.3t2\mathrm{G}_X(t)=0.2+0.5t+0.3t^2.

Key termsprobability generating function
Exam tip

Check GX(1)=1\mathrm{G}_X(1)=1 every time. If it does not, an error has crept in.

Section 2

Mean and variance from the PGF

Differentiate term by term: GX′(t)=∑r P(X=r)tr−1\mathrm{G}_X'(t)=\sum r\,\mathrm{P}(X=r)t^{r-1}, so at t=1t=1 E(X)=GX′(1).\mathrm{E}(X)=\mathrm{G}_X'(1). Differentiating again gives GX′′(1)=E[X(X−1)]=E(X2)−E(X)\mathrm{G}_X''(1)=\mathrm{E}[X(X-1)]=\mathrm{E}(X^2)-\mathrm{E}(X). Hence Var(X)=GX′′(1)+GX′(1)−[GX′(1)]2.\mathrm{Var}(X)=\mathrm{G}_X''(1)+\mathrm{G}_X'(1)-\left[\mathrm{G}_X'(1)\right]^2. Example: GX(t)=0.2+0.5t+0.3t2\mathrm{G}_X(t)=0.2+0.5t+0.3t^2 gives GX′(1)=0.5+0.6=1.1\mathrm{G}_X'(1)=0.5+0.6=1.1, GX′′(1)=0.6\mathrm{G}_X''(1)=0.6 and Var(X)=0.6+1.1−1.21=0.49\mathrm{Var}(X)=0.6+1.1-1.21=0.49. Proofs of these standard results may be asked for in the exam.

Key termsmeanvariance
Common mistake

Forgetting to subtract [G′(1)]2[\mathrm{G}'(1)]^2, or using G′′(1)\mathrm{G}''(1) as the variance. G′′(1)=E(X2)−E(X)\mathrm{G}''(1)=\mathrm{E}(X^2)-\mathrm{E}(X) only.

Section 3

PGFs of standard distributions

Binomial B(n,p)\mathrm{B}(n,p): G(t)=(q+pt)n\mathrm{G}(t)=(q+pt)^n, from the binomial theorem. Poisson Po(λ)\mathrm{Po}(\lambda): G(t)=eλ(t−1)\mathrm{G}(t)=\mathrm{e}^{\lambda(t-1)}, from ∑e−λ(λt)rr!\sum\mathrm{e}^{-\lambda}\frac{(\lambda t)^r}{r!}. Geometric (trials to first success), P(X=r)=qr−1p\mathrm{P}(X=r)=q^{r-1}p: G(t)=pt1−qt\mathrm{G}(t)=\dfrac{pt}{1-qt}, from a geometric series, valid for ∣qt∣<1|qt|<1. Then G′(t)=p(1−qt)2\mathrm{G}'(t)=\dfrac{p}{(1-qt)^2} and E(X)=1p\mathrm{E}(X)=\frac1p. Negative binomial (trials to the rrth success): G(t)=(pt1−qt)r\mathrm{G}(t)=\left(\dfrac{pt}{1-qt}\right)^r, because it is a sum of rr independent geometric variables. For each, you must be able to derive the PGF from the definition and use it for the mean and variance.

Key termsbinomialPoissongeometricnegative binomial
Common mistake

Using the geometric PGF for the number of failures. pt1−qt\frac{pt}{1-qt} counts trials including the success.

Section 4

Sums of independent random variables

If XX and YY are independent, then E(tX+Y)=E(tX)E(tY)\mathrm{E}(t^{X+Y})=\mathrm{E}(t^X)\mathrm{E}(t^Y), so GX+Y(t)=GX(t)×GY(t).\mathrm{G}_{X+Y}(t)=\mathrm{G}_X(t)\times\mathrm{G}_Y(t). Derivation of this result is not required. It extends to any number of independent variables: the sum of nn independent copies of XX has PGF [GX(t)]n[\mathrm{G}_X(t)]^n. Examples: B(3,0.4)\mathrm{B}(3,0.4) plus an independent Po(2)\mathrm{Po}(2) has PGF (0.6+0.4t)3e2(t−1)(0.6+0.4t)^3\mathrm{e}^{2(t-1)}; the sum of two independent Po(λ)\mathrm{Po}(\lambda) and Po(μ)\mathrm{Po}(\mu) has PGF e(λ+μ)(t−1)\mathrm{e}^{(\lambda+\mu)(t-1)}, so is Po(λ+μ)\mathrm{Po}(\lambda+\mu).

Key termsindependent
Common mistake

Adding PGFs instead of multiplying them for a sum of variables.

Section 5

Worked example: geometric mean and variance

X∼Geo(p)X\sim\mathrm{Geo}(p) with p=14p=\frac14: G(t)=t4−3t\mathrm{G}(t)=\dfrac{t}{4-3t}. Then G′(t)=4(4−3t)2\mathrm{G}'(t)=\dfrac{4}{(4-3t)^2}, so E(X)=G′(1)=4\mathrm{E}(X)=\mathrm{G}'(1)=4. Next G′′(t)=24(4−3t)3\mathrm{G}''(t)=\dfrac{24}{(4-3t)^3}, so G′′(1)=24\mathrm{G}''(1)=24 and Var(X)=24+4−42=12\mathrm{Var}(X)=24+4-4^2=12, which agrees with qp2=3/41/16=12\frac{q}{p^2}=\frac{3/4}{1/16}=12. For the third success: N=X1+X2+X3N=X_1+X_2+X_3 has GN=[GX]3\mathrm{G}_N=[\mathrm{G}_X]^3, with E(N)=3×4=12\mathrm{E}(N)=3\times4=12 in this case.

Exam tip

Simplify G′(t)\mathrm{G}'(t) before differentiating a second time; it saves quotient-rule errors.

That's the notes covered.

Carry on to the next subtopic.

Exam questions on Probability generating functions

  1. The discrete random variable XX has P(X=0)=0.2\mathrm{P}(X=0)=0.2, P(X=1)=0.5\mathrm{P}(X=1)=0.5 and P(X=2)=0.3\mathrm{P}(X=2)=0.3. The probability generating function of XX is GX(t)\mathrm{G}_X(t).
    Use GX(t)\mathrm{G}_X(t) to find Var(X)\mathrm{Var}(X).2 marks
  2. The random variable XX has a geometric distribution with parameter p=14p=\frac14, so XX is the number of trials up to and including the first success, with P(X=r)=qr−1p\mathrm{P}(X=r)=q^{r-1}p for r=1,2,3,…r=1,2,3,\ldots and q=1−pq=1-p.
    Given that GX′(t)=4(4−3t)2\mathrm{G}_X'(t)=\dfrac{4}{(4-3t)^2}, find GX′′(1)\mathrm{G}_X''(1) and hence show that Var(X)=12\mathrm{Var}(X)=12.2 marks
  3. The random variable XX has distribution B(3,0.4)\mathrm{B}(3,0.4) and the random variable YY has distribution Po(2)\mathrm{Po}(2). The variables XX and YY are independent and Z=X+YZ=X+Y.
    Show, from the definition of a probability generating function, that GX(t)=(0.6+0.4t)3\mathrm{G}_X(t)=(0.6+0.4t)^3.3 marks
See the full worksheet

Written by the Exaim team, led by Shaun Daswani (Head of Upper Secondary, Improve ME Institute; MSc Financial Mathematics, Imperial College London; BSc, UCL) and Jason Daswani (operational lead, Improve ME Institute; LSE).