FS1 7.3 - Generating functions of sums

FS1 7.3 - Generating functions of sums

Independent observations combine by multiplying their PGFs. You will extract sum probabilities, recognise distributions after simplifying products, and handle shifted or weighted totals with the correct independence assumptions (FS1 7.3).

Two fresh observations are different from twice one observation

Is tossing two fair coins equivalent to tossing one fair coin and recording its result twice?

Let a head score 11 and a tail score 00. Both procedures have expected total 11, but two fresh tosses can score 11, whereas doubling one score gives only 00 or 22. Their different PGFs will expose the distinction immediately.

You already know

For a PGF, probabilities are coefficients and G(1)=1G(1)=1. From value transformations, Ga+bX(t)=taGX(tb)G_{a+bX}(t)=t^aG_X(t^b). From the moment lesson, derivatives at one give the mean and variance. Equal means do not imply identical distributions.

Independence supplies a new operation: multiplication. It is a condition on the observations, not a consequence of giving them different names.

PGF of an independent sum

If XX and YY are independent, then for S=X+YS=X+Y,

GS(t)=GX(t)GY(t)G_S(t)=G_X(t)G_Y(t)

Use the product where both expectations exist. No derivation of this rule is required.

In the formula booklet.

Repeat the rule for mutually independent X1,…,XnX_1,\ldots,X_n: the sum's PGF is ∏i=1nGXi(t)\prod_{i=1}^nG_{X_i}(t). If they also have the same distribution, it is [GX(t)]n[G_X(t)]^n. These repeated-product forms are consequences of the supplied two-variable rule, rather than separate printed formulae.

For the opening fair scores, GX(t)=(1+t)/2G_X(t)=(1+t)/2. Independence gives

GX1+X2(t)=(1+t2)2=14+12t+14t2G_{X_1+X_2}(t)=\left(\frac{1+t}{2}\right)^2 =\frac14+\frac12t+\frac14t^2

By contrast, doubling one score gives

G2X(t)=GX(t2)=12+12t2G_{2X}(t)=G_X(t^2)=\frac12+\frac12t^2

The coefficient at tt is 1/21/2 in the first expression and zero in the second. The same average hides a completely different possibility of an odd score.

Key idea

Independent copies produce powers of the PGF. Multiplying one random value changes the PGF's argument. In general, [GX(t)]2≠GX(t2)[G_X(t)]^2\ne G_X(t^2).

The following check asks about one observation, so no independence assumption can turn it into two fresh observations.

Multiplication collects every way to make a total

The coefficient of a total's power collects all pairs of exponents that add to it.

In the product GX(t)GY(t)G_X(t)G_Y(t), a term with power xx times a term with power yy contributes to power x+yx+y. Under independence its coefficient is P(X=x)P(Y=y)\mathrm P(X=x)\mathrm P(Y=y), the probability of that pair. Add all contributions to the required total.

[DIAGRAM: asset_name: Coefficients collect pairs with the same total; asset_slug: 9fm0-fs115-sum-coefficients; caption: Three independent outcome pairs contribute to the coefficient for total two; description: A three-by-three grid has X and Y values zero, one and two along its headers. Each cell shows its pair and the term t to the power of the pair sum divided by nine. The three cells for pairs zero comma two, one comma one and two comma zero are sage; the other cells are pale blue. A side panel gives probability of total two as three ninths, and the expanded sum PGF below is one plus two t plus three t squared plus two t cubed plus t to the fourth, all divided by nine.]
Three independent outcome pairs contribute to the coefficient for total two

Worked example 1: adding two equal discrete distributions

Question. Independent variables X1,X2X_1,X_2 each take 0,1,20,1,2 with equal probabilities. Find the PGF of S=X1+X2S=X_1+X_2, and hence find P(S=2)\mathrm P(S=2).

Plan. Square the individual PGF and read its quadratic coefficient.

GX(t)=1+t+t23individual distributionGS(t)=(1+t+t23)2independent copies=1+2t+3t2+2t3+t49collect matching powersP(S=2)=39=13coefficient of t2\begin{aligned} G_X(t)&=\frac{1+t+t^2}{3} &&\text{individual distribution}\\ G_S(t)&=\left(\frac{1+t+t^2}{3}\right)^2 &&\text{independent copies}\\ &=\frac{1+2t+3t^2+2t^3+t^4}{9} &&\text{collect matching powers}\\ \mathrm P(S=2)&=\frac39=\frac13 &&\text{coefficient of }t^2 \end{aligned}

Why this step? The quadratic term comes from (0,2)(0,2), (1,1)(1,1) and (2,0)(2,0). Reading only the middle pair would lose two valid possibilities.

Answer. GS(t)=(1+2t+3t2+2t3+t4)/9G_S(t)=(1+2t+3t^2+2t^3+t^4)/9, with P(S=2)=1/3\mathrm P(S=2)=1/3.

Check. There are nine equally likely pairs and three have total 22; the expanded coefficients add to 9/9=19/9=1.

Worked example 2 changes only the number of independent observations, from 22 to 33. The individual distribution and target total remain unchanged.

Worked example 2: adding a third observation

Question. For three independent variables with the same distribution, find the PGF of S=X1+X2+X3S=X_1+X_2+X_3, and hence P(S=2)\mathrm P(S=2).

Plan. Cube the PGF, counting only products contributing to power 22.

GS(t)=(1+t+t2)327three independent copiescoefficient of t2=3+(32)=6one quadratic or two linear choicesP(S=2)=627=29include the outside factor\begin{aligned} G_S(t)&=\frac{(1+t+t^2)^3}{27} &&\text{three independent copies}\\ \text{coefficient of }t^2&=3+\binom32=6 &&\text{one quadratic or two linear choices}\\ \mathrm P(S=2)&=\frac6{27}=\frac29 &&\text{include the outside factor} \end{aligned}

Why this step? You can obtain total 22 with one value 22, in three positions, or two values 11, in three pairs of positions. A full cubic expansion is unnecessary.

Answer. GS(t)=(1+t+t2)3/27G_S(t)=(1+t+t^2)^3/27, and P(S=2)=2/9\mathrm P(S=2)=2/9.

Check. The six permitted ordered triples each have probability 1/271/27. Also GS(1)=33/27=1G_S(1)=3^3/27=1.

For unequal distributions, multiply their different PGFs instead. The powers still collect totals; identical distributions are needed only to replace a product by a power of one common function.

Watch out

Wrong: add the PGFs to obtain a sum distribution, or multiply coefficients at the same power only. Right: multiply the functions and collect every exponent pair making the required total. Adding PGFs gives value 22, rather than 11, at t=1t=1.

Increase the number of copies once more.

Recognising the distribution after multiplying

Simplify the product, then match it to a supplied standard PGF, accounting for any shift in the variable.

The following are consequences of the independent-product rule and the supplied standard PGFs; they are not separately printed sum-distribution rules.

Independent variablesProduct PGFDistribution of the sum
Po⁡(λ),Po⁡(μ)\operatorname{Po}(\lambda),\operatorname{Po}(\mu)e(λ+μ)(t−1)\mathrm e^{(\lambda+\mu)(t-1)}Po⁡(λ+μ)\operatorname{Po}(\lambda+\mu)
B(m,p),B(n,p)\mathrm B(m,p),\mathrm B(n,p)(1−p+pt)m+n(1-p+pt)^{m+n}B(m+n,p)\mathrm B(m+n,p)
NB⁡(r,p),NB⁡(s,p)\operatorname{NB}(r,p),\operatorname{NB}(s,p)[pt/(1−(1−p)t)]r+s[pt/(1-(1-p)t)]^{r+s}NB⁡(r+s,p)\operatorname{NB}(r+s,p)

The binomial and negative binomial rows require the same pp. In particular, a sum of rr independent geometric trial counts with common pp is NB⁡(r,p)\operatorname{NB}(r,p). Different probabilities do not generally combine into one of these standard forms.

For example, independent Bernoulli probabilities 0.20.2 and 0.80.8 give product 0.16+0.68t+0.16t20.16+0.68t+0.16t^2, not the 0.25+0.5t+0.25t20.25+0.5t+0.25t^2 of B(2,0.5)\mathrm B(2,0.5). Their mean alone does not identify the sum distribution.

Worked example 3: a nonstandard pair produces a standard wait

Question. Independent X1,X2X_1,X_2 each have PGF (3−2t)−1/2(3-2t)^{-1/2}. Let Y=X1+X2+1Y=X_1+X_2+1. Find GY(t)G_Y(t), identify the distribution of YY, and find P(Y>5)\mathrm P(Y>5) exactly.

Plan. Multiply the two PGFs, include the shift factor, then match the geometric trial-count form.

GY(t)=t(3−2t)−1/2(3−2t)−1/2independent sum and shift=t3−2tcombine equal bases=(1/3)t1−(2/3)tstandard geometric formY∼Geo⁡(1/3)include its parameterP(Y>5)=(23)5=32243five initial failures\begin{aligned} G_Y(t)&=t(3-2t)^{-1/2}(3-2t)^{-1/2} &&\text{independent sum and shift}\\ &=\frac t{3-2t} &&\text{combine equal bases}\\ &=\frac{(1/3)t}{1-(2/3)t} &&\text{standard geometric form}\\ Y&\sim\operatorname{Geo}(1/3) &&\text{include its parameter}\\ \mathrm P(Y>5)&=\left(\frac23\right)^5=\frac{32}{243} &&\text{five initial failures} \end{aligned}

Why this step? Both operations are needed: multiplication combines the independent values, while the outside tt adds 11. Omitting either changes the variable.

Answer. GY(t)=t/(3−2t)G_Y(t)=t/(3-2t), Y∼Geo⁡(1/3)Y\sim\operatorname{Geo}(1/3), and P(Y>5)=32/243\mathrm P(Y>5)=32/243.

Check. The PGF's first non-zero power is tt, with coefficient 1/31/3; that matches a first-success trial count starting at 11.

Continue this lesson.

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

Create free account