Iowa State University Capstones, Theses and Graduate Theses and Dissertations Dissertations
2011 Coalescence in Bellman-Harris and multi-type branching processes Jyy-i Joy Hong Iowa State University
Follow this and additional works at: https://lib.dr.iastate.edu/etd Part of the Mathematics Commons
Recommended Citation Hong, Jyy-i Joy, "Coalescence in Bellman-Harris and multi-type branching processes" (2011). Graduate Theses and Dissertations. 10103. https://lib.dr.iastate.edu/etd/10103
This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Graduate Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Coalescence in Bellman-Harris and multi-type branching processes
by
Jyy-I Hong
A dissertation submitted to the graduate faculty
in partial fulfillment of the requirements for the degree of
DOCTOR OF PHILOSOPHY
Major: Mathematics
Program of Study Committee:
Krishna B. Athreya, Major Professor
Clifford Bergman
Dan Nordman
Ananda Weerasinghe
Paul E. Sacks
Iowa State University
Ames, Iowa
2011
Copyright c Jyy-I Hong, 2011. All rights reserved. ii
DEDICATION
I would like to dedicate this thesis to my parents Wan-Fu Hong and Wen-Hsiang Tseng for their un- conditional love and support. Without them, the completion of this work would not have been possible. iii
TABLE OF CONTENTS
ACKNOWLEDGEMENTS ...... vii
ABSTRACT ...... viii
CHAPTER 1. PRELIMINARIES ...... 1
1.1 Introduction ...... 1
1.2 Discrete-time Single-type Galton-Watson Branching Processes ...... 3
1.2.1 Definitions and Notations ...... 3
1.2.2 Limit Theorems ...... 4
1.3 Discrete-time Multi-type Galton-Watson Branching Processes ...... 7
1.3.1 Definitions, Assumptions and Notations ...... 7
1.3.2 Limit Theorems ...... 10
1.4 Continuous-time Single-type Age-dependent Bellman-Harris Branching Processes . . 13
1.4.1 Definitions, Assumptions and Notations ...... 13
1.4.2 Limit Theorems ...... 14
1.4.3 Age Distribution in Bellman-Harris Processes ...... 16
1.5 Continuous-time Multi-type Age-dependent Bellman-Harris Branching Processes . . . 18
1.5.1 Definitions, Assumptions and Notations ...... 18
1.5.2 Limit Theorems ...... 19
1.5.3 Age Distribution in Multi-type Age-dependent Branching Processes ...... 19
CHAPTER 2. REVIEW OF THE COALESCENCE IN DISCRETE-TIME SINGLE-TYPE
GALTON-WATSON BRANCHING PROCESSES ...... 21
2.1 Introduction ...... 21
2.2 The Supercritical Case ...... 23 iv
2.3 The Critical Case ...... 23
2.4 The Subcritical Case ...... 24
2.5 The Explosive Case ...... 25
CHAPTER 3. COALESCENCE IN DISCRETE-TIME MULTI-TYPE GALTON-WATSON
BRANCHING PROCESSES ...... 27
3.1 Introduction ...... 27
3.2 Results in The Supercritical Case ...... 28
3.2.1 The Statements of Results ...... 29
3.2.2 The Proof of Theorem 3.1 ...... 30
3.2.3 The Proof of Theorem 3.2 ...... 32
3.2.4 The Proof of Theorem 3.3 ...... 33
3.3 Results in The Critical Case ...... 35
3.3.1 The statements of Results ...... 35
3.3.2 The Proof of Theorem 3.5 ...... 37
3.3.3 The Proof of Theorem 3.6 ...... 41
3.3.4 The Proof of Theorem 3.7 ...... 43
3.4 Results in The Subcritical Case ...... 46
3.4.1 The Statements of Results ...... 46
3.4.2 The Proof of Theorem 3.8 ...... 48
3.4.3 The Proof of Theorem 3.9 ...... 50
3.4.4 The Proof of Theorem 3.10 ...... 51
3.5 The Markov Property on Types ...... 54
3.5.1 The Statements of Results ...... 54
3.5.2 The Proof of Theorem 3.11 ...... 56
3.5.3 The Proof of Theorem 3.12 ...... 59
CHAPTER 4. COALESCENCE IN CONTINUOUS-TIME SINGLE-TYPE AGE-DEPEN-
DENT BELLMAN-HARRIS BRANCHING PROCESSES ...... 66
4.1 Introduction ...... 66 v
4.2 Results in The Supercritical Case ...... 67
4.2.1 The statement of Results ...... 67
4.2.2 The proof of Theorem 4.1 ...... 69
4.2.3 The proof of Theorem 4.2 ...... 73
4.2.4 The proof of Theorem 4.3 ...... 75
4.2.5 The proof of Theorem 4.4 ...... 85
4.3 Results in The Subcritical Case ...... 86
4.3.1 The statements of Results ...... 86
4.3.2 The proof of Theorem 4.5 ...... 87
4.3.3 The proof of Theorem 4.6 ...... 96
CHAPTER 5. COALESCENCE IN CONTINUOUS-TIME MULTI-TYPE AGE-DEPEND-
ENT BELLMAN-HARRIS BRANCHING PROCESSES ...... 101
5.1 Introduction ...... 101
5.2 Results in The Supercritical Case ...... 102
5.2.1 The statements of Results ...... 102
5.2.2 The proof of Theorem 5.1 ...... 104
5.2.3 The proof of Theorem 5.2 ...... 114
5.2.4 The proof of Theorem 5.3 ...... 119
5.3 The Generation Problem in Supercritical Case ...... 122
5.3.1 The statement of Result ...... 123
5.3.2 The proof of Theorem 5.4 ...... 123
CHAPTER 6. APPLICATION TO BRANCHING RANDOM WALKS ...... 128
6.1 Introduction ...... 128
6.2 Review of Results in The Supercritical Case ...... 129
6.3 Results in The Explosive Case ...... 130
6.3.1 The Statements of Theorems in The Explosive Case ...... 131
6.3.2 The Proof of Theorem 6.3 ...... 132
6.3.3 The Proof of Theorem 6.4 ...... 137 vi
CHAPTER 7. OPEN PROBLEMS ...... 139
7.1 Problems in Discrete-time Multi-type Galton-Watson Branching Processes ...... 139
7.2 Problems in Continuous-time Single-type Bellman-Harris Branching Processes . . . . 139
7.3 Problems in Continuous-time Multi-type Bellman-Harris Branching Processes . . . . . 140
BIBLIOGRAPHY ...... 141 vii
ACKNOWLEDGEMENTS
First of all, I would like to take this opportunity to express my appreciation to my advisor Dr.
Krishna B. Athreya for his excellent teaching and expert guidance. From the knowledge of mathematics to the advice of study, research and even life, I always feel so comfortable and encouraged to discuss with him and learn from him. Secondly, I would like to thank the members of my committee: Dr.
Ananda Weerasinghe, Dr. Paul E. Sacks, Dr. Clifford Bergman and Dr. Dan Nordman for their precious advices and suggestions. Also, my thanks are given to Dr. Jerold Mathews for the wonderful reading time he spent in improving my English and to Dr. Leslie Hogben for her helps in many ways.
Furthermore, I would like to give my thanks to Dr. Jhishen Tsay (National Sun Yat-sen University,
Taiwan) and his wife Chun-Tor Lee for their continuous encouragement when I was struggling and depressed. In addition, I would like to thank all the friends from LIFE at Campus Baptist Church for everything they have done for me.
Finally, I would like to give my deepest appreciation to my mother, my sister, Mou-Mou, My brother, Shih-Hsun, and my brother-in-low, Cohuahua, for their love, care, understanding and patience.
With their endless support in my life, I was able to pursue my dreams. I also give my special thanks to my husband, Chao-Chun, for his love and company during these years and it means a lot to me. viii
ABSTRACT
For branching processes, there are many well-known limit theorems regarding the evolution of the population in the future time. In this dissertation, we investigate the other direction of the evolution, that is, the past of the processes. We pick some individuals at random by simple random sampling without replacement and trace their lines of descent backward in time until they meet. We study the coalescence problem of the discrete-time multi-type Galton-Watson branching process and both the continuous- time single-type and multi-type Bellman-Harris branching processes including the generation number, the death time (in the continuous-time processes) and the type (in the multi-type processes) of the last common ancestor ( also called the most recent common ancestor) of the randomly chosen individuals for the different cases (supercritical, critical, subcritical and explosive). 1
CHAPTER 1. PRELIMINARIES
1.1 Introduction
The study of branching processes has a long history and was essentially motivated by the observa- tion of the extinction of certain family lines of the European aristocracy in contrast to the rapid exponen- tial growth of the whole population. Francis Galton formulated this extinction problem and originally posed it in the Educational Times in 1874 and the Reverend Henry William Watson replied with a so- lution (see Harris (1963)). Seneta and Heyde (1977) have pointed out that the French mathematician
Bienayme´ had formulated essentially the same model fifty years earlier.
The model of Galton and Watson (called the Galton-Watson branching process) appeared to have been neglected for many years after its creation. After 1940, interest in this model increased, partly because of the analogy between the growth of families and nuclear chain reactions and also partly because of the increased general interest in applications of probability theory. Since then, branching processes have been regarded as appropriate probability models for the description of the behavior of systems whose components (cells, particles, individuals in general) reproduce, are transformed, or die
(see Harris [20], Athreya and Ney [5], Jagers [22], Mode [28] and Sevastyanov [35]). Nowadays, this theory is an area of active and interesting research.
There are many generalizations of the single-type Galton-Watson branching process in discrete time. Of these the multi-type branching process model in discrete time is a natural one. The multi- type branching process is important because it is constructed in a way that closely matches real-life situations and hence can be used to study a wide variety of real-life problems, including those related to differences in types of ethnicities, types of genes, types of cosmic rays, etc. Another generalization is the continuous-time single-type case known as the Bellman-Harris branching process which is widely used by many fields. This device was suggested by Scott and Uhlenbeck (1942) in their treatment of 2 cosmic rays, where the continuous variable is energy, and was used by Bartlett (1946) and Leslie (1948) in dealing with human population, where the continuous variable is age.
In the rest of this chapter, we review basic definitions and results of single-type (Section 1.2) and multi-type (Section 1.3) discrete-time Galton-Watson branching processes and single-type (Section 1.4) and multi-type (Section 1.5) continuous-time Bellman-Harris age-dependent branching processes. We discuss results on the extinction probabilities, the growth rates of population and some other convergent properties. The results are fundamental and may be found in the books on branching processes men- tioned earlier. Here, we state the results which are needed in this thesis based on the books, Branching
Processes written by Athreya and Ney [5] and The Theory of Branching Processes written by T. E.
Harris [20].
In Chapter 2, we state the problem of coalescence in branching processes and review the results for all cases (supercritical, critical ,subcritical and explosive) of the discrete-time single-type Galton-
Watson branching processes.
In Chapter 3, we extend the results of the problem of coalescence to the discrete-time multi-type
Galton-Watson branching process including supercritical (Section 3.2), critical (Section 3.3) and sub- critical (Section 3.4) cases. Also, we present the Markov property on Types (Section 3.5) along the line of descent of an individual randomly chosen from the current generation by simple random sampling.
In Chapter 4, we consider the continuous-time single-type Bellman-Harris branching processes and give proofs to the problem of coalescence in the supercritical (Section 4.2) and subcritical (Section
4.3) cases. By the results on the problem of coalescence, we also are able to investigate the branching random walks (Section 4.4).
Although the research of branching processes has a long history, the study of the problem of co- alescence is still in its infancy. In Chapter 5, we state some interesting open questions related to this topic. 3
1.2 Discrete-time Single-type Galton-Watson Branching Processes
1.2.1 Definitions and Notations
A discrete-time single-type Galton-Watson branching process is the simplest type of branching process. This process can be thought as a population evolving in time. It starts at time 0 with Z0 individuals, each of which lives a unit of time and produces its offsprings upon death according to the probability distribution {p j} j≥0 independently of others. Let Z1 be the total number of children produced by the Z0 individuals, that is,
XZ0 ξ0,i i=1 where {ξ0,i}i≥1 are i.i.d. random variables with the probability distribution {p j} j≥0. It constitutes the first generation and then these individuals in the first generation go on to produce the second generation of population Z2 and so on. So, the total size of the population in the (n + 1)st generation, n = 0, 1, 2, ··· , is given by XZn ξn,i if Zn > 0 Z = n+1 i=1 0 if Zn = 0 where {ξn,i : i ≥ 1, n ≥ 0} are i.i.d. copies with the probability distribution {p j} j≥0.
Then {Zn}n≥0 is called a Galton-Watson branching process with initial population Z0 and offspring distribution {p j} j≥0. Here, ξn,i is the number of offspring of the ith individual of the nth generation. Let
X∞ m ≡ jp j j=1 be the mean of the offspring distribution {p j} j≥1. We shall refer to the Galton-Watson process as sub- critical, critical, supercritical or explosive according as 0 < m < 1, m = 1, 1 < m < ∞ or m = ∞, respectively.
Moreover, if T denotes the full family tree generated in this way, every individual in T can be identified by a finite string (i0, i1, ··· , in) meaning that this individual is in the nth generation and is the inth child of the individual (i0, i1, ··· , in−1) of the (n − 1)st generation. 4
1.2.2 Limit Theorems
In this section, we collect some well-known results for discrete-time single-type Galton-Watson branching processes.
Theorem 1.1. (Supercritical Case) Let p0 = 0 and 1 < m < ∞. Then
(a) P(Zn → ∞|Z0 > 0) = 1.
Z (b) (Harris, 1960) W ≡ n : n ≥ 0 is a nonnegative martingale and hence n mn
lim Wn ≡ W exists w.p.1. n→∞
(c) (Kesten and Stigum, 1966)
X∞ ( j log j)p j < ∞ if and only if E(W|Z0 = 1) = 1 j=1
and then W has an absolutely continuous distribution on (0, ∞) with a positive density.
(d) (Seneta and Heyde, 1970)
Cn+1 Zn ∃Cn s.t. → m and → W w.p.1 Cn Cn
X∞ n and P(0 < W < ∞) = 1. In particular, ( j log j)p j < ∞ if and only if Cn ∼ m . j=1 (e) (Athreya and Schuh [4])
E(W : W ≤ x) ≡ L(x)
L(cx) is slowly varying at ∞, i.e. ∀0 < c < ∞, → 1 as x → ∞. L(x)
Under the assumption p0 = 0, the population size Zn of a supercritical Galton-Watson branching process goes to infinity as n → ∞ with probability 1 and it grows like mn. This is the stochastic analogue of the so-called Malthusian law of geometric population growth.
In the next two theorems, we present the results for the critical and subcritical cases.
∞ 2 X 2 Theorem 1.2. (Critical Case) Let m = 1, p j , 1 for any j ≥ 1 and σ ≡ j p j − 1 < ∞. Then j=1 5
(a) P(Zn → 0|0 < Z0 < ∞) = 1.
(b) (Kolmogrov, 1938)
σ2 nP(Z > 0) → as n → ∞. n 2
(c) (Yaglom, 1947)
Z − 2 n 2 x P > x Zn > 0 → e σ , 0 < x < ∞. n
(d) (Athreya [12]) For 1 ≤ k ≺ n, let
(k) Zn−k,i ≡ (k) ≤ ≤ Vn,k I(Z >0) : 1 i Zk n − k n−k,i
(k) on the event {Zk > 0}, where {Z j,i : j ≥ 0} is the G-W process initiated by the ith individual in the kth generation.
k Let k → ∞, n → ∞ such that n → u, 0 < u < 1.
Then the sequence of point processes {Vn,k}n≥1 conditioned on {Zn ≥ 1} converges weakly to the point process
V ≡ {η j : j = 1, 2, ··· , Nu}
k−1 where {η j} j≥1 are i.i.d. exp(1),Nu is Geo(u), i.e., P(Nu = k) = (1 − u)u , k ≥ 1 and {η j} j≥1 and
Nu are independent. X∞ Theorem 1.3. (Subcritical Case) (Yaglom, 1947) Let 0 < m ≡ jp j < 1. Then j=1 ∞ ∞ X X j (a) For j ≥ 1, lim P(Zn = j|Zn > 0) ≡ b j exists, b j = 1 and B(s) ≡ b j s , 0 ≤ s ≤ 1 is the n→∞ j=0 j=0 unique solution of the functional equation
B( f (s)) = mB(s) + (1 − s) , 0 ≤ s ≤ 1
∞ X j where f (s) ≡ p j s , in the class of all probability generating functions vanishing at 0. j=0
X∞ X∞ (b) jb j < ∞ iff ( j log j)p j < ∞. j=1 j=1 6
P(Z > 0|Z = 1) 1 (c) lim n 0 = . n→∞ mn P∞ jb j j=1
(d) If Z0 is a random variable and EZ0 < ∞, then
lim P(Zn = j|Zn > 0) = b j , ∀ j ≥ 1 n→∞ X∞ and if, in addition, ( j log j)p j < ∞ then j=1
∞ X P(Zn > 0) EZ0 jb j < ∞ and lim = . n→∞ mn P∞ j=1 jb j j=1
In both of the critical and subcritical Galton-Watson branching processes, the population will die out eventually with probability 1. But, conditioned on the event of non-distinction, i.e. the set {Zn > 0},
Zn will go to infinity in distribution with the growth rate of n in the critical case while Zn will converge to a proper random variable in distribution in the subcritical case as n → ∞.
We present the results of P. L. Davies and D. R. Grey for the explosive Galton-Watson branching process as follows.
Theorem 1.4. (Explosive Case) Let p0 = 0, m = ∞ and, for some 0 < α < 1, P p j j>x → 1 as x → ∞ xαL(x) where L : (1, ∞) → (0, ∞) is a function slowly varying at ∞. Then
(a) (Davies [16])
n α log Zn → η w.p.1
and P(0 < η < ∞) = 1 and η has a continuous distribution.
(1) (2) (b) (Grey [19]) Let {Zn }n≥0 and {Zn }n≥0 be two i.i.d. copies of a Galton-Watson branching process X∞ { } ≡ ∞ (1) (2) with the offspring distribution p j j≥0 satisfying p0 = 0, m jp j = and Z0 = Z0 = 1. j=1 Then, w.p.1, 1 Z(1) 0, with prob. n → 2 (2) Zn 1 ∞, with prob. 2 . 7
It is easy to deduce b) from a) in the above theorem as
n (1) (2) α log Zn − log Zn → η1 − η2 ≡ η, say and
1 P(η > 0) = P(η < 0) = . 2
1.3 Discrete-time Multi-type Galton-Watson Branching Processes
1.3.1 Definitions, Assumptions and Notations
In a discrete-time single-type Galton-Watson branching process, we assume that each individual lives for a fixed unit time and then produces its children according to the same offspring distribution.
In this section, we allow a number of distinguishable types of individuals having different offspring distributions.
First, we consider a finite number d of individual types. Such processes arise in a variety of appli- cations in biology and physics and they could represent genetic or mutant types in the real populations such as animal population, bacterial population or photons, etc.
Through out this section and next chapter, we adopt the following conventions.
1. N0 is the set of all nonnegative integers.
n o d ≡ ≡ ··· ∈ ··· 2. N0 j ( j1, j2, , jd): ji N0, i = 1, 2, , d
··· ··· d 3. 0 = (0, 0, , 0) and 1 = (1, 1, , 1) in N0
··· ··· ∈ d 4. ei = (0, , 0, 1, 0, , 0) N0 with the 1 in the ith component.
5. u ≤ v means ui ≤ vi for i = 1, 2, ··· , d while u < v means ui ≤ vi for all i and ui < vi for at least one i.
6. The vector of absolute values is
|x| = |x1| + |x2| + ··· + |xd| 8
7. The sup norm is
kxk = max{|x1|, |x2|, ··· , |xd|}
8. The product notation is
Yd y yi x = xi i=1
9. For a matrix M, the super norm is
kMk = max{|mi j| : i, j = 1, 2, ··· , d}
Let Zn = (Zn,1, Zn,2, ··· , Zn,d) be the population vector in the nth generation, n = 0, 1, 2, ··· , where
Zn,i is the number of individuals of type i in the nth generation. We assume that each individual of type i, i = 1, 2, ··· , d, lives a unit of time and, upon death, produces children of all types and according to
(i) (i) the offspring distribution p (j) ≡ p ( j1, j2, ··· , jd) j∈Nd and independently of other individual, where (i) p ( j1, j2, ··· , jd) is the probability that a type i parent produces j1 children of type 1, j2 children of type 2, ··· , jd children of type d. Therefore, each component of the vector of the probability generating functions f = f (1), f (2), ··· , f (d) can be written as:
X (i) ··· (i) ··· j1 j2 ··· jd f (s1, s2, , sd) = p ( j1, j2, , jd)s1 s2 sd j1, j2,··· , jd≥0 where 0 ≤ sr ≤ 1, r = 1, 2, ··· , d, being the probability generating function of the number of various types produced by a type i individual,
d Thus, a discrete-time multi-type Galton-Watson branching process Zn n≥0 is a Markov chain on N0 with the transition function
d P(i, j) = P(Zn+1 = j|Zn = i) ∀i, j ∈ N0 ∞ X i such that P(i, j)sj = f(s) (see notation (8)). ∈ d j N0 When the process is initiated in state ei, we will denote the process {Zn}n≥0 by
(i) (i) (i) ··· (i) Zn = Zn,1, Zn,2, , Zn,d
(i) where Zn, j is the number of type j individuals in the nth generation for a process Z0 = ei. The generating (i) (i) function of Zn will be denoted by fn (s). 9
( j) Also, we let ξn,r be the vector of offsprings of the rth individual of type j in the nth generation then
( j) ( j) ( j) ( j) ξn,r ∼ {p (·)}, i.e., P(ξn,r = ·) = p (·). Then, the population in the (n + 1)th can be expressed as
d Zn, j X X ( j) Zn+1 = ξn,r. j=1 r=1
Let mi j = E(Z1, j|Z0 = ei) be the expected number of type j offspring of a single type i individual in one generation for any i, j = 1, 2, ··· , d. Then, we define the mean matrix
M = {mi j : i, j = 1, 2, ··· , d}.
n (n) n Clearly, we get E(Zn|Z0) = Z0M . We let mi j be the (i, j)th element of M . When the higher moments exist, we can denote them by the following notations. First, we let
(r) (r) (r) (r) qn (i, j) = E Zn,i Zn, j − δi, jZn,i i, j, r = 1, 2, ··· , d and define the matrix
(r) (r) Qn = qn (i, j): i, j = 1, 2, ··· , d ,