Why Duration Can Matter
Two people can occupy the same state at the same calendar time but have different transition risks because they have spent different amounts of time there. A recently diagnosed patient and a patient who has remained in a disease state for years, for example, need not have the same instantaneous risk of leaving that state. A model that assigns one fixed transition rate to everyone in the state cannot represent this difference.
A duration-dependent transition intensity makes the risk of moving from state \(i\) to state \(j\) depend on the elapsed time since entry into \(i\). This is a natural extension of the Markov model, in which transition intensities depend on the current state but not on how long it has been occupied. The earlier material on the Markov property and sojourn time establishes why elapsed time can be part of the information needed to describe future transitions.
The definition is conditional on still being in state \(i\) at duration \(u\). For small \(h\), the transition probability is approximately \(\alpha_{ij}(u)h\); it is not the probability of eventually entering \(j\), nor a probability that remains fixed over a long interval. If state \(i\) has several possible destinations, each destination has its own intensity. Under the usual assumption that simultaneous transitions have probability of smaller order than \(h\), the probability of any exit in \((t,t+h]\) is approximately the sum of the destination-specific intensities multiplied by \(h\).
From Intensities to Exit Risk
Suppose the possible exits from state \(i\) are to states in a set \(J_i\). The total exit intensity at duration \(u\) is $$ \lambda_i(u)=\sum_{j\in J_i}\alpha_{ij}(u). $$ It describes the instantaneous rate of leaving \(i\), conditional on still occupying \(i\) at duration \(u\). The destination-specific intensities describe how that instantaneous exit rate is allocated among destinations.
A useful consequence is the relationship between the total intensity and the probability of remaining in state \(i\). Assume the intensities are nonnegative and locally integrable, and that the model's transition mechanism is specified by these intensities conditional on the current state and duration. Write \(S_i(u)\) for the probability of no exit from \(i\) during the first \(u\) units after entry, for an individual entering at duration zero under the model. Over a short additional interval \(h\), conditional on remaining in \(i\) until duration \(u\), the chance of exit is \(\lambda_i(u)h+o(h)\). Thus $$ S_i(u+h)=S_i(u)\{1-\lambda_i(u)h+o(h)\}. $$ Subtract \(S_i(u)\), divide by \(h\), and let \(h\) decrease to zero. At durations where the derivative exists, this gives \(S_i'(u)=-\lambda_i(u)S_i(u)\). Since \(S_i(0)=1\), integrating the differential equation yields $$ S_i(u)=\exp\left\{-\int_0^u\lambda_i(v)\,dv\right\}. $$
This derivation also shows why an intensity and a fixed-interval probability should not be confused. The probability of at least one exit by duration \(u\) is \(1-S_i(u)\), which depends on the accumulated total intensity, not just its value at the endpoint. If there are multiple possible destinations, the same total exit probability applies, but the destination of an exit depends on the separate intensities.
Worked Example: An Increasing Exit Intensity
Consider a hypothetical state with one possible exit and duration \(u\) measured in months. Suppose the exit intensity is \(\alpha(u)=0.02+0.01u\) per month, with the numerical value of \(u\) expressed in months. This is nonnegative and increases with duration. At duration 3 months, the intensity is \(0.02+0.01(3)=0.05\) per month. For a short interval of 0.1 month, the conditional exit probability is approximately \(0.05(0.1)=0.005\), or 0.5%, to first order in the interval length.
To find the model probability of no exit during the first 4 months, integrate the intensity: $$ \int_0^4 (0.02+0.01u)\,du =0.02(4)+0.005(4^2) =0.08+0.08 =0.16. $$
Therefore, $$ S(4)=e^{-0.16}=0.8521, \qquad P(\text{exit by 4 months})=1-0.8521=0.1479. $$ Under this model, the probability of at least one exit by 4 months is about 14.79%. The increasing intensity means the instantaneous exit rate is higher at longer durations; it does not mean that the probability of exit in every fixed-length interval is itself the intensity.
What Different Duration Patterns Mean
The shape of \(\alpha_{ij}(u)\) gives a conditional, instantaneous description of how the transition risk changes as state duration increases. A rising curve indicates that, among individuals still in state \(i\) at each duration, the instantaneous rate of transition to \(j\) is higher at greater durations. A falling curve indicates a lower instantaneous rate at greater durations. Neither statement says that every individual follows the same path or that an individual’s risk necessarily changes causally because time has passed.
A duration effect can be monotone or nonmonotone. Risk may rise soon after entry, peak, and then decline; it may fall initially and later rise. Such patterns can arise when the mechanisms associated with leaving a state change over the course of the sojourn. A flexible model should allow the data and scientific context to support the shape rather than impose a trend by default.
Worked Example: A Decreasing Exit Intensity
Suppose, hypothetically, that a single exit has intensity \(\alpha(u)=0.12e^{-0.3u}\) per week, where \(u\) is measured in weeks. It is positive and decreases with duration because the exponential factor decreases. At duration 2 weeks, the intensity is \(0.12e^{-0.6}=0.06586\) per week. For a short interval of 0.1 week, the conditional exit probability is approximately \(0.06586(0.1)=0.006586\), or 0.659%.
The accumulated intensity over the first 5 weeks is $$ \int_0^5 0.12e^{-0.3u}\,du =\frac{0.12}{0.3}(1-e^{-0.3(5)}) =0.4(1-e^{-1.5}) =0.3108. $$
Hence the probability of no exit by week 5 is \(e^{-0.3108}=0.7329\), and the probability of exit by week 5 is \(1-0.7329=0.2671\), or 26.71%. The intensity falls over the sojourn, but cumulative exit probability still increases with time: it is the probability of having exited by a later time, not the instantaneous rate at that later time.
Worked Example: A Nonmonotone Destination-Specific Intensity
Consider a hypothetical state with two possible exits, to states 2 and 3. Let duration \(u\) be in months, with intensities, in transitions per month, $$ \alpha_{12}(u)=0.02+0.04u e^{-u}, \qquad \alpha_{13}(u)=0.03. $$ The first intensity is nonmonotone. Its derivative is \(0.04e^{-u}(1-u)\), so it increases for \(u<1\), reaches its maximum at \(u=1\), and decreases for \(u>1\). The second intensity is constant.
At duration 1 month, \(\alpha_{12}(1)=0.02+0.04e^{-1}=0.03472\) per month, while \(\alpha_{13}(1)=0.03000\) per month. The total exit intensity is therefore \(0.06472\) per month. For a short interval of 0.1 month, the approximate probabilities are \(0.03472(0.1)=0.003472\) for a transition to state 2 and \(0.03000(0.1)=0.003000\) for a transition to state 3. The approximate probability of any exit is their sum, \(0.006472\), or 0.6472%.
The probability of remaining in state 1 for 2 months uses the accumulated total intensity: $$ \int_0^2\{\alpha_{12}(u)+\alpha_{13}(u)\}\,du =0.04+0.04\{1-3e^{-2}\}+0.06 =0.1238. $$ Thus the probability of no exit by 2 months is \(e^{-0.1238}=0.8836\), or about 88.36%. This calculation gives the probability of remaining in state 1; the intensities to the two destinations are still needed to describe which exit occurs.
Assumptions, Model Choices, and Diagnostics
A duration-only model makes a substantive sufficiency assumption: once the current state and the elapsed time in it are known, the transition intensities do not additionally depend on calendar time or on the earlier route through the state system. It also assumes that entry into a state resets the duration clock for that state. The clock-reset and clock-forward conventions distinguish which time scale is used; this tutorial's duration-dependent formulation uses time since the relevant entry. If prior history or calendar time matters beyond duration, a duration-only specification may be inadequate.
A constant intensity is the special case \(\alpha_{ij}(u)=c_{ij}\) for each transition. It is simpler and corresponds to the familiar Markov specification when rates depend only on the current state. A duration-dependent model is more flexible, but flexibility has a cost: it needs enough observed transitions across the duration range to estimate the shape. A piecewise-constant intensity can be easier to communicate and estimate, while a smooth parametric curve can be compact but risks imposing a wrong shape. Flexible smoothers can represent nonmonotone patterns but may be unstable where data are sparse.
Useful checks include plotting estimated transition rates or smoothed transition counts against duration, and examining how many individuals remain at risk at long durations. Compare observed transitions with model-based expected transitions across duration bands. Large discrepancies, systematic residual patterns, or conclusions that change sharply when the duration grouping or smoothing choices change are warnings of poor fit. Pay particular attention to the longest durations: few individuals may remain there, so an apparently dramatic late increase or decrease can be driven by very little information.
A further limitation is that duration patterns in observed data need not identify a changing individual-level mechanism. Individuals who remain in a state for a long time may differ from those who leave earlier; unmeasured heterogeneity can therefore contribute to an observed duration pattern. The intensity is a conditional model quantity for those still in the state, not by itself a causal effect of waiting longer. Scientific interpretation should account for relevant covariates and study design.
Common Pitfalls
- Calling the intensity a probability. An intensity has units of inverse time and can exceed 1. Multiply it by a sufficiently short interval to obtain a first-order probability approximation.
- Interpreting a rising intensity as a rising cumulative probability. Cumulative exit probability rises over time whenever exit is possible; the duration effect describes whether the instantaneous rate rises or falls.
- Using one intensity for all destinations without checking the question. The sum describes any exit. Separate transition intensities are required to describe destination-specific risk.
- Ignoring duration reset. If duration is time since entry, it restarts when the process enters the state again. Using time since the study began instead describes a different time scale.
- Trusting a tail pattern with little support. Few individuals at risk or few observed transitions at long durations make estimates there uncertain and potentially unstable.
Check Your Understanding
Try each question, then compare your reasoning with the solutions below.
- A single exit has intensity \(0.04+0.02u\) per day, with \(u\) in days. What is the intensity at day 3, and what is the approximate exit probability over the next 0.1 day?
- For the intensity in Question 1, calculate the probability of no exit during the first 2 days.
- Two exits have intensities \(0.02\) and \(0.05\) per week at a given duration. What is the total exit intensity, and what is the first-order probability of either exit during the next 0.2 week?
- Why does a decreasing intensity not imply a decreasing probability of having exited by time \(u\)?
Solutions. (1) The intensity is \(0.04+0.02(3)=0.10\) per day. Over 0.1 day, the approximate exit probability is \(0.10(0.1)=0.010\), or 1%. (2) The accumulated intensity is \(\int_0^2(0.04+0.02u)\,du=0.08+0.04=0.12\). Therefore the probability of no exit is \(e^{-0.12}=0.8869\), about 88.69%. (3) The total intensity is \(0.02+0.05=0.07\) per week. The first-order probability of either exit over 0.2 week is \(0.07(0.2)=0.014\), or 1.4%. (4) The intensity is an instantaneous rate, whereas the probability of having exited by \(u\) accumulates risk over the entire interval from entry to \(u\). Even if the rate declines, that accumulated probability does not decline with time.