Tutorials › Competing Risks and Multi-State Models › Risk Differences Between Covariate Profiles

Absolute Risk Prediction · Tutorial 98 of 699

Risk Differences Between Covariate Profiles

Calculate a fixed-horizon difference between two predicted cumulative incidences, interpret its direction and scale, and recognize what the contrast does—and does not—say.

Advanced 10 min read

What You'll Learn

  • Define a cumulative-incidence risk difference for two specified covariate profiles and a fixed horizon.
  • Calculate the contrast from fitted profile-specific predictions.
  • Interpret positive, negative, and zero differences in probability and percentage-point units.
  • Explain why competing-event predictions affect the risk difference for the event of interest.
  • Distinguish a profile contrast from a standardized population contrast and a causal effect.
  • Check whether the model predictions and covariate profiles support a useful comparison.

Comparing Predicted Risks on an Absolute Scale

A predicted risk can be easier to interpret when it is compared with another predicted risk on the same probability scale. For example, a model might predict the probability of an event of interest by three years for two specified covariate profiles. Subtracting these predictions gives an absolute contrast: the difference in predicted cumulative incidence. Its size is expressed as a difference in probability, or equivalently in percentage points.

The contrast answers a focused prediction question: at the chosen horizon, how far apart are the model’s predicted cumulative incidences for these two profiles? It does not by itself describe an average difference in a population, and it does not show what would happen if a person’s covariates were changed. Those distinctions matter especially when the profiles differ in several characteristics.

Definition: Let \(F_k(t\mid z)\) be the cumulative incidence of event type \(k\) by time \(t\), conditional on covariate profile \(z\). For two specified profiles \(z_A\) and \(z_B\), the risk difference at fixed horizon \(t\) is $$ \Delta_k(t;z_A,z_B)=F_k(t\mid z_A)-F_k(t\mid z_B). $$ If fitted predictions are used, the reported contrast is $$ \widehat{\Delta}_k(t;z_A,z_B)=\widehat{F}_k(t\mid z_A)-\widehat{F}_k(t\mid z_B). $$ The fitted risk difference reported here is a point estimate of the corresponding contrast. Its precision and the uncertainty in the fitted predictions require additional methods.

The subtraction order is part of the definition: profile A minus profile B. A positive value means the predicted cumulative incidence is higher for profile A; a negative value means it is lower. A value of zero means the two fitted predictions are equal at that horizon, not that every feature of their risk processes is the same.

What the Contrast Represents

A cumulative incidence is the probability of experiencing event type \(k\) by time \(t\) before a competing event occurs. It is therefore not interchangeable with a cause-specific hazard, which concerns the instantaneous event rate among people still event-free. A risk difference compares two cumulative-incidence predictions, incorporating how both the event of interest and competing events affect the probability of remaining event-free.

The profiles must be specified clearly. They may differ in one covariate or in several. If they differ in several, the subtraction gives a contrast between the two complete profiles; it does not isolate the contribution of any one covariate. Nor is it automatically an average across people. A standardized risk difference instead compares risks averaged over a chosen target population. The profile contrast and standardized contrast answer different questions.

A model-based contrast is also not automatically causal. A causal interpretation would require a causal question and additional assumptions appropriate to that question. In ordinary prediction, the profiles are descriptions used to obtain conditional model predictions. Their difference is a model-based association on the absolute-risk scale, not an effect of intervening to move an individual from one profile to another.

Worked Example: Calculating a Difference From Discrete-Time Predictions

Consider a hypothetical fitted prediction model with two event types: event \(k\), the event of interest, and a competing event. For two successive intervals, suppose the model predicts the following conditional probabilities of each event among people who are event-free at the start of that interval. Within an interval, the event types are mutually exclusive.

ProfileIntervalEvent \(k\)Competing event
A10.120.08
A20.200.10
B10.080.12
B20.100.15

For a profile, the probability of remaining event-free through an interval is one minus the two event probabilities in that interval. Thus profile A has probability \(1-0.12-0.08=0.80\) of remaining event-free after interval 1. Its predicted cumulative incidence for event \(k\) through interval 2 is $$ \widehat{F}_k(2\mid z_A)=0.12+(0.80)(0.20)=0.28. $$ The first term is the probability of event \(k\) in interval 1. The second is the probability of remaining event-free through interval 1 and then experiencing event \(k\) in interval 2.

For profile B, the probability of remaining event-free after interval 1 is \(1-0.08-0.12=0.80\), giving $$ \widehat{F}_k(2\mid z_B)=0.08+(0.80)(0.10)=0.16. $$ With profile A first in the subtraction, the fitted risk difference is $$ \widehat{\Delta}_k(2;z_A,z_B)=0.28-0.16=0.12. $$ Thus, in this hypothetical model, profile A has a predicted cumulative incidence 0.12 higher than profile B through interval 2, or 12 percentage points higher. This is a comparison of the two predicted risks; it is not a causal effect of changing a person’s profile.

The competing-event predictions are part of the calculation. For example, they determine the probability of being event-free and therefore still at risk of event \(k\) in the next interval. Ignoring them would not produce these cumulative-incidence predictions.

Why Competing Events Matter to the Difference

The discrete-time recursion makes explicit why a difference in cumulative incidence is not simply a difference in event-\(k\) probabilities at each time. Let \(S_{j-1}(z)\) denote the probability of remaining free of all event types through the start of interval \(j\), and let \(p_{k,j}(z)\) be the conditional probability of event \(k\) in interval \(j\). Then $$ F_{k,j}(z)=F_{k,j-1}(z)+S_{j-1}(z)p_{k,j}(z). $$ The increment in cumulative incidence depends on both the conditional probability of event \(k\) and the probability of still being event-free.

For two profiles, subtract their recursions. Adding and subtracting \(S_{j-1}(z_A)p_{k,j}(z_B)\) gives $$ \begin{aligned} F_{k,j}(z_A)-F_{k,j}(z_B) ={}&F_{k,j-1}(z_A)-F_{k,j-1}(z_B)\\ &+S_{j-1}(z_A)\{p_{k,j}(z_A)-p_{k,j}(z_B)\}\\ &+\{S_{j-1}(z_A)-S_{j-1}(z_B)\}p_{k,j}(z_B). \end{aligned} $$

The second line reflects a difference in the interval-specific event-\(k\) probabilities, weighted by profile A’s event-free probability. The third reflects a difference in the probabilities of still being event-free. Because competing events affect event-free probability, differences in competing-event occurrence can contribute to the risk difference for event \(k\), even when the interval-specific event-\(k\) probabilities are similar.

Worked Example: The Risk Difference Can Change With the Horizon

Suppose a hypothetical fitted model predicts the following cumulative incidences for event \(k\) at two horizons. The competing-event probabilities are included to show that each profile’s predictions concern a full competing-risks process.

HorizonProfile A: event \(k\)Profile A: competing eventProfile B: event \(k\)Profile B: competing event
1 year0.080.060.110.09
3 years0.240.180.190.28

At one year, profile A minus profile B is $$ \widehat{\Delta}_k(1;z_A,z_B)=0.08-0.11=-0.03. $$ The model predicts a cumulative incidence 3 percentage points lower for profile A at one year. At three years, $$ \widehat{\Delta}_k(3;z_A,z_B)=0.24-0.19=0.05. $$ The model predicts a cumulative incidence 5 percentage points higher for profile A at three years. Both profiles’ event-\(k\) cumulative incidences increase over time, as cumulative incidence should, but the contrast changes direction. A conclusion about which profile has the higher predicted risk must therefore name the horizon.

These predictions are coherent with the displayed competing-event probabilities: at each horizon, the two cumulative incidences sum to no more than one for either profile. The values are hypothetical model predictions, not observed event counts. The example illustrates why a risk difference is tied to a specified time horizon rather than being a single time-independent feature of two profiles.

Scale, Direction, and Profile Choice

A risk difference is on an absolute probability scale. A difference of \(0.04\) means four percentage points, not a 4% relative increase. When useful, it can also be translated to a scale of 100 comparable profiles: a difference of \(0.04\) corresponds to four more predicted events per 100 profiles by the horizon, under the model’s conditional predictions. That wording describes predicted risks, not a promise about the number of events in a particular group of 100 people.

The profiles should be stated in enough detail that a reader can tell what is being compared. For continuous covariates, report the values and units. For categorical covariates, identify the category. If several characteristics differ, say so rather than suggesting that the risk difference is attributable to just one of them.

Worked Example: Interpreting a Contrast Between Multivariable Profiles

Suppose a hypothetical model predicts a 5-year cumulative incidence of \(0.31\) for event \(k\) for profile A and \(0.22\) for profile B. Profile A represents a 70-year-old person with covariate category X; profile B represents a 60-year-old person with category Y. Assume the profiles are within the range of covariate values represented in the data.

The estimated contrast, in the specified order, is $$ \widehat{\Delta}_k(5;z_A,z_B)=0.31-0.22=0.09. $$ The fitted model predicts a 5-year cumulative incidence 0.09 higher for profile A than for profile B, or 9 percentage points higher. Reversing the subtraction gives \(0.22-0.31=-0.09\), so the sign changes while the magnitude remains 9 percentage points.

Because age and covariate category both differ, the 9-percentage-point contrast describes the two complete profiles. It does not identify the contribution of age separately from category, and it does not establish that changing either characteristic would change risk by a particular amount. The fitted risk difference is a point estimate; its precision and the uncertainty in the fitted predictions require additional methods.

Assumptions, Diagnostics, and Alternatives

A profile risk difference is only as credible as the predictions being subtracted. Its interpretation as a model-based prediction relies on a clearly defined event of interest, correctly recorded competing-event status, and a fitted model that can provide meaningful cumulative-incidence predictions for both profiles at the chosen horizon. The model’s estimation must also handle censoring appropriately; for example, censoring assumptions must be suitable for the fitting and prediction procedure used.

Check that the two profiles are well defined and supported by the observed covariate data. A profile far outside the data’s range can require extrapolation, so a numerical contrast may be precise-looking but weakly supported. Also check that predicted cumulative incidences are nonnegative, do not decrease over time, and, across event types, do not sum to more than one. Investigate coding, model specification, and prediction implementation if these checks fail.

A large or small point contrast does not by itself show how stable the estimate is. Uncertainty from estimating model parameters and baseline quantities is not represented by the subtraction alone. Formal uncertainty intervals or tests require additional methods. Nor does a nonzero difference establish practical importance: its relevance depends on the event, time horizon, setting, and consequences of the outcome.

The main alternatives address different comparison questions. A standardized risk difference averages each profile’s or group’s predictions over a stated target population; it is not the same as comparing two single covariate profiles. A subgroup-specific contrast refers to a defined subgroup and should not be confused with an overall population contrast. Relative measures such as risk ratios express a contrast on a different scale, while hazard contrasts concern a different quantity from cumulative incidence. Choose the measure that matches the question, and do not substitute one interpretation for another.

Key takeaway: At a fixed horizon, subtract the predicted cumulative incidence for profile B from that for profile A. Report the subtraction order, horizon, event, profiles, and absolute units. Interpret the result as a model-based conditional contrast—not as an average population difference or a causal effect.

Check Your Understanding

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

  1. At 2 years, a model predicts event-\(k\) cumulative incidences of 0.18 for profile A and 0.25 for profile B. What is the risk difference for A minus B, and how should its sign be interpreted?
  2. Why can a difference in competing-event probabilities affect a risk difference for event \(k\)?
  3. Two profiles differ in age and disease category. What can and cannot be concluded from their predicted risk difference?
  4. At one horizon, profile A has a higher predicted cumulative incidence than profile B, but at a later horizon it has a lower prediction. Is that necessarily a calculation error?
  5. Name two checks that can reveal problems with a profile-based risk difference.

Solutions. (1) \(0.18-0.25=-0.07\): profile A’s predicted risk is 7 percentage points lower at 2 years. (2) Competing events change the probability of remaining event-free and therefore the probability of subsequently experiencing event \(k\). (3) The contrast compares the two complete profiles; it does not isolate age or disease category, and it is not automatically causal. (4) No. The contrast can change with the horizon, provided the underlying predictions are coherent. (5) For example, check that cumulative incidence is nonnegative and nondecreasing, that predictions across causes sum to no more than one, and that profiles are supported by the observed data.

Back to Tutorials
Clinical Biostats