Tutorials › Competing Risks and Multi-State Models › Influential Subjects In Transition Validation

Transition-Specific Prediction Validation · Tutorial 423 of 699

Influential Subjects In Transition Validation

Use subject-level deletion to quantify how much each person changes a transition-specific validation metric, then investigate large changes without overinterpreting them.

Advanced 10 min read

What You'll Learn

  • Define leave-one-subject-out changes for transition-specific validation metrics.
  • Calculate deletion effects for average Brier scores and mean calibration residuals.
  • Assess subject influence on a transition-specific area under the ROC curve.
  • Distinguish influential subjects from subjects with large residuals.
  • Diagnose instability caused by rare transitions, censoring, and metric definitions.
  • Report influence with its direction, scale, and validation context.

Why Examine Subject Influence?

A transition-specific validation result is usually a summary across subjects: a prediction score, a calibration summary, or a measure of discrimination for a specified transition and time horizon. A single summary can conceal how unevenly the subjects contribute. For example, a small number of observed transitions may determine much of a transition-specific metric. A leave-one-subject-out analysis makes this dependence visible by recalculating the metric after removing each subject in turn.

The purpose is sensitivity analysis, not a test of whether a subject is erroneous. A large deletion effect says that the reported metric depends materially on that subject under the chosen validation procedure. It does not establish that the subject is an outlier, that the prediction model is wrong, or that the subject should be excluded. Residual diagnostics examine prediction errors; deletion influence instead asks how much the overall metric changes when a subject is removed.

Definition: Let \(M\) be a transition-specific validation metric calculated from \(n\) subjects. Let \(M_{(-i)}\) be the same metric, calculated by the same procedure after deleting subject \(i\). Define the leave-one-subject-out change as \(\Delta_i=M_{(-i)}-M\). The signed change describes the direction of change on the metric's scale; its magnitude, \(|\Delta_i|\), describes how much the estimate changes.

The transition, prediction horizon, target population, outcome definition, and metric must be held fixed when comparing \(M_{(-i)}\) with \(M\). If a metric is better when larger, a positive \(\Delta_i\) indicates that deleting subject \(i\) increases the metric. If smaller is better, as for a Brier score, a negative \(\Delta_i\) indicates improvement after deletion. Always state which direction is desirable; the sign alone is not a judgement of quality.

Calculation and Interpretation

For a metric that is the mean of subject-level contributions \(L_i\), such as an individual Brier loss, the full-sample metric is \(M=\sum_{j=1}^{n}L_j/n\). Deleting subject \(i\) gives the mean of the remaining \(n-1\) contributions. Subtracting the full-sample mean yields a simple exact expression:

$$ M_{(-i)}=\frac{nM-L_i}{n-1}, \qquad \Delta_i=M_{(-i)}-M=\frac{M-L_i}{n-1}. $$

This identity shows that, for an average loss, a subject with a loss far above the sample average pulls the full-sample loss upward; removing that subject lowers the mean. A loss below the average has the opposite effect. The factor \(n-1\) also matters: the same difference between an individual's contribution and the mean tends to produce a larger deletion change in a smaller sample.

Not every validation metric is an average of independent subject contributions. An area under the ROC curve (AUC), for instance, compares pairs of subjects. A calibration slope is estimated from a fitted relationship. For such metrics, calculate \(M_{(-i)}\) by recomputing the metric on the reduced data, rather than applying the mean-loss identity. If the validation procedure estimates nuisance quantities, such as censoring weights, specify whether and how those quantities are re-estimated after deletion. The deletion calculation should reflect the procedure used to obtain the reported metric.

Worked Example: Influence on a Transition-Specific Brier Score

Suppose five hypothetical subjects are eligible for the same transition-specific prediction target at a fixed horizon. Their outcomes are fully observed, and their individual Brier losses are \(0.04,\ 0.09,\ 0.10,\ 0.11,\) and \(0.26\). This simple example has no censoring; in a censored validation sample, the loss contributions would need to be defined using the prespecified method for handling censoring.

The full-sample mean Brier score is

$$ M=\frac{0.04+0.09+0.10+0.11+0.26}{5} =\frac{0.60}{5}=0.120. $$

For the subject with loss \(0.26\), the leave-one-out score is

$$ M_{(-i)}=\frac{0.04+0.09+0.10+0.11}{4} =\frac{0.34}{4}=0.0850, \qquad \Delta_i=0.0850-0.120=-0.0350. $$

Removing this subject lowers the Brier score by \(0.0350\) score units, so the full-sample score is materially higher when this relatively large loss is included. For the subject with loss \(0.04\), the reduced score is \((0.09+0.10+0.11+0.26)/4=0.140\), giving \(\Delta_i=0.0200\). This subject's inclusion pulls the full-sample score down. These calculations identify sensitivity in this sample; they do not justify removing either subject.

Nonlinear Metrics and Transition-Specific Examples

For a non-additive metric, leave-one-out changes can depend on a subject's role in comparisons or denominators, not just on how unusual that subject's prediction is. The following example uses a binary outcome at a fixed horizon: whether the specified transition occurred by that horizon. For this illustration, outcomes are completely observed. If competing transitions define non-events for this target, that definition must be set in advance; censoring requires an appropriate outcome and metric procedure rather than simply treating an unknown outcome as a non-event.

Worked Example: Influence on Transition-Specific AUC

Consider three hypothetical subjects who experienced the target transition, with predicted probabilities \(0.90,\ 0.80,\) and \(0.40\), and three who did not, with predictions \(0.70,\ 0.30,\) and \(0.20\). The empirical AUC is the proportion of case-control pairs in which the case has a higher prediction, with ties conventionally counted as one half. There are \(3\times3=9\) pairs and no ties here.

The case with prediction \(0.90\) outranks all three non-cases, contributing 3 concordant pairs. The case with prediction \(0.80\) also contributes 3. The case with prediction \(0.40\) outranks \(0.30\) and \(0.20\), but not \(0.70\), contributing 2. Thus \(M=8/9=0.8889\).

Deleting the case with prediction \(0.40\) leaves two cases and three non-cases. Both remaining cases outrank every non-case, so \(M_{(-i)}=6/6=1.000\) and \(\Delta_i=1.000-0.8889=0.1111\). Deleting the case with prediction \(0.90\) leaves 5 concordant pairs among 6: the \(0.80\) case contributes 3 and the \(0.40\) case contributes 2. Therefore \(M_{(-i)}=5/6=0.8333\) and \(\Delta_i=-0.0556\). Deleting the non-case with prediction \(0.70\) leaves 6 concordant pairs among 6, so \(\Delta_i=1.000-0.8889=0.1111\). The result depends on the pair comparisons and the case-control composition, not only on an individual's prediction error.

The deletion effects are large relative to the AUC's possible range, but the sample is deliberately small. They indicate that this estimated discrimination depends strongly on particular subjects and should be interpreted with the limited number of cases and non-cases in view.

A calibration summary can also be sensitive to individual subjects. One useful descriptive quantity is the mean observed-minus-predicted transition frequency at a specified horizon. It is a mean residual, not a calibration slope or a complete assessment of calibration: positive values indicate average underprediction and negative values average overprediction on this scale.

Worked Example: Influence on Mean Calibration Residual

For five hypothetical subjects, suppose the predicted probabilities of the target transition are \(0.10,\ 0.20,\ 0.30,\ 0.40,\) and \(0.80\), and the observed binary outcomes are \(0,\ 0,\ 0,\ 1,\) and \(1\). Assume complete outcome ascertainment for this calculation. The observed-minus-predicted residuals are \(-0.10,\ -0.20,\ -0.30,\ 0.60,\) and \(0.20\).

The mean residual is

$$ M=\frac{-0.10-0.20-0.30+0.60+0.20}{5} =\frac{0.20}{5}=0.0400. $$

The positive value indicates that the observed event proportion, \(2/5=0.400\), exceeds the mean predicted probability, \((0.10+0.20+0.30+0.40+0.80)/5=0.360\), by \(0.0400\). Deleting the subject with residual \(0.60\) gives \(M_{(-i)}=(-0.10-0.20-0.30+0.20)/4=-0.100\), and hence \(\Delta_i=-0.100-0.0400=-0.140\). Removing that subject changes the summary from slight average underprediction to average overprediction in the remaining four. This is a substantial sensitivity for a small sample, not evidence that the subject should be discarded.

Assumptions, Diagnostics, and Comparisons

A deletion analysis is meaningful only when its target and calculation are well defined. In particular:

  • Use the subject as the deletion unit. A subject may contribute several transitions, time points, or records. Remove all information belonging to that subject together, rather than treating repeated records as independent subjects.
  • Keep the validation target fixed. Specify the transition, origin state or landmark, horizon, endpoint definition, and eligibility rules. Deleting a subject should not silently change these definitions.
  • Reapply the metric procedure consistently. If the metric uses censoring adjustment or other estimated quantities, document whether those quantities are recomputed after each deletion. A method designed for censoring should not be replaced by treating censored outcomes as observed non-events.
  • Keep the prediction model's role clear. In external validation, predictions are typically fixed and the deletion assesses their validation metric. If the validation procedure itself includes fitting or recalibration, determine whether that step is part of the estimand and repeat it consistently for each reduced sample.

A useful diagnostic is a table or plot of all \(\Delta_i\), with the subject's transition status and relevant prediction or outcome information available for review. Check whether large changes arise from a rare event, a small risk set, unusual follow-up, high censoring weight, an extreme prediction, or a small metric denominator. Also check computational stability: the full-sample and reduced-sample metrics must use the same definitions, and a reduced sample that cannot support a metric should be reported as undefined rather than assigned a misleading value.

There is no universal cutoff for a “large” deletion effect. Compare changes with the metric's scale, the scientific importance of the transition, the uncertainty of the estimate, and any prespecified tolerance for practical impact. A change can be numerically small but important when the metric is near a decision threshold. Conversely, a noticeable change in a very small validation sample may mainly reflect limited information. Confidence intervals describe sampling uncertainty in a metric; deletion changes describe its sensitivity to specific observed subjects. The two answer different questions.

Compared with examining residuals alone, deletion analysis directly measures the effect on the chosen aggregate metric. Compared with a bootstrap, it is a targeted diagnostic that removes one observed subject at a time; it does not approximate the full sampling distribution or account for all possible cohort variation. A subgroup analysis asks whether performance differs across defined groups, whereas leave-one-out influence asks whether the metric is sensitive to individual subjects. These analyses can complement one another, but none by itself establishes why performance differs.

1
Define the validation quantity.
State the transition, horizon, subject set, endpoint handling, and metric, including which direction represents better performance.
2
Calculate the full-sample result.
Apply the prespecified validation procedure to all eligible subjects and record \(M\).
3
Delete one subject at a time.
Remove all records for that subject, repeat the same procedure, and calculate \(\Delta_i=M_{(-i)}-M\).
4
Investigate and report sensitivity.
Review the largest changes in context, check for instability or calculation failure, and report the transition, metric, direction, and scale.

Common Pitfalls and Limits

  • Calling influence an error. A subject can have a large deletion effect without a data error or an implausible outcome. Verify records, but do not equate sensitivity with invalidity.
  • Removing a subject because the result improves. Deleting observations after seeing their effect can bias the reported validation result. Report the full-sample estimate as primary and use deletion results as diagnostics unless exclusion rules were specified independently.
  • Ignoring transition-specific denominators. A deletion may remove an event, a non-event, or a key comparison pair. With rare transitions, the metric can change sharply or become undefined.
  • Mixing deletion with model refitting. Refitting predictions changes the question from influence on a validation metric for fixed predictions to influence on a broader modeling procedure. State and follow the intended procedure.
  • Overgeneralizing one cohort. Leave-one-subject-out effects describe dependence on the observed validation sample. They do not show that the same subject would be influential in another population or prove a causal mechanism.

The central practical value is transparency: a transition-specific metric should not be presented as if every subject contributed equally when one or a few subjects may dominate its value. Reporting the full estimate alongside a concise summary of the largest deletion changes helps readers judge stability while preserving the prespecified validation analysis.

Check Your Understanding

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

  1. A mean transition-specific loss is \(0.18\) for \(n=6\) subjects. One subject's loss is \(0.43\). Calculate that subject's leave-one-out mean and change.
  2. An AUC rises when one subject is deleted. Does that establish that the subject is an outlier or should be excluded?
  3. Why must all records from a subject be removed together if each subject has repeated transition records?
  4. A subject is censored before the prediction horizon. Why is it generally inappropriate to classify that subject automatically as a non-event in a transition-specific Brier score?

Solutions. (1) The total loss is \(6(0.18)=1.08\). After deleting the subject, the remaining total is \(1.08-0.43=0.65\), so \(M_{(-i)}=0.65/5=0.130\) and \(\Delta_i=0.130-0.18=-0.0500\). (2) No. The deletion indicates sensitivity of the AUC to that subject under the chosen procedure. It does not establish an error, outlier status, or a valid reason for exclusion. (3) Records from one subject are not separate independent people; deleting only some records leaves that subject partly represented and can understate subject-level influence. (4) Censoring before the horizon leaves the target outcome unknown. Treating it as a known non-event misclassifies an unobserved outcome and can distort the score; use an appropriate prespecified censoring-adjusted procedure.

Clinical Biostats