Tutorials › Competing Risks and Multi-State Models › Conditional and Marginal Intensities

Counting Processes and Martingale Theory · Tutorial 471 of 699

Conditional and Marginal Intensities

You will derive the intensity under observed information from an enlarged-filtration model and see when averaging changes the model’s functional form.

Advanced 12 min read

What You'll Learn

  • Distinguish an intensity conditional on enlarged history from one conditional on observed history.
  • Derive the observed-filtration intensity by predictable projection.
  • Calculate how survival and event history update beliefs about an unobserved risk type.
  • Explain why a conditional model’s functional form may not survive marginalization.
  • Compare cause-specific marginal intensities with their enlarged-history counterparts.
  • Identify assumptions and practical warning signs when choosing an intensity model.

Why the Information Set Changes an Intensity

An event rate is conditional on what is known just before the event could occur. If a model includes a latent risk type, its intensity describes event occurrence conditional on that type and the rest of the enlarged history. An analyst who does not observe the type instead needs the rate conditional on the observed history. These rates generally differ, even when the latent-type model is correctly specified.

The distinction matters because a rate averaged over a fixed population distribution is not usually the rate for people who remain event-free or have experienced a particular event history. Those histories change the distribution of latent types among the observed risk set. The observed-history intensity must account for that change.

Definition: Let \((\mathcal F_t)\) be the observed-history filtration and \((\mathcal G_t)\) an enlarged filtration with \(\mathcal F_t\subseteq\mathcal G_t\). For an event-counting process \(N(t)\) adapted to both filtrations, an intensity \(\lambda^{\mathcal G}(t)\) under \((\mathcal G_t)\) satisfies that \(N(t)-\int_0^t\lambda^{\mathcal G}(u)\,du\) is a local martingale under the enlarged filtration. Under suitable integrability conditions, the observed-history intensity is a predictable version of \(\lambda^{\mathcal F}(t)=E\{\lambda^{\mathcal G}(t)\mid\mathcal F_{t-}\}\). Here \(\mathcal F_{t-}\) denotes the observed history strictly before time \(t\).

The conditional expectation is taken at each time using the information available just before that time. It is not, in general, an average over the latent variable using its baseline distribution. Conditional on observed event history, the relevant distribution of the latent variable may have changed.

Why Conditional Expectation Gives the Observed Intensity

Suppose \(N(t)\) is adapted to the observed filtration and has enlarged-filtration intensity \(\lambda^{\mathcal G}(t)\). Assume the compensator is absolutely continuous in time, the intensity is locally integrable, and the conditional expectation has a predictable version. For any bounded, nonnegative, \((\mathcal F_t)\)-predictable process \(H(t)\), the compensator property gives

$$ E\left\{\int_0^\tau H(t)\,dN(t)\right\} = E\left\{\int_0^\tau H(t)\lambda^{\mathcal G}(t)\,dt\right\}. $$

Because \(H(t)\) uses only the observed history just before \(t\), conditioning the integrand on that history yields

$$ E\left\{\int_0^\tau H(t)\lambda^{\mathcal G}(t)\,dt\right\} = E\left\{\int_0^\tau H(t)E\{\lambda^{\mathcal G}(t)\mid\mathcal F_{t-}\}\,dt\right\}. $$

Thus the observed-history compensator has rate \(E\{\lambda^{\mathcal G}(t)\mid\mathcal F_{t-}\}\); uniqueness of the compensator identifies this as the observed-filtration intensity, up to the usual equality almost everywhere. This argument requires that the same counting process is adapted to the smaller filtration. If the event itself is not observed, its observed-filtration intensity is not obtained by this formula without first specifying an observed event process.

Worked Example: Averaging Given the Current Observed History

Suppose a person is known to be at risk at time \(t\), but a binary severity type \(Z\) is unobserved. Under the enlarged filtration, the event intensity is \(0.04\) per day for \(Z=0\) and \(0.12\) per day for \(Z=1\). Given the observed history just before \(t\), suppose the posterior probability of \(Z=1\) is \(0.25\). Assume these are the correct conditional rates while the person is at risk, and that the event count is observed.

The observed-history intensity is the conditional mean of the two enlarged-history rates:

$$ \lambda^{\mathcal F}(t) =0.75(0.04)+0.25(0.12) =0.030+0.030 =0.060\ \text{per day}. $$

The rate is \(0.060\) events per day conditional on the observed history, not the rate for either known severity type. The calculation uses the posterior probability \(0.25\) at time \(t\); using a baseline type probability instead would be justified only if it still equalled the probability given the current observed history.

Survival and Event History Update the Average

A simple latent-class model shows why marginalization often changes the shape of an intensity model. Let \(Z\) be an unobserved type with \(P(Z=0)=P(Z=1)=0.5\). Suppose the first event time has an exponential distribution conditional on type, with rates \(0.5\) and \(1.5\) per month, respectively. The enlarged filtration includes \(Z\); the observed filtration includes the event history but not \(Z\). Assume everyone is initially at risk, the event is observed, and there is no other censoring or entry process in this calculation.

Among people who have had no event by time \(t\), Bayes’ rule gives a posterior probability of the higher-rate type proportional to its prior probability times its event-free probability:

$$ P(Z=1\mid T\ge t) = \frac{0.5e^{-1.5t}}{0.5e^{-0.5t}+0.5e^{-1.5t}} = \frac{e^{-1.5t}}{e^{-0.5t}+e^{-1.5t}}. $$

Consequently, while still event-free, the observed-history intensity is the posterior mean of the type-specific rates. At \(t=0\), it is the baseline average \(0.5(0.5)+0.5(1.5)=1.0\) event per month. At later times, the higher-rate type is less common among those who remain event-free.

Worked Example: The Marginal Hazard Falls Among Survivors

Use the two-type model above and consider a person known to be event-free at \(t=2\) months. The event-free probabilities for the two types are \(e^{-0.5(2)}=e^{-1}=0.3679\) and \(e^{-1.5(2)}=e^{-3}=0.04979\).

The posterior probability of the higher-rate type is

$$ P(Z=1\mid T\ge2) = \frac{0.04979}{0.3679+0.04979} =0.1192. $$

Therefore the observed-history intensity at two months, conditional on still being event-free, is

$$ \lambda^{\mathcal F}(2) = (1-0.1192)(0.5)+(0.1192)(1.5) =0.4404+0.1788 =0.6192\ \text{events per month}. $$

This rate is below the baseline average of \(1.0\) event per month because event-free survival is more likely for the lower-rate type. It is not the rate for a person whose type is known, and it is not a constant-rate exponential model in the observed filtration. The calculation conditions on the observed fact of remaining event-free through two months.

More generally, if \(S_z(t)\) is the event-free survival probability for type \(z\), the marginal survival function is \(S(t)=\sum_zP(Z=z)S_z(t)\). Its hazard is

$$ -\frac{d}{dt}\log S(t) = \sum_z P(Z=z\mid T\ge t)\lambda_z(t), \qquad P(Z=z\mid T\ge t) = \frac{P(Z=z)S_z(t)}{\sum_jP(Z=j)S_j(t)}. $$

The equality follows by differentiating \(S(t)\): since \(S'_z(t)=-\lambda_z(t)S_z(t)\), we have \(-S'(t)/S(t)=\sum_zP(Z=z)S_z(t)\lambda_z(t)/S(t)\), which is the stated posterior-weighted average. The weights depend on \(t\) whenever types have different survival distributions. Thus even constant type-specific hazards can produce a time-varying marginal hazard.

Cause-Specific Intensities Under Coarsened Information

For multiple event types, let \(N_k(t)\) count events of type \(k\), and let \(\lambda_k^{\mathcal G}(t)\) be its enlarged-history intensity. If the event types cannot jump simultaneously and each count is observed, the corresponding observed-history intensity is

$$ \lambda_k^{\mathcal F}(t) = E\{\lambda_k^{\mathcal G}(t)\mid\mathcal F_{t-}\}. $$

Linearity of conditional expectation means the total intensity is also the sum of the observed cause-specific intensities. But the cause-specific rates can respond differently to changing beliefs about unobserved type. A constant cause-specific form conditional on type therefore need not remain constant after type is hidden.

Worked Example: Two Causes and an Unobserved Type

Suppose a person is at risk for two mutually exclusive event types. Type \(H\) has posterior probability \(0.40\) given the observed history; type \(L\) has probability \(0.60\). Conditional on type, the cause-specific rates are constant while at risk: for type \(H\), cause 1 has rate \(0.10\) and cause 2 has rate \(0.02\) per month; for type \(L\), the respective rates are \(0.02\) and \(0.08\) per month. Assume the type-specific rates and posterior probabilities are valid for the same time and risk set.

Averaging each cause separately gives

$$ \lambda_1^{\mathcal F} =0.40(0.10)+0.60(0.02) =0.052\ \text{per month}, \qquad \lambda_2^{\mathcal F} =0.40(0.02)+0.60(0.08) =0.056\ \text{per month}. $$

The total rate is \(0.052+0.056=0.108\) events per month. Given an event at that instant, the cause 1 rate is \(0.052/0.108=0.4815\) of the total instantaneous rate. These are conditional rates for this observed history; a different history that changes the posterior probability of type \(H\) generally changes both cause-specific rates and their relative contribution.

What the Averaging Does—and Does Not—Preserve

Conditional expectation preserves sums and nonnegative weighted averages, but it does not generally preserve a model’s functional form. If an enlarged-history intensity is \(\lambda^{\mathcal G}(t)=Y(t)\lambda_0(t)\exp\{\beta Z\}\), where \(Y(t)\) is an observed at-risk indicator and \(Z\) is unobserved, then while at risk

$$ \lambda^{\mathcal F}(t) = Y(t)\lambda_0(t)E\{\exp(\beta Z)\mid\mathcal F_{t-}\}. $$

The final factor is a history-dependent conditional expectation. It need not be constant, and it need not be expressible as an exponential of the same covariates with the same coefficient. Treating it as \(E\{\exp(\beta Z)\}\) ignores the information contained in the event-free and event history. Conversely, a model specified directly for the observed filtration need not imply a simple model conditional on latent variables.

This distinction also separates an observed-history intensity from a population-standardized summary. The former conditions on each person’s observed history just before time \(t\). A standardized rate or risk may instead average predictions over a chosen covariate distribution. The target, conditioning information, and averaging distribution should be stated before comparing estimates.

ApproachConditioning informationWhat it representsMain caution
Enlarged-history modelObserved history plus latent or fully measured typeRate conditional on the more complete informationNot directly an observed-history rate when type is hidden
Observed-history intensityObserved history just before time \(t\)Rate among histories represented in the observed filtrationMay depend on past event-free time and observed events through changing posterior beliefs
Standardized summaryA specified distribution over covariates or historiesAn average over the chosen target distributionChanging the averaging distribution changes the summary

Assumptions, Diagnostics, and Common Pitfalls

The projection result is useful when the observed process is genuinely adapted to the observed filtration, the enlarged intensity is well defined, and the needed conditional expectation is integrable and predictably measurable. In applications, the observation scheme must also be described correctly: if event types, covariates, or censoring indicators are not observed as assumed, the filtration and resulting intensity change. If entry or censoring affects which histories are observed, those mechanisms must be handled in the model rather than silently treated as irrelevant.

A useful diagnostic question is whether the estimated rate is compatible with the information being conditioned on. If a fitted observed-history model assumes a constant rate, inspect whether rates or residual patterns vary with time since entry, event-free duration, or recorded event history. Such patterns may indicate model misspecification or unmodeled heterogeneity, but do not by themselves prove that a particular latent type exists. Standard residual checks can reveal some lack of fit; they cannot identify an unobserved mixing distribution from the observed data without additional assumptions.

  • Wrong: “The marginal intensity is the baseline average of the type-specific intensities.” Correction: It is the conditional average given \(\mathcal F_{t-}\); the relevant type probabilities can change with observed history.
  • Wrong: “A constant conditional rate implies a constant observed-history rate.” Correction: Different latent types can be selectively represented in the histories still at risk, making the observed rate time-dependent.
  • Wrong: “The projected intensity retains the same regression coefficients.” Correction: Conditional expectation generally changes the functional form; the new intensity must be derived or modeled for the observed filtration.

A conditional model is appropriate when the scientific target is a rate given the more complete information and the relevant variables are measured or modeled. An observed-history model is appropriate when prediction or inference conditions only on the available history. A joint latent-variable model can represent both levels, but its conclusions rely on assumptions about the latent distribution and its relation to event occurrence. None of these approaches makes the other rate universally preferable; they answer different conditional questions.

Key takeaway: Coarsening information replaces an enlarged-filtration intensity by its predictable conditional expectation given the observed past. The posterior distribution used in that average can evolve with event history, so a simple conditional rate model may become a time- or history-dependent marginal intensity.

Check Your Understanding

Answer each question using the distinction between enlarged and observed histories. Brief solutions follow the questions.

  1. A hidden type has rates \(0.02\) and \(0.08\) events per day. Given observed history, the probability of the higher-rate type is \(0.30\). What is the observed-history intensity?
  2. In a two-type first-event model, why can the hazard among event-free individuals change over time even when each type has a constant hazard?
  3. If \(N_1\) and \(N_2\) count mutually exclusive causes, what is the observed-history intensity for cause \(k\) in terms of its enlarged-history intensity?
  4. Does a residual pattern suggesting a time-varying observed rate prove that an unobserved frailty variable is present? Explain.

Solutions. (1) The rate is \(0.70(0.02)+0.30(0.08)=0.014+0.024=0.038\) events per day. (2) The higher-rate type is less likely to remain event-free, so its posterior weight among survivors falls; the observed hazard is the posterior-weighted average of the type-specific hazards. (3) It is \(\lambda_k^{\mathcal F}(t)=E\{\lambda_k^{\mathcal G}(t)\mid\mathcal F_{t-}\}\), assuming the event count is observed and the projection conditions hold. (4) No. The pattern may signal lack of fit, but residuals alone do not establish a particular latent mechanism; other model or observation features may also explain it.

Clinical Biostats