Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Information When Variance Is Estimated

Information Fraction, Sample Size and Events · Tutorial 85 of 699

Information When Variance Is Estimated

See how estimating the variance affects information and precision for a normal mean, and when a t statistic is needed instead of a normal approximation.

Intermediate 12 min read

What You'll Learn

  • Distinguish planned information based on an assumed variance from information estimated using a sample variance.
  • Calculate standard errors and plug-in information for a normal mean.
  • Explain why a t statistic, rather than a standard normal statistic, accounts for variance estimation in a finite sample.
  • Construct and interpret a t confidence interval and t test under the normal model.
  • Use a t-based planning calculation to choose a sample size for a target interval half-width.
  • Recognize when normal-theory precision calculations may be unreliable.

Why Estimating the Variance Changes Precision

For a normal mean, precision depends on the variance: with \(n\) independent observations, the variance of the sample mean is \(\sigma^2/n\). If \(\sigma\) is known, its standard error can be calculated directly. In most applications, however, \(\sigma\) is unknown and must be estimated from the data. Replacing \(\sigma\) by the sample standard deviation \(s\) gives a useful estimate of precision, but it also means the standard error itself is estimated and variable.

This distinction matters in information-based planning. A design calculation usually starts with a planning value for the standard deviation and uses it to project how much information a proposed sample size will provide. After data are collected, the observed \(s\) gives a data-based estimate of precision. The two quantities answer related but different questions: one is a planning calculation, and the other reflects the variability observed in the sample.

Definition: Under independent observations from a population with variance \(\sigma^2\), the variance of the sample mean is \(\sigma^2/n\), so its information for estimating \(\mu\) is \(I=n/\sigma^2\). If \(\sigma^2\) is unknown, a plug-in estimate of that information is \(\hat{I}=n/s^2\), and the estimated standard error of \(\bar{x}\) is \(SE=s/\sqrt{n}\). The plug-in value \(\hat{I}\) is an estimate, not the known information itself.

For a mean measured in milligrams, for example, information has units of inverse milligrams squared, and a standard error has units of milligrams. Larger information corresponds to a smaller standard error when the model and estimand are fixed. But because \(s^2\) varies across samples, \(\hat{I}\) and \(SE\) vary too.

Planning Information and Observed Information

A planning calculation generally specifies an assumed standard deviation \(\sigma_{\text{plan}}\). The planned information at sample size \(n\) is \(n/\sigma_{\text{plan}}^2\). This is a model-based projection used to compare sample sizes or set a target; it does not claim that the future sample variance will equal the planning variance. Once outcomes are observed, \(n/s^2\) summarizes the information estimated from that sample.

The usual normal-theory formulas require independent observations with a common variance and a normal outcome distribution when exact small-sample inference is intended. The sample size is fixed for the calculation, and there are no missing observations or other complications in this simple model. If observations are dependent, variances differ across observations, or the outcome model is not appropriate, the formulas may not describe the actual precision. Those issues cannot be repaired merely by substituting a different value for \(s\).

Worked Example: Comparing Planned and Estimated Information

Suppose a hypothetical study measures a continuous outcome in milligrams. Before data collection, investigators assume a standard deviation of \(10\) mg and plan to enroll \(n=25\) independent participants. In the completed sample, suppose the sample standard deviation is \(12\) mg. For this calculation, assume the observations are independent and come from a common normal distribution.

Using the planning standard deviation, planned information is

$$ I_{\text{plan}}=\frac{n}{\sigma_{\text{plan}}^2} =\frac{25}{10^2}=0.2500\ \text{mg}^{-2}. $$

The planned standard error is \(10/\sqrt{25}=2.00\) mg. Using the observed sample standard deviation instead gives

$$ \hat{I}=\frac{25}{12^2}=\frac{25}{144} =0.1736\ \text{mg}^{-2}, \qquad SE=\frac{12}{\sqrt{25}}=2.40\ \text{mg}. $$

The observed sample therefore implies less estimated information and a larger estimated standard error than the planning calculation. This difference does not establish that the population standard deviation is exactly \(12\) mg; \(s\) is itself a sample-based estimate.

Why the t Statistic Accounts for Variance Estimation

If \(\sigma\) is known and the normal model holds, standardizing \(\bar{x}\) by \(\sigma/\sqrt{n}\) gives a standard normal statistic. When \(\sigma\) is unknown, the natural statistic uses \(s/\sqrt{n}\) instead. The extra variation from estimating \(\sigma\) changes the reference distribution.

For independent observations \(X_1,\ldots,X_n\) from \(N(\mu,\sigma^2)\), with \(n\geq 2\), the sample mean is independent of the sample variance. In particular, \(Z=(\bar{x}-\mu)/(\sigma/\sqrt{n})\) has a standard normal distribution, and \(U=(n-1)s^2/\sigma^2\) has a chi-square distribution with \(n-1\) degrees of freedom. Their independence gives the following result.

Result: Under the independent, identically distributed normal model with unknown \(\sigma\),
$$ T=\frac{\bar{x}-\mu}{s/\sqrt{n}} =\frac{Z}{\sqrt{U/(n-1)}}\sim t_{n-1}. $$

The denominator includes the random sample standard deviation, so \(T\) has heavier tails than a standard normal variable for finite \(n\). As \(n\) grows, the \(t\) distribution approaches the standard normal distribution. The exact result depends on the normal model; for small samples, substantial skewness or outliers can make that model questionable.

1
Specify the model and target.
For a one-sample mean, state the normal, independent, common-variance model and the value \(\mu_0\) being tested, or the confidence level required.
2
Estimate the standard error.
Calculate \(SE=s/\sqrt{n}\), using the sample standard deviation rather than treating it as a known population value.
3
Use the t reference distribution.
For a one-sample mean, use \(n-1\) degrees of freedom to obtain a t critical value, p-value, or confidence interval.

For a two-sided test of \(H_0:\mu=\mu_0\), the statistic is \(t=(\bar{x}-\mu_0)/(s/\sqrt{n})\), and its p-value under the null is \(P(|T|\geq |t|)\) for \(T\sim t_{n-1}\). A \(100(1-\alpha)\%\) confidence interval is \(\bar{x}\pm t_{1-\alpha/2,n-1}s/\sqrt{n}\). The interpretation of the interval is about the procedure's long-run coverage under its assumptions, not a probability assigned to the fixed parameter after observing the interval.

Worked Example: A One-Sample t Test and Confidence Interval

Consider hypothetical independent normal measurements in a sample of \(n=16\), with \(\bar{x}=54\) mg and \(s=8\) mg. Test \(H_0:\mu=50\) mg against a two-sided alternative. Assume the observations are independent, have a common variance, and follow a normal distribution; these assumptions make the t reference distribution exact for this calculation.

The estimated standard error is \(8/\sqrt{16}=2\) mg. Thus,

$$ t=\frac{54-50}{2}=2.000,\qquad df=16-1=15. $$

The two-sided p-value is \(2P(T_{15}\geq 2.000)\approx 0.0639\). At significance level \(\alpha=0.05\), the result is not statistically significant: the statistic does not reach the two-sided 5% t critical value. This does not show that \(H_0\) is true; it means these data do not provide sufficient evidence to reject it at that level.

For a 95% confidence interval, \(t_{0.975,15}\approx 2.131\). The half-width is \(2.131(2)=4.262\) mg, so the interval is

$$ 54\pm 4.262\ \text{mg} =(49.738,\ 58.262)\ \text{mg}. $$

Under the stated model, this method has 95% coverage in repeated samples. The interval includes \(50\) mg, consistent with the two-sided test not rejecting at the 5% level.

Using Estimated Variance in Precision Planning

For known \(\sigma\), a normal-approximation calculation targets a half-width \(h\) by choosing \(n\) so that \(z_{1-\alpha/2}\sigma/\sqrt{n}\leq h\), or equivalently \(n\geq (z_{1-\alpha/2}\sigma/h)^2\). When variance must be estimated, a common planning adjustment uses a t critical value:

$$ t_{1-\alpha/2,n-1}\frac{\sigma_{\text{plan}}}{\sqrt{n}}\leq h. $$

Because the t critical value depends on \(n\), this is checked over integer sample sizes rather than solved by a single rearrangement. The calculation is still a planning approximation: it inserts an assumed \(\sigma_{\text{plan}}\) in place of the future random \(s\). It does not guarantee that the realized confidence interval will have half-width at most \(h\). The future sample standard deviation may be larger or smaller than the planning value.

Worked Example: Choosing n for a Target Half-Width

Suppose a hypothetical study plans to estimate a normal mean with a 95% confidence interval whose half-width is at most \(4\) mg. Use a planning standard deviation of \(\sigma_{\text{plan}}=12\) mg. Assume the future observations will be independent and normally distributed with a common variance. First, the normal approximation gives

$$ n\geq \left(\frac{1.96(12)}{4}\right)^2 =34.5744, $$

which suggests \(n=35\). To account approximately for estimating the variance, check the t-based half-width at candidate integer sample sizes. At \(n=37\), \(t_{0.975,36}\approx2.028094\), so

$$ 2.028094\frac{12}{\sqrt{37}}\approx 4.001\ \text{mg}. $$

This slightly exceeds the \(4\) mg target, so \(n=37\) does not meet it. At \(n=38\), \(t_{0.975,37}\approx2.026192\), giving

$$ 2.026192\frac{12}{\sqrt{38}}\approx 3.944\ \text{mg}. $$

Thus the t-based planning approximation selects \(n=38\). The half-width decreases as \(n\) increases here, since both \(1/\sqrt{n}\) and the t critical value decrease; therefore the failure at \(n=37\) rules out smaller sample sizes as well. The result is conditional on the chosen planning standard deviation and is not a guarantee about the realized interval.

Choosing Between Normal and t Calculations

SituationPrecision or reference distributionImportant qualification
Population standard deviation known; normal model\(SE=\sigma/\sqrt{n}\); standard normal referenceAppropriate only when \(\sigma\) is genuinely known for the inferential calculation.
Population standard deviation estimated; normal model\(SE=s/\sqrt{n}\); t reference with \(n-1\) degrees of freedomExact for independent normal observations; finite-sample uncertainty is reflected in heavier tails.
Large sample using a normal approximationOften close to the t calculationApproximation quality depends on the data and design; large \(n\) does not correct dependence or a badly unsuitable model.

In group sequential settings, estimating a variance can also affect the joint behavior of statistics across analyses. The fixed-sample t result above does not, by itself, justify using ordinary fixed-sample t critical values at repeated interim looks. Sequential error control requires a procedure calibrated for the planned repeated analyses; the point here is the distinction between known-variance information and precision estimated from data.

Diagnostics, Limitations, and a Common Mistake

Before trusting a t-based calculation, check whether its model is plausible. Look for unusual observations, marked skewness in a small sample, inconsistent measurement procedures, or dependence among observations. In a small sample, a single extreme value can substantially change \(s\), hence both the estimated information \(n/s^2\) and the t statistic. A zero sample standard deviation makes the usual t statistic undefined and is a warning to inspect the data and measurement process, not evidence of infinite precision.

  • Do not treat a planning standard deviation as known. It is an assumption used to plan; the realized standard deviation can differ.
  • Do not use a normal critical value automatically when \(s\) replaces \(\sigma\). For small samples under the normal model, use the t distribution with \(n-1\) degrees of freedom.
  • Do not interpret a larger observed \(s\) as a change in enrollment. At fixed \(n\), a larger \(s\) lowers estimated information and increases the estimated standard error.
  • Do not assume the t adjustment solves every design problem. It addresses variance estimation in the stated one-sample normal model, not unequal variances, clustering, missing outcomes, or repeated-look calibration.

The practical value of the distinction is that planning can be transparent about its assumed variability, while inference can reflect the uncertainty introduced by estimating that variability. Reporting \(n\), \(s\), the estimated standard error, and the inferential reference distribution makes the precision calculation easier to evaluate.

Check Your Understanding

Try each question, then compare your reasoning with the solutions.

  1. A hypothetical normal sample has \(n=25\) and \(s=5\) units. Calculate its estimated standard error and plug-in information.
  2. Why does a one-sample statistic using \(s/\sqrt{n}\) have a t distribution rather than a standard normal distribution under the normal model?
  3. For a 95% interval with \(n=16\), \(\bar{x}=20\) units, and \(s=4\) units, use \(t_{0.975,15}\approx2.131\) to calculate the interval.
  4. A planning calculation assumes \(\sigma_{\text{plan}}=10\) units, but the eventual sample has \(s=13\) units at the same \(n\). What happens to the estimated standard error and plug-in information relative to the planning values?

Solutions. (1) \(SE=5/\sqrt{25}=1\) unit, and \(\hat{I}=25/5^2=1\ \text{unit}^{-2}\). (2) The sample standard deviation is random. Under the normal model, the standardized mean using the true \(\sigma\) is independent of the sample variance, and dividing by the square root of the scaled sample variance gives a \(t_{n-1}\) statistic. (3) The half-width is \(2.131(4/\sqrt{16})=2.131\) units, so the interval is \((17.869,\ 22.131)\) units. (4) At fixed \(n\), the estimated standard error is \(13/\sqrt{n}\) rather than \(10/\sqrt{n}\), so it is larger; plug-in information is \(n/13^2\), which is smaller than the planned \(n/10^2\).

Clinical Biostats