Tutorials › Joint Modeling of Longitudinal and Time-to-Event Data › Harrell’s Concordance for Survival Predictions

Calibration and Prediction Validation · Tutorial 555 of 697

Harrell's Concordance for Survival Predictions

Compute Harrell’s concordance for joint-model risk predictions and understand which censored-outcome pairs contribute to the estimate.

Advanced 11 min read

What You'll Learn

  • Define comparable pairs using observed event times and censoring indicators
  • Calculate Harrell’s concordance, including half-credit for tied risk scores
  • Orient joint-model predictions so higher scores consistently mean higher risk
  • Interpret the index as pairwise discrimination, not calibration or causation
  • Recognize how censoring affects which pairs contribute and limits interpretation
  • Distinguish an overall pairwise concordance measure from horizon-specific discrimination

Why Pairwise Concordance Matters

A survival prediction can rank people by how soon they are expected to experience an event. Harrell’s concordance index, often written as Harrell’s C, evaluates that ranking by asking whether the person who experiences an event earlier was assigned the higher predicted risk. It is a measure of discrimination: how well predictions order individuals, rather than how close predicted probabilities are to observed probabilities.

For joint-model predictions, the risk score may summarize information from a participant’s longitudinal history and survival model. To evaluate rankings meaningfully, compare predictions made for the same prediction setting—for example, at the same landmark time and over the same prediction horizon. A score for a different horizon may encode a different prediction target.

Observed Follow-Up and Comparable Pairs

For individual \(i\), let \(T_i\) be the event time and \(C_i\) the censoring time, measured from the prediction origin. The observed follow-up time is \(Y_i=\min(T_i,C_i)\), and \(\delta_i=I(T_i\leq C_i)\) indicates whether the event was observed: \(\delta_i=1\) for an event and \(\delta_i=0\) for censoring. Let \(r_i\) be the predicted risk score, oriented so that larger values indicate greater risk and, in general, an earlier event.

Definition: A pair of individuals is comparable for Harrell’s C if one has an observed event before the other individual’s observed follow-up ends. Specifically, if \(Y_i<Y_j\) and \(\delta_i=1\), the pair is comparable and \(i\) is the earlier observed event. The pair is concordant if \(r_i>r_j\), discordant if \(r_i<r_j\), and receives half-credit if \(r_i=r_j\). A pair with tied observed follow-up times is excluded under the convention used here.

A censored individual does not automatically make every pair involving them unusable. If an event is observed before that individual’s censoring time, the ordering is known: that event occurred first. But if censoring happens before the other individual’s observed event, their eventual event ordering is unknown, so that pair does not contribute. This rule makes the calculation possible without pretending that a censored observation is an event.

If \(\mathcal{P}\) denotes the set of comparable pairs, Harrell’s estimate is the concordant-pair proportion, with ties in risk score receiving half-credit:

$$ \widehat C_{\mathrm{Harrell}} = \frac{ \displaystyle\sum_{(i,j)\in\mathcal{P}} \left[I(r_i>r_j)+\tfrac{1}{2}I(r_i=r_j)\right] }{ |\mathcal{P}| }. $$

In each pair \((i,j)\), the indexing is chosen so that \(i\) is the person with the earlier observed event and \(j\) is the person whose follow-up extends beyond that event. The index ranges from 0 to 1 when at least one comparable pair exists. A value near 1 indicates mostly correctly ordered pairs, a value near 0 indicates mostly reversed pairs, and a value near 0.5 indicates little pairwise ordering ability in that evaluation set. A value of 0.5 is not proof that the prediction contains no useful information.

Worked Example: Counting Concordant and Discordant Pairs

Consider five hypothetical participants with predictions made at the same landmark and for the same horizon. Their risk scores are larger for higher predicted risk.

ParticipantObserved follow-up \(Y_i\) (months)Event observed \(\delta_i\)Risk score \(r_i\)
A210.80
B410.60
C500.40
D710.30
E800.70

Participant A has an event at month 2, before each of B, C, D, and E ends follow-up. These are four comparable pairs. The score for A is higher than each of the four other scores, so all four are concordant. Participant B’s event at month 4 precedes the ends of follow-up for C, D, and E. Of these pairs, B’s score exceeds C’s and D’s, but not E’s: two concordant and one discordant. Participant C is censored at month 5, so C’s pairs with D and E cannot establish who would have had the earlier event. Participant D’s event at month 7 precedes the end of E’s follow-up, giving one more comparable pair; it is discordant because \(0.30<0.70\).

There are \(4+3+1=8\) comparable pairs and \(4+2+0=6\) concordant pairs, with no score ties. Therefore,

$$ \widehat C_{\mathrm{Harrell}} = \frac{6}{8} = 0.750. $$

In this hypothetical evaluation set, 75.0% of comparable pairs are correctly ranked. This describes discrimination among the pairs that can be ordered from the observed data; it does not mean that 75.0% of individuals receive accurate risk estimates.

Ties, Score Direction, and a Short Justification

Risk scores do not need to be probabilities. They can be linear predictors or another continuous model output, provided their direction is clear. For example, a model may output predicted survival rather than predicted event risk. At a fixed horizon, larger survival probabilities imply lower risk, so those values must be reversed in direction before applying the higher-risk-is-earlier-event rule. One option is to use \(r_i=1-\widehat S_i\), where \(\widehat S_i\) is predicted survival at that horizon.

Half-credit for a tied score is a convention with a useful symmetry: if the tied pair is ordered arbitrarily, reversing that arbitrary ordering would turn a concordant pair into a discordant one. Giving half-credit avoids favoring either direction. It does not make tied predictions informative; it records that the model did not distinguish the pair.

Justification: For every comparable pair, the observed data establish that \(i\) had an event before \(j\)’s follow-up ended. Under the intended ranking, the model should assign greater risk to \(i\) than to \(j\). The indicator \(I(r_i>r_j)\) contributes 1 when that ordering is correct and 0 otherwise; \(I(r_i=r_j)/2\) contributes 0.5 for a tie. Summing these contributions and dividing by the number of comparable pairs therefore gives the average pairwise credit. If every comparable pair is concordant, the ratio is 1; if every pair is discordant, it is 0.

Worked Example: Risk-Score Ties

Suppose four hypothetical participants have the following observed follow-up and scores: A has an event at month 3 with \(r_A=0.90\); B is censored at month 5 with \(r_B=0.50\); C has an event at month 6 with \(r_C=0.90\); and D is censored at month 9 with \(r_D=0.20\).

A’s event precedes the ends of follow-up for B, C, and D, so these three pairs are comparable. A’s score is higher than B’s and D’s, giving two concordant pairs. A and C have tied scores, giving half-credit for their pair. C’s event at month 6 precedes the end of D’s follow-up, creating one additional comparable pair; \(0.90>0.20\), so it is concordant. B’s censoring at month 5 occurs before C’s event, so the B–C pair is not comparable.

The total credit is \(1+1+0.5+1=3.5\), and there are four comparable pairs. Thus,

$$ \widehat C_{\mathrm{Harrell}} = \frac{3.5}{4} = 0.875. $$

The estimate is 0.875 for these data. The tied pair contributes half, not a full concordance, because the predictions do not order its members.

What Censoring Changes

Harrell’s C uses only pairs whose ordering is known from the observed follow-up. As censoring increases, fewer pairs may be comparable. More importantly, the set of retained pairs can become selective: it depends on when follow-up ends, and therefore on the censoring process as well as the event times and scores. The resulting estimate is directly a summary of the observed comparable pairs. It need not equal the concordance that would be calculated if every participant’s event time were known.

Independent censoring is a useful condition in survival analysis, but it does not make the set of comparable pairs identical to all pairs with known latent event ordering. Even when censoring is independent of event time under an appropriate conditioning set, the probability that a particular pair becomes comparable can vary with its event times. Heavy censoring can therefore make Harrell’s C unstable or make it reflect a different mix of pairs than a less-censored evaluation would. If censoring depends on participant characteristics or predicted risk, the observed comparisons can be especially selective.

Always report the number of comparable pairs alongside the estimate when possible. A high C based on very few pairs is not as informative as the same value supported by many pairs. Check follow-up and censoring across ranges of predicted risk, and investigate whether early censoring is concentrated in particular groups. Such checks do not repair selection, but they help reveal when the estimate may be difficult to interpret.

Worked Example: Why a Censored Pair Is Not Automatically Comparable

In a hypothetical validation set, participant F is censored at month 4 with risk score \(0.85\). Participant G has an observed event at month 6 with risk score \(0.30\). Since F’s follow-up ends before G’s event, the data do not establish whether F would have had an event before or after month 6. The F–G pair is not comparable, even though F’s score is higher.

Now suppose participant H has an observed event at month 3 with risk score \(0.40\), while F is still followed until month 4. H’s event occurs before F’s censoring time, so H–F is comparable. The score ordering is discordant because \(0.40<0.85\): the model assigned higher risk to F even though H had the earlier observed event. The distinction is determined by follow-up ordering, not by which score appears more plausible.

Using Concordance for Joint-Model Predictions

For a joint model, a prediction may be conditional on surviving to a landmark and on the longitudinal history observed up to that point. Harrell’s C should be computed using a clearly defined cohort and a common prediction target. For instance, if the target is event risk over a particular horizon after a landmark, generate the risk predictions at that landmark using information available by that time, and define follow-up from that landmark. Do not use later biomarker measurements to construct predictions that are meant to represent what could have been predicted at the landmark.

For a fixed horizon, higher predicted event probabilities should generally correspond to higher risk scores. Harrell’s C then evaluates how well those scores order observed event times among comparable pairs. It does not check whether a predicted risk of 0.30 corresponds to an observed event frequency of 30%; that is a calibration question. Nor is it the same as the time-dependent AUC, which evaluates discrimination using a specified time and case/control definitions. The measures can answer related but distinct questions.

Compared with a horizon-specific discrimination measure, Harrell’s C summarizes ordering across the comparable event-time pairs available in the data. That can be useful as a broad summary, but it can also obscure changes in prediction performance over time. Under heavy censoring, a censoring-adjusted concordance approach may be preferred; the choice of method should match the evaluation target and censoring conditions.

Common Mistakes and Limitations

  • Counting every pair. A pair where the earlier observed follow-up is censored does not reveal the event ordering. Include it only if the other person’s observed event occurred before that censoring time.
  • Reversing the score direction. If larger values mean better survival rather than greater event risk, the resulting C will be interpreted backwards unless the score is reoriented.
  • Calling C an accuracy percentage. A value of 0.75 means 75% pairwise credit among comparable pairs, not that 75% of individual predictions are correct or calibrated.
  • Ignoring censoring and sample size. The estimate depends on the observed comparable-pair set. Report its size and assess whether censoring patterns could make that set unrepresentative.
  • Comparing unlike targets. Predictions from different landmarks or horizons need not be comparable as if they were one score. Define the target before ranking.

Harrell’s C is most useful when its pairwise meaning is kept explicit. It reduces a complex prediction output to an interpretable ranking summary, but it does not by itself establish calibration, clinical usefulness, or causal effects. Its value should be considered alongside the prediction target, censoring pattern, and other performance summaries.

Check Your Understanding

Answer each question, then check the solutions below.

  1. Participant A has an event at month 2 and risk score 0.7. Participant B is censored at month 5 and has score 0.4. Is the pair comparable, and what credit does it receive?
  2. For a comparable pair, the earlier event has score 0.3 and the other participant has score 0.8. Is the pair concordant or discordant?
  3. Three comparable pairs have two concordant pairs and one tied-score pair. What is Harrell’s C?
  4. In one sentence, explain why a pair with censoring at month 4 and an observed event at month 6 is not comparable when the censored participant is the one followed only to month 4.
  5. What does a Harrell’s C of 0.80 mean, and what does it not mean?

Solutions. (1) Yes. A’s event is observed before B’s follow-up ends, so the pair is comparable; \(0.7>0.4\), giving one concordant credit. (2) It is discordant because the earlier-event participant has the lower risk score, \(0.3<0.8\). (3) The total credit is \(2+0.5=2.5\), so \(\widehat C=2.5/3=0.833\). (4) Censoring occurs before the other participant’s event, so the event ordering that would involve the censored person is unknown. (5) It means the predictions receive 80% pairwise credit among comparable pairs in that evaluation set; it does not mean 80% of individual risks are accurate, nor that predicted probabilities are calibrated.

Clinical Biostats