Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Cumulative Test Statistics at Interim Looks

The Canonical Joint Distribution · Tutorial 31 of 699

Cumulative Test Statistics at Interim Looks

Learn to calculate and interpret cumulative test statistics at successive interim looks, and to separate accumulated evidence from the contribution of a new data block.

Advanced 10 min read

What You'll Learn

  • Define a cumulative test statistic using all observations available at a scheduled look.
  • Distinguish the cumulative statistic from a statistic calculated on only the newest data block.
  • Derive the exact update identity for a normal-mean Z statistic with known variance.
  • Calculate cumulative and block statistics at successive looks using hypothetical data.
  • Recognize why cumulative statistics can rise or fall as observations accrue.
  • Diagnose common errors in monitoring calculations and interpretation.

What Is Cumulative at an Interim Look?

At an interim look, investigators have data from participants or units enrolled so far. The test statistic at that look is ordinarily calculated using all data accrued to date, not just the observations added since the preceding look. This distinction matters: the new observations contribute fresh information, but the cumulative statistic summarizes the evidence in both the old and new observations.

Suppose the monitoring schedule has looks at cumulative sample sizes \(n_1<n_2<\cdots<n_K\). Let \(n_k\) denote the total number of observations available at look \(k\), and let \(m_k=n_k-n_{k-1}\) denote the number added since look \(k-1\), with \(n_0=0\). The \(m_k\) observations form the new data block at look \(k\). A cumulative statistic uses observations \(1,\ldots,n_k\); a block statistic uses only observations \(n_{k-1}+1,\ldots,n_k\).

Definition: A cumulative test statistic at look \(k\) is the test statistic calculated from all observations available through that scheduled look. A new-block statistic at look \(k\) is a statistic calculated using only the observations accrued since the previous look. They answer different questions and are not interchangeable.

The cumulative statistic is the one associated with the test at that look. The new-block statistic describes the latest data alone. The cumulative statistic may incorporate the block statistic, but it is not generally equal to it, nor is it generally the unweighted sum of earlier and new statistics. The monitoring schedule and the rule for comparing look-specific statistics with thresholds are separate from the definition of the statistic itself.

A Normal-Mean Model and the Update Identity

To make the distinction precise, consider independent observations \(X_1,X_2,\ldots\) from a normal population with known standard deviation \(\sigma\). The monitoring question is whether the population mean exceeds a specified null value \(\mu_0\):

$$ H_0:\mu=\mu_0 \qquad\text{versus}\qquad H_a:\mu>\mu_0. $$

Assume the observation distribution, independence, known \(\sigma\), and the ordered look sizes are appropriate for the analysis. At look \(k\), the sample mean \(\bar{x}_k\) is computed from all \(n_k\) observations. The cumulative one-sided Z statistic is

$$ Z_k=\frac{\sqrt{n_k}(\bar{x}_k-\mu_0)}{\sigma}. $$

For the new block at look \(k\), let \(\bar{x}_{B,k}\) be its mean. Its null-centered standardized statistic is

$$ V_k=\frac{\sqrt{m_k}(\bar{x}_{B,k}-\mu_0)}{\sigma}. $$

The block statistic is defined here on the same scale and against the same null value as the cumulative statistic. It is not the cumulative test at look \(k\). Since the cumulative sample mean is a weighted average of the earlier and newly added observations, the two statistics satisfy an exact update identity.

1
Write the cumulative sum.
The centered sum through look \(k\) is the previous centered sum plus the new block's centered sum: \(n_k(\bar{x}_k-\mu_0)=n_{k-1}(\bar{x}_{k-1}-\mu_0)+m_k(\bar{x}_{B,k}-\mu_0)\).
2
Standardize by the cumulative standard error.
Divide the centered sum by \(\sigma\sqrt{n_k}\), the denominator corresponding to all \(n_k\) observations.
3
Express each part using its own standardized statistic.
Replace the previous and new-block centered sums by their definitions to obtain the update formula below.
$$ Z_k =\sqrt{\frac{n_{k-1}}{n_k}}\,Z_{k-1} +\sqrt{\frac{m_k}{n_k}}\,V_k. $$

For the first look, \(n_0=0\), so the statistic is calculated directly from the first block. At later looks, the first term carries forward the contribution of the observations already collected, while the second term adds the contribution of the new block. The weights depend on the relative block sizes. This identity is an algebraic consequence of combining the data; it does not say that a cumulative statistic is the same as a block statistic.

Worked Example: A Cumulative Statistic Can Fall After New Data

Suppose \(\sigma=12\), and the null value is \(\mu_0=0\). At the first look, \(n_1=4\) observations have mean \(9\) units above \(\mu_0\). Two more observations are then collected, and their mean is \(6\) units below \(\mu_0\). Calculate the first-look statistic, the new-block statistic, and the cumulative statistic at the second look.

At look 1, the cumulative sample size is 4, so

$$ Z_1=\frac{\sqrt{4}(9)}{12}=1.500. $$

The new block has \(m_2=2\) observations and mean difference \(-6\), giving

$$ V_2=\frac{\sqrt{2}(-6)}{12}=-0.7071. $$

At look 2, \(n_2=6\). The cumulative mean difference is the weighted average of the two block means:

$$ \bar{x}_2-\mu_0 =\frac{4(9)+2(-6)}{6} =\frac{24}{6} =4. $$

Thus the cumulative statistic is

$$ Z_2=\frac{\sqrt{6}(4)}{12}=0.8165. $$

The update identity gives the same result:

$$ Z_2 =\sqrt{\frac{4}{6}}(1.500)+\sqrt{\frac{2}{6}}(-0.7071) \approx 1.2247-0.4082 =0.8165. $$

The latest block points below the null, so it pulls the cumulative statistic down from 1.500 to about 0.8165. The look-2 statistic is still based on all six observations; it is not the block value \(-0.7071\).

Calculating from the Accumulated Data

The update identity is useful for seeing how a new block contributes, but the direct cumulative calculation is often the safest check: combine all observations available at the look, calculate their mean, and use the cumulative sample size in the standard error. Both methods should agree when the model and inputs are the same.

Worked Example: Checking an Update Against the Cumulative Mean

Suppose \(\sigma=12\), \(\mu_0=0\), and the first nine observations have mean difference \(4\). At the next look, four new observations have mean difference \(-2\). Calculate both the direct cumulative statistic and its update from the first-look and block statistics.

At the first look,

$$ Z_1=\frac{\sqrt{9}(4)}{12}=1.000. $$

For the new block,

$$ V_2=\frac{\sqrt{4}(-2)}{12}=-0.3333. $$

The second look has \(n_2=9+4=13\) observations. Their combined mean difference is

$$ \bar{x}_2-\mu_0 =\frac{9(4)+4(-2)}{13} =\frac{28}{13} \approx 2.1538. $$

Therefore, the direct cumulative calculation is

$$ Z_2 =\frac{\sqrt{13}(28/13)}{12} =\frac{28}{12\sqrt{13}} \approx \frac{28}{43.2666} =0.6472. $$

Using the update identity gives

$$ Z_2 =\sqrt{\frac{9}{13}}(1.000)+\sqrt{\frac{4}{13}}(-0.3333) =\frac{3}{\sqrt{13}}-\frac{8}{12\sqrt{13}} =\frac{28}{12\sqrt{13}} \approx 0.6472. $$

The new block is unfavorable to the positive direction of the alternative, so the cumulative statistic falls from 1.000 to 0.6472. The value \(0.6472\) is the statistic from all 13 observations; \(-0.3333\) is only the statistic from the four newly observed units.

Why Look Statistics Need Not Move in One Direction

A larger cumulative sample does not guarantee a larger test statistic. New observations may support the alternative, oppose it, or be close to the null. Consequently, the cumulative sample mean and its standardized statistic can move either way. Sample size affects the standard error, but the observed effect estimate also changes as data accrue.

Worked Example: A Larger Sample Does Not Ensure a Larger Z Statistic

Suppose \(\sigma=10\) and \(\mu_0=0\). At the first look, 100 observations have mean difference \(1.6\). A new block of 50 observations has mean difference \(-1.6\). Find the cumulative statistic at each look and verify the second look using the update identity.

The first-look statistic is

$$ Z_1=\frac{\sqrt{100}(1.6)}{10}=1.600. $$

For the new block,

$$ V_2=\frac{\sqrt{50}(-1.6)}{10}\approx -1.13137. $$

The cumulative mean difference at look 2 is

$$ \bar{x}_2-\mu_0 =\frac{100(1.6)+50(-1.6)}{150} =\frac{80}{150} \approx 0.5333. $$

So the direct cumulative statistic is

$$ Z_2=\frac{\sqrt{150}(0.5333\ldots)}{10} \approx 0.6532. $$

The update identity confirms the result:

$$ Z_2 =\sqrt{\frac{100}{150}}(1.600) +\sqrt{\frac{50}{150}}(-1.13137) \approx 1.3064-0.6532 =0.6532. $$

Although the cumulative sample size increases from 100 to 150, the statistic falls because the new block's mean difference is in the opposite direction. The conclusion at a look must be based on the prespecified decision rule and the cumulative statistic, not on the sample-size increase alone.

Diagnostics, Alternatives, and Limitations

A useful calculation check is to compute the cumulative statistic in two ways: directly from the combined data and by updating the previous statistic with the new-block contribution. A disagreement usually signals a bookkeeping or arithmetic error. Check that the cumulative sample size is the sum of all block sizes, that the cumulative mean weights each block by its size, and that the standard error uses the cumulative sample size.

The main tempting alternative is to use only the latest block's statistic as though it were the cumulative test. That discards the earlier observations and answers a different question. A block-only analysis may be useful for describing what the newest data show, but it should not silently replace the prespecified cumulative analysis. Conversely, simply repeating the first-look statistic at later looks ignores the new data. Neither shortcut produces the cumulative statistic defined here.

QuantityData usedInterpretation
Cumulative statistic \(Z_k\)All observations through look \(k\)Standardized evidence at the current look
New-block statistic \(V_k\)Only observations since look \(k-1\)Standardized contribution from the latest block alone
Cumulative mean \(\bar{x}_k\)All observations through look \(k\)Estimated mean at the current look

The simple identity shown here relies on a one-sample normal-mean model with known \(\sigma\), independent observations, and the same null value used for each block and cumulative calculation. With an estimated standard deviation, a different endpoint, or a more complex model, the test statistic and its updating calculation must be derived for that setting; this formula cannot be assumed to apply unchanged. The definition of cumulative data remains straightforward, but the exact arithmetic may differ.

Finally, a cumulative statistic is not itself a calibrated repeated-testing procedure. Repeatedly comparing it with an ordinary single-look cutoff can inflate the overall type I error rate. The earlier material on repeated-testing thresholds explains why decision boundaries require calibration. Here the central task is to identify correctly which data enter each look-specific statistic.

  • Wrong: “The look-2 statistic is the statistic from the 50 new observations.” Correction: That is a block statistic; the look-2 cumulative statistic uses all observations through look 2.
  • Wrong: “The cumulative statistic must increase because the sample size increased.” Correction: The new block can move the estimate toward or away from the null, so the statistic can rise or fall.
  • Wrong: “I can add the two Z statistics directly.” Correction: Their contributions have sample-size-dependent weights; use the derived update identity or recompute from the combined data.
Key takeaway: At look \(k\), calculate the test statistic from all \(n_k\) observations available. Keep the latest block statistic separate: it describes the new observations, while the cumulative statistic combines their contribution with the accumulated data.

Check Your Understanding

Try each question, then check the solutions below.

  1. At look 2, the cumulative sample size is 12 and the previous look had 8 observations. How many observations are in the new block, and which observations enter the cumulative statistic?
  2. With known \(\sigma\), suppose \(n_1=4\), \(n_2=6\), \(Z_1=1.2\), and the new-block statistic is \(V_2=-0.4\). Calculate \(Z_2\) using the update identity.
  3. Explain why a block statistic of \(-0.4\) does not mean that the cumulative statistic at that look is \(-0.4\).
  4. For the model in this tutorial, state two conditions needed for the displayed Z-statistic formula to be exact.

Solutions. (1) The new block has \(m_2=12-8=4\) observations. The cumulative statistic uses all 12 observations, including the previous 8 and the new 4. (2) The update is \(Z_2=\sqrt{8/12}(1.2)+\sqrt{4/12}(-0.4)\approx 0.9798-0.2309=0.7489\). (3) The block statistic uses only the newest observations; the cumulative statistic also includes the earlier observations, weighted as in the update identity. (4) For example, observations must be independent and normally distributed, and \(\sigma\) must be known. The specified null value and the sample size used in the standard error must also match the analysis.

Clinical Biostats