Tutorials › Stochastic Processes › Meaning of n-Step Transition Probabilities

Chapman-Kolmogorov Equations and n-Step Transitions · Tutorial 108 of 994

Meaning of n-Step Transition Probabilities

Learn to read an n-step transition probability as the chance of being in a specified state at a fixed future time, and compute and check it in finite chains.

Intermediate 12 min read

What You'll Learn

  • Define n-step transition probabilities for a time-homogeneous Markov chain
  • Distinguish the zero-step, one-step and n-step transition matrices
  • Interpret matrix powers as fixed-horizon state probabilities
  • Check that finite-state transition matrices and their computed powers are stochastic
  • Separate an exact-time probability from the probability of having hit a state by that time

From One Step to a Fixed Future Time

A one-step transition probability describes the next move. Many questions instead concern where a chain will be after several moves: what is the probability that tomorrow’s weather is sunny, given today’s state, or that a particle is in a particular location after four updates? The n-step transition probability answers this fixed-horizon question. Its meaning is simple, but keeping the time horizon precise prevents confusion with one-step probabilities and with the probability of having visited a state at any earlier time.

Let \((X_n)_{n\ge0}\) be a discrete-time Markov chain on a finite state space \(S\), defined on a probability space \((\Omega,\mathcal{F},\mathbb{P})\), and let \(\mathcal{F}_n=\sigma(X_0,\ldots,X_n)\) be its natural filtration. The chain is time-homogeneous when its one-step transition probabilities do not depend on \(n\): for every \(i,j\in S\), \(\mathbb{P}(X_{n+1}=j\mid\mathcal{F}_n)=p_{ij}\) almost surely on \(\{X_n=i\}\). Write \(P=(p_{ij})\) for its transition matrix, with rows indexed by the current state and columns by the next state. Each row sums to \(1\).

Definition: For \(i,j\in S\) and an integer \(n\ge0\), the n-step transition probability from \(i\) to \(j\) is the probability of being in \(j\) exactly \(n\) steps after starting in \(i\): \(p_{ij}^{(n)}=\mathbb{P}_i(X_n=j)\). Equivalently, when \(\mathbb{P}(X_m=i)>0\), time-homogeneity gives \(\mathbb{P}(X_{m+n}=j\mid X_m=i)=p_{ij}^{(n)}\), independent of \(m\). The matrix of these probabilities is \(P^{(n)}=(p_{ij}^{(n)})\). By convention \(P^{(0)}=I\) and \(p_{ij}^{(0)}=\mathbf{1}_{\{i=j\}}\). The standard result that \(P^{(n)}=P^n\), the nth matrix power, is derived later from the Chapman–Kolmogorov equations.

The defining event is \(X_n=j\), not “the chain moves from \(i\) to \(j\) in one step.” Thus \(p_{ij}^{(1)}=p_{ij}\), while \(p_{ij}^{(n)}\) for \(n>1\) concerns the state at a later time. The chain may visit \(j\) before time \(n\), leave it, and end in a different state; an n-step probability does not record those visits.

Proposition: For each starting state \(i\) and each \(n\ge0\), the n-step probabilities are nonnegative and satisfy \(\sum_{j\in S}p_{ij}^{(n)}=1\). Also \(P^{(0)}=I\) and \(P^{(1)}=P\).

Proof. Fix \(i\) and \(n\). The events \(\{X_n=j\}\), as \(j\) ranges over the finite state space \(S\), are disjoint and their union is the whole sample space under \(\mathbb{P}_i\), since \(X_n\) takes exactly one value in \(S\). Their probabilities are nonnegative and sum to \(1\). At time \(0\), a chain started in \(i\) is in \(i\) with probability \(1\), and in every \(j\ne i\) with probability \(0\), giving the identity matrix. At time \(1\), the probabilities are the one-step transition probabilities \(p_{ij}\), giving \(P^{(1)}=P\). \(\square\)

Consequently, every matrix of n-step probabilities is row-stochastic: it has nonnegative entries and each row sums to \(1\). In examples, we can calculate \(P^n\) by ordinary matrix multiplication and interpret its entries as n-step probabilities using the standard identity stated above. We will check each matrix row and the resulting probabilities directly.

Reading the Entries of a Matrix Power

Worked Example: Two-State Weather Chain

Suppose \(S=\{H,C\}\), where \(H\) means a hot day and \(C\) a cool day. From \(H\), the next day is hot with probability \(\tfrac34\) and cool with probability \(\tfrac14\). From \(C\), the next day is hot with probability \(\tfrac13\) and cool with probability \(\tfrac23\). Find the probability of being hot in two days, starting from each state.

The transition matrix, with rows and columns in the order \(H,C\), is

$$ P=\begin{bmatrix} \tfrac34&\tfrac14\\[2pt] \tfrac13&\tfrac23 \end{bmatrix}. $$

The first row sums to \(\tfrac34+\tfrac14=1\), and the second row sums to \(\tfrac13+\tfrac23=1\), so the matrix is a valid transition matrix. Using \(P^{(2)}=P^2\), compute each entry:

$$ P^2= \begin{bmatrix} \tfrac34\cdot\tfrac34+\tfrac14\cdot\tfrac13& \tfrac34\cdot\tfrac14+\tfrac14\cdot\tfrac23\\[3pt] \tfrac13\cdot\tfrac34+\tfrac23\cdot\tfrac13& \tfrac13\cdot\tfrac14+\tfrac23\cdot\tfrac23 \end{bmatrix} = \begin{bmatrix} \tfrac{31}{48}&\tfrac{17}{48}\\[2pt] \tfrac{17}{36}&\tfrac{19}{36} \end{bmatrix}. $$

For example, starting hot, the two-day probability of being hot is \(\tfrac{31}{48}\), not the one-day probability \(\tfrac34\). Starting cool, that probability is \(\tfrac{17}{36}\). Check the rows: \(\tfrac{31}{48}+\tfrac{17}{48}=1\), and \(\tfrac{17}{36}+\tfrac{19}{36}=1\). Each row therefore still describes a probability distribution over the state at the specified future time.

Worked Example: A Deterministic Three-State Cycle

Consider a chain on \(S=\{0,1,2\}\) that moves from \(0\) to \(1\), from \(1\) to \(2\), and from \(2\) to \(0\), each with probability \(1\). Determine its two-step and three-step transition probabilities.

In the state order \(0,1,2\), its transition matrix is

$$ P=\begin{bmatrix} 0&1&0\\ 0&0&1\\ 1&0&0 \end{bmatrix}. $$

Every row sums to \(1\): the sums are \(0+1+0=1\), \(0+0+1=1\), and \(1+0+0=1\). Multiplying gives

$$ P^2= \begin{bmatrix} 0&0&1\\ 1&0&0\\ 0&1&0 \end{bmatrix}, \qquad P^3= \begin{bmatrix} 1&0&0\\ 0&1&0\\ 0&0&1 \end{bmatrix}. $$

Thus, for example, \(p_{0,2}^{(2)}=1\): starting at \(0\), the chain is at \(2\) exactly two steps later. But \(p_{0,2}^{(3)}=0\), because after three steps the chain is back at \(0\). Each row of both displayed matrices sums to \(1\), as expected. This example also shows why the time horizon belongs in the notation: the same pair of states can have different transition probabilities at different numbers of steps.

Worked Example: A Walk with Holding Probabilities

A particle moves on \(S=\{0,1,2\}\). At the endpoints it stays put or moves inward, each with probability \(\tfrac12\). From state \(1\), it moves to either endpoint with probability \(\tfrac14\) each and stays at \(1\) with probability \(\tfrac12\). Find the two-step probability of being at state \(2\), starting from each state.

In state order \(0,1,2\),

$$ P=\begin{bmatrix} \tfrac12&\tfrac12&0\\ \tfrac14&\tfrac12&\tfrac14\\ 0&\tfrac12&\tfrac12 \end{bmatrix}. $$

The row sums are \(\tfrac12+\tfrac12+0=1\), \(\tfrac14+\tfrac12+\tfrac14=1\), and \(0+\tfrac12+\tfrac12=1\). The two-step matrix is obtained by multiplying \(P\) by \(P\). Its entries are

$$ P^2= \begin{bmatrix} \tfrac12\cdot\tfrac12+\tfrac12\cdot\tfrac14& \tfrac12\cdot\tfrac12+\tfrac12\cdot\tfrac12& \tfrac12\cdot0+\tfrac12\cdot\tfrac14\\[3pt] \tfrac14\cdot\tfrac12+\tfrac12\cdot\tfrac14+\tfrac14\cdot0& \tfrac14\cdot\tfrac12+\tfrac12\cdot\tfrac12+\tfrac14\cdot\tfrac12& \tfrac14\cdot0+\tfrac12\cdot\tfrac14+\tfrac14\cdot\tfrac12\\[3pt] 0\cdot\tfrac12+\tfrac12\cdot\tfrac14+\tfrac12\cdot0& 0\cdot\tfrac12+\tfrac12\cdot\tfrac12+\tfrac12\cdot\tfrac12& 0\cdot0+\tfrac12\cdot\tfrac14+\tfrac12\cdot\tfrac12 \end{bmatrix} = \begin{bmatrix} \tfrac38&\tfrac12&\tfrac18\\ \tfrac14&\tfrac12&\tfrac14\\ \tfrac18&\tfrac12&\tfrac38 \end{bmatrix}. $$

The requested probabilities are the third-column entries: \(\tfrac18\) from state \(0\), \(\tfrac14\) from state \(1\), and \(\tfrac38\) from state \(2\). Check the row sums: \(\tfrac38+\tfrac12+\tfrac18=1\), \(\tfrac14+\tfrac12+\tfrac14=1\), and \(\tfrac18+\tfrac12+\tfrac38=1\). In particular, starting at \(0\), being at \(2\) after two steps requires ending in that state at the specified time; the entry does not describe the chance of ever reaching \(2\).

Exact-Time Probabilities and Hitting Probabilities

For a set \(A\subseteq S\), the hitting time is \(T_A=\inf\{n\ge0:X_n\in A\}\), with \(\inf\varnothing=\infty\). The event \(T_A\le n\) means the chain has visited \(A\) at least once by time \(n\). By contrast, \(X_n\in A\) asks only whether it is in \(A\) at time \(n\). The events need not be equal: a chain can visit \(A\) and then leave.

Proposition: For any \(A\subseteq S\), starting state \(i\), and integer \(n\ge0\), \(\mathbb{P}_i(X_n\in A)\le\mathbb{P}_i(T_A\le n)\).

Proof. If \(X_n\in A\), then the chain has been in \(A\) at some time no later than \(n\), namely at time \(n\). Hence \(\{X_n\in A\}\subseteq\{T_A\le n\}\). Taking probabilities of the two events gives the inequality. The inclusion can be strict if the chain visits \(A\) by time \(n\) but is outside \(A\) at time \(n\). \(\square\)

For example, return to the weather chain and start from \(C\), with \(A=\{H\}\) and \(n=2\). The exact-time probability is \(p_{C,H}^{(2)}=\tfrac{17}{36}\). The probability of having hit \(H\) by time \(2\) is larger: the chain can go from \(C\) to \(H\) on the first step, with probability \(\tfrac13\), or stay at \(C\) first and then move to \(H\), with probability \(\tfrac23\cdot\tfrac13=\tfrac29\). These disjoint possibilities give \(\tfrac13+\tfrac29=\tfrac59=\tfrac{20}{36}\), which exceeds \(\tfrac{17}{36}\). The difference \(\tfrac{3}{36}\) is the probability of having reached \(H\) by time \(2\) but ending at \(C\) at time \(2\).

A Common Misreading

Pitfall: squaring a one-step entry. In the weather example, an incorrect calculation for the probability of being hot after two days starting hot is \((\tfrac34)^2=\tfrac{9}{16}\). This counts only the route that is hot after the first day and hot again after the second. It omits the possibility of being cool after the first day and hot after the second. The matrix-power calculation gives \(\tfrac34\cdot\tfrac34+\tfrac14\cdot\tfrac13=\tfrac{27}{48}+\tfrac{4}{48}=\tfrac{31}{48}\). Substituting into the two-step row gives \(\tfrac{31}{48}+\tfrac{17}{48}=1\), whereas the incorrect entry \(\tfrac{9}{16}\) alone is not the two-step probability.

When reading \(p_{ij}^{(n)}\), identify the starting state \(i\), the destination \(j\), and the exact time \(n\). Then check that the row of the n-step matrix gives probabilities over all possible destination states and sums to \(1\). If the question instead asks whether the chain has visited a state or set by time \(n\), it is a hitting-probability question, not merely an entry of \(P^{(n)}\).

Check Your Understanding

For each question, distinguish the stated time horizon from any event involving an earlier visit.

  1. What does \(p_{ij}^{(0)}\) equal when \(i=j\), and when \(i\ne j\)?
  2. In the weather example, what is the probability of being cool in two days when starting from \(H\)?
  3. For the deterministic cycle, what is \(p_{1,0}^{(2)}\), and what is \(p_{1,1}^{(3)}\)?
  4. Explain why the probability of having hit a set by time \(n\) can exceed the probability of being in the set exactly at time \(n\).
  5. In the holding-probability example, verify by addition that the two-step probabilities in the row for state \(1\) sum to \(1\).