Tutorials › Competing Risks and Multi-State Models › Rate Ratios Versus Event-Hazard Ratios

Recurrent-Event Estimands and Models · Tutorial 596 of 699

Rate Ratios Versus Event-Hazard Ratios

Learn what each ratio compares, how to calculate and interpret it, and why a marginal recurrent-event rate ratio need not equal a conditional event-intensity hazard ratio.

Advanced 12 min read

What You'll Learn

  • Define a marginal recurrent-event rate using event counts and person-time.
  • Define a conditional event intensity in terms of a participant’s event history and current risk status.
  • Interpret rate ratios and event-intensity hazard ratios with their distinct denominators and conditioning sets.
  • Derive why averaging conditional intensities need not preserve a common hazard ratio.
  • Identify assumptions and diagnostics relevant to each type of comparison.
  • Recognize when differing event-history compositions can make marginal and conditional ratios diverge.

Two Ratios, Two Questions

A rate ratio and an event-hazard ratio may both compare how often recurrent events occur in two groups, but they do not generally answer the same question. A marginal recurrent-event rate ratio compares event counts relative to accumulated person-time at the population level. A conditional event-intensity hazard ratio compares instantaneous event intensities for individuals at risk with specified observed histories. The denominator and the conditioning set determine which interpretation is justified.

This distinction matters when participants have different numbers or timings of earlier events. Their histories can affect later event intensity, and the groups being compared may spend different amounts of follow-up time in different histories. Consequently, a population-average rate ratio can differ from a history-specific intensity ratio, even when the latter is constant across histories.

Definition: Let \(N_i(t)\) be the number of recurrent events experienced by participant \(i\) by time \(t\). A marginal event rate over an interval is the expected number of events in the target population divided by the corresponding expected person-time under observation. Its units are events per person-time, such as events per person-year. A conditional event intensity, denoted \(\lambda_i(t\mid\mathcal{H}_i(t-))\), is the instantaneous rate of an event at time \(t\), given participant \(i\)'s observed event history \(\mathcal{H}_i(t-)\) just before \(t\) and that participant's being under observation and at risk.

Here, “at risk” means that another event can be counted at that time; it does not mean that a participant must be event-free. A recurrent-event intensity is not simply the hazard of a first event. A participant can remain at risk after an event, with their updated history included in the conditioning information.

What Is in Each Denominator?

For a group followed over an interval, the descriptive rate is total events divided by total observed person-time. Someone observed for two years contributes twice as much person-time as someone observed for one year. An individual can contribute multiple events, and each event contributes to the numerator. The resulting group rate describes the average event frequency per unit of observed time; it does not, by itself, specify how event intensity varies with an individual’s past history.

A conditional intensity has a different denominator in its interpretation: it describes the instantaneous event rate among participants under observation and at risk, conditional on their specified histories and covariates. A model might compare two covariate values for people with the same modeled history. For example, a proportional-intensity model can be written as follows, where \(X_i\) is a group indicator and \(\lambda_0\) may depend on the history:

$$ \lambda_i(t\mid\mathcal{H}_i(t-),X_i) = \lambda_0(t\mid\mathcal{H}_i(t-))\exp(\beta X_i). $$

Under this model, \(\exp(\beta)\) is the conditional event-intensity ratio for \(X_i=1\) versus \(X_i=0\), at the same time and modeled history. Its interpretation requires the model’s intensity structure to be appropriate; it is not automatically a ratio of total events per person-time.

A marginal model for recurrent-event counts instead relates a population-average count, or a mean count per unit of exposure, to covariates. For instance, a log-mean model with a person-time offset may use \(\log E[N_i]=\log(T_i)+\alpha+\beta X_i\), where \(T_i\) is observed person-time. In that model, \(\exp(\beta)\) compares modeled mean event rates for the covariate groups, conditional on any other covariates included in the mean model. It does not condition on each participant’s full event history in the way an event-intensity model does. GEE methods address the correlation among repeated observations when estimating a marginal mean model; that correlation adjustment does not turn its rate-ratio coefficient into a conditional intensity ratio.

Worked Example: Calculating a Marginal Recurrent-Event Rate Ratio

Suppose, hypothetically, that one group has 60 recurrent events over 300 person-years of observed follow-up, while a comparison group has 45 events over 300 person-years. Assume the event counts are complete for the recorded person-time and that both groups use the same event definition and follow-up rules.

The group rates are:

$$ \widehat r_1=\frac{60\text{ events}}{300\text{ person-years}} =0.200\text{ events per person-year}, \qquad \widehat r_0=\frac{45\text{ events}}{300\text{ person-years}} =0.150\text{ events per person-year}. $$

The marginal rate ratio is:

$$ \widehat{RR}_{\mathrm{rate}} = \frac{\widehat r_1}{\widehat r_0} = \frac{0.200}{0.150} = 1.333. $$

In this hypothetical data, the first group experienced 1.333 times as many events per observed person-year as the comparison group: its average event rate was 33.3% higher. This calculation does not establish that the instantaneous intensity was 1.333 times as high for every participant at the same event history. The groups could have spent different proportions of their person-time in histories associated with different event rates.

Why Averaging Can Change the Ratio

The distinction can be made precise using the defining property of an intensity. Let \(Y_i(t)\) indicate that participant \(i\) is under observation and at risk, and write \(\lambda_i(t\mid\mathcal{H}_i(t-),X_i)\) for the intensity when \(Y_i(t)=1\). For a sufficiently small time interval \(dt\), the expected event increment, conditional on the observed history, is approximately the intensity multiplied by \(dt\). Averaging over histories in a group gives:

$$ E[dN_i(t)\mid X_i=x] \approx E\!\left[ Y_i(t)\lambda_i(t\mid\mathcal{H}_i(t-),x) \mid X_i=x \right]dt. $$

Dividing by the expected person-time contribution, \(E[Y_i(t)\mid X_i=x]dt\), gives the group’s marginal rate at that time:

$$ r_x(t) = \frac{ E\!\left[Y_i(t)\lambda_i(t\mid\mathcal{H}_i(t-),x)\mid X_i=x\right] }{ E[Y_i(t)\mid X_i=x] }. $$

Thus, the marginal rate is an average of conditional intensities over the histories represented among the group’s observed person-time. If those history distributions differ between groups, the two averages need not retain a common history-specific ratio. In particular, a constant conditional ratio \(\exp(\beta)\) does not guarantee that \(r_1(t)/r_0(t)=\exp(\beta)\) when the groups average over different histories. If the histories are standardized to the same distribution and the conditional ratio is constant over those histories, then the standardized marginal rates have that same ratio. This is a special condition, not an automatic property of crude rates.

Worked Example: A Conditional Event-Intensity Ratio

Suppose, hypothetically, an event-intensity model gives a treatment coefficient \(\beta=0.5\). At a particular time, the reference intensity for a participant with a specified history is 0.120 events per person-year. Assume the model is appropriate for the compared covariate values, and that this participant is under observation and at risk in either group.

The model’s conditional intensity ratio is:

$$ \exp(\beta)=\exp(0.5)=1.648721\ldots. $$

The corresponding treatment-group intensity, calculated without rounding the multiplier, is:

$$ 0.120\exp(0.5) = 0.1978465\ldots \text{ events per person-year} \approx 0.1978 \text{ events per person-year}. $$

The unrounded ratio is \(\exp(0.5)=1.648721\ldots\), or 1.649 to four significant figures. This is an instantaneous comparison for individuals with the same modeled history and at-risk status. It is not, solely on the basis of this model calculation, a claim that the treatment group has 1.649 times the total event count per person-year in the population.

Worked Example: Different History Distributions, Different Marginal Ratio

Consider hypothetical person-time divided into two history categories: a lower-intensity history and a higher-intensity history. Suppose that within each category the treatment-group intensity is exactly twice the comparison-group intensity. Assume each reported category rate applies to the person-time in that category, and that event counts and person-time are measured consistently.

Group and historyPerson-timeRate (events per person-year)Events
Comparison, lower-intensity history800 person-years0.1080
Comparison, higher-intensity history200 person-years0.50100
Treatment, lower-intensity history200 person-years0.2040
Treatment, higher-intensity history800 person-years1.00800

Within the lower-intensity history, the rate ratio is \(0.20/0.10=2\); within the higher-intensity history, it is \(1.00/0.50=2\). But the overall rates weight those history-specific rates by the person-time distribution in each group. For the comparison group:

$$ \widehat r_0 = \frac{80+100}{800+200} = \frac{180}{1000} = 0.180 \text{ events per person-year}. $$

For the treatment group:

$$ \widehat r_1 = \frac{40+800}{200+800} = \frac{840}{1000} = 0.840 \text{ events per person-year}. $$

The marginal rate ratio is \(0.840/0.180=4.667\), not 2. The treatment group contributes a larger proportion of its person-time in the higher-intensity history category. This example demonstrates a difference in weighting, not evidence that the within-history ratio changed. In an actual study, history distributions may be related to treatment, follow-up, and participants’ underlying event processes; interpreting a conditional ratio causally requires additional design and modeling assumptions.

Assumptions, Diagnostics, and Alternatives

A rate comparison is interpretable only if its numerator and denominator refer to compatible populations and periods. Check that recurrent events are defined and counted consistently, person-time is calculated consistently, and the rules for starting and ending follow-up are clear. If censoring or loss to follow-up depends on event risk, a simple events-per-observed-person-time calculation may not represent the intended population rate; appropriate modeling or adjustment may be needed. A model-based rate ratio also relies on an adequate mean model for the covariates and exposure time. Robust standard errors can address some misspecification of within-person correlation, but they do not repair a misspecified mean or biased person-time.

An event-intensity ratio requires a clearly defined risk process and a history set that is available just before each event time. Covariates used to predict the intensity must be appropriately timed, and the model must describe the conditional intensity adequately. If a constant ratio is claimed, the model must support that constancy over the relevant times and histories. A large difference between a crude rate ratio and a conditional ratio can be a prompt to inspect history distributions, follow-up patterns, and model definitions; it is not, by itself, proof that either estimate is wrong.

The main alternatives target different summaries. A marginal rate model, including a GEE count model with a person-time offset, targets a population-average mean rate and accommodates repeated observations through its estimating approach. A conditional event-intensity model targets instantaneous rates given specified histories. Neither is universally preferable: choose according to whether the scientific question concerns average event burden per unit follow-up or event intensity among people with comparable histories. They may also differ in covariates, time scales, risk sets, and assumptions, so comparing their coefficients as though they were interchangeable can be misleading.

1
Name the estimand.
State whether the target is events per person-time in a population or an instantaneous event intensity conditional on history.
2
Write the denominator and conditioning set.
For a marginal rate, identify the person-time contributing to the rate. For a conditional intensity, specify the at-risk status and the history being held fixed.
3
Check who contributes follow-up.
Compare follow-up rules and, where relevant, the distributions of histories across the groups being averaged.
4
Report the ratio in context.
State the event type, comparison groups, time unit, and whether the interpretation is marginal or conditional.

Common Mistakes and Limitations

  • Calling every event-rate ratio a hazard ratio. Events divided by person-time give a rate comparison; a conditional event-intensity ratio additionally specifies an instantaneous risk set and history.
  • Assuming equal conditional ratios force equal marginal ratios. Different history distributions can give the groups different weights when conditional intensities are averaged.
  • Ignoring the person-time definition. Different follow-up or censoring rules can change both the denominator and the population represented by the rate.
  • Interpreting an intensity ratio as a causal effect without justification. A conditional model coefficient is an association under its model unless the study design and assumptions support a causal interpretation.
  • Using a robust variance estimator as a cure-all. Robust variance estimation may protect uncertainty estimates against some correlation misspecification; it does not make an inappropriate estimand, biased follow-up, or incorrect mean model valid.

These ratios are summaries, not full descriptions of recurrent-event experience. A single rate ratio can hide when events occur, how event histories evolve, and how much participants differ. Conversely, a conditional intensity ratio is tied to its conditioning set and may not summarize the population’s average burden. Good reporting makes those boundaries explicit.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. A group has 36 events over 240 person-years and a comparison group has 30 events over 250 person-years. Calculate the marginal rate ratio and state its interpretation.
  2. In an intensity model, \(\beta=\log(1.5)\). What does \(\exp(\beta)\) compare, and what information must be held fixed for the conditional interpretation?
  3. Two groups have the same conditional event-intensity ratio within each history category. Give one reason their crude marginal rate ratios could differ from that common conditional ratio.
  4. Why does a robust variance estimator in a GEE count model not make its rate-ratio coefficient a conditional event-intensity ratio?

Solutions. (1) The rates are \(36/240=0.150\) and \(30/250=0.120\) events per person-year. Their ratio is \(0.150/0.120=1.25\). The first group’s average rate is 1.25 times the comparison group’s rate over the stated observed person-time. (2) \(\exp(\beta)=1.5\) is the modeled instantaneous intensity ratio for the compared covariate values, conditional on the same specified event history, time, and at-risk status. (3) The groups may contribute different proportions of person-time from history categories with different baseline intensities, so the marginal averages use different weights. (4) Robust variance estimation concerns uncertainty under within-person correlation; the GEE mean coefficient still targets the specified population-average count or rate relationship and does not condition on each person’s event history as an intensity model does.

Clinical Biostats