Why quantify departure from the Neyman share?
The Neyman allocation identifies the two-arm sample-size share that minimizes the variance of a difference in sample means when the total sample size is fixed. In practice, a design may use a different share: equal allocation is common, a target may be chosen for operational reasons, or estimated variances may lead to a share that differs from the variance-minimizing one. Knowing only that an allocation is “not optimal” is not enough to judge its consequences. The exact variance penalty shows how much precision is lost under the stated model.
This tutorial studies a fixed total of \(n\) participants assigned between treatment 1 and treatment 0. It holds the outcome variances fixed and compares allocations for estimating the contrast \(\mu_1-\mu_0\). It does not assign a value to other design considerations: the calculation is about contrast variance alone.
Model and assumptions
For the calculation, assume that the outcome in arm \(j\) has mean \(\mu_j\) and finite variance \(\sigma_j^2\), for \(j=0,1\). Outcomes are independent within and between arms, and the arm sample sizes are fixed in advance. Each arm uses the sample mean to estimate its mean. The total sample size \(n\) is fixed, and every participant has the same cost, so fixing \(n\) is also fixing the total cost. Both arms must have positive sample sizes: \(0<q<1\).
Under these assumptions, the two sample means are independent, with variances \(\sigma_1^2/n_1\) and \(\sigma_0^2/n_0\). Therefore the variance of their difference is the sum of those variances. For an interior Neyman share and the relative-penalty formula below, assume in addition that \(\sigma_1>0\) and \(\sigma_0>0\). This positivity condition matters: zero-variance cases need separate treatment.
The formula is an exact sampling variance under the stated model; it does not require normally distributed outcomes. Normality may matter for some inferential procedures, but it is not needed to add the variances of independent sample means. In a response-adaptive design, arm sizes may depend on previous outcomes. The expression here describes the fixed-allocation comparison; it is not, without further justification, a complete variance calculation for every outcome-dependent allocation process.
Deriving the minimizing share and the penalty identity
When both standard deviations are strictly positive, differentiate \(V(q)\) over \(0<q<1\). The second derivative is positive throughout the interval, so any stationary point is the unique minimum. The first derivative vanishes when the two positive terms in the resulting equality are equal.
Because \(V(q)\) tends to infinity as \(q\) approaches either endpoint when both variances are positive, the unique stationary point is the global minimum. Its allocation ratio is \(n_1/n_0=q^*/(1-q^*)=\sigma_1/\sigma_0\), and its variance is \((\sigma_1+\sigma_0)^2/n\).
A useful result is an exact expression for the variance away from this share. Divide \(V(q)\) by \(V(q^*)\), substitute \(q^*=\sigma_1/(\sigma_1+\sigma_0)\), and combine terms over the common denominator \(q(1-q)\). Subtracting 1 leaves a squared difference, so the excess is nonnegative.
Thus the relative variance penalty is exactly \((q-q^*)^2/[q(1-q)]\), and the variance itself is the minimum variance multiplied by one plus that penalty. The penalty is zero at the Neyman share and positive at every other interior share. Its denominator also explains why a fixed absolute departure can be especially costly near an allocation boundary.
Worked numerical comparisons
Worked Example: The Cost of Equal Allocation
Suppose, hypothetically, that a fixed total of \(n=180\) participants is split between two arms with outcome standard deviations \(\sigma_1=3\) units and \(\sigma_0=6\) units. Compare equal allocation with the Neyman share.
Both standard deviations are positive, the total is fixed, and the assumed arm outcomes are independent with finite variances. The Neyman share is \(q^*=3/(3+6)=1/3\), giving \(n_1=60\) and \(n_0=120\). Its contrast variance is \(V(q^*)=(3+6)^2/180=81/180=0.4500\) square units.
Equal allocation has \(q=0.5\), so each arm has 90 participants. Substitution in the variance formula gives \(V(0.5)=3^2/90+6^2/90=9/90+36/90=0.5000\) square units.
The variance ratio is \(0.5000/0.4500=1.1111\): equal allocation has an 11.11% larger contrast variance. The exact penalty identity gives the same result: \(1+(0.5-1/3)^2/[0.5(0.5)]=1+1/9=1.1111\). The standard-error ratio is \(\sqrt{1.1111}=\sqrt{10/9}\approx1.0541\), or about 5.409% greater. The standard error is measured in the outcome units.
Worked Example: A Near-Neyman Alternative
Suppose, hypothetically, \(n=200\), \(\sigma_1=4\) units, and \(\sigma_0=2\) units. Compare the Neyman share with an alternative share \(q=0.60\).
The assumptions for the fixed-allocation variance formula hold, and both standard deviations are positive. Here \(q^*=4/(4+2)=2/3\), so the Neyman allocation is \(n_1=133.33\) and \(n_0=66.67\) before any integer rounding. Its theoretical minimum variance is \(V(q^*)=6^2/200=0.1800\) square units.
At \(q=0.60\), the continuous allocation is \(n_1=120\) and \(n_0=80\). Its variance is \(V(0.60)=4^2/120+2^2/80=16/120+4/80=0.133333+0.050000=0.183333\) square units.
The variance ratio is \(0.183333/0.1800=1.01852\), a 1.852% variance increase. From the identity, the ratio is \(1+(0.60-2/3)^2/[0.60(0.40)]=1.01852\). The standard-error ratio is \(\sqrt{1.01852}\approx1.00922\), about 0.922% greater. The alternative is close to Neyman, so its penalty is small in this example.
Worked Example: An Unequal-Variance Case With Equal Allocation
Suppose, hypothetically, \(n=140\), \(\sigma_1=5\) units, and \(\sigma_0=2\) units. For illustration, compare the equal share with the Neyman share.
The fixed sample size, independent arm outcomes, finite variances, and positive standard deviations meet the conditions for the result. The Neyman share is \(q^*=5/7\), with minimum variance \(V(q^*)=(5+2)^2/140=49/140=0.3500\) square units.
At equal allocation, \(n_1=n_0=70\). Direct substitution gives \(V(0.5)=25/70+4/70=29/70=0.414286\) square units. The variance ratio is \((29/70)/(49/140)=58/49=1.18367\), a 18.367% increase. The identity gives \(1+(1/2-5/7)^2/(1/4)=1+9/49=58/49\).
The standard-error ratio is \(\sqrt{58/49}\approx1.08889\), so equal allocation has about 8.889% greater standard error. This example shows why the variance penalty and the standard-error penalty should not be reported as though they were the same percentage.
Reading the penalty and diagnosing a poor comparison
The penalty identity separates the effect of allocation mismatch from the scale of the outcome variances. The minimum variance \(V(q^*)\) sets the scale; the dimensionless factor \(1+(q-q^*)^2/[q(1-q)]\) tells how much the chosen share inflates it. Because the mismatch is squared, shares equally far above and below \(q^*\) have the same penalty only when their denominators \(q(1-q)\) are also equal. In general, departures on opposite sides need not cost the same amount.
A sound numerical comparison should report the target and alternative shares, the arm counts or their continuous equivalents, the variance in squared outcome units, and the standard-error ratio if discussing precision on the original outcome scale. Check that shares sum to one, that the sample sizes add to \(n\), and that the direct variance calculation agrees with the penalty identity. A disagreement usually signals a switched arm label, an incorrect denominator, or confusion between a variance ratio and a standard-error ratio.
The comparison is only as informative as its inputs. If the variances used to calculate \(q^*\) are uncertain or misspecified, the share called “Neyman” may not minimize the true variance. Likewise, if arm sizes are rounded to integers, the realized share may differ slightly from the theoretical continuous target; use the actual counts in the variance formula to assess that allocation. The displayed exact identity applies to the continuous share formulation. It does not replace evaluating feasible integer allocations or a design with additional constraints.
Boundary cases and limits
The assumption \(\sigma_1,\sigma_0>0\) is essential to the interior minimization argument. If \(\sigma_1=0\) and \(\sigma_0>0\), then \(V(q)=\sigma_0^2/[n(1-q)]\), which increases with \(q\). Its infimum over \(0<q<1\) is approached as \(q\) tends to zero, but is not attained at any share in that interval. At \(q=0\), the treatment-1 mean cannot be estimated from a sample, so that endpoint is not a valid two-arm contrast estimator. If \(\sigma_0=0\) and \(\sigma_1>0\), the corresponding infimum is approached as \(q\) tends to one and is likewise not attained in the interior. If both variances are zero, \(V(q)=0\) for every interior share, and the positive-variance formula for \(q^*\) does not determine a unique target.
In applications, minimum arm sizes, integer counts, unequal costs, or other design objectives can alter the relevant comparison. Those cases require their own constraints or objective functions; the unconstrained, equal-cost result here should not be applied as if those features were absent. In particular, a variance-optimal share is not automatically the preferred allocation for every scientific or operational purpose.
Common mistake
A second useful check is not to substitute a zero standard deviation into the interior formula and claim an interior optimum at an endpoint. The derivation requires both standard deviations to be positive; the zero-variance cases have the boundary behavior described above.
Check Your Understanding
Use the fixed-total, independent-arm variance model unless a question states otherwise. The solutions are included below each question.
- For \(\sigma_1=2\) and \(\sigma_0=3\), what is the Neyman share \(q^*\) and the corresponding allocation ratio \(n_1/n_0\)?
- Suppose \(q^*=0.40\) and an alternative uses \(q=0.50\). Calculate its relative variance and standard-error ratios using the penalty identity.
- Why does an 11.11% variance increase correspond to about a 5.409% standard-error increase rather than an 11.11% increase?
- If \(\sigma_1=0\) and \(\sigma_0>0\), what happens to \(V(q)\) as \(q\) increases, and why is the limiting share not an interior minimizing share?
Solutions. (1) \(q^*=2/(2+3)=0.40\), and \(n_1/n_0=2/3\). (2) The relative variance is \(1+(0.50-0.40)^2/[0.50(0.50)]=1+0.04=1.04\). The standard-error ratio is \(\sqrt{1.04}\approx1.0198\), about 1.98% greater. (3) Standard error is the square root of variance, so the ratio is \(\sqrt{1.1111}\approx1.0541\), which is about 5.409% above 1. (4) \(V(q)=\sigma_0^2/[n(1-q)]\) increases as \(q\) increases. Its infimum is approached as \(q\) tends to zero, but zero is outside the required interior shares and would leave no treatment-1 sample for estimating the contrast.