Tutorials › Stochastic Processes › Deriving the Mean Population Recursion

Galton-Watson Branching Processes · Tutorial 348 of 994

Deriving the Mean Population Recursion

See why the expected next generation is the offspring mean times the expected current population, and when this identity is valid.

Advanced 9 min read

What You'll Learn

  • Express the next generation as a random sum of offspring counts.
  • State the independence and finite-mean assumptions needed for conditioning.
  • Prove the conditional mean identity given the current generation.
  • Derive the unconditional population-mean recursion.
  • Distinguish an expectation identity from a pathwise population equation.
  • Account for a fixed or random initial population.

From Offspring Counts to a Mean Recursion

A Galton–Watson process describes a population in generations: each individual produces a random number of children, and different individuals reproduce independently according to the same offspring law. The population in the next generation is therefore a random sum. Conditioning on the current population fixes how many terms appear in that sum, which is the key to computing its mean.

Work on a probability space \((\Omega,\mathcal{F},\mathbb{P})\). Let \(Z_n\) be the number of individuals in generation \(n\), with state space \(S=\{0,1,2,\ldots\}\) and discrete time \(n\geq 0\). Let \(\mathcal{F}_n\) contain the information about the process through generation \(n\), before the individuals in generation \(n\) reproduce. In particular, \(Z_n\) is \(\mathcal{F}_n\)-measurable.

For \(r\geq1\), write \(\xi_{n,r}\) for the number of children produced by the \(r\)th individual in generation \(n\), when that individual exists. Assume these offspring variables are independent and identically distributed, are independent of \(\mathcal{F}_n\) at generation \(n\), and have a common finite mean

$$ m=\mathbb{E}[\xi_{n,r}]<\infty. $$

Here each \(\xi_{n,r}\) is a nonnegative integer-valued random variable. We assume \(Z_n<\infty\) almost surely, as is the case when the initial population is finite almost surely and every individual has a finite number of children. The offspring-sum representation is

$$ Z_{n+1}=\sum_{r=1}^{Z_n}\xi_{n,r}, $$

where a sum with no terms is defined to be \(0\). Thus, if \(Z_n=0\), then \(Z_{n+1}=0\): there are no individuals available to reproduce.

Proposition (Conditional mean of the next generation): Under the offspring assumptions above, for any \(n\) such that \(Z_n<\infty\) almost surely, \(\mathbb{E}[Z_{n+1}\mid\mathcal{F}_n]=mZ_n\) almost surely, where the conditional expectation is understood as a nonnegative extended conditional expectation if its expectation is not finite.

Proof. For each \(r\geq1\), define \(Y_r=\mathbf{1}_{\{Z_n\geq r\}}\xi_{n,r}\). The event \(\{Z_n\geq r\}\) is in \(\mathcal{F}_n\), and \(\xi_{n,r}\) is independent of \(\mathcal{F}_n\). Therefore \(\mathbb{E}[Y_r\mid\mathcal{F}_n]=\mathbf{1}_{\{Z_n\geq r\}}\mathbb{E}[\xi_{n,r}\mid\mathcal{F}_n]=m\mathbf{1}_{\{Z_n\geq r\}}\) almost surely. The offspring-sum representation gives \(Z_{n+1}=\sum_{r\geq1}Y_r\): for each outcome, the indicators are one exactly for the first \(Z_n\) terms. Since all \(Y_r\) are nonnegative, conditional monotone convergence gives \(\mathbb{E}[Z_{n+1}\mid\mathcal{F}_n]=\sum_{r\geq1}\mathbb{E}[Y_r\mid\mathcal{F}_n]=m\sum_{r\geq1}\mathbf{1}_{\{Z_n\geq r\}}=mZ_n\) almost surely. The last sum has exactly \(Z_n\) terms, including zero terms when \(Z_n=0\). \(\square\)

This conditional statement says that after the generation-\(n\) population is known, each of its \(Z_n\) members contributes an expected \(m\) children. It is an identity of conditional expectations, not an assertion that the realized population is \(mZ_n\). The offspring counts remain random even after \(Z_n\) is known.

Theorem (Mean population recursion): Suppose the offspring assumptions above hold and \(\mathbb{E}[Z_n]<\infty\). Then \(\mathbb{E}[Z_{n+1}]=m\mathbb{E}[Z_n]\). In particular, the theorem applies at \(n=0\) when the initial population is either a fixed finite number or a random, finite-almost-surely population with finite mean, independent of the future offspring variables.

Proof. By the conditional-mean proposition and the defining property of conditional expectation, \(\mathbb{E}[Z_{n+1}]=\mathbb{E}[\mathbb{E}[Z_{n+1}\mid\mathcal{F}_n]]=\mathbb{E}[mZ_n]=m\mathbb{E}[Z_n]\). The assumption \(\mathbb{E}[Z_n]<\infty\), together with \(m<\infty\), makes the right-hand side finite, so the unconditional expectation is finite as well. If \(Z_0\) is a fixed finite number, it has finite expectation; if \(Z_0\) is random with finite mean, the same calculation applies at generation zero. In either case the offspring variables must have the stated law independently of the information about the initial population. \(\square\)

Worked Examples

Worked Example: Three Individuals with Bernoulli Offspring

Suppose the current population is \(Z_n=3\), and each individual has either \(0\) or \(1\) child, with probability \(\tfrac13\) of \(0\) and \(\tfrac23\) of \(1\). Find the expected next population and verify it directly from the distribution of \(Z_{n+1}\).

The offspring mean is \(m=0\cdot\tfrac13+1\cdot\tfrac23=\tfrac23\). Since \(Z_n=3\), the conditional-mean identity gives \(\mathbb{E}[Z_{n+1}\mid\mathcal{F}_n]=3\cdot\tfrac23=2\).

For a direct check, \(Z_{n+1}\) is the sum of three independent Bernoulli variables with success probability \(\tfrac23\). Its probabilities are \(\mathbb{P}(Z_{n+1}=0)=\tfrac1{27}\), \(\mathbb{P}(Z_{n+1}=1)=\tfrac6{27}\), \(\mathbb{P}(Z_{n+1}=2)=\tfrac{12}{27}\), and \(\mathbb{P}(Z_{n+1}=3)=\tfrac8{27}\). These sum to \(\tfrac{1+6+12+8}{27}=1\). Therefore \(\mathbb{E}[Z_{n+1}]=0\cdot\tfrac1{27}+1\cdot\tfrac6{27}+2\cdot\tfrac{12}{27}+3\cdot\tfrac8{27}=\tfrac{6+24+24}{27}=\tfrac{54}{27}=2\), as required.

Worked Example: A Random Initial Population

Let \(Z_0\) equal \(0\), \(1\), or \(3\) with probabilities \(\tfrac14\), \(\tfrac12\), and \(\tfrac14\), respectively. Independently, let the offspring count \(\xi\) have probabilities \(\mathbb{P}(\xi=0)=\tfrac12\), \(\mathbb{P}(\xi=1)=\tfrac14\), and \(\mathbb{P}(\xi=3)=\tfrac14\). Find \(\mathbb{E}[Z_1]\), both by conditioning on \(Z_0\) and by using the recursion.

First compute the offspring mean: \(m=0\cdot\tfrac12+1\cdot\tfrac14+3\cdot\tfrac14=0+\tfrac14+\tfrac34=1\). The initial mean is \(\mathbb{E}[Z_0]=0\cdot\tfrac14+1\cdot\tfrac12+3\cdot\tfrac14=0+\tfrac12+\tfrac34=\tfrac54\). Thus the mean recursion gives \(\mathbb{E}[Z_1]=1\cdot\tfrac54=\tfrac54\).

For the conditional calculation, if \(Z_0=0\), then \(Z_1=0\). If \(Z_0=1\), the conditional mean of \(Z_1\) is \(1\cdot m=1\). If \(Z_0=3\), the conditional mean is \(3\cdot m=3\). Averaging over the three possible initial populations gives \(\mathbb{E}[Z_1]=0\cdot\tfrac14+1\cdot\tfrac12+3\cdot\tfrac14=\tfrac12+\tfrac34=\tfrac54\). This agrees with the recursion and shows explicitly how the initial-population distribution enters through its mean at this step.

Worked Example: The Mean Is Not the Realized Population

Suppose \(Z_n=2\) and each individual has \(0\) children with probability \(\tfrac12\) and \(2\) children with probability \(\tfrac12\), independently. Find the law and mean of \(Z_{n+1}\), and explain why the recursion is not a pathwise equality.

The offspring mean is \(m=0\cdot\tfrac12+2\cdot\tfrac12=1\), so the conditional-mean identity predicts \(\mathbb{E}[Z_{n+1}\mid\mathcal{F}_n]=1\cdot2=2\). The two offspring counts can be \((0,0)\), \((0,2)\), \((2,0)\), or \((2,2)\). By independence, each ordered pair has probability \(\tfrac12\cdot\tfrac12=\tfrac14\). Hence \(Z_{n+1}\) equals \(0\) with probability \(\tfrac14\), \(2\) with probability \(\tfrac12\), and \(4\) with probability \(\tfrac14\). The probabilities sum to \(\tfrac14+\tfrac12+\tfrac14=1\), and \(\mathbb{E}[Z_{n+1}]=0\cdot\tfrac14+2\cdot\tfrac12+4\cdot\tfrac14=1+1=2\).

In particular, the realized next population can be \(0\), \(2\), or \(4\), not always \(2\). The value \(2\) is its conditional mean given the current population, not a prediction that each reproduction outcome must equal the mean.

Assumptions, Initial Conditions, and a Common Error

The finite-mean condition has a precise role. To conclude a finite numerical identity from the theorem, we need both \(m<\infty\) and \(\mathbb{E}[Z_n]<\infty\). A finite population almost surely is not, by itself, enough to guarantee finite expectation: a random initial population can take finite values on every outcome and still have an infinite mean. In that case the ordinary finite-expectation calculation in the theorem is not available. The conditional identity remains a nonnegative conditional-expectation statement, but it should not be presented as an equality between finite numbers unless the integrability assumptions hold.

The initial population need not be deterministic. If \(Z_0=z\) for a fixed nonnegative integer \(z\), the first expected population is \(mz\). If \(Z_0\) is random with finite mean and is independent of the offspring variables, the first expected population is \(m\mathbb{E}[Z_0]\). The case \(z=0\) is included: the offspring sum is empty, so the next population is zero and the formula gives \(m\cdot0=0\). The case \(m=0\) is also included. Since offspring counts are nonnegative and have mean zero, they equal zero almost surely, and the next population is zero whenever the current population is finite.

Common mistake. A frequent incorrect calculation is to set \(\mathbb{E}[Z_{n+1}\mid\mathcal{F}_n]=m\), as if the next generation consisted of the offspring of just one individual. For example, with \(Z_n=3\) and offspring mean \(m=\tfrac23\), this gives the wrong answer \(\tfrac23\). The offspring-sum representation has three terms, so the correct calculation is \(\mathbb{E}[Z_{n+1}\mid\mathcal{F}_n]=\tfrac23+\tfrac23+\tfrac23=2\). Equivalently, it is \(mZ_n=\tfrac23\cdot3=2\). The factor \(Z_n\) counts how many individuals reproduce.

The conditional form is the most informative version of the result: it states what the next-generation mean is once the current population is known. Taking expectations then produces the unconditional recursion. Keeping those two steps separate helps prevent confusion between the random population \(Z_{n+1}\), its conditional mean, and its unconditional mean.

Check Your Understanding

Use the offspring-sum representation and check each numerical mean by substituting the stated parameters.

  1. If \(Z_n=4\) and the offspring mean is \(\tfrac32\), what is \(\mathbb{E}[Z_{n+1}\mid\mathcal{F}_n]\)?
  2. If \(Z_0\) is random with finite mean \(3\), and the offspring mean is \(\tfrac23\), what is \(\mathbb{E}[Z_1]\), assuming independence?
  3. Why does \(Z_n=0\) imply \(Z_{n+1}=0\) in the offspring-sum representation?
  4. For offspring probabilities \(\mathbb{P}(\xi=0)=\tfrac34\) and \(\mathbb{P}(\xi=4)=\tfrac14\), compute \(m\). If \(Z_n=2\), compute the conditional mean of \(Z_{n+1}\).
  5. Explain why \(Z_{n+1}=mZ_n\) is generally not a pathwise identity, even though the corresponding conditional-expectation identity holds.