3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Section 3.1: Definitions and Basic Concepts STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 2 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Stochastic Process Definition: {X (t), t 2 T } is called a stochastic process if X (t) is a rv (or possibly a random vector) for any given t 2 T . T is referred to as the index set and is often interpreted in the context of time. As such, X (t) is often called the state of the process at time t. We note that: Index set T 8><>: can be a continuum of values such as T = {t : t 0}, can be a set of discrete points such as T = {t0, t1, t2, . . .}. Since there is a one-to-one correspondence between the sets T = {t0, t1, t2, . . .} and N = {0, 1, 2, . . .}, we will use T = N as the general index set for a discrete-time stochastic process (unless otherwise stated). In other words, {X (n), n 2 N} or {Xn, n 2 N} will represent a general discrete-time stochastic process. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 3 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Discrete-time Stochastic Process Some examples of discrete-time stochastic processes {Xn, n 2 N} might include: (1) Xn represents the outcome of the nth toss of a die, (2) Xn represents the price of a stock at the end of day n trading, (3) Xn represents the maximum temperature in Waterloo during the nth month, (4) Xn represents the number of goals scored in game n by the varsity hockey team, (5) Xn represents the number of STAT 333 students in class for the nth lecture. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 4 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Discrete-time Markov Chain Definition: A stochastic process {Xn, n 2 N} is said to be a discrete-time Markov chain (DTMC) if the following two conditions hold true: (1) For n 2 N, Xn is a discrete rv (i.e., the state space S of Xn is of discrete type). (2) For n 2 N and all states x0, x1, . . . , xn+1 2 S, the Markov property must hold: P(Xn+1 = xn+1|Xn = xn,Xn 1 = xn 1, . . . ,X1 = x1,X0 = x0) = P(Xn+1 = xn+1|Xn = xn). In mathematical terms, this property states that the conditional distribution of any future state Xn+1 given the past states X0,X1, . . . ,Xn 1 and the present state Xn is independent of the past states. In a more informal way , the Markov property tells us, for a random process, that if we know the value taken by the process at a given time, we will not get any additional information about the future behaviour of the process by gathering more knowledge about the past. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 5 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Discrete-time Markov Chain Remarks: (1) Unless otherwise stated, the state space S of a DTMC {Xn, n 2 N} will be assumed to be N. (2) In general, the sequence of rvs {Xn}1n=0 are neither independent nor identically distributed. (3) The Markov property does not require “full” information on the past to ensure independence. For example, consider the following conditional probability: P(Xn+1 = xn+1|Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) = P(Xn+1 = xn+1,Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) P(Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) , which is “missing” the information for Xn 1. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 6 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Discrete-time Markov Chain However, note that P(Xn+1 = xn+1,Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) = 1X xn 1=0 P(Xn+1 = xn+1,Xn = xn,Xn 1 = xn 1,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) = 1X xn 1=0 P(Xn+1 = xn+1|Xn = xn,Xn 1 = xn 1,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) ⇥ P(Xn = xn,Xn 1 = xn 1,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) = P(Xn+1 = xn+1|Xn = xn) ⇥ 1X xn 1=0 P(Xn = xn,Xn 1 = xn 1,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) by Markov property = P(Xn+1 = xn+1|Xn = xn)P(Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0). STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 7 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Discrete-time Markov Chain We have: P(Xn+1 = xn+1|Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) = P(Xn+1 = xn+1,Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) P(Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) and P(Xn+1 = xn+1,Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) = P(Xn+1 = xn+1|Xn = xn)P(Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0). Substituting this latter expression into the numerator of the top equation yields P(Xn+1 = xn+1|Xn = xn,Xn 2 = xn 2, . . . ,X1 = x1,X0 = x0) = P(Xn+1 = xn+1|Xn = xn). It is straightforward to extend the above result to any number of past time points from 0, 1, . . . , n 1 with “missing” information. This is the essence of the Markov property. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 8 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts One-step Transition Probability Matrix Definition: For any pair of states i and j , the transition probability from state i at time n to state j at time n + 1 is given by Pn,i ,j = P(Xn+1 = j |Xn = i), n 2 N. Let Pn be the associated matrix containing all these transition probabilities, referred to as the one-step transition probability matrix (TPM) from time n to time n + 1. It looks like Pn = [Pn,i ,j ] = 26666664 0 1 2 · · · j · · · 0 Pn,0,0 Pn,0,1 Pn,0,2 · · · Pn,0,j · · · 1 Pn,1,0 Pn,1,1 Pn,1,2 · · · Pn,1,j · · · ... ... ... ... ... i Pn,i ,0 Pn,i ,1 Pn,i ,2 · · · Pn,i ,j · · · ... ... ... ... ... 37777775, where, for convenience, the states of the DTMC are labelled along the margins of the matrix. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 9 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts One-step Transition Probability Matrix For each pair of states i and j , if Pn,i ,j = Pi ,j 8 n 2 N, then we say that the DTMC is stationary or homogeneous. In this situation, the one-step TPM becomes: P = [Pi ,j ] = 26666664 0 1 2 · · · j · · · 0 P0,0 P0,1 P0,2 · · · P0,j · · · 1 P1,0 P1,1 P1,2 · · · P1,j · · · ... ... ... ... ... i Pi ,0 Pi ,1 Pi ,2 · · · Pi ,j · · · ... ... ... ... ... 37777775. Remark: In STAT 333, we only consider stationary DTMCs. Moreover, from the construction of the TPM P , it is clear that Pi ,j 0 8 i , j 2 N and P1 j=0 Pi ,j = 1 8 i 2 N (i.e., each row sum of P must be 1). Such a matrix whose elements are non-negative and whose row sums are equal to 1 is said to be stochastic. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 10 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts One-step Transition Probability Matrix Example 3.1. On a given day, the weather is either clear, overcast, or raining. If the weather is clear today, then it will be clear, overcast, or raining tomorrow with respective probabilities 0.6, 0.3, and 0.1. If the weather is overcast today, then it will be clear, overcast, or raining tomorrow with respective probabilities 0.2, 0.5, and 0.3. If the weather is raining today, then it will be clear, overcast, or raining tomorrow with respective probabilities 0.4, 0.2, and 0.4. Construct the underlying DTMC and determine its TPM. Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 11 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts n-step Transition Probability Matrix Definition: For any pair of states i and j , the n-step transition probability is given by P(n)i ,j = P(Xm+n = j |Xm = i), m, n 2 N. Due to the stationary assumption, this quantity is actually independent of m (which is why we do not include m in its notation). Thus, P(n)i ,j = P(Xn = j |X0 = i), n 2 N. Furthermore, it is evident that P(0)i ,j = P(X0 = j |X0 = i) = ( 0 , if i 6= j , 1 , if i = j . In a similar manner, let P(n) = h P(n)i ,j i represent the associated n-step TPM. Clearly, when n = 1, P(1) = P . When n = 0, P(0) = I , where I represents the identity matrix. Just as with the one-step TPM, it follows that the row sums of P(n) must equal 1 as well. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 12 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Chapman-Kolmogorov Equations For n 2 Z+, let us consider P(n)i ,j = P(Xn = j |X0 = i) = 1X k=0 P(Xn = j |Xn 1 = k ,X0 = i)P(Xn 1 = k |X0 = i) = 1X k=0 P(n 1)i ,k P(Xn = j |Xn 1 = k ,X0 = i) = 1X k=0 P(n 1)i ,k P(Xn = j |Xn 1 = k) due to the Markov property = 1X k=0 P(n 1)i ,k Pk,j . (3.1) STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 13 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Chapman-Kolmogorov Equations We have: P(n)i ,j = P1 k=0 P (n 1) i ,k Pk,j (3.1) Recall: If A = [ai ,j ], B = [bi ,j ], and C = AB where C = [ci ,j ], then ci ,j = P k ai ,kbk,j . As a result, note that (3.1) implies that P(n) = P(n 1)P , n 2 Z+. More generally, P(n)i ,j can be expressed as P(n)i ,j = 1X k=0 P(m)i ,k P (n m) k,j 8 i , j 2 N and 0 m n, which are referred to as the Chapman-Kolmogorov equations for a DTMC. In matrix form, this translates to P(n) = P(m)P(n m), 0 m n. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 14 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts n-step Transition Probability Matrix Coming back to P(n) = P(n 1)P , n 2 Z+, let us look at a few values of n: Take n = 2 =) P(2) = P(1)P = PP = P2, Take n = 3 =) P(3) = P(2)P = P2P = P3, Take n = 4 =) P(4) = P(3)P = P3P = P4. Clearly, we see that P(n) = Pn, and so the n-step TPM is simply the one-step TPM multiplied by itself n times. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 15 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Marginal pmf of Xn For n 2 N, let us now introduce a particular row vector, which we will denote as either ↵n = (↵n,0,↵n,1, . . . ,↵n,k , . . .), or ↵n = ⇥ ↵n,0 ↵n,1 . . . ↵n,k . . . ⇤ , where ↵n,k = P(Xn = k) 8 k 2 N. In other words, ↵n,k represents the marginal pmf of Xn, n 2 N. As a consequence, it follows that P1k=0 ↵n,k = 1 8 n 2 N. If we focus on the case when n = 0, then ↵0 is referred to as the initial probability row vector of the DTMC, or simply the initial conditions of the DTMC. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 16 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Marginal pmf of Xn For n 2 Z+, note that ↵n,k = P(Xn = k) = 1X i=0 P(Xn = k |Xm = i)P(Xm = i) where 0 m n = 1X i=0 ↵m,iP(Xn m = k |X0 = i) due to the stationary assumption = 1X i=0 ↵m,iP (n m) i ,k . In matrix form, the above relation implies that ↵n = ↵mP (n m) = ↵mP n m, 0 m n, which subsequently leads to ↵n = ↵0P (n) = ↵0P n, n 2 N. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 17 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Probabilities of Interest Having knowledge of the initial conditions and the one-step transition probabilities, one can readily calculate various probabilities of possible interest such as P(Xn = xn,Xn 1 = xn 1, . . . ,X1 = x1,X0 = x0) = P(X0 = x0)P(X1 = x1|X0 = x0)P(X2 = x2|X1 = x1,X0 = x0) · · · ⇥ P(Xn = xn|Xn 1 = xn 1,Xn 2 = xn 2, . . . ,X0 = x0) = P(X0 = x0)P(X1 = x1|X0 = x0)P(X2 = x2|X1 = x1) · · ·P(Xn = xn|Xn 1 = xn 1) = ↵0,x0Px0,x1Px1,x2 · · ·Pxn 1,xn , where the second last equality follows due to the Markov property. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 18 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Probabilities of Interest Similarly, P(Xn+m = xn+m,Xn+m 1 = xn+m 1, . . . ,Xn+1 = xn+1|Xn = xn) = P(Xn+m=xn+m,Xn+m 1=xn+m 1,...,Xn+1=xn+1,Xn=xn)P(Xn=xn) = P(Xn=xn)P(Xn+1=xn+1|Xn=xn)···P(Xn+m=xn+m|Xn+m 1=xn+m 1,Xn+m 2=xn+m 2,...,Xn=xn)P(Xn=xn) = Pxn,xn+1Pxn+1,xn+2 · · ·Pxn+m 1,xn+m . The key observation here is that the DTMC is completely characterized by its one-step TPM P and the initial conditions ↵0. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 19 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Probabilities of Interest Example 3.2. A particle moves along the states {0, 1, 2} according to a DTMC whose TPM is given by P = 24 0 1 2 0 0.7 0.2 0.1 1 0 0.6 0.4 2 0.5 0 0.5 35. Let Xn denote the position of the particle after the nth move. Suppose that the particle is equally likely to start in any of the three positions. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 20 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Probabilities of Interest (a) Calculate P(X3 = 1|X0 = 0). Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 21 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Probabilities of Interest (b) Calculate P(X4 = 2). Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 22 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Probabilities of Interest (c) Calculate P(X6 = 0,X4 = 2). Solution:... (d) Calculate P(X9 = 0,X7 = 2|X5 = 1,X2 = 0). Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 23 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Accessibility and Communication With these basic results in place, we next consider the classification of states in a DTMC. Definition: State j is accessible from state i (denoted by i ! j) if 9 n 2 N such that P(n)i ,j > 0. If states i and j are accessible from each other, then states i and j communicate (denoted by i $ j). In other words, i $ j i↵ 9 m, n 2 N such that P(n)i ,j > 0 and P(m)j ,i > 0. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 24 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Accessibility and Communication In terms of accessibility, note that the size of the components of P do not matter. All that matters is which are positive and which are 0. In particular, if state j is not accessible from state i , then P(n)i ,j = 0 8 n 2 N and P(DTMC ever visits state j |X0 = i) = P ([1n=0{Xn = j}|X0 = i) 1X n=0 P(Xn = j |X0 = i) due to Boole’s inequality (see Exercise 1.1.1) = 1X n=0 P(n)i ,j = 0, implying that P(DTMC ever visits state j |X0 = i) = 0. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 25 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Equivalence Relation The concept of communication forms what is known as an equivalence relation, satisfying the following criteria: (i) Reflexivity: i $ i . Clearly true since P(0)i ,i = 1 > 0. (ii) Symmetry: i $ j =) j $ i . This is obviously true by definition. (iii) Transitivity: i $ j and j $ k =) i $ k . To see this holds formally, we know that 9 n 2 N such that P(n)i ,j > 0. Also, 9 m 2 N such that P(m)j ,k > 0. Using the Chapman-Kolmogorov equations, we have that P(n+m)i ,k = 1X l=0 P(n)i ,l P (m) l ,k P(n)i ,j P(m)j ,k > 0. Therefore, P(n+m)i ,k > 0, implying that i ! k . Using precisely the same logic, it is straightforward to show that k ! i . Thus, by definition, i $ k . STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 26 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Communication Classes The fact that communication forms an equivalence relation allows us to partition all the states of a DTMC into various communication classes, so that within each class, all states communicate. However, if states i and j belong to di↵erent classes, then i $ j is not true (i.e., at most one of i ! j or j ! i can be true). Definition: A DTMC that has only one communication class is said to be irreducible. On the other hand, a DTMC is called reducible if there are two or more communication classes. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 27 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Communication Classes Example 3.2. (continued) What are the communication classes of the DTMC? P = 24 0 1 2 0 0.7 0.2 0.1 1 0 0.6 0.4 2 0.5 0 0.5 35. Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 28 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Communication Classes Example 3.3. Consider a DTMC with TPM P = 2664 0 1 2 3 0 0 1 0 0 1 0 0 1 0 2 0 0 0 1 3 0.5 0 0.5 0 3775. What are the communication classes of this DTMC? Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 29 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Communication Classes Example 3.4. Consider a DTMC with TPM P = 266664 0 1 2 3 0 13 2 3 0 0 1 12 1 4 1 8 1 8 2 0 0 1 0 3 34 1 4 0 0 377775. What are the communication classes of this DTMC? Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 30 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Periodicity Definition: The period of state i is given by d(i) = gcd{n 2 Z+ : P(n)i ,i > 0}, where gcd{·} denotes the greatest common divisor of a set of positive integers. Remark: If d(i) = 1, then state i is said to be aperiodic. In fact, a DTMC is said to be aperiodic if d(i) = 1 8 i 2 N. Furthermore, if P(n)i ,i = 0 8 n 2 Z+, then we set d(i) =1. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 31 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Periodicity Example 3.5. Consider a DTMC with TPM P = 2664 0 1 2 3 0 13 0 0 2 3 1 12 1 4 1 8 1 8 2 0 0 1 0 3 34 0 0 1 4 3775. Determine the communication classes of this DTMC and the period of each state. Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 32 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Periodicity Example 3.3. (continued) Recall that {0, 1, 2, 3} is the only communication class for the DTMC with TPM P = 2664 0 1 2 3 0 0 1 0 0 1 0 0 1 0 2 0 0 0 1 3 0.5 0 0.5 0 3775. Determine the period of each state. Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 33 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Periodicity Example 3.6. Consider the DTMC with TPM P = 2664 0 1 2 3 0 12 1 2 0 0 1 23 1 3 0 0 2 0 0 0 1 3 0 0 1 0 3775. Find the communication classes of this DTMC and determine the period of each state. Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 34 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Periodicity The above examples illustrate some of the kinds of periodic behaviour that can be exhibited by DTMCs. However, we do observe that among the states within a given communication class, it seems as though the periodic behaviour is consistent. This is not a coincidence, as the next theorem indicates. Theorem 3.1. If i $ j , then d(i) = d(j). Proof:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 35 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Periodicity Example 3.7. Consider a DTMC with TPM P = 264 0 1 2 0 0 12 1 2 1 12 0 1 2 2 12 1 2 0 375. Find the communication classes of this DTMC and determine the period of each state. Solution:... STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 36 / 37 3. Discrete-time Markov Chains 3.1 Definitions and Basic Concepts Periodicity Remark: As the previous example demonstrates, it is still possible to observe aperiodic behaviour even though the main diagonal components of P are all zero. More generally, if d(i) = k , then this does not necessarily imply that P(k)i ,i > 0. Instead, it implies that if the DTMC is in state i at time 0, then it is impossible to observe the DTMC in state i at time n 2 Z+ if n is not a multiple of k (i.e., P(n)i ,i = 0 for such n). STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.1) Spring 2021 37 / 37 3. Discrete-time Markov Chains 3.2 Transience and Recurrence Section 3.2: Transience and Recurrence STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.2) Spring 2021 2 / 21 3. Discrete-time Markov Chains 3.2 Transience and Recurrence First Visit Probability We now wish to take a closer look at the likelihood of a DTMC beginning in some state of returning to that particular state. To that end, let us consider the probability that, starting from state i , the first visit of the DTMC to state j occurs at time n 2 Z+, to be denoted by f (n)i ,j = P(Xn = j ,Xn 1 6= j , . . . ,X2 6= j ,X1 6= j |X0 = i) 8 i , j 2 N. Clearly, we see that f (1)i ,j = Pi ,j . STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.2) Spring 2021 3 / 21 3. Discrete-time Markov Chains 3.2 Transience and Recurrence First Visit Probability For n 2, however, the determination of f (n)i ,j becomes more complicated, and so we wish to construct a procedure which will enable us to compute f (n)i ,j for such n. To do so, we consider the quantity P(n)i ,j , n 2 Z+, and condition on the time that the first visit to state j is made: P(n)i ,j = P(Xn = j |X0 = i) = nX k=1 P(Xn = j , first visit to state j occurs at time k |X0 = i) = nX k=1 P(Xn = j ,Xk = j ,Xk 1 6= j , . . . ,X2 6= j ,X1 6= j |X0 = i) = nX k=1 P(Xk = j ,Xk 1 6= j , . . . ,X2 6= j ,X1 6= j |X0 = i)P(Xn = j |Xk = j) = nX k=1 f (k)i ,j P (n k) j ,j , (3.2) where we applied the Markov property in the second last equality. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.2) Spring 2021 4 / 21 3. Discrete-time Markov Chains 3.2 Transience and Recurrence First Visit Probability We have: P(n)i ,j = Pn k=1 f (k) i ,j P (n k) j ,j (3.2) From (3.2), we can write P(n)i ,j = f (n) i ,j P (0) j ,j + n 1X k=1 f (k)i ,j P (n k) j ,j = f (n) i ,j + n 1X k=1 f (k)i ,j P (n k) j ,j , implying that f (n)i ,j = P (n) i ,j n 1X k=1 f (k)i ,j P (n k) j ,j , n 2 Z+. (3.3) For n 2, note that (3.3) yields a recursive means to compute f (n)i ,j . STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.2) Spring 2021 5 / 21 3. Discrete-time Markov Chains 3.2 Transience and Recurrence Transience and Recurrence Define the related quantity: fi ,j = P(DTMC ever makes a future visit to state j |X0 = i). Note that fi,j = 1X k=1 P(DTMC ever makes a future visit to state j ,DTMC visits state j for 1st time at time k|X0= i) = 1X k=1 P(DTMC visits state j for 1st time at time k|X0 = i) = 1X k=1 f (k)i,j 1. This leads to the following important concept in the study of Markov chains. Definition: State i is said to be transient if fi ,i < 1. On the other hand, state i is said to be recurrent if fi ,i = 1. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.2) Spring 2021 6 / 21 3. Discrete-time Markov Chains 3.2 Transience and Recurrence Transience and Recurrence Xn Time 0 i fi,iz }| { visit 1 fi,iz }| { visit 2 · · · · · · fi,iz }| { visit k 1 1 fi,i visit k · · · In what follows, we proceed to look at alternative ways of characterizing the notions of transience and recurrence. As such, let us first define Mi to be a rv which counts the number of (future) times the DTMC visits state i (disregarding the possibility of starting in state i at time 0). If we assume that fi ,i < 1, then the Markov property and stationary assumption jointly imply that P(Mi = k |X0 = i) = kY n=1 fi ,i ! (1 fi ,i ) = f ki ,i (1 fi ,i ), k = 0, 1, 2, . . . , (3.4) as the DTMC will return to state i k times with probability fi ,i , and then never return with probability 1 fi ,i . STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.2) Spring 2021 7 / 21 3. Discrete-time Markov Chains 3.2 Transience and Recurrence Transience and Recurrence We have: P(Mi = k |X0 = i) = f ki ,i (1 fi ,i ), k = 0, 1, 2, . . . (3.4) We recognize (3.4) as the pmf of a GEOf (1 fi ,i ) rv, thereby implying that E[Mi |X0 = i ] = 1 (1 fi ,i )1 fi ,i = fi ,i 1 fi ,i <1 since fi ,i < 1. However, if fi ,i = 1, then P(Mi =1|X0 = i) = 1, immediately implying that E[Mi |X0 = i ] =1. Therefore, an equivalent way of viewing transience/recurrence is as follows: E[Mi |X0 = i ] ( <1 , i↵ state i is transient, =1 , i↵ state i is recurrent. STAT 333 Ch 3. Discrete-time Markov Chains (Section 3.2) Spring 2021 8 / 21 3. Discrete-time Markov Chains 3.2 Transience and Recurrence Transience and Recurrence Following further on the notion of Mi , define a sequence of indicator rvs {An}1n=1 such that An = ( 0 , if Xn 6= i , 1 , if Xn = i . With this definition, note that Mi = P1 n=1 An. Now, E[Mi |X0 = i ] = E " 1X n=1 An