Tutorials › Multi-Arm Multi-Stage (MAMS) Designs › Continuous Biomarkers for Arm Selection

Time-to-Event MAMS With Intermediate Outcomes · Tutorial 308 of 499

Continuous Biomarkers for Arm Selection

Learn to choose and interpret standardized biomarker contrasts for interim arm selection, accounting for within-arm variability and dependence from a shared control.

Advanced 10 min read

What You'll Learn

  • Define a biomarker contrast with its favorable direction stated explicitly.
  • Distinguish a standardized mean difference from a precision-weighted Z statistic.
  • Calculate how arm-specific variances change interim rankings.
  • Derive the covariance between treatment-control contrasts that share a control.
  • Assess uncertainty in the ordering of two active arms.
  • Recognize when missingness, outliers, or inconsistent scales undermine a ranking.

Why Standardize a Continuous Biomarker Contrast?

At an interim analysis, a continuous biomarker can help identify which treatment arms warrant further study while the final time-to-event outcome is not yet available. A raw difference in biomarker means is easy to interpret in the biomarker’s units, but raw differences may be hard to compare if arms have different variability or if the biomarker’s scale is not clinically intuitive. Standardization can put contrasts on a common scale. It does not, by itself, make a biomarker a reliable surrogate for the final endpoint or establish that an arm will succeed on that endpoint.

The direction of benefit must be fixed before comparing arms. For example, if lower biomarker values are favorable, define a sign-reversed biomarker so that higher values always mean a more favorable response. Let arm \(0\) denote the shared control and arm \(a\) an active treatment. At a given interim look, write \(n_j\), \(\bar{x}_j\), and \(s_j\) for the number of observed participants, the sample mean, and the sample standard deviation in arm \(j\), respectively.

Definition: The estimated biomarker contrast for active arm \(a\) is \(\hat{\Delta}_a=\bar{x}_a-\bar{x}_0\), in biomarker units. A positive value favors arm \(a\) under the convention that larger values are favorable. A standardized mean difference divides this contrast by a prespecified reference standard deviation. A precision-standardized contrast divides it by its estimated standard error.

Those last two quantities answer different questions. A standardized mean difference describes the estimated effect relative to a variability scale. A Z statistic describes the estimated contrast relative to its sampling uncertainty. Neither is automatically the correct selection criterion: the design should specify the statistic and favorable direction in advance.

Assumptions and the Main Statistics

For the formulas below, assume participants are randomized to one of the arms, each participant contributes one interim biomarker measurement, and different participants’ measurements are independent. Within each arm, the biomarker has a finite variance; approximate normality within arms is a convenient condition for small samples, while sufficiently large samples may support normal approximations for means. The sample variances should estimate the variability of the measurements available at this interim look. If biomarker values are missing, the analysis also needs a defensible missing-data assumption and a prespecified method; simply analyzing observed values can be biased when availability depends on response or treatment.

With independent arms, an estimated standard error for the treatment-control contrast is

$$ SE(\hat{\Delta}_a)=\sqrt{\frac{s_a^2}{n_a}+\frac{s_0^2}{n_0}}, \qquad Z_a=\frac{\hat{\Delta}_a}{SE(\hat{\Delta}_a)}. $$

Here \(Z_a\) is a precision-standardized contrast. Its magnitude can be large because the estimated difference is large, its standard error is small, or both. It is not a standardized effect size in the sense of “number of within-arm standard deviations.” In small samples, inference for a single contrast may use a \(t\) statistic with an appropriate degrees-of-freedom approximation; an interim selection rule should use the method specified by the design.

For an effect-size scale, one possible pairwise measure is the pooled-SD standardized mean difference:

$$ \hat{d}_a= \frac{\bar{x}_a-\bar{x}_0}{s_{p,a}}, \qquad s_{p,a}= \sqrt{ \frac{(n_a-1)s_a^2+(n_0-1)s_0^2} {n_a+n_0-2} }. $$

This standardizes each contrast by a pooled estimate of the standard deviation in the two arms being compared. If variability differs across arms, then \(s_{p,a}\) can differ across active-control comparisons. Consequently, rankings by \(\hat{d}_a\) can differ from rankings by raw contrasts or by \(Z_a\). If a common standard deviation is used instead, it must be defined consistently and justified for the arms being compared.

Worked Example: A Larger Raw Contrast Does Not Guarantee a Higher Precision Statistic

Suppose a hypothetical interim analysis has \(n_0=n_A=n_B=40\), control mean \(\bar{x}_0=10\), and control standard deviation \(s_0=4\). Arm A has \(\bar{x}_A=14\) and \(s_A=2\); arm B has \(\bar{x}_B=16\) and \(s_B=8\). Higher biomarker values are favorable, and the assumptions above are taken to hold for these illustrative data.

The raw contrasts are \(\hat{\Delta}_A=14-10=4\) units and \(\hat{\Delta}_B=16-10=6\) units, so B ranks higher by observed raw contrast. For A,

$$ SE(\hat{\Delta}_A)=\sqrt{\frac{2^2}{40}+\frac{4^2}{40}} =\sqrt{0.5}=0.7071, \qquad Z_A=\frac{4}{0.7071}=5.657. $$

For B,

$$ SE(\hat{\Delta}_B)=\sqrt{\frac{8^2}{40}+\frac{4^2}{40}} =\sqrt{2}=1.414, \qquad Z_B=\frac{6}{1.414}=4.243. $$

Thus A has the larger precision-standardized statistic, despite its smaller raw contrast. B’s much greater within-arm variability makes its estimated contrast less precise. The pooled standard deviations are \(s_{p,A}=\sqrt{10}=3.162\) and \(s_{p,B}=\sqrt{40}=6.325\), giving \(\hat{d}_A=4/3.162=1.265\) and \(\hat{d}_B=6/6.325=0.9487\). The effect-size ranking also favors A here, but for a different reason: each contrast is divided by its pair-specific pooled variability. These statistics should not be treated as interchangeable selection rules.

Shared-Control Dependence

When two active arms share the same control participants, their estimated treatment-control contrasts are not independent. The same random control mean is subtracted from both. Under independent participants across arms, the covariance is

$$ \operatorname{Cov}(\hat{\Delta}_a,\hat{\Delta}_b) = \operatorname{Cov}(\bar{x}_a-\bar{x}_0,\bar{x}_b-\bar{x}_0) = \operatorname{Var}(\bar{x}_0) = \frac{\sigma_0^2}{n_0}, \qquad a\ne b. $$

The positive covariance applies when the contrasts use the same control observations and the same favorable-direction convention. If the control mean happens to be unusually low, both treatment-control contrasts tend to be higher; if it is unusually high, both tend to be lower. The correlation also depends on the variances and sample sizes in the active arms:

$$ \operatorname{Corr}(\hat{\Delta}_a,\hat{\Delta}_b) = \frac{\sigma_0^2/n_0} {\sqrt{(\sigma_a^2/n_a+\sigma_0^2/n_0) (\sigma_b^2/n_b+\sigma_0^2/n_0)}}. $$

This dependence matters when describing the uncertainty in a set of arm rankings or calculating the probability that one arm ranks above another. Treating shared-control contrasts as independent generally misrepresents their joint uncertainty. By contrast, in the difference between two active-control contrasts, the control mean cancels:

$$ \hat{\Delta}_a-\hat{\Delta}_b =(\bar{x}_a-\bar{x}_0)-(\bar{x}_b-\bar{x}_0) =\bar{x}_a-\bar{x}_b. $$

Under independent active-arm samples, its variance is \(\sigma_a^2/n_a+\sigma_b^2/n_b\). The shared control therefore affects the precision of each arm-versus-control contrast, but not the direct comparison between two active-arm sample means.

Worked Example: Quantifying Dependence From a Shared Control

Consider hypothetical arms with \(n_0=50\), \(s_0=2\), \(n_A=40\), \(s_A=3\), \(n_B=30\), and \(s_B=4\). The arm samples are independent, but A and B each use the same control participants. For this calculation, sample variances estimate the corresponding population variances.

The estimated covariance between the two contrasts is the control contribution:

$$ \widehat{\operatorname{Cov}}(\hat{\Delta}_A,\hat{\Delta}_B) =\frac{s_0^2}{n_0} =\frac{4}{50}=0.0800. $$

Their estimated standard errors are

$$ SE(\hat{\Delta}_A)=\sqrt{\frac{9}{40}+\frac{4}{50}} =\sqrt{0.305}=0.5523, \qquad SE(\hat{\Delta}_B)=\sqrt{\frac{16}{30}+\frac{4}{50}} =\sqrt{0.6133}=0.7832. $$

The estimated correlation is therefore \(0.0800/(0.5523\times0.7832)=0.1849\). The shared control produces positive dependence, though it is not overwhelming in this example. A calculation that assumes the two contrasts are independent would omit this covariance and give an incorrect joint description of their uncertainty.

Interpreting an Interim Ranking

A ranking is a summary of observed data, not proof that the top-ranked arm has the largest true biomarker effect. In particular, two arms can have similar estimated contrasts but substantial uncertainty about their ordering. A direct comparison of their estimated contrasts uses \(\bar{x}_a-\bar{x}_b\); for independent active-arm samples, its standard error is \(\sqrt{s_a^2/n_a+s_b^2/n_b}\). A confidence interval for this difference helps distinguish a clear ordering from a ranking that is mostly numerical.

Worked Example: Checking Whether Two Arms Are Clearly Ordered

Use the hypothetical data from the first example: A has mean 14 and standard deviation 2, while B has mean 16 and standard deviation 8, with 40 participants per arm. The estimated A-minus-B difference is \(14-16=-2\) biomarker units. The control arm is not part of this direct comparison.

The estimated standard error is

$$ SE(\bar{x}_A-\bar{x}_B)= \sqrt{\frac{2^2}{40}+\frac{8^2}{40}} =\sqrt{1.7}=1.304. $$

Using 1.96 as a large-sample normal multiplier, an approximate 95% interval for the A-minus-B difference is

$$ -2 \mathbin{\pm} 1.96(1.304) =-2 \mathbin{\pm} 2.556 =(-4.556,\ 0.556). $$

This interval includes zero. The data therefore do not provide a clear ordering of the two active-arm mean responses at this precision, even though B has the larger observed raw contrast and A the larger \(Z\) statistic and pooled-SD effect estimate in the earlier calculations. Selection may still require choosing an arm under a prespecified rule, but the ranking should not be reported as a certain difference in true effects.

Diagnostics, Alternatives, and Limitations

Before using a standardized biomarker contrast for selection, check the analysis scale and the data that determine its numerator and denominator. Useful checks include:

  • Direction and units: verify that larger values mean the same favorable direction in every arm, and that raw contrasts are reported in the biomarker’s units.
  • Variability: inspect arm-specific spreads and influential observations. A few extreme measurements can change both a sample mean and a sample standard deviation, especially at a small interim sample size.
  • Availability: compare the number and pattern of observed biomarker measurements across arms. Differential missingness can make an observed-case contrast misleading.
  • Control sharing: retain the covariance induced by the shared control whenever joint uncertainty in multiple contrasts is needed.
  • Rule stability: check whether the selected arm changes under reasonable, prespecified choices of standardization or robust summaries. A highly unstable ordering signals weak information, not a license to choose the most favorable analysis after seeing the data.

The main alternatives serve different purposes. Raw mean differences preserve the biomarker’s units and are often easiest to interpret clinically. A common-reference standardized mean difference aids comparison on a common scale if the reference standard deviation is meaningful and fixed consistently. Pairwise pooled-SD effects adapt to each comparison but may reorder arms when variability differs. Precision-standardized statistics account for both variability and sample size and can be useful for evidence-based selection, but are not direct measures of clinical effect magnitude. Robust or transformed analyses may be needed for strong skewness or influential values, but their estimands and selection rules must be stated explicitly.

All these approaches have limits. A biomarker contrast is an interim measure; it does not guarantee a corresponding benefit on the final event endpoint. Very small or unequal interim samples can yield unstable variances and rankings. Standardization also does not repair bias from informative missingness, measurement differences, or departures from the randomized comparison. As with the earlier material on binary early response as a selection variable, the selection criterion should be prespecified and its relationship to the final objective considered separately from the arithmetic of the interim ranking.

A Short Derivation for the Shared-Control Formula

Write each sample mean as its arm’s population mean plus a sampling error. Under independent sampling across arms, the active-arm sampling errors are independent of the control sampling error and of one another. Thus the covariance between contrasts is obtained by expanding their product:

$$ \begin{aligned} \operatorname{Cov}(\bar{x}_a-\bar{x}_0,\bar{x}_b-\bar{x}_0) &=\operatorname{Cov}(\bar{x}_a,\bar{x}_b) -\operatorname{Cov}(\bar{x}_a,\bar{x}_0)\\ &\quad-\operatorname{Cov}(\bar{x}_0,\bar{x}_b) +\operatorname{Var}(\bar{x}_0)\\ &=0-0-0+\frac{\sigma_0^2}{n_0} =\frac{\sigma_0^2}{n_0}. \end{aligned} $$

The zero terms follow from independence of distinct arm samples. The remaining term is the variance of the control mean. This derivation also shows why the covariance disappears when subtracting one treatment-control contrast from another: the two occurrences of the same control mean cancel algebraically.

Key takeaway: State whether an interim ranking uses raw contrasts, standardized mean differences, or precision-standardized statistics. Arm-specific variability can change the ranking, and shared-control contrasts are correlated even when the active-arm samples are independent.

Check Your Understanding

Use the assumptions and definitions in this tutorial. Answers follow each question so the calculations can be checked directly.

  1. Two arms each have a treatment-control contrast of 5 biomarker units. Arm A has standard error 1.0 and arm B has standard error 2.0. Which has the larger precision-standardized statistic? Solution: \(Z_A=5/1=5\) and \(Z_B=5/2=2.5\), so A has the larger statistic. The estimated effects in biomarker units are tied.
  2. Two treatment-control contrasts use the same control with estimated control variance 9 and \(n_0=45\). What is their estimated covariance? Solution: \(9/45=0.2000\), in squared biomarker units.
  3. Explain why dividing each contrast by its pair-specific pooled standard deviation can change an arm ranking. Solution: The pooled standard deviations can differ across active-control comparisons when arm variability differs. Since each effect estimate is divided by its own denominator, the ratios need not preserve the ordering of the raw contrasts.
  4. Two independent active arms have standard deviations 3 and 4 with 36 participants in each. Find the standard error of their difference in means. Solution: \(\sqrt{3^2/36+4^2/36}=\sqrt{25/36}=0.8333\) biomarker units. The common control is not included in this direct active-arm comparison.
  5. A biomarker has lower values when patients improve, but the selection rule treats positive treatment-control differences as favorable. What must be done before applying the rule? Solution: Reverse the scale or otherwise define the contrast so the favorable direction is consistently positive, then apply that convention to every arm and interpret the results on the stated scale.
Clinical Biostats