Next Article in Journal
Investigation of the Fork–Join System with Markovian Arrival Process Arrivals and Phase-Type Service Time Distribution Using Machine Learning Methods
Next Article in Special Issue
Estimation–Calibration of Continuous-Time Non-Homogeneous Markov Chains with Finite State Space
Previous Article in Journal
Assembly Function Recognition in Embedded Systems as an Optimization Problem
Previous Article in Special Issue
A Bilevel DEA Model for Efficiency Evaluation and Target Setting with Stochastic Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Strong Ergodicity in Nonhomogeneous Markov Systems with Chronological Order

by
P.-C.G. Vassiliou
Department of Statistical Sciences, University College London, Gower St., London WC1E 6BT, UK
Mathematics 2024, 12(5), 660; https://doi.org/10.3390/math12050660
Submission received: 11 January 2024 / Revised: 5 February 2024 / Accepted: 18 February 2024 / Published: 23 February 2024

Abstract

:
In the present, we study the problem of strong ergodicity in nonhomogeneous Markov systems. In the first basic theorem, we relax the fundamental assumption present in all studies of asymptotic behavior. That is, the assumption that the inherent inhomogeneous Markov chain converges to a homogeneous Markov chain with a regular transition probability matrix. In addition, we study the practically important problem of the rate of convergence to strong ergodicity for a nonhomogeneous Markov system (NHMS). In a second basic theorem, we provide conditions under which the rate of convergence to strong ergodicity is geometric. With these conditions, we in fact relax the basic assumption present in all previous studies, that is, that the inherent inhomogeneous Markov chain converges to a homogeneous Markov chain with a regular transition probability matrix geometrically fast. Finally, we provide an illustrative application from the area of manpower planning.

1. Introductory Notes

Consider a stochastic system that has a population of members categorized in different states. Three types of movements are possible in the system. Firstly, movements of members are probable among the states of the system; secondly, members are leaving the system from the various states; and thirdly, new members are entering the system to replace leavers and to expand the population. When the various movements of the system are modeled by a nonhomogeneous Markov chain, we call such a system a nonhomogeneous Markov system (NHMS).
An NHMS is actually a generalization of the classical Markov chain where we have one particle moving among the states without the possibility of leaving the system and probably being replaced by another with possibly different characteristics. Of great importance is the vector of absolute probabilities which consists of the probabilities of the particle to be in any state of the Markov chain. On the other hand, in an NHMS, we have a population of particles categorized according to their characteristics in the various states. Particles are leaving the population from all the states, and new particles are entering the population to replace them and to expand the population. Of great importance is the vector of the expected relative population structure. Hence, the problems to be solved are a lot harder, and new strategies and tools are used other than the simple Markov chain. The roots of the motive for the development of the theory of NHMS, which was first introduced in Vassiliou [1], could be summarized in the use of Markov models in manpower systems. This started with the work of Young and Almond [2], Young [3,4], and Bartholomew [5,6], and it was extended in the works of Young and Vassiliou [7], Vassiliou [8,9], and McClean [10,11]. In the book by Vassiliou [12], one can find the evolution of the theory of NHMS and the large diversity of its developments in various directions, that is, NHMS in discrete and continuous time, stochastic control in NHMS, Laws of Large Numbers for NHMS, Perturbations theory, NHMS in a stochastic environment, Markov systems, and others. In Section 5.4 of [12], there is a synopsis of real and potential applications of NHMS which illustrates the breadth of applications and some of the reasons why the entire theory is central to these processes. Work on manpower planning using results from NHMS and new areas has continued throughout the years up to nowadays, for example, Garg, et al. [13,14,15], Ugwugo and McClean [16], Vassiliadis [17,18], Georgiou, et al. [19], Guerry [20,21], Pollard [22], Esquivel, et al. [23,24,25].
In Section 2 of the present study, we define and describe the NHMS in discrete time and space in a compact but hopefully readable way for the reader who comes in contact for the first time with these processes. We also state the expected relative population structure in the various states as a function of the parameters of the population that could be estimated from the available data. In Section 3, we start with some basic definitions of concepts and mathematical tools, as well as useful known results that will be used in what follows. The novel part of this section is Theorem 1, where we study strong ergodicity for NHMS by relaxing the basic assumption present in all studies of strong ergodicity for NHMS. That is, we will not assume that the inherent inhomogeneous Markov chain converges as time goes to infinity to a homogeneous Markov chain with a regular (it consists of one communicating class of states, which is aperiodic) transition probability matrix or, equivalently, that the inhomogeneous Markov chain is strongly ergodic. In Section 3, we start with some basic definitions of concepts and mathematical tools, as well as useful known results that will be used in what follows. The novel part of this section is Theorem 4, where we prove under what conditions the rate of convergence of strong ergodicity in a NHMS is geometrically fast. This is an important question in NHMS due to its large practical value. In Theorem 4, we relax the basic assumption present in all studies of the rate of convergence to its asymptotic behavior for a NHMS. That is, we will not assume that the inherent inhomogeneous Markov chain converge as time goes to infinity to a homogeneous Markov chain with a regular transition probability matrix geometrically fast. Finally, in Section 5, we provide an illustrative application from the area of manpower planning.

2. The NHMS in Discrete Time and Space

Let a population consist of any kind of entities and let us denote by T ( t ) for t = 1 , 2 , , the total number of memberships at time t, that is, at the end of the interval t 1 , t , which are being held by its members. At every point of time that a member leaves the population, the membership is being transferred to a new member. For example, members could be patients in a hospital and memberships the beds they occupy. It is assumed that the total number of memberships are known or the sequence T t t = 0 is a realization of a known stochastic process depending on the application. The memberships are distributed in a finite number of states and let S = 1 , 2 , , k be the state space. Important aspect of an NHMS (see Vassiliou [12], Section 5.2) is the population structure, that is, the vector of random variables
N t = N 1 t , N 2 t , , N k t ,
where N i t is the number of memberships in state i at time t. Also, very important is the relative population structure, which is the vector of random variables q t = N t / T t . We denote by P t the transition probability matrix of the internal transitions of the members of the population during the interval t 1 , t , that is, the t-th interval. Also, we have probable leavers from the states in S in every time interval t and let us denote by
p k + 1 t = p 1 , k + 1 t , p 2 , k + 1 t , , p k , k + 1 t ,
where the state k + 1 represents the external environment. Finally, we have new entrants of memberships to the population in order to replace leavers and to expand the population. Let us collect the probabilities of allocation of the memberships to the various states in the t-th interval in the following stochastic vector
p 0 t = p 01 t , p 02 t , , p 0 k t ,
where the state 0 represents the new entrants with their memberships waiting to be allocated in the various states. Note that in what follows we assume that Δ T t = T t T t 1 0 .
The transition probability matrix of the memberships Q t in the t-th interval can be shown (Vassiliou [12], p. 193), and it is given by
Q t = P t + p k + 1 t p 0 t .
We call the inhomogeneous Markov chain defined uniquely by the sequence of transition probability matrices Q t t = 0 the embedded or inherent nonhomogeneous Markov chain of NHMS. A population or any physical phenomenon that could be modeled in the above described way is defined to be a nonhomogeneous Markov system.
Now, as previously, we define the relative population structure for a population which started at time s and is at time t to be q s , t = N s , t / T t , where N s , t is the population structure for the population. It could be proved that (Georgiou and Vassiliou [26] and Vassiliou [12], p. 195)
E [ q s , t ] = E [ q s , t 1 ] α t 1 Q t + b t 1 p 0 t ,
where s is the initial time and, therefore, q s is the initial relative population structure which is known and
α t 1 = T t 1 T t and b t 1 = T t T t 1 T t ,
and where q s , t is the relative population structure at time t for the system that started with initial relative population structure q s , which apparently is a random variable, and we denote by E [ q s , t ] its expected value.
From (2), recursively, we obtain
E [ q s , t ] = q s T s T t Q s , t + 1 T t τ = 1 t Δ T s + t p 0 s + t Q s + τ , t ,
where
Q s , t = Q s + 1 Q s + 2 Q t .

3. Strong Ergodicity in NHMS with Chronological Order

The asymptotic behavior of NHMS and of nonhomogeneous Markov chains, as well as of homogeneous Markov chains, has been one of the central problems for many years, as can be seen in Refs. [12,27,28,29,30,31]. The asymptotic behavior of NHMS started with Vassiliou [1,32], and an updated evolution of these theorems and their variants could be found in Vassiliou [12]. In the present section, we provide and prove a basic theorem for strong ergodicity in NHMS when the transition probabilities matrices of the inherent Markov chain Q t t = 0 are given in chronological order; that is, it is assumed that the time order of the elements of the sequence Q t t = 0 is given and will not be changed. In Theorem 1 we relax the basic assumption present in all studies of asymptotic behavior for NHMS. That is, we will not assume that the inherent inhomogeneous Markov chain converge as time goes to infinity to a homogeneous Markov chain with a regular (it consists of one communicating class of states which is aperiodic) transition probability matrix or, equivalently, that the inhomogeneous Markov chain is strongly ergodic. We start with some basic definitions and results, which will be useful in what follows.
For what follows, we assume a complete probability space Ω , F , P and consider an NHMS in discrete time and space. Therefore, we will not repeat this in every Definition, Lemma, Proposition, or Theorem.
Definition 1.
We say that the NHMS is strongly ergodic if and only if there exists a stochastic vector ψ such that
lim ν E q t , t + ν ψ = 0 . for t = 0 , 1 , 2 , ,
where . from now on is any vector norm, except if it is otherwise stated.
Definition 2.
If Q is any finite stochastic matrix with state space S = 1 , 2 , , k , then the Dobrushin ergodicity coefficient that is known is given by
τ 1 Q = 1 2 max i , r j = 1 k q i j q k j .
We now define another class of ergodicity coefficients which will be generated by different norms.
Definition 3
(Seneta [27]). Let the set
D n = x : x R n , x 0 , x 1 = 1 ,
and by d any metric on this set. Then, the quantity
τ P = sup x , y D n d xP , yP d x , y with x y ,
for any stochastic matrix, P is called a coefficient of ergodicity.
Remark 1.
All matrix norms on R n provide an appropriate metric on D n via d x , y = x y . Then for any stochastic matrix P we get
τ . P = sup x , y D n x y P x y with x y .
We call τ . P the coefficient of ergodicity induced by the norm . . It is proved that (8) can be written equivalently (Vassiliou [12], p. 118) as follows:
τ . P = sup z = 1 z 1 = 1 zP over all z R n .
When the L 1 norm is used, that is, for an n × n matrix A with elements from C , we obtain that
A 1 = i , j = 1 n a i j ,
then
τ 1 P = sup z 1 = 1 z 1 = 1 zP 1 over all z R n .
We now give the definitions of strong and weak ergodicity for a nonhomogeneous Markov chain:
Definition 4.
Consider an inhomogeneous Markov chain X t t = 0 in discrete time and space. We say that X t t = 0 with a sequence of transition probability matrices Q t t = 0 is strongly ergodic if there exists a stable stochastic matrix Q such that for every t
lim ν Q t , t + ν Q = 0 .
If the limit is zero uniformly in t we say that X t t = 0 or equivalently Q t t = 0 is uniformly strongly ergodic.
Definition 5.
Consider an inhomogeneous Markov chain X t t = 0 in discrete time and space. We say that X t t = 0 with sequence of transition probability matrices Q t t = 0 is weakly ergodic if for all states i , j , r , t
q i r t , t + ν q j r t , t + ν 0 as ν .
Note that in (13), lim ν q i j t , t + ν is not actually necessary to exist.
Remark 2.
Note that equivalently a non-homogeneous Markov chain is weakly ergodic if τ 1 Q t , t + ν < 1 for every t.
We now state the following Lemma the proof of which exists in Vassiliou ([12], p. 119).
Lemma 1.
The coefficient of ergodicity generated by any metric as in (7) or induced by any vector norm on R n has the following properties:
(1)
τ P 1 P 2 τ P 1 τ P 2 for any P 1 , P 2 stochastic matrices .
(2)
For any stochastic matrix P , τ P = 0 if and only if r a n k P = 1 .
The following Lemma proof, which can be found in Paz [33], is useful in what follows:
Lemma 2.
If P is a stochastic matrix and if R is any real matrix such that R 1 = 0 and in addition R < , then
RP R τ 1 P .
We now state and prove one of the basic theorems of the present paper.
Theorem 1.
Let there be the complete probability space Ω , F , P and consider an NHMS in discrete time and space. We assume that
lim t T t = T with Δ T t 0 .
Let Y t t = 0 be the inherent nonhomogeneous Markov chain of the movement of memberships. If Y t t = 0 is weakly ergodic with Q t t = 0 , the sequence of transition probabilities, and in addition, there exists a stochastic vector
ψ t = ψ 1 t , ψ 2 t , , ψ k t ,
which is the left eigenvector of Q t for t = 0 , 1 , 2 , , that is
ψ t = ψ t Q t for t = 0 , 1 , 2 , ,
and
t = 0 ψ t + 1 ψ t < ,
then the NHMS is strongly ergodic.
Proof. 
From (17) we obtain that there exists a stochastic vector ψ such that
lim t ψ t ψ = 0 .
From (16) and (18), we obtain that there exists a stochastic matrix Q such that
lim t Q t Q = 0 with Q 1 = 1 .
Denote by Ψ t to be the stable matrix with row the stochastic vector ψ t and Ψ the stable matrix with row ψ . Then, we have
Ψ t + 1 Ψ t = ψ t + 1 ψ t and Ψ t Ψ = ψ t ψ .
We now show that
lim ν Q t , t + ν Ψ = 0 .
We have that
Q t , t + ν Ψ = Q t , t + r Q t + r , t + ν Ψ t + r Q t + r , t + ν +
Ψ t + r Q t + r , t + ν Ψ t + ν 1 + Ψ t + ν 1 Ψ .
Taking norms on (22) we obtain that
Q t , t + ν Ψ Q t , t + r Q t + r , t + ν Ψ t + r Q t + r , t + ν
+ Ψ t + r Q t + r , t + ν Ψ t + ν 1 + Ψ t + ν 1 Ψ .
We now have
N 1 t , ν = Q t , t + r Q t + r , t + ν Ψ t + r Q t + r , t + ν
[ Q t , t + r Ψ t + r ] Q t + r , t + ν ( Due to Lemma 2 )
Q t , t + r Ψ t + r τ 1 Q t + r , t + ν 2 τ 1 Q t + r , t + ν ,
where τ 1 Q t + r , t + ν is less than one due to weak ergodicity of Y t t = 0 the inherent nonhomogeneous Markov chain. Also for fixed r, we can always choose ν such that
N 1 t , ν ε 3 , with ε > 0 a small number .
We now have that
Ψ t + r Q t + r , t + ν = ( due to Equation 16 ) =
Ψ t + r Ψ t + r + 1 Q t + r + 1 , t + ν +
+ Ψ t + r + 1 Q t + r + 1 , t + ν .
Now, similarly, we obtain
Ψ t + r + 1 Q t + r + 1 , t + ν =
Ψ t + r + 1 Ψ t + r + 2 Q t + r + 2 , t + ν
+ Ψ t + r + 2 Q t + r + 2 , t + ν ,
and using this equation recursively, we obtain
Ψ t + r Q t + r , t + ν = j = t + r + 1 t + ν 1 Ψ j 1 Ψ j Q j , t + ν
+ Ψ t + ν 1 Q t + ν 1 , t + ν .
From condition (16), we obtain that
Ψ t + ν 1 = Ψ t + ν 1 Q t + ν 1 ,
and
ψ t + ν 1 = ψ Ψ t + ν 1 .
Therefore, from (28)–(30), as well as, Lemma 2, we have that
N 2 t , ν = Ψ t + r Q t + r , t + ν Ψ Ψ t + ν 1
j = t + r + 1 t + ν 1 Ψ j 1 Ψ j τ 1 Q j , t + ν
j = t + r + 1 t + ν 1 ψ j 1 ψ j τ 1 Q j , t + ν
( τ 1 Q t + r , t + ν < 1 due to weak ergodicity of Y t t = 0 ; see also Remark 2 ) .
j = t + r + 1 t + ν 1 ψ j 1 ψ j < ,
due to condition (17). Since N 2 t , ν < for every r and ν and it is a sum of positive numbers, we have that its tail goes to zero. Hence, we could, for every ε > 0 , fix t + r such that N 2 t , ν < ε / 3 for ν 1 r and we can always take a ν large enough so that N 1 t , ν < ε / 3 and
N 3 t , ν = Ψ t + ν 1 Ψ ε 3 ,
therefore, we obtain that
lim ν Q t , t + ν Ψ = 0 .
Now, from (4), we have that
E q t , t + ν ψ =
q t T t T t + ν Q t , t + ν + 1 T t + ν τ = t t + ν Δ T τ Q τ , t + ν ψ
1 T t + ν q t T t Q t , t + ν + τ = t t + ν Δ T τ Q τ , t + ν T t + ν ψ
q t T t Q t , t + ν T t ψ +
τ = t t + ν Δ T τ p 0 τ Q τ , t + ν T t + ν T t ψ .
Now, we have
A t , t + ν = q t T t Q t , t + ν T t ψ
T t q t Q t , t + ν ψ .
By fixing ε > 0 , we can always find ν 0 such that
A t , t + ν ε 2 for ν ν 0 .
B t , t + ν = τ = t t + ν Δ T τ p 0 τ Q τ , t + ν T t + ν T t ψ =
= sin ce p 0 τ Ψ = p 0 Ψ = ψ =
τ = t t + ν Δ T τ p 0 τ Q τ , t + ν τ = t t + ν Δ T τ p 0 Ψ =
τ = t t + ν Δ T τ p 0 τ Q τ , t + ν p 0 Ψ
τ = t t + ν Δ T τ p 0 τ Q τ , t + ν p 0 Ψ
τ = t t + ν Δ T τ p 0 τ Q τ , t + ν Ψ + p 0 τ Ψ p 0 Ψ
τ = t t + ν Δ T τ Q τ , t + ν Ψ
Now, for ν > ν 0 from (37), we have that
B t , t + ν τ = t t + ν ν 0 Δ T τ Q τ , t + ν Ψ + τ = t + ν ν 0 t + ν Δ T τ Q τ , t + ν Ψ
τ = t t + ν ν 0 Δ T τ ε 2 + 2 T t + ν T t + ν ν 0 + 1 ,
where the second part of (38) for ν > > ν 0 is less than ε / 2 due to condition (14). The first part is also for, ν > > ν 0 , less than ε / 2 since Δ T t t 0 due to condition (14). From (34) and (36), and a ν > ν 0 large enough ν > > ν 0 , we easily see that
E q t , t + ν ψ ε for every ν > > ν 0 ,
hence, the NHMS is strongly ergodic. □

4. Rate of Convergence in NHMS with Chronological Order

An important question in nonhomogeneous Markov chains and NHMS, due to its very large practical value, is the rate of convergence to their asymptotic structure. In fact, it is important to find the necessary conditions under which the rate of convergence is geometric because then the value of the asymptotic result is greater. The roots of the study of finding conditions under which the rate of convergence is geometric for nonhomogeneous Markov chains are in Huang, et al. [34,35] and Seneta [31]. The study of the geometric rate of convergence in NHMS started in Vassiliou and Tsaklidis [36] and Georgiou and Vassiliou [26], and an updated evolution of these theorems and their variants could be found in Vassiliou [12]. The importance of answering this problem for practical purposes is apparent in Bartholomew [37] for the homogeneous Markov system, which is a very special case of an NHMS. In the present section, we provide and prove a basic theorem for the rate of convergence to strong ergodicity in NHMS when the transition probabilities matrices of the inherent Markov chain Q t t = 0 are given in chronological order. In Theorem 4, we relax the basic assumption present in all studies of the rate of convergence to its asymptotic behavior for an NHMS. That is, we will not assume that the inherent inhomogeneous Markov chain converges as time goes to infinity to a homogeneous Markov chain with a regular transition probability matrix geometrically fast. We start with some basic definitions and results which will be useful in what follows.
Definition 6.
We say that a sequence of matrices A n n = 0 converges with geometrical rate to a matrix A if there exists constants c > 0 and 0 < b < 1 such that
A n A c b n for n = 1 , 2 , .
Lemma 3
(Vassiliou and Tsaklidis [36]). The following statements are equivalent:
(i)
The sequence  Δ T t t = 0  converges to zero geometrically fast.
(ii)
The sequence  T t t = 0  converges to T geometrically fast.
Definition 7.
i A stochastic matrix P is called Markov if at least one column of P is entirely positive. Let M be the set of all Markov matrices. i i We say that the stochastic matrix P G 2 if a P G 1 the set of all regular matrices; b QP G 1 for any Q G 1 . i i i We say that the stochastic matrix P G 3 the set of all scrambling matrices if τ . 1 P < 1 .
Remark 3.
The distinction of the set G 2 from all stochastic regular matrices is due to the fact that the product of two regular matrices is not always regular. In addition, the product of two nonregular stochastic matrices could be regular. A practical way to check if a stochastic matrix of small dimension is scrambling is the following: given any two rows i , j , there is at least one column k such that p i k > 0 and p j k > 0 .
Definition 8.
The incidence matrix of a stochastic matrix P is a matrix where in the positions of positive elements we put the number 1. Therefore, two stochastic matrices P and Q of the same dimension have the same incidence matrix if they have the positive elements in the same positions. Then, we write P Q .
We now state some known Lemmas and Theorems, the proofs of which can be found in Vassiliou ([12], p. 143).
Theorem 2.
For all stochastic matrices, we have M G 3 G 2 G 1 .
Lemma 4.
If P t , t + ν G 1 with t 0 , n 1 , then P t , t + ν M for t + ν μ the number of distinct incidence matrices corresponding to G 1 with the same dimension as P t , t + ν .
Theorem 3.
Let there be a complete probability space Ω , F , P , and consider a nonhomogeneous Markov chain X t t = 0 in discrete time and space with a sequence of transition probabilities matrices P t t = 0 . If P t G 2 for every t = 1 , 2 , and
min i , j p i j t , 0 + γ > 0 ,
uniformly for all t 1 , then weak ergodicity obtains at a uniform geometric rate.
We now define the geometrically strongly ergodic NHMS, which is a central concept in the present section.
Definition 9.
Let there be a complete probability space Ω , F , P , and consider an NHMS in discrete time and space. We say that the NHMS is strongly ergodic if there exists a stochastic vector ψ and constants c > 0 and 0 < b < 1 such that
E q t , t + ν ψ c b ν for c > 0 and 0 < b < 1 .
We now state and prove the basic theorem of this section.
Theorem 4.
Let there be a complete probability space Ω , F , P , and consider an NHMS in discrete time and space. We assume that the total number of memberships is increasing ( Δ T t 0 ), and it converges geometrically fast to T. That is,
lim t T t = T geometrically fast with Δ T t 0 .
Let Y t t = 0 be the inherent nonhomogeneous Markov chain of the movement of memberships. If Y t t = 0 is weakly ergodic with Q t t = 0 the sequence of transition probabilities to be such that
i Q t G 2 for every t = 1 , 2 , ,
i i min i , j q i j t , 0 + γ > 0 ,
and if in addition there exists a stochastic vector
ψ t = ψ 1 t , ψ 2 t , , ψ k t ,
which is the left eigenvector for t = 1 , 2 , of Q t , that is,
ψ t = ψ t Q t for every t ,
and in addition
lim ν ψ t + ν ψ t + ν 1 = 0 geometrically fast ,
then the NHMS is geometrically strongly ergodic.
Proof. 
In order to prove that the NHMS is geometrically strongly ergodic, we must show that the expected relative structure satisfies Definition 9.
From (46), we have that there exists a vector ψ and constants c 1 > 0 and 0 < b 1 < 1 such that
ψ t + ν ψ c 1 b 1 t + ν .
We define Ψ t and Ψ as the stable matrices with rows ψ t and ψ , respectively.
We now show that the inherent nonhomogeneous Markov chain Y t t = 0 is geometrically strongly ergodic. That is, we need to show that there exists constants c 2 > 0 and 0 < b 2 < 1 such that
Q t , t + ν Ψ c 2 b 2 ν for every t .
Let us denote by μ the number of distinct incidence matrices corresponding to G 1 with the same dimension as P t , t + ν . Then, for ν > μ , we have that
Q t , t + ν Ψ Q t , t + μ Q t + μ , t + ν Ψ t + μ Q t + μ , t + ν
+ Ψ t + μ Q t + μ , t + ν Ψ t + ν 1 + Ψ t + ν 1 Ψ .
Now, we have that
D 1 t , ν , μ = Q t , t + μ Q t + μ , t + ν Ψ t + μ Q t + μ , t + ν
Q t , t + μ Ψ t + μ Q t + μ , t + ν
Q t , t + μ Ψ t + μ τ 1 Q t + μ , t + ν
2 τ 1 Q t + μ , t + ν .
For arbitrary but fixed t and ν with m μ ν with m, the largest such integer from (50), we have that
D 1 t , ν , μ 2 τ 1 Q t + μ , t + 2 μ τ 1 Q t + μ , t + 2 μ
τ 1 Q t + m 1 μ , t + m μ τ 1 Q t + m μ , t + ν .
Since Q t G 2 and by Lemma 4 we have that τ 1 Q t + i 1 μ , t + i μ M and
τ 1 Q t + i 1 μ , t + i μ 1 γ μ with 0 < γ < 1 for i = 2 , 3 , , m .
From the weak ergodicity of the inherent nonhomogeneous Markov chain Y t t = 0 , we have that
τ 1 Q t + m μ , t + ν 1 .
From (51)–(53), we arrive at
D 1 t , ν , μ 1 γ μ m ,
which as ν goes to zero at a uniform geometric rate.
Following the steps of arriving at relation (31), we straightforwardly obtain that
D 2 t , ν , μ = Ψ t + μ Q t + μ , t + ν Ψ t + ν 1
Ψ t + μ Q t + μ , t + ν Ψ Ψ t + ν 1
j = t + μ + 1 t + ν 1 ψ j 1 ψ j τ 1 Q j , t + ν .
We now have the largest integer for m, such that m μ ν
τ 1 Q t + μ , t + ν = τ 1 ( Q t + μ , t + 2 μ Q t + 2 μ , t + 3 μ
Q t + m 1 μ , t + m μ Q t + m μ , t + ν μ )
= ( by Lemma 1 )
τ 1 Q t + μ , t + 2 μ τ 1 Q t + 2 μ , t + 3 μ
τ 1 t + m 1 μ , t + m μ τ 1 t + m μ , t + ν
using 52
1 γ μ 1 γ μ 1 γ μ = 1 γ μ m .
From (56), we have that
τ 1 Q t + μ , t + ν ν 0 uniformly geometrically fast .
From (55), we have that
D 2 t , ν , μ ψ t + μ 1 ψ t + μ 1 γ μ m +
ψ t + μ ψ t + μ + 1 1 γ μ m 1 + ψ t + μ ψ t + μ + 1 1 γ μ m 1
+ + ψ t + μ m 1 1 ψ t + μ m 1 1 γ μ
+ ψ t + μ m 1 ψ t + μ m .
Now, we have that
1 γ μ m m 0 geometrically fast ,
since 1 γ μ < 1 . Also, from condition (46), we have that
ψ t + μ m 1 ψ t + μ m m 0 geometrically fast ,
Therefore,
D 2 t , ν , μ ν 0 geometrically fast .
Now, from (49), (54) and (59), and condition (46), we easily have that
Q t , t + ν Ψ ν 0 geometrically fast ;
hence, there exist constants c > 0 and 0 < b < 1 such that
Q t , t + ν Ψ c b ν for every t , and c > 0 , 0 < b < 1 .
What remains, according to Definition 9, is to show that the expected relative population structure E q t , t + ν converges for every t to the vector ψ as ν goes to infinity geometrically fast.
From (35), we have that
A t , t + ν = q t T t Q t , t + ν T t ψ
T t q t Q t , t + ν ψ .
From (61), we have that
A t , t + ν T t q t Q t , t + ν ψ c 1 b ν ,
with c 1 > 0 and 0 < b < 1 .
From Lemma 3, we have that since T t converges geometrically fast to T, then Δ T t converges geometrically fast to zero. Hence,
Δ T t c 2 b 2 t with c 2 > 0 and 0 < b 2 < 1 .
Now, from (37), we have that
B t , t + ν τ = t t + ν Δ T τ Q τ , t + ν Ψ .
With no loss of generality, we may assume that b > b 2 , and then, from (63) and (64), we have that
B t , t + ν c 1 c 2 τ = t t + ν b 2 b τ b t + ν = c 3 b t + ν τ = t t + ν b 2 b τ
c 3 b t + ν 1 b 2 b ν 1 b 2 b 1 c 4 b 3 ν .
From (63) and (66), we easily arrive at the proof of the theorem. □

5. An Illustrative Application

We will illustrate the results in the previous section with an example from a population of manpower. To possibly better visualize, the reader may have in mind a University system with three grades. That is, grade one is those with the level of Professors, grade two belongs to the Associate Professors, and finally, in grade three, there are the Assistant Professors. The University has a plan for funding T t for t = 0 , 1 , 2 , memberships for the next few years. When a member of staff is leaving, their membership remains with the University, that is, the funding of their position is not lost but remains, and the University could go on to appoint someone at any grade. The external environment to which the leavers go and their membership is retained by the University and it is the population of members of academic staff from Universities almost all over the world, as practices have shown. Hence, it is from this external environment that the new members will obtain the memberships available from the organization of the University.
A fundamental question for the practitioner is the estimation of the transition probabilities from the historical records. The way to estimate the transition probabilities of the memberships is a small extension of the way it is performed in manpower planning, which is well documented in Bartholomew [38]. Another quite similar problem is that of the competing risk model in the medical literature, as it was presented by Kalbfleisch and Prentice ([39], Chapter 8), Lee [40], and Cox and Oakes [41]. Also, similar problems exist in the study of reliability models and various actuarial studies as discussed in Elandt–Johnson and Johnson ([42], Chapter 7).
Now define by
N i j t : the number of memberships moving from grade i to grade j in the t-th interval.
N i t 1 : the number of memberships in state i at the beginning of the t-th interval.
Assume that the number of years available in the historical data of the University is n. Then, the maximum likelihood estimate of the probability q i j t is the following:
q ^ i j t = N i j t N i t 1 for any t .
It is an apparent advantage that the probabilities q ^ i j t are separately estimated for every i , j . In this way, the number of years of historical records necessary are significantly reduced. At this point, it is useful to test the hypothesis that the probabilities q ^ i j t are indeed functions of time. That is,
H 0 : q ^ i j t = q ^ i j for every t .
Considering the flow of memberships which move from grade i to grade j as a multinomial random variable, then (see Andersen and Goodman [43]) hypothesis (68) is tested by the statistic:
χ 2 i , j = t = 1 n N i t 1 q ^ i j t q ^ i j q ^ i j ,
where
q ^ i j = t = 0 n N i j t t = 1 n N i t 1 ,
is the maximum likelihood estimate under the null hypothesis and is chi-square distributed with n 1 degrees of freedom.
Now, let that the χ 2 i , j showed that the probabilities q ^ i j t are functions of time. Then, there is a need to predict their values as functions of time. For a specific pair i , j , let that
x 1 i j t , x 2 i j t , , x m i j t ,
are probable covariates for the specific application. Then, logistic stepwise regression is an appropriate model for these probabilities. Let us define by
log i t x = log x 1 x ,
then we obtain
log i t q ^ i j t = a 0 + a 1 x 1 i j t + a 2 x 2 i j t + + a m x m i j t .
Now, it is obvious that stepwise regression will show what are the important covariates to predict q ^ i j t (see also Vassiliou [44]).
Let that following the above-described steps, we ended with the matrix:
Q t = 0.2 + 1 4 + t 3 0.8 1 4 + t 3 0 0.3 0.5 1 8 + t 4 0.2 + 1 8 + t 4 0 0.2 + 2 10 + t 8 0.8 2 10 + t 8 for t = 1 , 2 , .
The total population of memberships was planned according to the following sequence
T 0 = 400 , T 1 = 430 , T 2 = 450 ,
T 3 = 475 , T 4 = 500 and T t = 500 for t = 5 , 6 , ,
that is, the total number of memberships converge geometrically fast with Δ T t 0 , satisfying the condition (41) in Theorem 4.
It is not difficult to check that the sequence (74) satisfies condition (42), that is, Q t G 2 for t = 1 , 2 , . Also, it satisfies condition (43) since
min i , j q i j t , 0 + γ = 0.2 > 0 .
For the sequence of transition probability matrices for the memberships (74), we find for the condition (45) that the vectors ψ t that satisfy it were the following:
For t = 1 ψ 1 = 0.216 , 0.432 , 0.352 ,
for t = 2 ψ 2 = 0.166 , 0.390 , 0.444 ,
for t = 3 ψ 3 = 0.164 , 0.416 , 0.420 ,
for t = 4 ψ 4 = 0.164 , 0.424 , 0.412 ,
for t = 5 ψ 5 = 0.163 , 0.426 , 0.411 ,
for t = 6 ψ 6 = 0.163 , 0.427 , 0.410 ,
for t = 7 ψ 7 = 0.163 , 0.427 , 0.410 .
We observe that already for t = 6 in (82) and (83) we have convergence of ψ t which satisfies condition (45) of Theorem 4, that is the convergence is geometrically fast. Hence, we conclude that
ψ = 0.16 , 0.43 , 0.41 .
Calculating using the transition probability matrices given in (74) the matrix product Q t , t + ν we find that for ν 9 it converges to ψ = 0.16 , 0.43 , 0.41 . Hence, we conclude that Q t , t + ν as ν converges geometrically fast as was expected from Theorem 4 relation (48).
Now, given the convergence of Q t , t + ν as ν it is straight forward to find from Equation (2) applying it recursively or equivalently from Theorem 1 that E q s , t converges geometrically fast (in fact in 9 time steps) to
ψ = 0.16 , 0.43 , 0.41 .
Hence, the NHMS is geometrically strongly ergodic.

6. Conclusions and Further Research

Two fundamental theorems have been founded that relax previous assumptions and provide conditions for ergodicity and for the convergence rate of the relative population structure in an NHMS. More specifically, in Theorem 1 the strong ergodicity of an NHMS is studied without assuming the convergence of the inherent inhomogeneous Markov chain to infinity to a homogeneous Markov chain with a regular transition probability matrix. In Theorem 4, it is proved under what conditions the rate of convergence of strong ergodicity in a NHMS is geometrically fast. This is done by departing from the basic assumption in which the inherent inhomogeneous Markov chain converges after a large time to a homogeneous Markov chain with a regular transition probability matrix. The proved theorems are expanding the understanding of the dynamics and behavior of an NHMS. The paper concludes with an illustrative application from the field of Manpower Planning, showcasing the vital practical relevance of the discussed concepts. Further research paths may include the relaxation of relative assumptions in the many variant models of the NHMS in diverse populations. Of particular interest may be the theorems of Laws of Large numbers in an NHMS.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The author declare no conflict of interest.

References

  1. Vassiliou, P.-C.G. Asymptotic behavior of Markov systems. J. Appl. Prob. 1982, 19, 433–438. [Google Scholar] [CrossRef]
  2. Young, A.; Almond, G. Predicting distributions of staff. Comput. J. 1961, 3, 144–153. [Google Scholar] [CrossRef]
  3. Young, A. Models for planning recruitment and promotion of staff. Brit. J. Indust. Rel. 1965, 3, 301–310. [Google Scholar] [CrossRef]
  4. Young, A. Demographic and ecological models for manpower planning. In Aspects of Manpower Planning; Bartholomew, D.J., Morris, B.R., Eds.; English University Press: London, UK, 1971. [Google Scholar]
  5. Bartholomew, D.J. A multistage renewal processes. J. R. Stat. Soc. B 1963, 25, 150–168. [Google Scholar]
  6. Bartholomew, D.J. Stochastic Models for Social Processes, 1st ed.; Wiley: New York, NY, USA, 1967. [Google Scholar]
  7. Young, A.; Vassiliou, P.-C.G. A non-linear model on the promotion of staff. J. R. Stat. Soc. A 1974, 138, 584–595. [Google Scholar] [CrossRef]
  8. Vassiliou, P.-C.G. A Markov model for wastage in manpower systems. Oper. Res. Quart. 1976, 27, 57–70. [Google Scholar] [CrossRef]
  9. Vassiliou, P.-C.G. A high order non-linear Markovian model for promotion in manpower systems. J. R. Stat. Soc. A 1978, 141, 86–94. [Google Scholar] [CrossRef]
  10. McClean, S.I. A continuous-time population model with Poisson recruitment. J. Appl. Prob. 1976, 13, 348–354. [Google Scholar] [CrossRef]
  11. McClean, S.I. Continuous-time stochastic models for multigrade population. J. Appl. Prob. 1978, 15, 26–37. [Google Scholar] [CrossRef]
  12. Vassiliou, P.-C.G. Non-Homogeneous Markov Chains and Systems, Theory and Applications; Chapman and Hall: London, UK; CRC Press: Boca Raton, FL, USA, 2023. [Google Scholar]
  13. Garg, L.; McClean, S.I.; Meenan, B.; Millard, P. A non-homogeneous discrete time Markov model for admission scheduling and resource planning in a cost capacity constaint healthcare system. Health Care Manag. Sci. 2010, 13, 155–169. [Google Scholar] [CrossRef]
  14. Garg, L.; McClean, S.I.; Meenan, B.; Millard, P. Non-homogeneous Markov models for sequential pattern mining of healthcare data. Ima J. Manag. Math. 2009, 20, 327–344. [Google Scholar] [CrossRef]
  15. Garg, L.; McClean, S.I.; Meenan, B.; Millard, P. Phase-Type survival trees and mixed distribution syrvival trees for clustering patient’s hospital length of stay. Informatika 2011, 22, 57–72. [Google Scholar]
  16. Ugwuogo, F.I.; McClean, S.I. Modelling heterogeneity in manpower systems: A review. Appl. Stoch. Models Bus. Ind. 2000, 2, 99–110. [Google Scholar] [CrossRef]
  17. Vassiliadis, G. Transient analysis of the M/M/k/N/N queue using a continuous time homogeneous Markov chain system with finite state capacity. Commun. Stat. Theory Methods 2014, 43, 1548–1562. [Google Scholar] [CrossRef]
  18. Vassiliadis, G. Transient analysis of a finite source discrete-time queueing system using homogeneous Markov system with state size capacities. Commun. Stat. Theory Methods 2014, 45, 1403–1423. [Google Scholar] [CrossRef]
  19. Georgiou, A.C.; Thanassoulis, E.; Papdopoulou, A. Using data envelopment analysis in markovian decision making. Eur. J. Oper. Res. 2022, 298, 276–292. [Google Scholar] [CrossRef]
  20. Guerry, M.A. On the evolution of stock vectors in a deterministic integer-valued Markov system. Linear Algebra Its Appl. 2008, 429, 1944–1953. [Google Scholar] [CrossRef]
  21. Guerry, M.A. Some results on the embeddable problem for discrete time Markov models in manpower planning. Commun. Stat. Theory Methods 2014, 43, 1575–1584. [Google Scholar] [CrossRef]
  22. Pollard, B.S. Open Markov processes: A compositional perspective on a non-equilibrium steady state in biology. Entropy 2016, 18, 140. [Google Scholar] [CrossRef]
  23. Esquivel, M.L.; Fernandes, J.M.; Guerriero, G.R. On the evolution and asymptotic analysis of open Markov populations: Application to consumption credit. Stoch. Model. 2014, 30, 365–389. [Google Scholar] [CrossRef]
  24. Esquivel, M.L.; Guerriero, G.R.; Fernandes, J.M. Open Markov chain scheme models fed by second order stationary and non stationary processes. Revstat-Stat. J. 2017, 15, 277. [Google Scholar]
  25. Esquivel, M.L.; Krasil, N.P.; Guerriero, G.R. Open type population models: From discrete to continuous time. Mathematics 2021, 9, 1496. [Google Scholar] [CrossRef]
  26. Georgiou, A.C.; Vassiliou, P.-C.G. Periodicity of asymptotically attainable stuctures in non-homogeneous Markov systems. Linear Algebra Its Appl. 1992, 176, 137–174. [Google Scholar] [CrossRef]
  27. Seneta, E. Non-Negative Matrices and Markov Chains, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 1981. [Google Scholar]
  28. Brémaud, P. Markov Chains. In Gibbs Fields, Monte Carlo Simulation and Queues, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  29. Tweedie, R.L. Sufficient conditions for ergodicity and recurrence of Markov chains on a general state space. Stoch. Proc. Appl. 1975, 3, 385–403. [Google Scholar] [CrossRef]
  30. Mitrophanov, A.Y. Sensitivity and convergence of uniformly ergodic Markov chains. J. Appl. Prob. 2005, 42, 1003–1014. [Google Scholar] [CrossRef]
  31. Seneta, E. Inhomogeneous Markov chains and Ergodicity Coefficients: John Hajnal (1924–2008). Commun. Stat. Theory Methods 2014, 43, 1575–1584. [Google Scholar] [CrossRef]
  32. Vassiliou, P.-C.G. On the limiting behavior of a non-homogeneous Markov model in manpower systems. Biometrika 1981, 68, 557–561. [Google Scholar]
  33. Paz, A. Introduction to Probabilistic Automata; Academic Press: Cambridge, MA, USA, 1971. [Google Scholar]
  34. Huang, C.; Isaacson, D.; Vinograde, B. The rate of convergence of certain non-homogeneous Markov chains. Zeitsh. Wahrsch. Geb. 1976, 35, 141–146. [Google Scholar] [CrossRef]
  35. Huang, C.; Isaacson, D. Ergodicity using mean visit times. J. London Math. Soc. 1976, 14, 570–576. [Google Scholar] [CrossRef]
  36. Vassiliou, P.-C.G.; Tsaklidis, G. The rate of convergence of the vector of variances and covariances in non-homogeneous Markov systems. J. Appl. Prob. 1989, 27, 776–783. [Google Scholar] [CrossRef]
  37. Bartholomew, D.J. Stochastic Models for Social Processes, 3rd ed.; Wiley: New York, NY, USA, 1981. [Google Scholar]
  38. Bartholomew, D.J.; Forbes, A.F.; McClean, S. Statistical Techniques for Manpower Planning; John Wiley: Chichester, UK, 1991. [Google Scholar]
  39. Kalbfleisch, J.D.; Prentice, R.L. The Statistical Analysis of Failure Time Data, 2nd ed.; John Wiley: New York, NY, USA, 2002. [Google Scholar]
  40. Lee, E.T. Statistical Methods for Survival Data Analysis, 2nd ed.; John Wiley: New York, NY, USA, 1992. [Google Scholar]
  41. Cox, D.R.; Oakes, D. Analysis of Survival Data; Chapman and Hall: London, UK, 1984. [Google Scholar]
  42. Elandt-Johnson, R.C.; Johnson, N.L. Survival Models and Data Analysis; John Wiley: New York, NY, USA, 1980. [Google Scholar]
  43. Anderson, T.W.; Goodman, L.A. Statistical Inference about Markov chains. Ann. Math. Statist. 1957, 28, 89–110. [Google Scholar] [CrossRef]
  44. Vasileiou, A.; Vassiliou, P.-C.G. An inhomogeneous semi-Markov model for the term structure of credit risk spreads. Adv. Appl. Prob. 2006, 38, 171–198. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vassiliou, P.-C.G. Strong Ergodicity in Nonhomogeneous Markov Systems with Chronological Order. Mathematics 2024, 12, 660. https://doi.org/10.3390/math12050660

AMA Style

Vassiliou P-CG. Strong Ergodicity in Nonhomogeneous Markov Systems with Chronological Order. Mathematics. 2024; 12(5):660. https://doi.org/10.3390/math12050660

Chicago/Turabian Style

Vassiliou, P.-C.G. 2024. "Strong Ergodicity in Nonhomogeneous Markov Systems with Chronological Order" Mathematics 12, no. 5: 660. https://doi.org/10.3390/math12050660

APA Style

Vassiliou, P. -C. G. (2024). Strong Ergodicity in Nonhomogeneous Markov Systems with Chronological Order. Mathematics, 12(5), 660. https://doi.org/10.3390/math12050660

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop