Computing conditional extinction probabilities for pollen cells

  • Thread starter Thread starter psie
  • Start date Start date
  • Tags Tags
    Probability theory
Click For Summary
SUMMARY

This discussion focuses on computing conditional extinction probabilities for pollen cells using Galton-Watson processes. The extinction probability, denoted as ##\eta##, is derived from the equation ##t=g(t)##, where ##g(t)=q+pt^2##, leading to the conclusion that ##\eta=q/p<1##. The computed extinction probabilities for parts (a) and (b) are 7 and approximately 0.0206, respectively. The confusion arises regarding the interpretation of the number of ancestors and the application of the formulas under the binary splitting assumption.

PREREQUISITES
  • Understanding of Galton-Watson processes
  • Familiarity with extinction probabilities in branching processes
  • Knowledge of generating functions, specifically the function ##g(t)##
  • Basic probability theory, including conditional probabilities
NEXT STEPS
  • Study the properties of Galton-Watson processes in detail
  • Learn about generating functions and their applications in probability
  • Explore advanced topics in branching processes, including multi-type branching
  • Investigate the implications of binary splitting in population dynamics
USEFUL FOR

Researchers in mathematical biology, statisticians working with stochastic processes, and students studying probability theory and branching processes.

psie
Messages
315
Reaction score
40
Homework Statement
The growth dynamics of pollen cells can be modeled by binary splitting as follows: After one unit of time, a cell either splits into two or dies. The new cells develop according to the same law independently of each other. The probabilities of dying and splitting are 0.46 and 0.54, respectively.

(a) Determine the maximal initial size of the population in order for the probability of extinction to be at least 0.3.

(b) What is the probability that the population is extinct after two generations if the initial population is the maximal number obtained in (a)?
Relevant Equations
My previous post collects some useful assumptions about the Galton-Watson process: https://www.physicsforums.com/threads/understanding-extinction-probability-in-simple-gwp.1064663/
I'm studying branching processes, and I am bit confused about this exercise. The assumption of the processes I am studying (Galton-Watson processes) is that ##X(0)=1##, i.e. the population starts with one single ancestor. Then ##X(1)## denotes the number of children obtained by the ancestor. The answer I get for (a) is 7, and this is what the book gets. The way I computed this was using the fact that the extinction probability ##\eta## is the smallest nonnegative root to ##t=g(t)##, where ##g(t)=q+pt^2##, so ##\eta=q/p<1## (here ##q=1-p##). Now, from my previous post above, we've got $$P(\text{extinction}\mid A_k)=\eta^k,$$where ##A_k=\{\text{the founding member produces } k\text{ children}\}, k\geq 0##. So we want to solve for ##k## in $$\eta^k\geq0.3,$$which gives approximately 7.5, and we round off to 7, since the inequality does not hold if ##k=8##.

What confuses me about part (a) and how I solved it is that we are studying binary splitting, so the ancestor can only get 2 children. I am not sure if the number 7 relates to the ancestor's children or not. If 7 relates to the number of ancestors, then I'm not sure my formulas hold because they were derived under the assumption that ##X(0)=1##. Also, I am confused about the probability one wants to find in (b). I can compute ##P(X(n)=0)## for ##n\geq 1##, but I don't think this is the probability they are after. The answer here is 0.0206.
 
Physics news on Phys.org
Both parts are solved.

(a) Let there be ##n## cells to start with; each cell defines a Galton-Watson process and the extinction probability for each branch is ##\eta##. So the probability for the entire population to go extinct is ##\eta^n##.

(b) The probability we are looking for in (b) is that each branch goes extinct by generation ##2##. From my previous post, we have ##q_n=P(X(n)=0)## and ##q_n=\sum_{k=0}^\infty p_kq_{n-1}^k=p_0+p_2q_{n-1}^2##, where ##q_0=0##. So ##q_2\approx 0.5743##. For the entire population to go extinct by generation ##2##, this probability is just ##(0.5743)^7\approx 0.0206##.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
Replies
2
Views
1K
  • · Replies 8 ·
Replies
8
Views
1K
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 3 ·
Replies
3
Views
1K
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 5 ·
Replies
5
Views
1K