Tutorials › Exam ALTAM - Advanced Long-Term Actuarial Mathematics › Why Same-State Probabilities Can Fall Then Rise

Exam ALTAM · Transition Probabilities · Tutorial 49 of 999

Why Same-State Probabilities Can Fall Then Rise

See why return paths can make a same-state probability rise after its initial decline, and calculate that effect in a three-state continuous-time model.

Beginner 10 min read

What You'll Learn

  • Distinguish the probability of ending in the starting state from the probability of making no transition.
  • Set up a three-state continuous-time Markov chain with a directed cycle.
  • Derive the same-state transition probability for equal transition intensities.
  • Identify intervals on which the same-state probability falls or rises.
  • Apply the formula with explicit rates and time units.
  • Recognize when a model has no paths that permit a return.

Why “Still in the Same State” Does Not Always Mean “No Change”

A person can begin in a state, leave it, and later return. Consequently, the probability of occupying the starting state at time \(t\) includes paths with transitions as well as paths with no transitions. It need not behave like a holding-time survival probability, which counts only the chance of making no transition during the interval.

This distinction matters whenever a model allows transitions in more than one direction. A same-state probability can initially decrease as probability leaves the starting state, then rise as probability returns. The three-state example below makes this effect explicit using a hypothetical model with constant transition intensities.

Definition: In a continuous-time Markov model, \(p_{ij}(t)\) is the probability that the process is in state \(j\) after elapsed time \(t\), given that it started in state \(i\). The same-state probability is \(p_{ii}(t)\). It includes every possible path from \(i\) to \(i\) over the interval, including paths that leave and return.

A Three-State Cycle

Consider three states, \(A\), \(B\), and \(C\), with transitions only in the directed cycle \(A\to B\to C\to A\). Suppose each transition has constant intensity \(\lambda>0\) per year. The model is time-homogeneous: its transition intensities do not vary with calendar time or elapsed time. Each time the process enters a state, its holding time there is exponential with rate \(\lambda\), and the next state is determined by the stated cycle.

Using rows for starting states and columns for ending states, the intensity matrix is

$$ Q= \begin{pmatrix} -\lambda & \lambda & 0\\ 0 & -\lambda & \lambda\\ \lambda & 0 & -\lambda \end{pmatrix}. $$

The diagonal entry in each row is the negative of the total rate of leaving that state. The transition matrix after elapsed time \(t\) is \(P(t)=e^{Qt}\). We focus on its \(A,A\) entry, \(p_{AA}(t)\).

Worked Example: A Same-State Probability Falls and Then Rises

For a hypothetical three-state cycle, let \(\lambda=1\) per year. Calculate the probability of being in state \(A\) at two selected times, and determine whether it falls or rises between those times. Use the times corresponding to phases \(2\pi/3\) and \(\pi\), where the phase is \(\sqrt{3}\lambda t/2\).

Method. Use the derived expression \(p_{AA}(t)=\frac13+\frac23e^{-3\lambda t/2}\cos(\sqrt{3}\lambda t/2)\). The specified phase gives each time directly:

$$ t_1=\frac{2(2\pi/3)}{\sqrt{3}\lambda} =\frac{4\pi}{3\sqrt{3}} \approx 2.4184\text{ years}, \qquad t_2=\frac{2\pi}{\sqrt{3}\lambda} =\frac{2\pi}{\sqrt{3}} \approx 3.6276\text{ years}. $$

At \(t_1\), the phase is \(2\pi/3\), so the cosine is \(-1/2\). The exponential exponent is \(-3t_1/2=-2\pi/\sqrt{3}\approx-3.6276\), giving \(e^{-3.6276}\approx0.0266\). Therefore

$$ p_{AA}(t_1)=\frac13+\frac23(0.0266)\left(-\frac12\right) =\frac13-\frac13(0.0266) \approx 0.3245. $$

At \(t_2\), the phase is \(\pi\), so the cosine is \(-1\). The exponent is \(-3t_2/2=-\sqrt{3}\pi\approx-5.4414\), and \(e^{-5.4414}\approx0.0043\). Thus

$$ p_{AA}(t_2)=\frac13+\frac23(0.0043)(-1) \approx0.3304. $$

Rounded to four decimal places, the probability rises from \(0.3245\) to \(0.3304\). This is a return effect: by these times some paths have gone from \(A\) to \(B\) to \(C\) and back to \(A\). The calculation is for occupying \(A\) at the endpoint, not for remaining in \(A\) throughout.

Why the Formula Has This Form

For this equal-rate cycle, the same-state probability can be written in a compact form:

Result: For the time-homogeneous three-state cycle \(A\to B\to C\to A\), with all transition intensities equal to \(\lambda>0\),
$$ p_{AA}(t)=\frac13+\frac23 e^{-3\lambda t/2} \cos\left(\frac{\sqrt{3}\lambda t}{2}\right),\qquad t\geq0. $$

One justification uses the matrix exponential and a scalar differential equation. The characteristic polynomial of \(Q\) is \(r(r^2+3\lambda r+3\lambda^2)\), so \(Q\) satisfies \(Q^3+3\lambda Q^2+3\lambda^2Q=0\). Since \(P(t)=e^{Qt}\), differentiating \(P\) gives \(P'=PQ\). Applying the polynomial identity to \(P\) shows that \(f(t)=p_{AA}(t)\) satisfies

$$ f'''(t)+3\lambda f''(t)+3\lambda^2f'(t)=0. $$

The initial values follow from \(P(0)=I\), \(P'(0)=Q\), and \(P''(0)=Q^2\). In particular, \(Q_{AA}=-\lambda\) and \((Q^2)_{AA}=\lambda^2\), so \(f(0)=1\), \(f'(0)=-\lambda\), and \(f''(0)=\lambda^2\). The characteristic roots of the differential equation are \(0\) and \(-3\lambda/2\pm i\sqrt{3}\lambda/2\). Therefore its real-valued solution has the form

$$ f(t)=c+e^{-3\lambda t/2} \left(C\cos\left(\frac{\sqrt{3}\lambda t}{2}\right) +D\sin\left(\frac{\sqrt{3}\lambda t}{2}\right)\right). $$

Substitution of \(f(0)=1\) and the two derivative conditions gives \(c=1/3\), \(C=2/3\), and \(D=0\). This yields the stated expression. The constant \(1/3\) is also the long-run probability of each state in this symmetric cycle; the exponential factor makes the oscillations around that value decay.

What Determines Whether the Probability Rises?

Differentiating the formula gives

$$ p'_{AA}(t)=\frac23e^{-3\lambda t/2} \left[ -\frac{3\lambda}{2}\cos\left(\frac{\sqrt{3}\lambda t}{2}\right) -\frac{\sqrt{3}\lambda}{2}\sin\left(\frac{\sqrt{3}\lambda t}{2}\right) \right]. $$

At \(t=0\), this derivative is \(-\lambda\), so the probability initially decreases. It reaches its first minimum when the phase is \(2\pi/3\). For phases strictly between \(2\pi/3\) and \(\pi\), the bracket is positive, so the probability rises over that interval. Later it can fall and rise again in smaller oscillations, because the oscillating term is damped by \(e^{-3\lambda t/2}\).

The same conclusion can be understood as a balance of probability leaving and returning. For a time-homogeneous model, \(P'(t)=P(t)Q\). The \(A,A\) entry of that equation gives

$$ p'_{AA}(t) =-\lambda p_{AA}(t)+\lambda p_{CA}(t). $$

The first term is the rate at which probability mass currently in \(A\) leaves \(A\). The second is the rate at which mass in \(C\) moves into \(A\). Once return flow is large enough relative to outflow, the net derivative can be positive. At the initial instant, no mass is in \(C\), so there is outflow but no return flow.

Worked Example: When There Is No Route Back

Now consider a different hypothetical model with states \(A\), \(B\), and \(C\): transitions are \(A\to B\) and \(B\to C\), each at a constant rate of \(1\) per year, and \(C\) is absorbing. Starting in \(A\), calculate the probability of being in \(A\) after one year and after three years.

Method. There is no transition into \(A\). To be in \(A\) at time \(t\), the process must therefore have made no transition out of \(A\) by that time. The holding time in \(A\) is exponential with rate \(1\) per year, so its survival probability is \(e^{-t}\).

$$ p_{AA}(1)=e^{-1}\approx0.3679, \qquad p_{AA}(3)=e^{-3}\approx0.0498. $$

The probability decreases, as expected when there is no return route. The absorbing state \(C\) does not itself make \(A\) absorbing; it is simply the state from which the process cannot leave. This example distinguishes a model with no possible return to \(A\) from the cycle, where a return path exists.

Worked Example: Respecting the Time Unit

Suppose a hypothetical three-state cycle has equal intensities of \(0.25\) per month. Starting in \(A\), calculate the same-state probability after four months and after twelve months, and compare the results. Round probabilities to four decimal places.

Method. Keep the intensity and elapsed time in matching units. Here the products \(\lambda t\) are \(0.25(4)=1\) and \(0.25(12)=3\). For four months, the exponent is \(-3/2\) and the phase is \(\sqrt{3}/2\):

$$ p_{AA}(4)=\frac13+\frac23e^{-1.5}\cos(0.8660) \approx\frac13+\frac23(0.2231)(0.6479) \approx0.4297. $$

For twelve months, the exponent is \(-4.5\) and the phase is \(3\sqrt{3}/2\):

$$ p_{AA}(12)=\frac13+\frac23e^{-4.5}\cos(2.5981) \approx\frac13+\frac23(0.0111)(-0.8559) \approx0.3270. $$

The probability has moved closer to its long-run value of \(1/3\). The result also illustrates why an intensity stated per month must be paired with time in months, not years. Using four as a number of years would describe a different elapsed time and give a different probability.

Common Mistake and Practical Interpretation

A frequent error is to equate \(p_{AA}(t)\) with the probability of no transition. In the cycle, the no-transition probability from \(A\) is \(e^{-\lambda t}\), but this counts only paths that stayed in \(A\) continuously. It omits paths that left \(A\) and returned. For example, at the second time in the first worked example, the no-transition probability is \(e^{-5.4414}\approx0.0043\), while the same-state probability is about \(0.3304\). The much larger endpoint probability includes return paths.

The distinction is useful in actuarial models where a state can be entered again, such as hypothetical models of temporary disability, recovery, or care status. A same-state transition probability describes status at the endpoint; it does not say whether the person remained in that state at every time in between. If the product or question depends on whether a transition occurred at all, an endpoint transition probability alone is not enough.

Key takeaway: A same-state probability measures where the process is at the endpoint, not whether it stayed put. In a model with return paths, outflow can dominate initially and return flow can later make the probability rise.

Check Your Understanding

Use the stated models and units in each question. For cycle-model calculations, use the derived equal-rate formula.

  1. A directed three-state cycle has equal intensities of \(1\) per year. Calculate \(p_{AA}(2)\), rounding to four decimal places. (Answer: \(0.3280\))
  2. A directed three-state cycle has equal intensities of \(0.5\) per year. Calculate \(p_{AA}(2)\), rounding to four decimal places. (Answer: \(0.4297\))
  3. In the equal-rate cycle, determine whether \(p'_{AA}(t)\) is positive, negative, or zero when the phase \(\sqrt{3}\lambda t/2\) is \(5\pi/6\). (Answer: Positive)
  4. A hypothetical model has states \(A,B,C\), transitions only \(A\to B\) and \(B\to C\), and an \(A\to B\) intensity of \(0.4\) per month. Starting in \(A\), calculate \(p_{AA}(5)\), rounding to four decimal places. (Answer: \(e^{-2}\approx0.1353\))
  5. For the equal-rate directed cycle, what is the limiting value of \(p_{AA}(t)\) as \(t\) tends to infinity, and what feature of the formula gives this limit? (Answer: \(1/3\); the exponential term tends to zero.)
Clinical Biostats