Why Examine the Null Case?
A response-adaptive design is often motivated by the possibility of directing more patients to a treatment that appears to perform better. But when treatments have equal outcomes, apparent differences arise from random variation. The allocation rule can react to those chance differences, changing who is assigned to each arm and how later estimates behave. The null case is therefore a stress test: it asks whether the full design and analysis procedure behave properly when there is no treatment advantage to discover.
A null evaluation is not a claim that treatments are equal in practice. It is a controlled operating-characteristic assessment. It should preserve the proposed assignment rule, outcome model, enrollment size, and analysis procedure, and repeat the entire trial many times. Looking at one simulated trial, or at average allocation alone, cannot establish how often the analysis rejects a true null or how variable its estimates are.
For two treatments, let \(n_A\) and \(n_B\) be their final sample sizes, with total \(N=n_A+n_B\). One simple measure of allocation imbalance is \(|n_A-n_B|\). A related descriptive quantity is the equal-allocation effective sample size, \(N_{\mathrm{eff}}=4n_An_B/N\). It equals \(N\) for equal allocation and decreases as the allocation becomes more uneven. It describes the allocation's balance; it is not, by itself, a correction for adaptive sampling or a guarantee about estimator variance.
Specify the Null and the Whole Procedure
A useful null stress test starts by specifying exactly what is null. In a global null, each patient's potential outcomes have the same distribution under A and B. In the examples below, outcomes are independent Bernoulli variables with common success probability \(p=0.30\), regardless of assigned treatment. Assignment may depend on outcomes already observed, but the next patient's outcome is independent of that history. We also assume a fixed enrollment total, complete outcome ascertainment, no interference between patients, and a fully specified assignment rule.
These assumptions make the simulated trials reproducible and give a clear target: the null distribution of the complete procedure. If the real design includes delayed outcomes, varying enrollment, clustering, or a different null (for example, equality of means but not of outcome distributions), those features must be represented in its own evaluation. A result under one null model does not establish operating characteristics under every plausible data-generating process.
Each replicate must rerun the adaptive decisions from the beginning. For each trial, record the final allocation, the treatment-effect estimate, its uncertainty measure, and whether the prespecified test rejects \(H_0\). Across replicates, summarize the rejection fraction, the mean estimate, its standard deviation, and allocation summaries. This separates four related but distinct questions: Was allocation uneven? Was the estimate centered at the null value? How much did estimates vary? Did the test reject too often?
An Analytic Illustration of Adaptive Bias
The following deliberately small example isolates one mechanism by which an adaptive rule can affect a familiar estimator. It is not a recommended trial design: with only three patients, conventional inference would not be meaningful. Under the global null, let all outcomes be independent \(N(0,1)\). Randomly assign patient 1 to A or B with equal probability, then assign patient 2 to the opposite arm. Let \(X_A\) denote the outcome of the patient assigned to A. Assign patient 3 to A if \(X_A>0\), and to B otherwise. This specifies the rule for every possible observed history.
Worked Example: A Difference in Means Under a Complete Adaptive Rule
Write \(U\) for the initial B patient's outcome and \(V\) for patient 3's outcome. Under the null, \(X_A\), \(U\), and \(V\) are independent standard normal variables. The estimator is \(D=\bar{X}_A-\bar{X}_B\), in standard-deviation units.
If \(X_A>0\), the final sizes are \((n_A,n_B)=(2,1)\), so \(D=(X_A+V)/2-U\). Since \(E(X_A\mid X_A>0)=\sqrt{2/\pi}\), while the conditional means of \(U\) and \(V\) remain zero,
If \(X_A\leq 0\), the final sizes are \((1,2)\), and \(D=X_A-(U+V)/2\). By symmetry, \(E(X_A\mid X_A\leq 0)=-\sqrt{2/\pi}\). Therefore,
Each branch has probability one-half. Applying the law of total expectation gives
The true mean difference is zero, so this estimator has negative bias in this specified toy design. The result is not caused by unequal outcome distributions: the outcomes have the same distribution in both arms. It arises because the assignment decision uses the A patient's observed outcome, and that outcome also enters the final estimate. This derivation also illustrates why it is essential to define every branch of an adaptive rule before calculating its consequences.
Repeated-Trial Evaluation of Type I Error
For a larger illustrative evaluation, consider this fully specified response-adaptive rule with \(N=100\) binary outcomes. Randomize the first 20 patients independently, with probability \(0.5\) to each arm. After that, let \(n_j\) and \(s_j\) be the number assigned and the number of successes observed in arm \(j\), for \(j\in\{A,B\}\). Calculate the smoothed response rate \(r_j=(s_j+1)/(n_j+2)\), and assign the next patient to A with probability \(q_A=0.1+0.8r_A/(r_A+r_B)\); assign to B with probability \(1-q_A\). Under the stated Bernoulli null, each new outcome is generated independently with success probability \(0.30\), after its assignment. This keeps both assignment probabilities between 0.1 and 0.9.
For illustration, suppose a hypothetical simulation of 20,000 complete trials gives the summary below. These are hypothetical Monte Carlo results for this specified scenario, not empirical findings or universal properties of response-adaptive randomization. The prespecified analysis is the ordinary two-proportion pooled score test: calculate the pooled proportion \(\hat{p}\), compute \(Z=(\hat{p}_A-\hat{p}_B)/\sqrt{\hat{p}(1-\hat{p})(1/n_A+1/n_B)}\), and reject at the two-sided 0.05 level when \(|Z|>1.96\). The null outcomes are binary; if either arm has zero patients or the pooled proportion is 0 or 1, define the test decision as non-rejection, and otherwise use the stated rule.
| Summary across 20,000 trials | Hypothetical result |
|---|---|
| Mean number assigned to A | 50.0 |
| Mean absolute allocation difference, \(|n_A-n_B|\) | 14.8 patients |
| Mean equal-allocation effective sample size | 96.7 patients |
| Mean estimated difference, \(\hat{p}_A-\hat{p}_B\) | 0.0004 |
| Across-trial standard deviation of the estimated difference | 0.101 |
| Trials rejecting at the nominal 0.05 level | 1,180 of 20,000 |
Worked Example: Estimating Type I Error From Replicates
The Monte Carlo estimate of the type I error rate is the fraction of null trials that rejected:
Thus this hypothetical evaluation estimates a 5.90% rejection rate, above the nominal 5%. The approximate Monte Carlo standard error is \(\sqrt{\widehat{\alpha}(1-\widehat{\alpha})/R}\), where \(R=20{,}000\) is the number of simulated trials:
An approximate 95% Monte Carlo interval is \(0.0590\pm1.96(0.00167)\), or about 0.0557 to 0.0623. This interval describes simulation error in estimating the rejection probability for this specified scenario. It is not a confidence interval for a treatment effect. The result suggests that this particular test-and-design combination may be liberal under this null; additional checks and a calibrated analysis may be warranted.
Why Arm Sizes Alone Do Not Determine Variability
For independent outcomes with common variance \(\sigma^2\), the familiar variance formula for a difference in sample means is \(\sigma^2(1/n_A+1/n_B)\) when arm sizes are fixed, or are otherwise noninformative about the outcomes entering the means. In response-adaptive allocation, final arm sizes can depend on previously observed outcomes. Conditioning on final sizes can therefore select outcomes and change the conditional distribution of the estimator. One must not infer its conditional variance from the two final counts alone.
The appropriate null stress test estimates variability across complete repeated trials under the specified rule. In the hypothetical evaluation, the reported standard deviation of 0.101 is the standard deviation of the 20,000 trial-specific estimates. It summarizes their across-trial spread under that simulation model; it is not obtained by substituting a typical pair of arm sizes into a fixed-allocation formula. If a conditional variance for a particular allocation history is needed, it must be derived for that assignment rule or estimated by simulation conditional on that history.
Worked Example: Comparing Null Rejection Rates
Suppose a separate hypothetical evaluation uses the same Bernoulli null, total sample size, and pooled score test, but allocates exactly 50 patients to each arm. Assume 1,024 of 20,000 fixed-allocation trials reject. The estimated type I error is \(1,024/20{,}000=0.0512\). Comparing it with the adaptive estimate gives a difference of \(0.0590-0.0512=0.0078\), or 0.78 percentage points.
If the two sets of simulation runs are independent, the Monte Carlo standard error of this difference is approximately
An approximate 95% Monte Carlo interval for the difference is \(0.0078\pm1.96(0.00228)\), or about 0.0033 to 0.0123. In this hypothetical comparison, the simulated adaptive procedure has a higher rejection rate than the fixed-allocation procedure, with Monte Carlo uncertainty alone not accounting for the whole difference. This comparison is informative only for the stated outcome model, assignment rules, and test. It does not show that every adaptive design inflates type I error.
Diagnostics, Comparisons, and Limitations
A null stress test is most useful when it reports more than a single rejection fraction. Examine allocation distributions, not just their mean: a mean of 50 patients per arm can coexist with substantial trial-to-trial imbalance. Check the frequency of very small arm sizes, the estimator's mean and spread, and whether standard errors track empirical variability. Also inspect the assignment probabilities and the timing of updates. Unexpectedly extreme allocations, unstable estimates, or rejection rates that change markedly with simulation settings are reasons to investigate the implementation and assumptions.
The main comparison is with fixed allocation under the same null outcome model and analysis. Fixed allocation gives a useful benchmark because it does not respond to observed outcomes; it is not automatically a perfect standard if the chosen test has poor small-sample calibration. A second comparison is between the usual analysis and an analysis calibrated to the complete adaptive design. For example, a simulation-based null calibration can reproduce the assignment and analysis procedures under the stated null. Any calibration must preserve the dependence created by the actual adaptive rule rather than treating the final arm counts as if they had been fixed in advance.
The null case does not answer every design question. It does not show how often a genuinely effective treatment will be identified, nor how many patients will benefit when treatments differ. It also cannot establish robustness to alternative outcome distributions, missingness, or implementation errors unless those features are included. A favorable null result is evidence about one specified operating condition, not a general guarantee. Conversely, a poor result identifies a problem with the combination of assignment rule and analysis under the tested conditions; it does not by itself identify which component caused it.
Common Mistakes
- Calling balanced average allocation proof of good behavior. A mean of 50 per arm does not rule out large imbalances in individual trials. Inspect the distribution of allocations and estimator variability.
- Using the fixed-size variance formula after outcome-dependent assignment. Final arm sizes may select observed outcomes. Estimate variability across full simulated trials or derive it under the complete rule.
- Interpreting a Monte Carlo rejection fraction as exact. The fraction has simulation error. Report the number of replicates and its Monte Carlo standard error.
- Leaving branches of the allocation rule unspecified. Conditional expectations and overall bias cannot be derived reliably until the rule is defined for every reachable history.
- Treating the null as the only operating condition that matters. It tests false-positive behavior under a stated model, but says nothing by itself about performance when treatment outcomes differ.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- In the three-patient example, what is the bias of \(D\) under the stated null?
- In a simulation with 20,000 null trials and 1,160 rejections, calculate the estimated type I error rate.
- For \(N=100\), what is \(N_{\mathrm{eff}}\) when 50 patients are assigned to each arm?
- Why can conditioning on final arm sizes fail to justify the usual fixed-allocation variance formula?
- What does a null-case rejection fraction estimate, and what does it not establish?
Solutions. (1) The bias is \(E(D)-0=-\frac{1}{4}\sqrt{2/\pi}\approx-0.1995\) standard-deviation units. (2) The estimate is \(1,160/20{,}000=0.0580\), or 5.80%. (3) \(N_{\mathrm{eff}}=4(50)(50)/100=100\) patients. (4) In outcome-adaptive allocation, final sizes depend on observed outcomes, so conditioning on those sizes can select outcomes and alter the estimator's conditional variance. (5) It estimates the probability that the specified test-and-design procedure rejects under the specified null model. It does not establish type I error under other null models or performance when treatments differ.