1 Introduction
Let be a sequence of
independent and identically distributed (i.i.d.), non-negative and integer-valued random variables defined on a probability space . Let also , for , be a sequence of stochastic processes which consist of independent non-negative
integer-valued random variables on with the same one-dimensional distributions. Furthermore, let us assume that and
are independent.
A controlled branching process (CBP) is defined recursively as
|
|
|
(1) |
where is defined as 0 and is a non-negative, integer-valued,
square-integrable random variable which is independent of
and .
Here, denotes the size of the -th generation of a
population and is the offspring size of the -th
individual in the -th generation. We will assume that
the mean and variance are both finite.
The class of CBPs is a very general family of stochastic processes that collect as particular cases the simplest branching model, the standard Bienaymé–Galton–Watson (BGW) process, by considering a.s. for each , or a branching processes
with immigration, by setting , where
are i.i.d. random variables (writing in this way the immigrants give rise to offspring at the same generation as their arrival and with the same probability law as ),
among others. The monograph [7] provides an extensive description of its probabilistic theory.
The research of functional weak limit theorems for branching processes arises
a lot of interest since many years ago. It was firstly formulated for a BGW process by [3] and proved by [10] and [12]. These results have been extended to another classes of branching processes. For instance, a wide literature exists around weak convergence results for branching processes with immigration (BPI) since the pioneer work by [15], see also [1] and references therein. In this paper we focus our attention on a weak convergence theorem for a critical CBP with a random initial number of individuals and assuming finite second order moment on the this initial value. A similar result was already established for a single CBP in [14], and for an array of CBPs in [6], by assuming fixed initial numbers of progenitors using infinitesimal generators results for their proofs. Inspired in the paper [1] on BPI we will use limit theorems for random step processes
towards a diffusion process provided in [9] to obtain an alternative proof. The scheme of it follows similar steps to the ones in [1]. An important feature of a CBP is that the value of conditioned on the knowledge of the previous generation, , is a random sum of random variables, namely , instead of a non-random sum as in the case of a BPI. This leads to handle the proofs of each steps using conditioning arguments different from those used in [1].
Apart from this introduction, the paper is organized as follows. In Section 2 we provide
the notation and some auxiliary results about the behaviour of the first and second moments of the process.
Section 3 gathers the main theorem. For the ease of reading the paper, additional results are presented in the Appendix.
2 Notation and auxiliary results
We denote, for
|
|
|
|
|
|
|
|
|
|
and assume all finite.
It is easy to obtain that for ,
|
|
|
|
|
(2) |
|
|
|
|
|
(3) |
where is the algebra generated by the random variables , (see Proposition 3.5 in [7]).
We introduce the quantities
|
|
|
(4) |
The quantity represents a mean growth
rate. Intuitively, it can be interpreted as an average offspring
per individual for a generation of size .
Assuming that exists, the process can be classified as:
|
|
|
We are interested in critical CBPs that satisfy the following hypotheses:
-
A1)
, , ,
-
A2)
, , as .
The behavior of critical CBPs was studied in
[5]. Assuming that , i.e. 0 is an
absorbing state, and it is verified or , , it was established that under A1) and A2), if and an assumption on conditional moments holds, then In the present paper we will consider critical CBPs, , satisfying the above conditions, but with a reflecting barrier at zero, namely, . Thus will have a finite number of returns to the sate zero till the explosion to infinity, i.e. .
Notice that under A1), , , and, for simplicity in the posterior calculations, we will also assume throughout the paper that .
Remark 2.1.
The controlled branching process we are considering is such that migration may take place in the next generation no matter the size of the current generation (when there are no individuals in the populations only immigration is possible). BGW processes with immigration at 0 were considered firstly in [4] and [13].
In next result we calculate the first and second moments of a CBP which verifies A1) and A2).
Proposition 2.1.
Let be a CBP
with and satisfying hypotheses A1) and A2). It is verified as that
|
|
|
Proof.
From (2) and A1) it follows that
|
|
|
(5) |
Using (3) we have
|
|
|
|
|
Now, from A2), we have that there exists such that for all , so that .
Consequently, letting and , we have
|
|
|
(6) |
Hence,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
The latter inequality proves that .
Next Lemma presents certain relationships among the random variables which can be easily verified.
Lemma 2.1.
Let be a sequence of
i.i.d. random variables
with mean and , assumed finite.
Let denote , , , and for , and .
It is verified that
|
|
|
and
|
|
|
3 Main result
We introduce for each
, a stochastic process , for
, , denoting the integer part. It is
easy to see that is a sequence of random
functions that take values in , which is
the space of non-negative functions on that are right
continuous and have left limits. We also denote by
the space of infinitely differentiable
functions on which have a compact support. Throughout the paper “” denotes the convergence
of random functions in the Skorokhod topology.
Theorem 3.1.
Let be a CBP
with , satisfying hypotheses A1) and A2).
Then,
,
as being a non-negative diffusion process, with
generator , for .
The process is the pathwise unique solution of the stochastic
differential equation
|
|
|
(7) |
with initial value , denoting , , and where is a standard Wiener process.
Remark 3.1.
Taking into account Theorem A2 in Appendix, the stochastic differential equation (SDE) (7) has a pathwise unique solution for all initial values . Moreover if , then almost surely for all .
In order to prove Theorem 3.1, we will establish previously the weak convergence of random step processes defined from a martingale difference created from the CBP.
We introduce the following sequence of martingale differences with respect the filtration as:
|
|
|
Consider the random step processes:
|
|
|
(8) |
Theorem 3.2.
Let be a CBP
with , satisfying hypotheses A1) and A2).
It is verified that
|
|
|
where the limit process is the pathwise unique solution of
|
|
|
(9) |
Proof As was done in [1], we prove the result by applying Theorem A1 in Appendix with , , ,
where (yielding , as well), and with coefficient functions
and given by
|
|
|
Firstly, we check that the SDE (9) has a pathwise unique strong solution for all initial values . In fact, notice that if is a strong solution of the SDE
(9) with initial value , then, by Itô’s formula, the process
, is a solution of the SDE
|
|
|
(10) |
Conversely, if is a strong solution of the SDE (10) with initial value , then,
by Itô’s formula, the process
, is a strong solution of the SDE (9) with initial value . Notice that SDE (10) is the same as SDE (7), consequently the SDE (10) and therefore the SDE (9) as well admit a pathwise unique strong solution with arbitrary initial value, and
|
|
|
(11) |
Let us see that for all and . Indeed, taking into account (6) in Proposition 2.1,
|
|
|
(12) |
and, by the assumption in the statement of the theorem, for
Moreover, as , especially as .
For conditions (i), (ii) and (iii) of Theorem A1 in Appendix, we have to check that for each , , as :
-
a)
-
b)
-
c)
For all ,
Since , , , a) holds.
For each and , and for all , and , we have:
|
|
|
thus . Now, we have, for all and ,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
It is verified that, for and ,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Consequently,
|
|
|
|
|
|
|
|
|
|
Since for each ,
|
|
|
in order to show , it suffices to prove that for each ,
|
|
|
(13) |
and
|
|
|
(14) |
First we check (13). For each , we have , thus
|
|
|
and hence, for each and , we get
|
|
|
Consequently, in order to prove , it suffices to show
|
|
|
By (12), as , and therefore by Jensen’s inequality, as , and hence
|
|
|
Thus we obtain as implying .
Now, taking into account hypothesis A2)
|
|
|
(15) |
and hence (14).
We write
|
|
|
Let denote . It is verified for each , and that
|
|
|
and
|
|
|
Hence
|
|
|
In consequence, to check c) we will prove, as ,
-
c.1)
for all .
-
c.2)
for all .
-
c.3)
In what follows let be fixed.
Let us see c.3). It is verified that
|
|
|
|
|
|
|
|
|
|
This latter was proved by considering (15).
Now, we check c.1).
By the properties of conditional expectation with respect to a -algebra, we get for all ,
|
|
|
where on , with
|
|
|
Consider the decomposition with
|
|
|
|
|
|
|
|
|
|
Now, let denote , , . It is verified the inequality, for
|
|
|
We have, using Lemma 2.1,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Therefore
|
|
|
with
|
|
|
|
|
|
|
|
|
|
Using (5) in Proposition 2.1, it is verified that for
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
By applying the dominated convergence theorem we have
|
|
|
(16) |
It is also verified by using again (5) in Proposition 2.1 and A2) that
|
|
|
|
|
(17) |
|
|
|
|
|
Taking into account (16) and (17)
we have that, as ,
|
|
|
Let us now dealt with . It is verified that, using Cauchy-Schwarz’s inequality:
|
|
|
|
|
|
|
|
|
|
Now, using Markov’s inequality
|
|
|
and using
Lemma 2.1, we have
|
|
|
By Lyapunov’s inequality, . Consequently, in order to prove, as ,
|
|
|
is enough to check
that
|
|
|
In fact, using hypotheses A1) and A2) and Proposition 2.1, we have
|
|
|
Finally, we check c.2). We have that
|
|
|
where on , with ,
|
|
|
Now, again by Cauchy-Schwarz’s inequality and Markov’s inequality
|
|
|
|
|
|
|
|
|
|
In consequence from
|
|
|
c.2) follows.
Finally, using the weak convergence of , we will obtain weak convergence of .
Proof of Theorem 3.1.
A version of the continuous mapping theorem is applied (see Lemma 1 in Appendix). Let be the space of the real functions on that are right
continuous and have left limits. For each , by (8), , where the mapping
is given by
|
|
|
for and . Indeed, for each and :
|
|
|
Further, taking into account (11), , where the mapping is given by
|
|
|
The measurability of , and likewise the conditions for applying Lemma 1 are proved in [1].
Remark 3.2.
1) Notice that the result in Theorem 3.1 is also valid as , and even when the hypothesis A1) is replaced with the more general condition , as , , being the calculation in this latter scenario a little more cumbersome. In the case and the result provides an alternative proof of the weak convergence result for the BGW process (see [2], p. 388) for a non-array version.
2) As noted in the introduction, a BPI can be written as a special case of a CBP. For this case, by considering ( and ), we obtain an analogous result to that in [1]. Recall that the difference between both models is in which generation the immigrants will give rise to their offspring.
3) As was pointed out in the Introduction, the proof of the main result follows similar steps as those given in [1]. One can check that similar formulas often appear being the roles of the immigration mean and the offspring variance in the BPI case played in the CBP by and , respectively. However, new approaches by considering conditioning arguments are needed to dealt with b)- c) in Theorem 3.1, as a consequence that random sums of i.i.d. random variables arise in the proofs, see for instance the definition of and , in p.3. An extra work is required to calculate the mathematical expectation of these quantities.
Acknowledgements: The authors thank Professor M. Barczy for his constructive suggestions which have improved this paper.
This research has been supported by the Ministerio de Ciencia e Innovación of Spain (grant PID2019-108211GB-I00/AEI/10.13039/501100011033).