Power of Test Statistic
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:
- Pick Ξ± (e.g. 0.05) β gives you $c$.
- Pick a minimum effect size you care about detecting (this is ΞΈβ β ΞΈβ, or Ξ΄) β this is a judgment call, not derived from data.
- Plug into the power formula above (or software like
power.t.testin R, orstatsmodels.stats.powerin Python) to get a number like “80% power.” - 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.