Why Integrate AUC Across Horizons?
A transition-specific AUC describes how well a prediction score ranks cases above controls at a particular horizon. But a model may rank well at short horizons and less well at longer ones. Reporting a single horizon can therefore hide variation across the period in which predictions will be used. An integrated AUC summarizes discrimination across a prespecified range of horizons; its interpretation depends on both that range and the weights assigned to different times.
This summary is not a universal property of a model. It is a weighted average of horizon-specific discrimination for a chosen prediction origin, transition, target population, horizon interval, and definition of cases and controls. Changing any of those choices can change the estimand, not just its numerical estimate.
The case-control definition must be held fixed across horizons and follow the chosen competing-transition convention. For example, competing transitions may be treated as controls or handled another prespecified way; those choices define different discrimination targets. The scores must also refer to the same target transition and prediction origin. AUC values based on different populations or incompatible case-control definitions should not be combined as though they described one curve.
Choosing the Horizon Weights
The weight function expresses which parts of the horizon range should contribute most to the summary. With \(w(h)=1\), every unit of elapsed time receives equal weight. This is a natural choice when the entire interval is of equal interest, but it does not give equal weight to a set of irregularly spaced evaluation points unless the integration accounts for their spacing.
Other choices can emphasize horizons that matter more for a clinical decision or represent a target distribution of use. A weight based on the frequency of decisions at different horizons answers a different question from uniform time weighting. Weights can also be specified to emphasize early or late prediction, but should be chosen before examining model performance where possible. They must be nonnegative for the summary to remain a weighted average of AUC values.
| Weighting choice | What receives emphasis | Interpretation |
|---|---|---|
| Uniform over time | Each unit of time in the selected interval | Average discrimination over elapsed horizon time |
| Application-based | Horizons important for the intended use | Average discrimination under the specified use pattern |
| Event- or population-based | Horizons receiving more weight under a chosen distribution | Average discrimination weighted by that distribution |
A weighted summary remains between the smallest and largest AUC values over the interval, provided the weights are nonnegative and the denominator is positive. It does not, however, reveal whether performance is steady or whether a high summary is driven by a narrow part of the range. Always report the horizon interval and weighting scheme alongside the integrated value.
Calculating an Integrated AUC From a Horizon Grid
In practice, AUC values are usually estimated at a finite set of horizons. The integral can then be approximated numerically. With uniform weights and a grid \(h_1<\cdots<h_m\), the trapezoidal rule uses the width of each interval: $$ \widehat{\operatorname{iAUC}}_{\mathrm{uniform}} = \frac{ \displaystyle\sum_{j=1}^{m-1} (h_{j+1}-h_j) \frac{\widehat{\operatorname{AUC}}(h_j)+\widehat{\operatorname{AUC}}(h_{j+1})}{2} } {h_m-h_1}. $$ If \(w(h)\) varies, approximate the numerator by integrating \(w(h)\widehat{\operatorname{AUC}}(h)\) and the denominator by integrating \(w(h)\), using the same numerical scheme. With a grid of prespecified interval weights \(q_j\), a discrete weighted average is \(\sum_j q_j\widehat{\operatorname{AUC}}(h_j)/\sum_j q_j\). These are related summaries, but the discrete weights and the integration rule should be stated explicitly.
A short justification clarifies what the calculation does. When \(w(h)\geq0\), define the normalized weight density \(g(h)=w(h)/\int_{h_L}^{h_U}w(u)\,du\). Then \(g(h)\geq0\), its integral is one, and $$ \operatorname{iAUC}_{w}=\int_{h_L}^{h_U}g(h)\operatorname{AUC}(h)\,dh. $$ Because this is an average with nonnegative weights, if \(a\leq\operatorname{AUC}(h)\leq b\) throughout the interval, then \(a\leq\operatorname{iAUC}_{w}\leq b\). Thus the integrated value summarizes the curve; it cannot improve the underlying discrimination or erase poor performance at a particular horizon.
Worked Example: Uniform Integration on an Irregular Horizon Grid
Suppose a hypothetical model predicts a specified transition from a fixed origin state. The estimated AUCs at horizons 1, 2, and 4 years are 0.68, 0.74, and 0.80. Assume the same case-control convention and target population were used at all three horizons, and use uniform time weighting from year 1 through year 4. For this calculation, assume linear interpolation between adjacent AUC estimates.
The interval from year 1 to year 2 has width 1 year, and its trapezoidal area is \(1(0.68+0.74)/2=0.71\) AUC-years. The interval from year 2 to year 4 has width 2 years, and its area is \(2(0.74+0.80)/2=1.54\) AUC-years.
The total area is \(0.71+1.54=2.25\) AUC-years. The full interval is \(4-1=3\) years, so $$ \widehat{\operatorname{iAUC}}_{\mathrm{uniform}} = \frac{2.25\ \text{AUC-years}}{3\ \text{years}} =0.750. $$ The result is an average discrimination of 0.750 over years 1 through 4 under uniform time weighting and the stated interpolation approximation. It is not the AUC at year 4.
The assumptions for this summary are that the horizon-specific estimates target a compatible prediction problem and that the trapezoidal approximation is adequate for the curve between observed horizons. With only three estimates, the interpolation assumption is consequential; additional horizons may be needed if the curve changes sharply.
How Weighting Changes the Answer
An integrated AUC is easiest to interpret when its weights correspond to the question being asked. A uniform time average and an average weighted toward near-term decisions are both legitimate summaries, but they are not interchangeable. Differences between them reveal where the model performs better or worse relative to the intended use.
Worked Example: Emphasizing Near-Term Prediction
Consider hypothetical AUC estimates of 0.82 at 1 year, 0.76 at 2 years, and 0.64 at 3 years. Suppose the intended application values the 1-, 2-, and 3-year horizons with prespecified discrete weights 0.50, 0.30, and 0.20. These weights sum to one, so the weighted average is $$ \widehat{\operatorname{iAUC}}_{q} = 0.50(0.82)+0.30(0.76)+0.20(0.64) =0.410+0.228+0.128 =0.766. $$
For comparison, equal discrete weights at those same three points would give \((0.82+0.76+0.64)/3=2.22/3=0.740\). The near-term-weighted summary is higher because the largest AUC receives the largest weight. The value 0.766 therefore describes discrimination under the specified horizon-use weights, not equal weighting over elapsed time. The example assumes the three estimates are comparable and that the weights were chosen to represent use, not selected after seeing which weighting gives the most favorable result.
For continuous horizon weighting, the same principle applies, but weights are assigned across time rather than only at selected points. In particular, a weight proportional to the frequency of clinical decisions at each horizon yields an average for that decision distribution, if that distribution is defined and normalized over the analysis interval. The weight function should not be confused with the estimator's censoring weights: the former defines which horizons matter in the summary, whereas censoring weights address incomplete outcome observation when estimating horizon-specific discrimination.
Worked Example: Restricting the Horizon Range
Suppose a hypothetical transition model has estimated AUCs of 0.84, 0.78, 0.70, and 0.62 at years 1, 2, 3, and 4. Use uniform time weighting and linear interpolation. Over years 1 to 4, the trapezoidal area is $$ \frac{0.84+0.78}{2}(1) +\frac{0.78+0.70}{2}(1) +\frac{0.70+0.62}{2}(1) =0.81+0.74+0.66=2.21. $$ Dividing by 3 years gives \(2.21/3=0.737\).
If the decision is specifically about years 1 to 2, the same uniform rule gives \((0.84+0.78)/2=0.810\). The two summaries differ because they answer different questions: discrimination averaged over the first year after year 1 is 0.810, whereas average discrimination across years 1 to 4 is 0.737. The later decline is not a computational defect; it is part of the longer-range performance being summarized. This comparison assumes that the prediction origin, transition target, case-control definition, and evaluation population remain aligned across all horizons.
Estimation, Diagnostics, and Limitations
The integration step does not by itself solve the problem of estimating AUC under right censoring. At each horizon, use an estimator appropriate to the outcome-observation process and the chosen competing-transition definition. If inverse probability of censoring weighting (IPCW) is used, its assumptions and diagnostics apply to each horizon-specific estimate. In particular, censoring must be adequately addressed under the estimator's assumptions, and estimated censoring probabilities must not be so small that weights become unstable. Integrating unstable estimates can produce a precise-looking number that is still unreliable.
Before reporting an integrated value, check the following:
- Case and control support: At each included horizon, verify that the data provide cases and controls under the chosen definition. If either group is absent, the horizon-specific AUC is undefined; a numerical integration cannot repair that.
- Comparable targets: Confirm that the origin, target transition, target population, and competing-transition convention are consistent across horizons.
- Censoring and weights: Inspect the horizon-specific estimates and any censoring weights for instability. Sparse risk sets, very large weights, or abrupt changes in the estimates warrant investigation.
- Grid and interpolation: Include enough horizons to represent plausible changes in discrimination. A coarse grid may miss a sharp rise or fall between evaluation points.
- Weight sensitivity: Compare the prespecified summary with a small number of substantively reasonable alternatives. Large changes indicate that the single integrated value depends strongly on the weighting choice.
A common mistake is to take an unweighted arithmetic mean of AUC estimates at irregularly spaced horizons and call it a uniform time-integrated AUC. For example, the grid in the first worked example has intervals of 1 and 2 years, so giving all three points equal weight overrepresents the short interval relative to the long one. The correction is to integrate over time, accounting for interval widths, or to state explicitly that the target is a discrete average over equally weighted evaluation points rather than a uniform time average.
Another limitation is that AUC concerns ranking. A model can have a high AUC while systematically overpredicting transition probabilities; it can also have a lower AUC while providing useful calibrated probabilities. Integrated AUC does not measure calibration or overall prediction error. The integrated Brier score is a distinct summary based on prediction error across horizons, while a single-horizon AUC focuses on ranking at one time. The appropriate measure depends on whether the validation question concerns discrimination, calibration, or probability accuracy.
Check Your Understanding
Answer each question, then compare with the solutions below.
- A model has AUCs 0.70 and 0.80 at horizons 1 and 3 years. Under uniform time weighting and linear interpolation, what is the integrated AUC over years 1 to 3?
- Why can an integrated AUC based on a near-term-heavy weight be larger than a uniform time-integrated AUC?
- What should be done if the horizon-specific AUC is undefined at one horizon because there are no controls under the chosen definition?
- Does an integrated AUC of 0.80 show that predicted transition probabilities are well calibrated? Explain.
Solutions. (1) With linear interpolation, the trapezoidal average is \((0.70+0.80)/2=0.750\); the interval width cancels when dividing the area by 2 years. (2) The near-term-heavy summary gives greater influence to near-term AUCs. If those values are higher than the longer-horizon values, the weighted average can be larger; it describes a different target weighting. (3) Do not treat the undefined value as zero or silently interpolate across it. Investigate why support is absent, restrict or redefine the horizon range if substantively justified, and report the change. (4) No. AUC measures ranking of cases relative to controls, not agreement between predicted probabilities and observed transition frequencies; assess calibration separately.