Power of Test Statistic

Last modified: July 21, 2026

Parent: Probability and Statistics

Power is Ξ²(ΞΈ) evaluated at a specific ΞΈ in Ξ˜β‚

Recall the power function:

$$\beta(\theta) = P_\theta\big((X_1,\ldots,X_n) \in R\big)$$

This is defined for every ΞΈ, but it means different things depending on which region ΞΈ sits in:

  • When $\theta \in \Theta_0$: $\beta(\theta)$ = probability of a false positive at that ΞΈ. You want this small (≀ Ξ±).
  • When $\theta \in \Theta_1$: $\beta(\theta)$ = probability of correctly rejecting Hβ‚€ β€” this is the power at that ΞΈ.

So “power” isn’t a single number in general β€” it’s a function of ΞΈ, because how easy it is to detect an effect depends on how far that particular ΞΈ is from Hβ‚€. A coin biased at P(heads)=0.9 is way easier to detect than one biased at P(heads)=0.51. That’s why you’ll often see “power at ΞΈ = ___” rather than just “the power.”

Worked example: normal mean test

Setup: $X_1,\ldots,X_n \sim N(\theta, 1)$, testing $H_0: \theta = \theta_0$ vs $H_1: \theta \neq \theta_0$, rejecting when $|T_n| > c$ where $T_n = \sqrt{n}(\bar X_n - \theta_0)$.

Step 1 β€” find c from Ξ±. Under $H_0$, $T_n \sim N(0,1)$, so you pick $c = z_{\alpha/2}$ (the standard normal critical value, e.g. 1.96 for Ξ± = 0.05).

Step 2 β€” ask: what if the true mean is actually some $\theta_1 \neq \theta_0$? Under this true θ₁, $T_n$ is no longer standard normal β€” it’s shifted:

$$T_n = \sqrt{n}(\bar X_n - \theta_0) \sim N\big(\sqrt{n}(\theta_1 - \theta_0),\ 1\big)$$

Why? Because $\bar X_n \sim N(\theta_1, 1/n)$ when θ₁ is the true mean, so $\sqrt n(\bar X_n - \theta_0)$ has mean $\sqrt n (\theta_1-\theta_0)$ and variance 1. Call this shift $\delta = \sqrt{n}(\theta_1 - \theta_0)$ β€” it’s essentially “the true effect size, scaled up by sample size.”

Step 3 β€” compute the power directly.

$$\text{Power}(\theta_1) = P_{\theta_1}(|T_n| > c) = P(|Z + \delta| > c)$$

where $Z \sim N(0,1)$. Expanding this out:

$$\text{Power}(\theta_1) = 1 - \Phi(c - \delta) + \Phi(-c-\delta)$$

That messy-looking expression is just: “the probability the shifted normal curve still lands past your cutoff $c$ in either tail.”


The intuition, stripped of formulas

Picture two bell curves on the same axis:

  • Your rejection region is “past the cutoff $c$” (shaded tail).
  • Power = how much of the H₁ curve’s area falls in that shaded rejection zone.

This picture makes every dependency obvious:

  • Bigger effect size (θ₁ far from ΞΈβ‚€, so Ξ΄ is large) β†’ H₁ curve shifts far away from Hβ‚€ curve β†’ more of it clears the cutoff β†’ higher power.
  • Smaller Ξ± β†’ cutoff $c$ moves further out β†’ less area under H₁ curve clears it β†’ lower power (this is the Ξ±/power tradeoff).
  • Bigger sample size $n$ β†’ Ξ΄ = √n(ΞΈβ‚βˆ’ΞΈβ‚€) grows, and the curves get narrower (variance shrinks) β†’ separation between the two curves increases faster than the cutoff moves β†’ higher power without needing to loosen Ξ±. This is why more data breaks the tradeoff.

In practice (how it’s actually calculated)

For power analysis in real study design:

  1. Pick Ξ± (e.g. 0.05) β†’ gives you $c$.
  2. Pick a minimum effect size you care about detecting (this is θ₁ βˆ’ ΞΈβ‚€, or Ξ΄) β€” this is a judgment call, not derived from data.
  3. Plug into the power formula above (or software like power.t.test in R, or statsmodels.stats.power in Python) to get a number like “80% power.”
  4. Often flipped around: fix your desired power (commonly 80%) and solve for the sample size $n$ needed to hit it β€” this is the standard “sample size calculation” you see in papers.

One-line takeaway: power is the area of the H₁-true sampling distribution that falls inside your rejection region β€” driven by effect size, sample size, and how strict Ξ± forces the cutoff to be.