Tutorials › Multi-Omics Data Integration › Selecting Metrics for Integration Comparisons

Early, Intermediate and Late Integration · Tutorial 386 of 597

Selecting Metrics for Integration Comparisons

Choose metrics that match the outcome and prediction goal, and recognize when different valid metrics favor different integration strategies.

Advanced 12 min read

What You'll Learn

  • Match continuous-outcome metrics to the scale and consequences of prediction errors.
  • Distinguish threshold-based binary metrics from metrics that assess predicted probabilities or rankings.
  • Calculate MAE, RMSE, test-set \(R^2\), sensitivity, specificity, precision, and balanced accuracy.
  • Explain why class imbalance can make accuracy misleading.
  • Compare probability predictions using Brier score and log loss.
  • Specify the evaluation sample and decision goal needed to interpret a metric.

Why Metric Choice Matters

An integration benchmark can produce several reasonable summaries of prediction performance, and those summaries need not rank strategies in the same order. A method may make small errors for most individuals but a few very large errors; another may have lower average absolute error but worse extreme misses. For a binary outcome, a strategy may classify most observations correctly by favoring the common class while failing to identify the uncommon but important class. A metric is not a neutral decoration: it defines which aspects of prediction count as good.

The comparison is meaningful only when strategies are evaluated on the same held-out individuals, with the same outcome definition and compatible prediction targets. The held-out data should represent the population to which performance is meant to generalize, and all preprocessing and model fitting must be confined to the training data. A metric cannot repair an evaluation design that does not support the intended comparison.

Definition: A prediction metric is a numerical summary computed from observed outcomes and predictions, where its direction and scale express a specified aspect of predictive performance. For a loss, smaller values are better; for a score such as sensitivity or area under a curve, larger values are better. The target outcome, prediction type, evaluation population, and practical objective determine which metric is relevant.

Continuous Outcomes: Error Size and Error Shape

Let \(y_i\) be the observed continuous outcome and \(\hat{y}_i\) its prediction for evaluation individual \(i\), for \(i=1,\ldots,n\). The prediction error is \(e_i=\hat{y}_i-y_i\). Two common summaries are the mean absolute error (MAE) and root mean squared error (RMSE):

$$ \operatorname{MAE}=\frac{1}{n}\sum_{i=1}^{n}|e_i|, \qquad \operatorname{RMSE}=\sqrt{\frac{1}{n}\sum_{i=1}^{n}e_i^2}. $$

Both are expressed in the outcome's units. MAE gives the average absolute prediction error and is less dominated by a small number of large errors. RMSE gives greater weight to large errors because errors are squared before averaging. Neither is universally preferable: if unusually large errors have especially serious consequences, RMSE may reflect the objective better; if the desired summary is a typical absolute miss, MAE may be easier to interpret.

A related summary is the test-set coefficient of determination:

$$ R^2_{\mathrm{test}}= 1-\frac{\sum_{i=1}^{n}(y_i-\hat{y}_i)^2} {\sum_{i=1}^{n}(y_i-\bar{y})^2}, \qquad \bar{y}=\frac{1}{n}\sum_{i=1}^{n}y_i. $$

This compares squared prediction error with the squared error from predicting the evaluation-set mean \(\bar{y}\) for everyone. It is dimensionless; higher values indicate less squared error relative to that reference. On held-out data, \(R^2_{\mathrm{test}}\) can be negative when predictions are worse than the reference. It is not a percentage of individual outcomes correctly explained, and its value depends on the outcome variation in the evaluation sample.

Worked Example: Comparing Continuous-Outcome Summaries

Suppose a hypothetical test set contains four measured biomarker concentrations, in micromoles per litre (\(\mu\)mol/L). The observed values are \(10,12,15,17\), and one integration strategy predicts \(11,10,16,20\). The errors \(\hat{y}_i-y_i\) are \(1,-2,1,3\) \(\mu\)mol/L.

The absolute errors sum to \(1+2+1+3=7\), so

$$ \operatorname{MAE}=\frac{7}{4}=1.75\ \mu\text{mol/L}. $$

The squared errors sum to \(1+4+1+9=15\), giving

$$ \operatorname{RMSE}=\sqrt{\frac{15}{4}} =\sqrt{3.75}\approx1.936\ \mu\text{mol/L}. $$

The observed mean is \(\bar{y}=(10+12+15+17)/4=13.5\ \mu\)mol/L. The squared deviations from this mean sum to \(12.25+2.25+2.25+12.25=29\). Thus

$$ R^2_{\mathrm{test}}=1-\frac{15}{29}\approx0.4828. $$

On these four hypothetical observations, the strategy's squared error is about 48.3% lower than that of predicting the test-set mean for every individual. That statement is specific to this evaluation set and squared-error reference; it does not imply that 48.3% of variation in every future sample is explained.

The relationship between MAE and RMSE helps explain why they can disagree in comparisons. Squaring makes a large error count disproportionately: an error of 6 contributes 36 to the squared-error sum, whereas an error of 2 contributes 4. Consequently, a strategy with one extreme miss can have a competitive MAE but a poor RMSE. If all strategies are evaluated on the same observations, minimizing RMSE and minimizing mean squared error give identical rankings because the square root is increasing. MAE need not give that ranking.

Binary Outcomes: Thresholds, Probabilities, and Rankings

For a binary outcome, write \(y_i=1\) for the event and \(y_i=0\) otherwise. A model may provide a predicted probability \(\hat{p}_i\), or a predicted class obtained by applying a threshold \(c\): predict 1 when \(\hat{p}_i\ge c\), and 0 otherwise. The threshold is part of a threshold-based metric's definition. The resulting counts are true positives (TP), false negatives (FN), false positives (FP), and true negatives (TN).

$$ \operatorname{Accuracy}=\frac{\mathrm{TP}+\mathrm{TN}}{n}, \quad \operatorname{Sensitivity}=\frac{\mathrm{TP}}{\mathrm{TP}+\mathrm{FN}}, \quad \operatorname{Specificity}=\frac{\mathrm{TN}}{\mathrm{TN}+\mathrm{FP}}, \quad \operatorname{Precision}=\frac{\mathrm{TP}}{\mathrm{TP}+\mathrm{FP}}. $$

Sensitivity is the proportion of actual events detected; specificity is the proportion of non-events correctly ruled out; precision is the proportion of predicted events that are actual events. Balanced accuracy is the mean of sensitivity and specificity. These proportions are undefined if their denominators are zero, so report the relevant counts and handle such cases explicitly rather than silently assigning a value.

Accuracy uses one threshold and gives equal weight to each correctly classified observation. When the event is uncommon, predicting the non-event for everyone can yield high accuracy while sensitivity is zero. In contrast, sensitivity and specificity show performance for each actual class, and precision speaks directly to how often a positive prediction is correct. Precision depends on event prevalence in the evaluation population, so it may change when prevalence changes, even if sensitivity and specificity do not.

Worked Example: Accuracy Can Hide Missed Events

In a hypothetical evaluation of 100 patients, 10 have the event and 90 do not. At a prespecified probability threshold, an integration strategy predicts 20 patients as positive: 8 are true positives, 12 are false positives, 2 events are missed, and 78 non-events are correctly classified.

The accuracy is \((8+78)/100=0.860\), or 86.0%. Sensitivity is \(8/(8+2)=0.800\), or 80.0%; specificity is \(78/(78+12)=0.8667\), or 86.7%; and precision is \(8/(8+12)=0.400\), or 40.0%. Balanced accuracy is \((0.800+0.8667)/2=0.8333\), or 83.3%.

A rule predicting every patient to be negative would have accuracy \(90/100=0.900\), or 90.0%, but sensitivity \(0/10=0\). Thus accuracy alone would rank that rule above the strategy even though it detects none of the events. If detecting events is central to the intended use, sensitivity should be reported and the threshold selected in light of that goal; the figures do not establish that any threshold is clinically appropriate.

When probability predictions matter, rather than only classifications at one threshold, two standard losses are the Brier score and log loss. For binary outcomes, they are

$$ \operatorname{Brier}=\frac{1}{n}\sum_{i=1}^{n}(\hat{p}_i-y_i)^2, \qquad \operatorname{LogLoss}=-\frac{1}{n}\sum_{i=1}^{n} \left[y_i\log(\hat{p}_i)+(1-y_i)\log(1-\hat{p}_i)\right]. $$

Both require probabilities between 0 and 1, and smaller values are better. Brier score is mean squared probability error. Log loss penalizes confident predictions that assign very low probability to the outcome that occurs; a probability of exactly 0 for an event or exactly 1 for a non-event gives infinite loss. In computation, probabilities may be bounded away from 0 and 1 using a prespecified rule, which should be reported.

Area under the receiver operating characteristic curve (ROC AUC) summarizes ranking: it is the probability that a randomly selected event receives a higher score than a randomly selected non-event, with ties counted as one-half. It does not assess performance at a particular threshold. Precision-recall summaries focus on precision and sensitivity across thresholds and can be informative when events are rare; the precise area calculation can depend on the implementation, so the convention should be stated. Neither kind of curve summary replaces a decision-relevant threshold analysis.

Worked Example: Probability Losses and Ranking Can Tell Different Stories

Consider four hypothetical patients, two with events and two without. Their outcomes are \(1,1,0,0\). Strategy A predicts probabilities \(0.9,0.6,0.4,0.1\), while strategy B predicts \(0.7,0.7,0.3,0.3\), in the same order.

For A, squared probability errors are \(0.01,0.16,0.16,0.01\), so its Brier score is \((0.01+0.16+0.16+0.01)/4=0.085\). For B, all four squared errors are \(0.09\), giving Brier score \(0.090\). A's log loss is \[ -\frac{\log(0.9)+\log(0.6)+\log(0.6)+\log(0.9)}{4} =-\frac{\log(0.2916)}{4}\approx0.3083. \] For B it is \[ -\frac{2\log(0.7)+2\log(0.7)}{4} =-\log(0.7)\approx0.3567. \] Both probability losses favor A on this small set.

Both strategies also rank both events above both non-events, so each has ROC AUC \(1\). AUC therefore ties the strategies, while the probability losses distinguish them. This is not a contradiction: AUC summarizes ranking, whereas Brier score and log loss evaluate the probabilities themselves. Four observations are only an illustration, not evidence that the difference will persist in a target population.

A Practical Selection Framework

Before computing metrics, state what is predicted and what decision the prediction is meant to support. Then choose a small set of complementary summaries rather than selecting whichever metric makes a preferred strategy look best. The following distinctions help make that choice explicit.

Prediction target or priorityUseful metric(s)Interpretive caution
Continuous outcome; typical absolute missMAEDoes not emphasize a few very large errors as strongly as RMSE.
Continuous outcome; large misses especially costlyRMSE, with MAE for contextSquared errors give outliers substantial influence.
Binary class at a specified operating thresholdSensitivity, specificity, precision, and the confusion countsState the threshold; precision depends on prevalence.
Binary probabilities used as probabilitiesBrier score or log lossReport how exact 0 or 1 predictions are handled for log loss.
Binary ordering across thresholdsROC AUC; precision-recall summary when event rarity is importantA ranking summary alone does not select an operating threshold.

If the benchmark includes several simulated scenarios, define the primary metric before comparing strategies and report scenario-specific performance rather than relying only on a pooled value. A metric's meaning also depends on the outcome scale, class prevalence, threshold, and evaluation population. A single overall score can conceal which scenarios or types of errors drive the result.

Diagnostics, Alternatives, and Limitations

Check the ingredients of a metric before interpreting it. For continuous outcomes, inspect the distribution of errors and identify whether a few extreme errors dominate RMSE. Confirm that outcome units are consistent across strategies and scenarios; MAE and RMSE are not directly comparable across differently scaled outcomes without a deliberate standardization plan. Test-set \(R^2\) can be unstable when the observed outcomes have little variation and can be negative, so report an error metric as well.

For binary outcomes, report event and non-event counts, the probability threshold when using class metrics, and the confusion counts behind the proportions. With rare events, accuracy can be deceptively high, while precision can change substantially as prevalence changes. ROC AUC can look strong even when precision at useful thresholds is low. Brier score and log loss are defined for probabilities, not merely class labels; converting probabilities to classes discards information these losses use.

No metric identifies the best integration strategy independently of the use case. If false negatives and false positives have different consequences, sensitivity and specificity at a prespecified threshold may be more actionable than accuracy, but the costs or acceptable trade-offs must come from the application. If a strategy is better on one metric and worse on another, report that trade-off instead of treating the metrics as interchangeable or choosing a winner after seeing the results. These summaries also do not by themselves establish uncertainty in a ranking or explain which omics block contributes the performance; those require separate analyses.

Common pitfall: “Strategy A has higher accuracy, so it is better for the binary outcome.” This conclusion can fail under class imbalance, as a rule that predicts only the majority class may have high accuracy and zero sensitivity. The correction is to report class-specific performance and confusion counts, and to select a threshold and primary metric that match the intended prediction task.

Check Your Understanding

Answer each question, then compare with the brief solutions.

  1. For continuous predictions, what feature of RMSE makes it more sensitive than MAE to a few large errors?
  2. A binary classifier has 18 true positives, 2 false negatives, 30 false positives, and 50 true negatives. Calculate sensitivity, specificity, and precision.
  3. Why can two probability models have the same ROC AUC but different Brier scores?
  4. When would accuracy alone be especially misleading for a binary outcome?
  5. What information must accompany a reported sensitivity so its interpretation is clear?

Solutions. (1) RMSE squares errors before averaging, so a large error contributes disproportionately. (2) Sensitivity is \(18/(18+2)=0.900\); specificity is \(50/(50+30)=0.625\); precision is \(18/(18+30)=0.375\). (3) AUC evaluates ranking, while Brier score evaluates the numerical probability errors; models can rank cases identically but assign different probabilities. (4) It is especially misleading when one class is much more common, since predicting the majority class can achieve high accuracy while failing to detect the minority class. (5) Give the threshold used to define a positive prediction and the number of actual events, along with the evaluation context; sensitivity is the fraction of actual events classified positive at that threshold.

Clinical Biostats