Tutorials › Competing Risks and Multi-State Models › Comparing Predictions From Different Model Strategies

Absolute Risk Prediction · Tutorial 112 of 699

Comparing Predictions From Different Model Strategies

Learn how cause-specific hazard and Fine-Gray models produce absolute-risk predictions, why their results can differ, and how to compare them fairly.

Advanced 12 min read

What You'll Learn

  • Distinguish cause-specific hazards from subdistribution hazards as prediction strategies.
  • Calculate a cause-specific cumulative incidence using hazards for all competing causes.
  • Convert a Fine-Gray subdistribution hazard prediction into absolute risk.
  • Compare predictions only at matched covariate profiles and time horizons.
  • Diagnose differences, assess prediction quality, and avoid causal interpretations of model choice.

Why Compare Prediction Strategies?

A competing-risk prediction asks for the probability of experiencing a particular event by a stated time, before a different event prevents it from occurring. Two common strategies start from different hazard models: cause-specific hazard models describe each event cause’s instantaneous rate among people still free of every event, while a Fine-Gray model describes the subdistribution hazard associated with a particular cause. Both can be used to predict absolute risk, but they do not generally give identical predictions.

A fair comparison holds the prediction target, covariate profile, time horizon, and target population fixed. The comparison is about the predictive implications of two modeling strategies—not evidence that one model or hazard is a cause of the outcome. A difference in predicted risk can arise from the models’ structures, fitted coefficients, baseline estimates, or assumptions.

Definition: For event cause \(k\), the cumulative incidence function (CIF) \(F_k(t\mid\mathbf{x})\) is the probability of experiencing cause \(k\) by time \(t\), conditional on covariates \(\mathbf{x}\), in the presence of the competing event causes. A cause-specific hazard \(\lambda_j(t\mid\mathbf{x})\) is the instantaneous rate of cause \(j\) at time \(t\), among individuals who have experienced no event by that time. The subdistribution hazard for cause \(k\), \(\alpha_k(t\mid\mathbf{x})\), is defined by the rate of increase in \(F_k\) relative to the probability \(1-F_k\) of not yet having experienced cause \(k\).

How Each Strategy Produces Absolute Risk

A cause-specific strategy models the cause-specific hazard for every relevant cause. Let \(\Lambda_j(t\mid\mathbf{x})=\int_0^t\lambda_j(u\mid\mathbf{x})\,du\) be the cumulative cause-specific hazard for cause \(j\). Under continuous event-time distributions, the probability of remaining free of every event through time \(t\) is

$$ S(t\mid\mathbf{x})= \exp\left\{-\sum_j\Lambda_j(t\mid\mathbf{x})\right\}. $$

The cause-\(k\) CIF is then obtained by accumulating cause-\(k\) hazard only while a person remains event-free:

$$ F_k(t\mid\mathbf{x})= \int_0^t S(u\mid\mathbf{x})\lambda_k(u\mid\mathbf{x})\,du. $$

This formula shows why a cause-\(k\) hazard model alone is not enough to obtain the cause-\(k\) absolute risk: the hazards of the competing causes also affect \(S(u\mid\mathbf{x})\), and hence the number who remain at risk of cause \(k\).

A Fine-Gray strategy models the subdistribution hazard for the target cause. In its proportional-hazards form, \(\alpha_k(t\mid\mathbf{x})=\alpha_{0k}(t)\exp(\boldsymbol{\gamma}_k^\mathsf{T}\mathbf{x})\). If \(A_k(t\mid\mathbf{x})=\int_0^t\alpha_k(u\mid\mathbf{x})\,du\), the corresponding risk prediction is

$$ F_k(t\mid\mathbf{x})=1-\exp\{-A_k(t\mid\mathbf{x})\} =1-\exp\left\{-A_{0k}(t)\exp(\boldsymbol{\gamma}_k^\mathsf{T}\mathbf{x})\right\}. $$

The subdistribution hazard is tied directly to the target cause’s CIF. It is not interchangeable with a cause-specific hazard, and a Fine-Gray model for one cause does not by itself provide a complete model for every other cause. In particular, independently fitted Fine-Gray models for multiple causes need not produce risks that sum to at most one; check probability coherence rather than assuming it.

Assumptions Before Comparing Predictions

A model comparison is interpretable only if both predictions refer to the same well-defined quantity. Specify the event cause, time origin, covariate values, prediction horizon, and population of interest. For example, “the risk of cause 1 by five years for a person with \(\mathbf{x}=(1,0)\)” is more precise than “the risk for the treated group.” If the goal is an average over a population rather than risk at a single profile, calculate each strategy’s predictions for the same individuals and average them using the same target-population weights.

Both strategies also rely on sound outcome and follow-up information. Event causes must be classified correctly, and censoring must be independent of the event process conditional on the variables used in the analysis, or otherwise handled by an appropriate method. Prediction at a horizon requires adequate follow-up support there; extending either model past the observed support is extrapolation, not a direct comparison of observed-range prediction performance.

A proportional cause-specific hazards model assumes that each modeled cause-specific hazard has the stated covariate effect over time, unless time-varying effects are included. A proportional Fine-Gray model makes the corresponding proportionality assumption for the subdistribution hazard. Each also requires suitable baseline hazard estimation and an adequate covariate functional form. These assumptions are distinct: satisfying one does not imply that the other holds.

Worked Comparisons

Worked Example: Constant Cause-Specific Hazards Versus Fine-Gray Risk

Consider a hypothetical fitted cause-specific strategy for a fixed covariate profile \(B\). Suppose two causes have constant hazards \(\lambda_1=0.08\) per year and \(\lambda_2=0.12\) per year. Compare the predicted risk of cause 1 by five years with a hypothetical Fine-Gray prediction for the same profile and horizon, whose cumulative subdistribution hazard is \(A_1(5\mid B)=0.36\).

The total cause-specific hazard is \(\lambda_{\mathrm{tot}}=0.08+0.12=0.20\) per year. Since this total is positive, the event-free survival is \(S(u\mid B)=\exp(-0.20u)\), and the cause-1 CIF is

$$ F_{1,\mathrm{CS}}(5\mid B) =\int_0^5 0.08\exp(-0.20u)\,du =\frac{0.08}{0.20}\{1-\exp(-0.20\times 5)\} =0.40\{1-\exp(-1)\} \approx 0.2528. $$

The cause-specific strategy predicts about 25.28% risk by five years. The Fine-Gray prediction is

$$ F_{1,\mathrm{FG}}(5\mid B)=1-\exp(-0.36)\approx 0.3023. $$

For this same profile and horizon, the Fine-Gray prediction is about 30.23%, or 4.95 percentage points higher. This numerical difference alone does not identify which prediction is better: that requires assessing predictive performance and the plausibility of the respective model assumptions.

Worked Example: Comparing Predictions When Hazards Change Over Time

Suppose, hypothetically, that for profile \(C\), cause-specific hazards are constant within each of two intervals. From year 0 to year 2, let \(\lambda_1=0.05\) and \(\lambda_2=0.10\) per year. From year 2 to year 4, let \(\lambda_1=0.10\) and \(\lambda_2=0.05\) per year. Compare the cause-1 risk by year 4 with a Fine-Gray prediction having \(A_1(4\mid C)=0.25\).

In the first two years the total hazard is 0.15 per year, so event-free survival at year 2 is \(\exp(-0.15\times2)=\exp(-0.30)\approx0.7408\). The cause-1 incidence accumulated in that interval is

$$ \frac{0.05}{0.15}\{1-\exp(-0.30)\}\approx0.08639. $$

In the second interval the cause-1 incidence contribution must also account for surviving event-free to year 2. It is

$$ \exp(-0.30)\frac{0.10}{0.15}\{1-\exp(-0.30)\} \approx0.1280. $$

Adding the two contributions gives \(F_{1,\mathrm{CS}}(4\mid C)\approx0.08639+0.1280=0.2144\), or 21.44%. The Fine-Gray prediction is \(F_{1,\mathrm{FG}}(4\mid C)=1-\exp(-0.25)\approx0.2212\), or 22.12%. The predictions differ by about 0.680 percentage points. The calculation illustrates that a cause-specific prediction must account for the timing of all causes, not merely compare hazards at the horizon.

Worked Example: The Zero-Total-Hazard Boundary Case

For a hypothetical profile \(D\), suppose every cause-specific hazard is zero throughout the first three years. Then \(\lambda_{\mathrm{tot}}=0\), event-free survival is \(S(u\mid D)=1\), and the cause-1 hazard is also zero. Directly using the CIF integral gives

$$ F_1(3\mid D)=\int_0^3 1\times0\,du=0. $$

The positive-total-hazard shortcut \(\frac{\lambda_1}{\lambda_{\mathrm{tot}}}\{1-\exp(-\lambda_{\mathrm{tot}}t)\}\) must not be used here: it would divide by zero. More generally, with constant hazards, that shortcut applies when \(\lambda_{\mathrm{tot}}>0\); when \(\lambda_{\mathrm{tot}}=0\), all cause-specific hazards are zero and the CIF is zero. A comparison with another strategy should report this boundary case as zero risk, not as an undefined ratio.

Why the Predictions Can Differ

The cause-specific approach builds the target CIF from the joint pattern of all cause-specific hazards. The Fine-Gray approach models the target CIF through its subdistribution hazard. Even when both are fitted to the same data, their regression structures and assumptions differ. Different covariate effects, baseline estimates, or departures from proportionality can therefore lead to different risks for the same profile and horizon. Neither strategy is guaranteed to be numerically closer to the truth merely because it models a particular type of hazard.

One useful mathematical link makes the distinction clear. For a continuous subdistribution hazard, its definition gives \(dF_k(t\mid\mathbf{x})=\{1-F_k(t\mid\mathbf{x})\}\alpha_k(t\mid\mathbf{x})\,dt\). Dividing by \(1-F_k(t\mid\mathbf{x})\) and integrating from 0 to \(t\), using \(F_k(0\mid\mathbf{x})=0\), yields

$$ -\log\{1-F_k(t\mid\mathbf{x})\} =\int_0^t\alpha_k(u\mid\mathbf{x})\,du =A_k(t\mid\mathbf{x}). $$

Rearranging gives \(F_k(t\mid\mathbf{x})=1-\exp\{-A_k(t\mid\mathbf{x})\}\). This identity explains how a Fine-Gray model’s cumulative subdistribution hazard becomes an absolute-risk prediction. It does not make that hazard the same as a cause-specific hazard.

Diagnostics, Alternatives, and Limitations

Before interpreting a gap between predictions, check that it is not caused by a mismatched target. Confirm the cause label, time origin, horizon, covariate values, and population weights. Then inspect whether each model’s assumptions are credible and whether predictions behave sensibly over time. For example, a CIF should not decrease as the horizon increases, and a predicted probability must lie between zero and one. Check multiple causes jointly when predictions are meant to describe a coherent set of competing outcomes.

Where data permit, evaluate each strategy’s calibration and predictive error at the same horizon in the same target population, using appropriate methods for censoring. A comparison of apparent fit in the data used to fit the models is not, by itself, evidence of better prediction in new individuals. Uncertainty in the difference matters too: predictions from models fitted to the same sample are dependent, so a paired resampling procedure that refits both strategies can preserve that dependence when estimating uncertainty.

A nonparametric competing-risk estimate can be a useful reference for an unadjusted population or group, whereas regression strategies provide covariate-specific predictions under their model assumptions. It does not automatically answer an adjusted prediction question. Conversely, regression predictions depend on correct model specification and adequate data support. No single strategy is uniformly best; the appropriate comparison depends on the prediction target, assumptions, and validation evidence.

A common mistake is to say, “The Fine-Gray model gives a higher risk, so the covariate causes more events under that model,” or to infer that one hazard ratio is the causal effect because its risk prediction is larger. This reasoning is wrong: model choice is not a causal intervention, and hazard parameters from the two strategies describe different quantities. The correction is to report the two absolute-risk predictions at the matched profile and horizon, quantify their difference and uncertainty, then assess predictive accuracy and assumptions without assigning a causal interpretation to the model contrast.

Key takeaway: Compare cause-specific and Fine-Gray predictions only after matching the event, covariate profile, horizon, and target population. The two strategies construct absolute risk differently; a discrepancy is a prompt for diagnosis and validation, not a causal conclusion.

Check Your Understanding

Answer the questions, then check the brief solutions below.

  1. Why can a cause-\(k\) cause-specific hazard model alone be insufficient to predict the cause-\(k\) CIF?
  2. For a Fine-Gray prediction with cumulative subdistribution hazard \(A_k(t\mid\mathbf{x})=0.20\), calculate the predicted risk.
  3. With constant cause-specific hazards \(\lambda_1=0.06\) and \(\lambda_2=0.09\) per year, calculate the cause-1 risk by two years.
  4. What does a higher Fine-Gray prediction than a cause-specific prediction establish about causation?

Solutions. (1) The other causes affect event-free survival, which weights the cause-\(k\) hazard in the CIF integral. (2) \(1-\exp(-0.20)\approx0.1813\), or 18.13%. (3) The total hazard is \(0.15\) per year, so \(F_1(2)=\frac{0.06}{0.15}\{1-\exp(-0.15\times2)\}=0.40\{1-\exp(-0.30)\}\approx0.1037\), or 10.37%. (4) Nothing about causation follows from that difference; it is a predictive contrast between strategies that requires model assessment and validation.

Clinical Biostats