Journal of Statistical Theory and Applications

Volume 20, Issue 1, March 2021, Pages 111 - 131

A Note on Generalization of the Simplest Time-Dependent Discrete Markov Process: Linear Growth Process With Immigration-Emigration

Authors
Bijoy Kumar Pradhan1, *, Priyaranjan Dash1, Upasana
1Department of Statistics, Utkal University, Bhubaneswar, Odisha, 751004, India
*Corresponding author. Email: prdashjsp@gmail.com
Corresponding Author
Bijoy Kumar Pradhan
Received 24 January 2019, Accepted 27 April 2020, Available Online 5 February 2021.
DOI
10.2991/jsta.d.210126.003How to use a DOI?
Keywords
Birth; Death; Immigration and Emigration Process; Linear growth process
Abstract

Lineargrowth process with immigration and emigration is the general model in the study of population in biological and ecological systems, and their transient analysis is the most important factor in the understanding of the structural behavior of such systems. The probability-generating function π(z, t) of the probability distribution {pn(t)} of the random variable N(t) in many queuing situations concerning Birth, Death, Immigration, and Emigration were studied and we find the generalization of the transient solutions of the queuing systems and also studied its particular cases.

Copyright
© 2021 The Authors. Published by Atlantis Press B.V.
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

1. INTRODUCTION

Transient queuing systems became the focus of the researchers since early 1940s working in the area of queuing theory and also modeling different biological populations, disease modeling specifically in the area of epidemiology but subsequently it became more popular in several areas like modeling the number of animals belonging to a particular forest, the fish population in a lake, spreading of bacteria and viruses in a specific area, logistic services like modeling operation efficiency in railway container terminal, analyzing buffer arrivals in modeling telecommunication systems, signal processing, and many more. The present paper generalizes a time-dependent Markovian process in presence of births, deaths, immigration, and emigration. The characterizing principle of such type of processes can be highlighted on the basis of the generating function of probabilities of the number of instances/occurrences. So, we focus our attention on the generating function of the process.

The generating function

π(z,t)=n=0pn(t)zn(1.1)
for many transient queuing situations satisfy a partial differential difference equation of the form
Q(z)π(z,t)z+R(z)π(z,t)t=1,(1.2)
with pn(t)=Pr{N(t)=n}, where {pn(t)} represents probability distribution of the random variable N(t) for every value of t [?].

In order to solve this partial differential difference equation subject to some appropriate boundary conditions, we are required to form the subsidiary equations given by

dπ(z,t)1=dzQ(z)=dtR(z).(1.3)

Taking into considerations of the first and second parts of (1.3) and then from second and third parts of (1.3), we find two independent solutions of the subsidiary equations writing them in the form

u(z,t,π)=aconstant,ν(z,t)=aconstant.(1.4)

From (1.4), we find the general solution of (1.2), i.e., the probability-generating function π(z,t) of the probability distribution {pn(t)}. We may verify that this solution satisfies the partial differential difference equation given by (1.2).

Also, one may verify that the mean and variance of the probability distribution {pn(t)} found by considering the expression for pn(t) are same as the mean and variance of the probability distribution by considering the probability-generating function π(z,t) on the above procedure.

The present paper deals with linear growth process with immigration and emigration, where the probability-generating function is derived and its mean and variance are calculated directly by considering the expression for pn(t) and also from the the probability-generating function π(z,t) for the probability distribution {pn(t)}. We found that both procedures give the same result. From the probability-generating function π(z,t) of the probability distribution {pn(t)} of this process, we can derive the probability-generating functions of Linear Birth Process, Linear Death Process, Birth Immigration Process, Birth Emigration Process, Death Immigration Process, Death Emigration Process, Immigration-Emigration Process, Linear Growth Process with Immigration, Linear Growth Process with Emigration, Birth process with Immigration-Emigration, Death process with Immigration-Emigration as the particular cases and also we can find the mean and variance of all these transient queuing processes from the mean and variance of linear growth process with immigration and emigration as particular cases.

Getz [1] has also found the probability-generating function by considering the wrong procedure, i.e., by considering the first and third parts of (1.3) and the second and third parts of (1.3) in his paper. The probability-generating function found by him does not satisfy the partial differential equation of this process. He found the mean and variance by considering the expressions for pn(t) of the probability distribution {pn(t)} of this process under “Direct derivation of Mean and Variance.” He found the correct mean but incorrect variance. We would like to mention here that the probability-generating function given by Getz [1] is not correct.

Zheng et al. [2] have also tried to find the probability-generating function of this process. But they could not able to find the correct differential difference equation of this process because they considered additional term θ(112)P(t,0) in the differential difference equation of this process (here, θ stands for Emigration rate). In this regard, we would like to mention here that the term θP(t,0) does not exist in the process because when the population size reduces to zero at time t, there will be no emigration and hence we set θP(t,0)=0.

2. ASSUMPTIONS

  1. If n individuals are present at the instant from which the interval commences, the probability of one birth will occur in any short interval of length h is nλh+o(h) and the probability of more than one birth occurring in that small interval is o(h). Birth occurring in (t,t+h) are independent of time since the last occurrence.

  2. If n individuals are present at the instant from which the interval commences, the probability of one death will occur in any short interval of length h is nμh+o(h) and the probability of more than one death occurring in that small interval is o(h). Death occurring in (t,t+h) are independent of time since the last occurrence.

  3. During the small interval of time (t,t+h), the probability that a new member being added to the population by immigration is νh+o(h) and the probability that more than one individual is added to the population in that small interval of length h is o(h). Immigration occurring in (t,t+h) are independent of time since the last occurrence.

  4. During the small interval of time (t,t+h), the probability that a member being left by emigration from the population is ξh+o(h) and the probability that more than one emigration occurring in that small interval of length h is o(h). Emigration occurring in (t,t+h) are independent of time since the last occurrence.

  5. For the same population, there is no interaction among the birth, death, immigration, and emigration in small interval (t,t+h) of time.

3. DERIVATION OF PROBABILITY-GENERATING FUNCTION

Let pn(t) be the probability that the process starts with n individuals at instant t, i.e.,

pn(t)=Pr[N(t)=n].

To calculate pn(t+h) at the next time point (t+h), the system can be instate En only if one of the following conditions are satisfied:

  1. At time t, the population consist of n individuals and no birth, no death, no immigration, no emigration occur during the time interval (t,t+h).

  2. At time t, the population consist of (n1) individuals and a birth occurs during the time interval (t,t+h) but no death, no immigration, no emigration occur during that interval of time.

  3. At time t, the population consist of (n1) individuals and a new member is added to the population by immigration during the time interval (t,t+h) but no birth, no death, no emigration occur during that interval of time.

  4. At time t, the population consist of (n+1) individuals and a death occurs during the next time interval (t,t+h) but no birth, no immigration, no emigration occur during that interval of time.

  5. At time t, the population consist of (n+1) individuals and an individual left the process during (t,t+h) but no birth, no death, and no immigration occur during that interval of time.

  6. During (t,t+h) two or more transitions occur with probability o(h), i.e., the occurrence of more than one transition in (t,t+h) is negligible.

All the above contingencies are mutually exclusive and hence for n1, we have

pn(t+h)=pn(t)[1nλh+o(h)][1nμh+o(h)][1νh+o(h)][1ξh+o(h)]+pn1(t)[n1¯λh+o(h)][1n1¯μh+o(h)][1νh+o(h)][1ξh+o(h)]+pn1(t)[1n1¯λh+o(h)][1n1¯μh+o(h)][νh+o(h)][1ξh+o(h)]+pn+1(t)[1n+1¯λh+o(h)][n+1¯μh+o(h)][1νh+o(h)][1ξh+o(h)]+pn+1(t)[1n+1¯λh+o(h)][1n+1¯μh+o(h)][1νh+o(h)][ξh+o(h)]+o(h),  forλμandνξ.(3.1)

Hence,

pn(t)=limh0pn(t+h)pn(t)h=nλpn(t)nμpn(t)νpn(t)ξpn(t)+n1¯λpn1(t)+νpn1(t)+n+1¯μpn+1(t)+ξpn+1(t)aslimh0o(h)h=0,forn1.(3.2)

For n=0, in the similar manner we can find

p0(t)=ξp1(t)+μp1(t)νp0(t).(3.3)

We can write (3.2) and (3.3) in a single form by letting p1(t)=0,ξp0(t)=0,μp0(t)=0. Similarly, if the population size reduces to zero at any time point t, there will be no emigrations and no death and hence ξp0(t)=0,μp0(t)=0. Hence,

pn(t)=nλpn(t)nμpn(t)νpn(t)ξpn(t)+n1¯λpn1(t)+n+1¯μpn+1(t)+νpn1(t)+ξpn+1(t),n0.(3.4)

Define the probability-generating function of pn(t) as

π(z,t)=n=0pn(t)zn(3.5)

Hence,

π(z,t)t=n=0dpn(t)dtzn=(λzμ)(z1)π(z,t)z+(z1)(νzξ)zπ(z,t).(3.6)

We write this linear partial differential equation given by (3.6), in the form

(λzμ)zνzξlogπ(z,t)z+z(z1)(νzξ)logπ(z,t)t=1(3.7)

To solve this partial differential equation, we first form the subsidiary equations given by

dlogπ(z,t)1=dzz(λzμ)νzξ=dtz(z1)(νzξ).(3.8)

From the first and second parts of (3.8) we find

π(z,t)(λzμz)ξμ(λzμ)νλ=a constant.(3.9)

From the second and third part of (3.8) we find

λzμz1e(λμ)t=a constant(3.10)

Hence,

π(z,t)(λzμz)ξμ(λzμ)νλ=ψ[λzμz1e(λμ)t],(3.11)
where ψ() is an arbitrary function to be determined by initial conditions.

If the process starts with n0 individuals at t=0, we have

pn0(0)=1andpk(0)=0,forkn0.

And hence π(z,0)=zn0. Putting t=0 in (3.11), we find

ψ[λzμ1z]=zn0(λzμz)ξμ(λzμ)νλ(3.12)

Letting u=λzμz1, we find z=uμuλ and hence from (3.12)

ψ(u)=(uμuλ)n0[u(λμ)uμ]ξμ[u(λμ)uλ]νλ.(3.13)

From (3.11) and (3.13), after simplification, we find

π(z,t)=[(λzμz1μ)e(λμ)t]ξμ[(λzμz1λ)e(λμ)t]νλ×[λzμz1e(λμ)tμ]n0+ξμ[λzμz1e(λμ)tλ]n0+νλ,forλμandνξ.(3.14)

Hence, π(z,t) given by (3.14) is the probability-generating function of linear growth process with immigration and emigration.

Substituting z=eθ in the probability-generating function of linear growth process with immigration and emigration given by (3.14) and then find

M(t)=E[N(t)]=π(eθ,t)θ|θ=0,andM2(t)=E[N(t)]2=2π(eθ,t)θ2|θ=0
and considering V[N(t)]=M2(t)[M(t)]2, we find the mean and variance of this process given by
E[N(t)]=(νξλμ)[e(λμ)t1]+n0e(λμ)t,  andV[N(t)]=νμ(λμ)2[1e(λμ)t]+νλ(λμ)2e(λμ)t[e(λμ)t1]+ξλ(λμ)2[e(λμ)t1]+ξμ(λμ)2e(λμ)t[1e(λμ)t]+n0(λ+μλμ)e(λμ)t[e(λμ)t1].(3.15)

Corollary 1.

When λ=μ and ν=ξ, from the probability generating function of linear growth process with immigration and emigration given by (3.14), we find

π(z,t)=[z(1λt)+λt(1+λt)λtz]n0+νλzνλ.(3.16)

Also from (3.15), we find

E[N(t)]=n0andV[N(t)]=2n0λt+2νt.(3.17)

Corollary 2.

If we consider λn=nλ+ν and μn=nμ+ξ, we get also linear growth process with emigration and emigration with same probability-generating function given by (3.14).

4. DIRECT METHOD TO FIND MEAN AND VARIANCE

For linear growth process with immigration and emigration, we have

pn(t)=nλpn(t)nμpn(t)νpn(t)ξpn(t)+n1¯λpn1(t)+νpn1(t)+ξpn+1(t)+n+1¯μpn+1(t).(4.1)

Multiplying both sides of (4.1) by n and summing over all values of n, we find

M(t)=(λμ)M(t)+(νξ)(4.2)
M(t)(λμ)M(t)+(νξ)=1.(4.3)

Hence,

(λμ)M(t)+(νξ)=C1e(λμ)t,(4.4)
where C1 is a constant to be determined by the appropriate initial conditions.

When the process starts with n0 individuals at time t=0, we have

pn0(0)=1andpk(0)=0,forkn0.

And hence M(0)=n0. Substituting t=0 in (4.4), we find

C1=n0(λμ)+(νξ).(4.5)

Hence, from (4.4) and (4.5), we find after simplification

M(t)=(νξλμ)[e(λμ)t1]+n0e(λμ)t.(4.6)

Similarly, multiplying both sides of (4.1) by n2 and summing over all values of n, we find

M2(t)=2(λμ)M2(t)+(λ+μ+2ν2ξ)M(t)+(ν+ξ)=2(λμ)M2(t)+(λ+μ+2ν2ξ)[νξλμ{e(λμ)t1}+n0e(λμ)t]+(ν+ξ)(4.7), (4.8)

Hence, the general solution of (4.8) is given by

M2(t)=C2e(λμ)tνλ(λμ)2e(λμ)t+ξλ(λμ)2e(λμ)t+νλ2(λμ)2ξλ2(λμ)2n0λ(λμ)e(λμ)tνμ(λμ)2e(λμ)t+ξμ(λμ)2e(λμ)t+νμ2(λμ)2ξμ2(λμ)2n0μ(λμ)e(λμ)t2ν2(λμ)2e(λμ)t+2νξ(λμ)2e(λμ)t+ν2(λμ)2νξ(λμ)22n0ν(λμ)e(λμ)t+2νξ(λμ)2e(λμ)t2ξ2(λμ)2e(λμ)tνξ(λμ)2+ξ2(λμ)2+2n0ξ(λμ)e(λμ)tν2(λμ)ξ2(λμ)(4.9)

If the process starts with n0 individuals at time t=0, we have M2(0)=n02 and hence substituting t=0 in (4.9), we find after simplification

C2=n02+νλ(λμ)2ξλ(λμ)2+n0λ(λμ)+n0μ(λμ)+ξ(λμ)+2n0ν(λμ)2n0ξ(λμ)+ν2(λμ)2+ξ2(λμ)22νξ(λμ)2(4.10)

Substituting the value of C2 given by (4.10) in (4.9), we get the value of M2(t).

Since, V[N(t)]=M2(t)[M(t)]2, we find

V[N(t)]=νμ(λμ)2[1e(λμ)t]+νλ(λμ)2e(λμ)t[e(λμ)t1]+ξλ(λμ)2[e(λμ)t1]+ξμ(λμ)2e(λμ)t[1e(λμ)t]+n0(λ+μλμ)e(λμ)t[e(λμ)t1].(4.11)

5. SOME PARTICULAR CASES

The suggested linear growth process with immigration and emigration can be reduced to a number of time- -dependent Markov processes for suitable choice of values of λ,μ,ν, and ξ. These are discussed as follows.

5.1. Linear Birth Process

Substituting μ=0,ν=0, and ξ=0, in (3.14), we find the probability generating function of Linear Birth Process given by

π(z,t)=[zeλt1z(1eλt)]n0.(5.1)

Further letting μ=0,ν=0 and ξ=0, in (3.15), we find directly the mean and variance of Linear Birth Process given by

E[N(t)]=n0eλtandV[N(t)]=n0eλt(eλt1).(5.2)

5.2. Pure Death Process

Substituting λ=0,ν=0, and ξ=0, in (3.14), we find the probability generating function of Pure Death Process given by

π(z,t)=[1+(z1)eμt]n0(5.3)
and from (3.15) the mean and variance of the Pure Death Process are given by
E[N(t)]=n0eμtandV[N(t)]=n0eμt(1eμt).(5.4)

5.3. Immigration Process

Letting λ0,μ0, and ξ=0, in (3.14), we find the probability generating function of Immigration Process given by

π(z,t)=zn0e(1z)νt(5.5)
and from (3.15), we find the mean and variance of the Immigration Process given by
E[N(t)]=n0+νtandV[N(t)]=νt.(5.6)

5.4. Emigration Process

Substituting λ0,μ0, and ν=0, in (3.14), we find the probability generating function of the Emigration Process given by

π(z,t)=zn0e(11z)ξt(5.7)
and from (3.15), we find the mean and variance of the Immigration Process given by
E[N(t)]=n0ξtandV[N(t)]=ξt.(5.8)

5.5. Linear Growth Process

Substituting ν=0 and ξ=0, in (3.14), the probability-generating function of Linear Growth Process is given by

π(z,t)=[λzμz1e(λμ)tμλzμz1e(λμ)tλ]n0,ifλμ(5.9)
and from (3.15) we find the mean and variance of Linear Growth Process given by
E[N(t)]=n0e(λμ)tandV[N(t)]=n0λ+μλμe(λμ)t[e(λμ)t1].(5.10)

5.6. Immigration-Emigration Process

Letting λ0 and μ0, in (3.14), we find the probability generating function of the Immigration-Emigration Process is given by

π(z,t)=zn0e(1z)(ξνz)zt,(5.11)
and from (3.15), we find the mean and variance of the Immigration-Emigration Process as
E[N(t)]=n0+(νξ)tandV[N(t)]=(ν+ξ)t.(5.12)

5.7. Birth-Immigration Process

Substituting μ=0 and ξ=0, in (3.14), we find the probability-generating function of the Birth-Immigration Process given by

π(z,t)=zn0en0λteνt[1z(1eλt)]n0+νλ,(5.13)
and from (3.15), we find the mean and variance of the Birth-Immigration Process as
E[N(t)]=νλ(eλt1)+n0eλtandV[N(t)]=(n0+νλ)eλt(eλt1).(5.14)

5.8. Birth-Emigration Process

Substituting μ=0 and ν=0, in (3.14), we find the probability generating function of the Birth-Emigration Process is given by

π(z,t)=ez1z(1eλt)ξλ[zeλt1z(1eλt)]n0,(5.15)
and from (3.15), we find the mean and variance of the Birth-Emigration Process as
E[N(t)]=ξλ(1eλt)+n0eλtandV[N(t)]=(n0eλt+ξλ)(eλt1).(5.16)

5.9. Death-Immigration Process

Substituting λ=0 and ξ=0, in (3.14), we find the probability-generating function of the Death-Immigration Process given by

π(z,t)=eνμ(1z)(1eμt)[1(1z)eμt]n0,(5.17)
and from (3.15), we get the mean and variance of the Death-Immigration Process as
E[N(t)]=νμ(1eμt)+n0eμtandV[N(t)]=(n0eμt+νμ)(1eμt).(5.18)

5.10. Death-Emigration Process

Substituting λ=0 and ν=0, in (3.14), we find the probability generating function of the Death-Emigration Process given by

π(z,t)=zξμ[1+(z1)eμt]n0+ξμ,(5.19)
and from (3.15), we find the mean and variance of the Death-Emigration Process as
E[N(t)]=ξμ(eμt1)+n0eμtandV[N(t)]=(n0+ξμ)eμt(1eμt).(5.20)

5.11. Linear Growth Process with Immigration

Substituting ξ=0, in (3.14), the probability generating function of Linear Growth Process with Immigration given by

π(z,t)=[(λzμz1λ)e(λμ)tλzμz1e(λμ)tλ]νλ×[λzμz1e(λμ)tμλzμz1e(λμ)tλ]n0(5.21)
and from (3.15), we find the mean and variance of Linear Growth Process with Immigration given by
E[N(t)]=νλμ(e(λμ)t1)+n0e(λμ)t   and=νμ(λμ)2(1e(λμ)t)+νλ(λμ)2e(λμ)t[e(λμ)t1]+n0(λ+μλμ)e(λμ)t[e(λμ)t1].(5.22)

5.12. Linear Growth Process with Emigration

Substituting ν=0, in (3.14), we find the probability generating function of Linear Growth Process with Emigration given by

π(z,t)=[λzμz1e(λμ)tμ(λzμz1μ)e(λμ)t]ξμ×[λzμz1e(λμ)tμλzμz1e(λμ)tλ]n0(5.23)
and from (3.15), we find the mean and variance of Linear Growth Process with Emigration given by
E[N(t)]=ξλμ(1e(λμ)t)+n0e(λμ)t   andV[N(t)]=λξ(λμ)2[e(λμ)t1]+μξ(λμ)2e(λμ)t[1e(λμ)t]+n0(λ+μλμ)e(λμ)t[e(λμ)t1].(5.24)

5.13. Birth Process with Immigration-Emigration

Substituting μ=0 in (3.14), we find the probability generating function of Linear Birth Process with Immigration-Emigration given by

π(z,t)=zn0en0λteνt[1z(1eλt)]n0+νλeξλz[1{1z(1eλt)}eλt].(5.25)

Further, substituting μ=0 in (3.15), we find the mean and variance of Birth Process with Immigration-Emigration given by

E[N(t)]=n0eλt+(νξλ)(eλt1),   andV[N(t)]=n0eλt(eλt1)+νλeλt(eλt1)+ξλ(eλt1).(5.26)

5.14. Death Process with Immigration-Emigration

Substituting λ=0 in (3.14), we find the probability generating function of Linear Death Process with Immigration-Emigration given by

π(z,t)=zξμeνμ(1z)(1eμt)[1+(z1)eμt]n0+ξμ.(5.27)

Further, substituting λ=0 in (3.15), we find the mean and variance of Death Process with Immigration-Emigration given by

E[N(t)]=n0eμt+(νξμ)(1eμt),   andV[N(t)]=n0eμt(1eμt)+νμ(1eμt)+ξμeμt(1eμt).(5.28)

6. IMPACT OF UNRESTRICTED CONSTANTS ON THE MEAN AND VARIANCE OF LINEAR GROWTH PROCESS WITH IMMIGRATION-EMIGRATION

For Linear Growth Process with Immigration-Emigration, the mean M(t) is

M(t)=E[N(t)]=(νξλμ)[e(λμ)t1]+n0e(λμ)t.(6.1)

When ν=ξ, the immigration and emigration have no control over the population and the mean M(t) of the population distribution tends to a desired value N (say) at a given time tj, we have

N=M(t)=E[N(t)]=n0e(λμ)tj,(6.2)
so that
Nn0e(λμ)tj(6.3)
is the difference between the means when control and no control respectively are applied.

From (6.1), we have

νξ=(λμ)[Nn0e(λμ)tj]e(λμ)tj1.(6.4)

For all tj>0, we have

λμe(λμ)tj1>0,λandμ.(6.5)

If N>n0e(λμ)tj, we have νξ>0, and hence ν>ξ, i.e., immigration rate is greater than emigration rate.

If the emigration rate reduces to zero, we have

ν=(λμ)[Nn0e(λμ)tj]e(λμ)tj1.(6.6)

If N<n0e(λμ)tj, we have νξ<0, and hence ν<ξ, i.e., immigration rate is less than emigration rate.

If the immigration rate reduces to zero, we have

ξ=(μλ)[Nn0e(λμ)tj]e(λμ)tj1.(6.7)

Once we have reached a desired mean N at time tj, we may need to maintain the mean at a constant level for all time in future, i.e., M(t)=N for all t>tj.

Hence, from (6.4), we have

νξ=(λμ)[NNe(λμ)(ttj)][e(λμ)(ttj)1],(6.8)
for all t>tj i.e.,
νξ=(μλ)N.(6.9)

If N>Ne(λμ)(ttj), we have νξ>0 and hence λ<μ.

Hence, if the immigration rate reduces to zero, we have

ν=(μλ)N.(6.10)

If N<Ne(λμ)(ttj), we have νξ<0 and hence λ>μ.

Hence, if the immigration rate reduces to zero, we have

ξ=(λμ)N.(6.11)

For linear growth process with immigration and emigration, the variance, V[N(t)] is

V[N(t)]=νμ(λμ)2[1e(λμ)t]+νλ(λμ)2e(λμ)t[e(λμ)t1]+ξλ(λμ)2[e(λμ)t1]+ξμ(λμ)2e(λμ)t[1e(λμ)t]+n0(λ+μλμ)e(λμ)t[e(λμ)t1].(6.12)

Since e(λμ)t1λμ>0 for all nonnegative λ and μ, all the terms in (6.12) are nonnegative except

νμ(λμ)2[1e(λμ)t]+μξ(λμ)2e(λμ)t[1e(λμ)t],(6.13)
which is nonpositive.

So, V[N(t)] is always greater than or equal to zero when ν=0 and ξ=0.

However, when ν is nonzero, the greatest value ν can have over the interval [t0,tj] is such as to make N=0 is

νmax=n0e(λμ)tj(μλe(λμ)tj1).(6.14)

Similarly, when ξ is nonzero, the greatest value ξ can have over the interval [t0,tj] is such as to make N=0 is

ξmax=n0e(λμ)tj(λμe(λμ)tj1).(6.15)

Consider the terms

νμ(λμ)2[1e(λμ)tj]+μξ(λμ)2e(λμ)tj[1e(λμ)tj]+n0(λ+μλμ)e(λμ)tj[e(λμ)tj1](6.16)
in V[N(t)].

Substituting the values of ν and ξ found by (6.14) and (6.15) respectively in (6.16), we left with

n0λe(λμ)tj[e(λμ)tj1λμ],(6.17)
which by (6.5) is nonnegative for all nonnegative λ and μ.

Hence V[N(t)] is positive for all t and for all possible constant control parameters.

7. A SIMULATION STUDY

In the following we performed a simulation study of Birth, Death, Immigration, and Emigration process and some of its particular cases. The behavior of these properties were also discussed except birth, death, and immigration processes because these are most common processes found in various text books.

7.1. Birth-Emigration Process

The process depends on the initial population size. Let the process starts with n0 individuals. From (5.1) and (5.7), probability generation function of Birth Process and Emigration Process is given by

Φ(z,t)=zn0e(11z)ξt×[zeλt1z(1eλt)]n0.

From (5.15), the p.g.f. of Birth-Emigration Process is

Ψ(z,t)=ez1z(1eλt)ξλ[zeλt1z(1eλt)]n0.

From Figure 1, for the particular values of ξ,λ and n0, we can get,

ϕ(z,t)ψ(z,t)
i.e., the p.g.f. of two independent process is less than or equal to the p.g.f. of Birth-Emigration Process. From (5.2) and (5.7), sum of expectation of Birth Process and Emigration Process is given by
M(t)=n0ξt+n0eλt.

Figure 1

P.G.F. of birth-emigrationprocess.

The expectation of Birth-Emigration Process from (5.16) is

M1(t)=ξλ(1eλt)+n0eλt

From Figure 2, we can get for the particular values of ξ, λ, t, and n0

M(t)=n0ξt+n0eλtξλ(1eλt)+n0eλt=M1(t).

Figure 2

Mean of birth-emigration process.

Hence we can conclude that the sum of expectation of Pure Birth Process and Emigration Process is larger than the expectation of Birth-Emigration Process, where birth and emigration occurs simultaneously.

7.2. Chance of Extinction

The expectation of Birth-Emigration Process is

M1(t)=ξλ(1eλt)+n0eλt.

As t

M1(t)=ξλ(1eλt)+n0eλt>1.

Hence we can conclude that, if λ>ξ, there is a chance of population explosion. If λ<ξ, the process will extinct after certain transition and the chance of extinction is nearly unity.

7.3. Death-Emigration Process

If the process starts with n0 individuals, from (5.3) and (5.7), the probability-generating function of Death Process and Emigration Process is

Φ1(z,t)=zn0e(11z)ξt×[1+(z1)eμt]n0

From (5.19), the p.g.f. of Death-Emigration Process is

Ψ1(z,t)=zξμ[1+(z1)eμt]n0+ξμ.

These are the p.g.f. within the range 0z1. Let z0<1 and 0t<. From Figure 3, we can see that, in the range 0zz0

Φ1(z,t)Ψ1(z,t)
and for the range z0<z1,
Φ1(z,t)Ψ1(z,t).

Figure 3

P.G.F. of death-emigration process.

From (5.4) and (5.8), the sum of expectation of Emigration Process and Pure Death Process is

M(t)=n0ξt+n0eμt
and from (5.18) the expectation Death-Emigration Process is
M1(t)=ξμ(eμt1)+n0eμt.

From Figure 4, we can get,

M(t)=n0ξt+n0eμt>ξμ(eμt1)+n0eμt=M1(t)
i.e., the sum of expectation of Emigration Process and Pure Death Process is larger than the expectation of Death-Emigration Process.

Figure 4

Mean of death-emigration process.

7.4. Chance of Extinction

The expectation of Death-Emigration Process is given by

M1(t)=ξμ(eμt1)+n0eμt.

As t

M1(t)=ξμ(eμt1)+n0eμt0
i.e., if the process starts with certain individuals, then it will extinct after some transition.

Hence the chance of extinction is nearly unity.

7.5. Birth-Death-Emigration Process (λμ)

Let the process starts with n0 individuals. From (5.7) and (5.9) the p.g.f. of Emigration Process and Birth-Death Process is

ϕ2(z,t)=zn0e(11z)ξt×[λzμz1e(λμ)tμλzμz1e(λμ)tλ]n0,ifλμ.

The p.g.f. of Birth-Death-Emigration (BDE) Process from (5.23), is

ψ2(z,t)=[λzμz1e(λμ)tμ(λzμz1μ)e(λμ)t]ξμ×[λzμz1e(λμ)tμλzμz1e(λμ)tλ]n0.

Let z0<1 and 0t<. From Figures 5 and 6), we can see that, in the range 0zz0

Φ2(z,t)>Ψ2(z,t).
and for the range z0<z1,
Φ2(z,t)Ψ2(z,t)

Figure 5

P.G.F. of Birth-Death-Emigration (BDE) process (λ > μ).

Figure 6

P.G.F. of Birth-Death-Emigration (BDE) process (λ < μ).

From (5.8) and (5.10), the expectation of Emigration Process and Birth-Death Process is given by

M(t)=n0ξt+n0e(λμ)t.

From (5.24), the expectation of BDE Process is given by

M1(t)=ξλμ(1e(λμ)t)+n0e(λμ)t.

From Figures 7 and 8, we can see if the process starts with n0 individuals, for particular values of λ, μ and ξ when (λ>μ),

M(t)=n0ξt+n0e(λμ)tξλμ(1e(λμ)t)+n0e(λμ)t=M1(t)
and when (λ<μ),
M(t)=n0ξt+n0e(λμ)t>ξλμ(1e(λμ)t)+n0e(λμ)t=M1(t)

Figure 7

Mean of Birth-Death-Emigration (BDE) process (λ > μ).

Figure 8

Mean of Birth-Death-Emigration (BDE) process (λ < μ).

7.5.1. Chance of extinction

The expectation of BDE process is given by

M1(t)=ξ(λμ)[1e(λμ)t]+n0e(λμ)t

As t

M1(t)n0e(λμ)tξ(λμ)[e(λμ)t](λ>μ)M1(t)=ξ(λμ)(λ<μ).

For (λ>μ), As t M1(t)>1 but for (λ<μ) As t M1(t)1

Hence we can conclude that for (λ>μ), there is a chance of population explosion but if the emigration rate will be larger than the birth rate and death rate, then the process will extinct after certain transition. For (λ<μ), the chance of extinction is nearly unity, i.e., the process will extinct after finite number of transition with probability one.

7.6. Birth-Death-Immigration-Emigration Process (λμ)

From (5.5), (5.7), and (5.9), the p.g.f. of Immigration Process, Emigration Process, and Birth-Death Process is given by

ϕ3(z,t)=zn0e(1z)νt×zn0e(11z)ξt×[λzμz1e(λμ)tμλzμz1e(λμ)tλ]n0

From (3.14), the p.g.f. Birth-Death-Immigration-Emigration (BDIE) Process is given by

ψ3(z,t)=[(λzμz1μ)e(λμ)t]ξμ[(λzμz1λ)e(λμ)t]νλ×[λzμz1e(λμ)tμ]n0+ξμ[λzμz1e(λμ)tλ]n0+νλ

From Figure 9, We can see

ϕ3(z,t)ψ3(z,t)

Figure 9

P.G.F. of Birth-Death-Immigration-Emigration (BDIE) process.

From (5.6), (5.8) and (5.10), the expectation of Immigration Process, Emigration Process, and Birth-Death Process is given by

M(t)=n0+νt+n0ξt+n0e(λμ)t
and the expectation of BDIE Process (λμ) is given by
M1(t)=(νξλμ)[e(λμ)t1]+n0e(λμ)t

Case-I [λ>μ,ξ>ν]

From Figure 10, we can see, if the process starts with n0 individuals for different values of ξ, ν, λ, and μ,

M(t)=n0+νt+n0ξt+n0e(λμ)t(νξλμ)[e(λμ)t1]+n0e(λμ)t=M1(t)
i.e., the average growth of BDIE Process is more than the sum of average growth of individual process.

Figure 10

Mean of BDIE process (λ > μ, ξ > ν).

Case-II [λ>μ,ξ<ν]

From Figure 11, we can see, if the process starts with n0 individuals for different values of ξ, ν, λ, and μ,

M(t)=n0+νt+n0ξt+n0e(λμ)t(νξλμ)[e(λμ)t1]+n0e(λμ)t=M1(t)
i.e., the average growth of BDIE Process is less than the sum of average growth of individual process.

Figure 11

Mean of Birth-Death-Immigration-Emigration (BDIE) process (λ > μ, ξ < ν).

Case-III [λ>μ,ξ=ν]

If the process starts with certain individual then for (λ>μ,ξ=ν), from Figure 12, we can see

M(t)=n0+νt+n0ξt+n0e(λμ)t=(νξλμ)[e(λμ)t1]+n0e(λμ)t=M1(t)

Figure 12

Mean of Birth-Death-Immigration-Emigration (BDIE) process (λ > μ; ξ = ν).

Since ξ=ν, the expectation of BDIE Process and sum of the expectation of individual process are equal to the expectation of Birth-Death process. As (λ>μ), if the process starts with certain individual, then it shows a positive growth and for large t there is a chance of population explosion.

Case-IV [λ<μ,ξ<ν]

If the process starts with n0 individuals, for different values of λ, μ, ν and ξ, from Figure 13, we can see

M(t)=n0+νt+n0ξt+n0e(λμ)t(νξλμ)[e(λμ)t1]+n0e(λμ)t=M1(t)

Figure 13

Mean of Birth-Death-Immigration-Emigration (BDIE) process (λ < μ; ξ < ν).

Since birth rate is less than the death rate and emigration rate is less than the immigration rate, the sum of expectation of individual process is lager than the expectation of BDIE process the initial population as well as every immigrant leave the system either due to death or emigration. As (λ<μ) here we can expect that the process will extinct after certain transition. Therefore there is a chance of extinction. But in case when we add the expectation of individual process, it grows indefinitely with additional immigrant. So there is chance of population explosion.

Case-V [λ<μ,ξ>ν]

From Figure 14, if the process starts with n0 individuals for different values of λ, μ, ξ, and ν we can see

M(t)=n0+νt+n0ξt+n0e(λμ)t(νξλμ)[e(λμ)t1]+n0e(λμ)t=M1(t)
i.e., the expectation of the process where birth, death, immigration, and emigration occurs at a time is larger than the sum of expectation of individual process. Since (λ<μ) and (ξ>ν), both are shows a negative growth. We can say each process will extinct after certain transition. Hence in both the cases the chance of extinction is nearly unity.

Figure 14

Mean of Birth-Death-Immigration-Emigration (BDIE) process (λ < μ; ξ > ν).

Case-VI [λ<μ,ξ=ν]

For (λ<μ) and (ξ=ν) Figure 15, indicates that the process will become the expectation of Birth-Death Process and

M(t)=n0+νt+n0ξt+n0e(λμ)t=(νξλμ)[e(λμ)t1]+n0e(λμ)t=M1(t)

Figure 15

Mean of Birth-Death-Immigration-Emigration (BDIE) process (λ < μ; ξ = ν).

Since in both the case it becomes the expectation of Birth-Death Process (λ<μ), it shows a negative growth and process will extinct after certain transition.

7.7. Chance of Extinction

The expectation of BDIE process is given by

M1(t)=ν(λμ)[e(λμ)t1]ξ(λμ)[e(λμ)t1]+n0e(λμ)t

As t

M1(t)={ν(λμ)[e(λμ)t1]ξ(λμ)[e(λμ)t1]+n0e(λμ)t,λ>μν(λμ)+ξ(λμ),ifλ<μ

Hence we can conclude that as t, the chance of extinction and chance of Population explosion is depends on (λ>μ) and (λ<μ). If (λ>μ), there is a chance of population explosion and if (λ<μ), the chance of extinction is nearly unity.

8. CONCLUSION

The process discussed in this paper is the linear growth process with immigration and emigration, which is a generalization to several population processes. The probability generating function of this process is also derived. Mean and variance are obtained by direct method and also from the probability generating function. Several other processes are also derived as particular cases with their probability generating functions. Also, the impact of unrestricted constants on the mean and variance of this process are discussed. The behavior of different processes are also analyzed using a simulation study. The chance of extinction of these processes were also studied.

CONFLICTS OF INTEREST

The authors declare that there is no conflict of interest.

AUTHORS' CONTRIBUTIONS

All the authors are contributed for the derivation of expressions in the research article and Upasana has contributed towards the simulation study in the article.

Funding Statement

The authors have not received ant fund for the preparation of this manuscript from any agency or organization.

ACKNOWLEDGMENTS

The authors are very much thankful to the reviewers for their valuable inputs to bring this research paper to this form.

Journal
Journal of Statistical Theory and Applications
Volume-Issue
20 - 1
Pages
111 - 131
Publication Date
2021/02/05
ISSN (Online)
2214-1766
ISSN (Print)
1538-7887
DOI
10.2991/jsta.d.210126.003How to use a DOI?
Copyright
© 2021 The Authors. Published by Atlantis Press B.V.
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

Cite this article

TY  - JOUR
AU  - Bijoy Kumar Pradhan
AU  - Priyaranjan Dash
AU  - Upasana
PY  - 2021
DA  - 2021/02/05
TI  - A Note on Generalization of the Simplest Time-Dependent Discrete Markov Process: Linear Growth Process With Immigration-Emigration
JO  - Journal of Statistical Theory and Applications
SP  - 111
EP  - 131
VL  - 20
IS  - 1
SN  - 2214-1766
UR  - https://doi.org/10.2991/jsta.d.210126.003
DO  - 10.2991/jsta.d.210126.003
ID  - Pradhan2021
ER  -