4.5 - PGFs of independent sums

4.5 - PGFs of independent sums

Adding independent counts corresponds to multiplying their probability generating functions. You will explain that connection, recognise binomial and Poisson sum distributions, and use coefficients and derivatives to find probabilities and moments of a total.

The product rule

If two counts each have a short probability distribution, can you find their total without listing every outcome separately? Suppose independent variables X,YX,Y each take values 0,10,1, with probabilities 2/3,1/32/3,1/3. A total of one can occur as (0,1)(0,1) or (1,0)(1,0), so its probability is 2/9+2/9=4/92/9+2/9=4/9.

Multiplication of their PGFs makes exactly the same calculation:

(23+13t)(23+13t)=49+49t+19t2.\left(\frac23+\frac13t\right) \left(\frac23+\frac13t\right) =\frac49+\frac49t+\frac19t^2.

An exponent in a product is the sum of the original exponents, so collecting equal powers collects all pairs giving the same total.

[DIAGRAM: asset_name: Collecting probabilities for a total; asset_slug: cam9231-pgf-sum-coefficients; caption: The two outcome pairs with total one contribute to the same coefficient; description: Exact probability grid for independent X and Y, each taking zero with probability two thirds and one with probability one third. The four cells show totals zero, one, one and two with probabilities four ninths, two ninths, two ninths and one ninth. A neighbouring bar chart shows the sum distribution four ninths, four ninths, one ninth on zero, one, two, using matching pale colours.]
The two outcome pairs with total one contribute to the same coefficient

For independent X,YX,Y, the general argument is

GX+Y(t)=E(tX+Y)=E(tXtY)=E(tX) E(tY)=GX(t)GY(t)\begin{aligned} G_{X+Y}(t) &=\mathrm E(t^{X+Y}) \\ &=\mathrm E(t^Xt^Y) \\ &=\mathrm E(t^X)\,\mathrm E(t^Y) \\ &=G_X(t)G_Y(t) \end{aligned}

Independence justifies splitting the expectation of a product. Having the same distribution does not imply independence.

PGF of an independent sum

For independent nonnegative integer-valued variables X1,…,XnX_1,\ldots,X_n, where the PGFs converge (in particular, 0≤t≤10\le t\le1),

GX1+⋯+Xn(t)=∏i=1nGXi(t).G_{X_1+\cdots+X_n}(t)=\prod_{i=1}^{n}G_{X_i}(t).

If the variables also have the same distribution, with PGF G(t)G(t), the result is {G(t)}n\{G(t)\}^n.

The product rule is not listed in the formula booklet: learn the result and its independence condition.

The variables need not have identical distributions. Each factor belongs to the corresponding independent count.

Worked example 1: two independent observations

Question. Independent variables X,YX,Y each satisfy P(X=0)=P(Y=0)=2/3\mathrm P(X=0)=\mathrm P(Y=0)=2/3 and P(X=1)=P(Y=1)=1/3\mathrm P(X=1)=\mathrm P(Y=1)=1/3. Find the PGF and distribution of S=X+YS=X+Y.

Plan. Multiply the PGFs and read the three coefficients.

GS(t)=(23+13t)2=49+49t+19t2\begin{aligned} G_S(t)&=\left(\frac23+\frac13t\right)^2 \\ &=\frac49+\frac49t+\frac19t^2 \end{aligned}

Answer. The PGF is GS(t)=(4+4t+t2)/9G_S(t)=(4+4t+t^2)/9. The probabilities of S=0,1,2S=0,1,2 are 4/9,4/9,1/94/9,4/9,1/9, respectively.

Check. Directly, both variables are zero with probability (2/3)2=4/9(2/3)^2=4/9, exactly one is one with probability 2(2/3)(1/3)=4/92(2/3)(1/3)=4/9, and both are one with probability (1/3)2=1/9(1/3)^2=1/9.

Change only the dependence: use the same observation twice instead of taking a second independent observation. The individual distributions stay the same, but the total changes.

Worked example 2: repeating the same observation

Question. XX takes values 0,10,1, with probabilities 2/3,1/32/3,1/3. Define Y=XY=X and T=X+YT=X+Y. Find GT(t)G_T(t), and explain why GX(t)GY(t)G_X(t)G_Y(t) is unsuitable.

Plan. Express TT as 2X2X, then encode its actual possible values.

T=2XGT(t)=23+13t2=GX(t2)\begin{aligned} T&=2X \\ G_T(t)&=\frac23+\frac13t^2 \\ &=G_X(t^2) \end{aligned}

There is no outcome T=1T=1: either both summands are zero or both are one.

Answer. GT(t)=(2+t2)/3G_T(t)=(2+t^2)/3. The product of the two individual PGFs would give a positive probability for T=1T=1, because it would incorrectly treat XX and Y=XY=X as independent.

Check. Directly, P(T=0)=P(X=0)=2/3\mathrm P(T=0)=\mathrm P(X=0)=2/3 and P(T=2)=P(X=1)=1/3\mathrm P(T=2)=\mathrm P(X=1)=1/3. These agree with the two coefficients and sum to one.

More generally, for nonnegative integers a,ba,b,

GaX+b(t)=E(taX+b)=tbGX(ta).G_{aX+b}(t)=\mathrm E(t^{aX+b})=t^bG_X(t^a).

This is a transformation of one observation. In particular, G2X(t)=GX(t2)G_{2X}(t)=G_X(t^2), whereas two independent copies of XX have sum PGF {GX(t)}2\{G_X(t)\}^2.

Watch out

Squaring the PGF describes a sum of two independent observations. Replacing tt by t2t^2 describes twice one observation.

Use multiplication for the following independent counts, whose success probabilities differ.

Recognising a sum distribution

A product of PGFs sometimes has a standard form. Since the coefficients determine the whole distribution, recognising that form identifies the sum without a long expansion.

For independent X∼B(m,p)X\sim\mathrm B(m,p) and Y∼B(n,p)Y\sim\mathrm B(n,p),

GX+Y(t)=(1−p+pt)m(1−p+pt)n=(1−p+pt)m+n\begin{aligned} G_{X+Y}(t)&=(1-p+pt)^m(1-p+pt)^n \\ &=(1-p+pt)^{m+n} \end{aligned}

Hence X+Y∼B(m+n,p)X+Y\sim\mathrm B(m+n,p). The common success probability is essential to this conclusion: different linear factors generally cannot be combined into a power of one common factor.

For independent U∼Po⁡(λ)U\sim\operatorname{Po}(\lambda) and V∼Po⁡(μ)V\sim\operatorname{Po}(\mu),

GU+V(t)=eλ(t−1)eμ(t−1)=e(λ+μ)(t−1)\begin{aligned} G_{U+V}(t) &=\mathrm e^{\lambda(t-1)}\mathrm e^{\mu(t-1)} \\ &=\mathrm e^{(\lambda+\mu)(t-1)} \end{aligned}

Thus U+V∼Po⁡(λ+μ)U+V\sim\operatorname{Po}(\lambda+\mu). This derives the Poisson sum result from the PGF product rule.

Key idea

Match the entire PGF to a standard form, including its normalisation. A resemblance in one factor is not enough to identify the distribution.

First recognise a sum within the binomial family.

Worked example 3: combining binomial counts

Question. X∼B(5,14)X\sim\mathrm B(5,\tfrac14) and Y∼B(7,14)Y\sim\mathrm B(7,\tfrac14) are independent. Use PGFs to identify the distribution of S=X+YS=X+Y.

Plan. Multiply the two binomial PGFs and combine equal bases.

GS(t)=(34+14t)5(34+14t)7=(34+14t)12\begin{aligned} G_S(t) &=\left(\frac34+\frac14t\right)^5 \left(\frac34+\frac14t\right)^7 \\ &=\left(\frac34+\frac14t\right)^{12} \end{aligned}

Answer. S∼B(12,14)S\sim\mathrm B(12,\tfrac14).

Check. The means add to 5/4+7/4=35/4+7/4=3, matching 12(1/4)12(1/4). Independence makes the variances add to 5(1/4)(3/4)+7(1/4)(3/4)=9/45(1/4)(3/4)+7(1/4)(3/4)=9/4, matching the identified binomial variance.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account