Why a Reduction in Residual Error Needs a Benchmark
Adding predictors to a regression model cannot increase its residual sum of squares when the models are fitted to the same outcome and observations and the smaller model is nested within the larger one. The larger model has more freedom to fit the data, so some reduction in residual error can occur even when the added predictors do not capture a meaningful signal in the population. The question is whether the improvement is large relative to the residual variation that remains.
A partial F-test addresses that question for a specified set of added predictors. Its logic has two parts: it averages the improvement in fit across the restrictions being tested, then compares that average with the full model’s residual variation per residual degree of freedom. This makes the test sensitive not just to how much the residual sum of squares falls, but also to how many added coefficient restrictions are being assessed.
The numerator is the reduction in residual sum of squares per restriction. The denominator is the full model’s residual mean square, an estimate of the error variance under the full model. The statistic therefore compares an average improvement in fit with the scale of the residual noise left after fitting the full model. A larger \(F\) means the average improvement is large relative to that residual scale.
What Is Being Tested?
Suppose the full model contains all the terms in the reduced model plus a specified group of additional terms. The null hypothesis says that the added terms satisfy the restrictions defining the reduced model. For example, if the full model adds two predictors with coefficients \(\beta_5\) and \(\beta_6\), the null hypothesis may be \(H_0:\beta_5=0\) and \(\beta_6=0\). The alternative is that at least one of those restrictions does not hold.
The number \(q\) counts independent restrictions, not necessarily the number of words in the hypothesis or the number of variables in a data table. In the example above, \(q=2\). A categorical predictor represented by several indicator coefficients may contribute multiple restrictions, as can a group of polynomial terms. The comparison must be genuinely nested: every mean pattern allowed by the reduced model must also be allowed by the full model.
The division by \(q\) matters. A reduction of 20 units from adding one coefficient and a reduction of 20 units from adding five coefficients do not represent the same average improvement per restriction. The division by \(df_F\) matters as well: a residual sum of squares has no direct scale-free interpretation without accounting for how many residual degrees of freedom remain.
Conditions for the Usual F Reference Distribution
Before interpreting a p-value from the usual partial F-test, state the assumptions that support its reference distribution. The reduced and full models must use the same outcome and the same observations, and the reduced model must be nested in the full model. The restrictions must be independent, so that \(q\) correctly counts their number. The full model’s estimated coefficients must be identifiable; in particular, its design matrix must have full column rank.
For the exact classical F reference distribution, the regression model is assumed to have a correct conditional mean, with independent normally distributed errors of mean zero and constant variance \(\sigma^2\), conditional on the predictor values. Equivalently, the error terms have the same variance and are normal, with independence across observations. In applications, the sampling or study design and residual diagnostics can inform whether these conditions are plausible. If the assumptions are seriously questionable, the usual p-value may not be reliable.
Under the null hypothesis, the reduced model’s restrictions are true. The reduction in fit then reflects random variation rather than a real contribution from the added terms. In that setting, both mean squares in the ratio estimate the same error-variance scale, which explains why their ratio is a natural benchmark. Under the alternative, the added terms can account for systematic variation, increasing the numerator relative to the full model’s residual variation. This is the reasoning behind the test; the exact distributional result requires the stated assumptions.
Worked Examples
Worked Example: Why the Number of Restrictions Matters
Suppose two nested models are fitted to the same hypothetical outcome and 32 observations. The reduced model has \(RSS_R=240\). The full model estimates eight coefficients, including the intercept, and has \(RSS_F=210\). It adds three independent coefficients, so \(q=3\), and its residual degrees of freedom are \(32-8=24\). Assume the nested-model comparison and regression assumptions described above are appropriate.
The residual sum of squares falls by $$ RSS_R-RSS_F=240-210=30. $$
That improvement is spread across three restrictions, giving a numerator mean square of $$ \frac{RSS_R-RSS_F}{q}=\frac{30}{3}=10. $$
The full model’s residual mean square is $$ \frac{RSS_F}{df_F}=\frac{210}{24}=8.75. $$
Thus $$ F=\frac{10}{8.75}=1.143. $$
The statistic says that the average reduction in residual sum of squares per restriction is about 1.143 times the full model’s residual mean square. The reduction of 30 alone would not tell us that: dividing by the three restrictions gives the relevant average improvement. Whether this \(F\) value is unusually large is judged using an \(F\) distribution with \(3\) and \(24\) degrees of freedom.
Worked Example: Comparing Fit Improvement With Residual Variation
In a hypothetical analysis of 36 independent observations, a full model adds two predictors to a reduced model. Suppose the full model has six estimated coefficients, including the intercept, so \(df_F=36-6=30\). The reduced and full residual sums of squares are \(190\) and \(150\), respectively. Assume the added predictors correspond to two independent restrictions, that the full design matrix has full rank, and that the error assumptions are plausible for this illustration.
The reduction is $$ 190-150=40, $$ so the reduction per restriction is \(40/2=20\). The full model’s residual mean square is \(150/30=5\). Therefore, $$ F=\frac{(190-150)/2}{150/30}=\frac{20}{5}=4.000. $$
The test uses an \(F\) distribution with \(2\) numerator and \(30\) denominator degrees of freedom under \(H_0\). For these degrees of freedom, the upper-tail p-value is approximately \(0.02884\). At \(\alpha=0.05\), the p-value is smaller than the significance level, so we reject the null hypothesis that both added coefficients are zero. In context, the data provide evidence that the two added predictors, considered together, improve the model’s fit beyond what would be expected from residual variation under that null. This conclusion does not establish that either predictor causes changes in the outcome.
Worked Example: A Small Statistic Is Not Proof of No Improvement
Consider hypothetical nested models fitted to the same 32 observations. The full model estimates five coefficients and has \(RSS_F=108\), giving \(df_F=32-5=27\). It adds one coefficient to a reduced model, so \(q=1\). Suppose \(RSS_R=122.4\). Assume the model is identifiable and that the study design and residual checks make the exact-test assumptions reasonable.
The reduction in residual sum of squares is $$ 122.4-108=14.4. $$
Because there is one restriction, the reduction per restriction is \(14.4/1=14.4\). The full model’s residual mean square is $$ \frac{108}{27}=4. $$
Thus $$ F=\frac{14.4}{4}=3.600. $$
With \(1\) and \(27\) degrees of freedom, the upper-tail p-value is approximately \(0.06852\). Since \(0.06852>0.05\), we fail to reject the null hypothesis at the 5% significance level. The result does not prove that the added coefficient is zero or that the predictor has no value. It says that this comparison does not provide sufficiently strong evidence against the specified null at the chosen significance level.
A Short Justification of the Mean-Square Comparison
The form of the statistic follows from asking what quantities should be compared. The raw reduction \(RSS_R-RSS_F\) tends to be larger when more restrictions are released, so divide it by \(q\) to express improvement per restriction. The residual sum of squares in the full model also tends to scale with the number of residual degrees of freedom, so divide it by \(df_F\) to obtain the full-model residual mean square. Taking their ratio gives $$ \frac{(RSS_R-RSS_F)/q}{RSS_F/df_F}. $$
Under \(H_0\), neither mean square is expected to reflect a systematic contribution from the added terms; both are on the error-variance scale. If the added terms explain real systematic variation, the numerator tends to be larger relative to the denominator. This reasoning motivates the statistic and its upper-tail test. It does not mean the ratio itself is the probability that \(H_0\) is true; the p-value is calculated from the null distribution of the statistic.
Common Pitfalls and Why the Test Matters
- Using only the raw reduction. A decrease in residual sum of squares is expected when predictors are added. Compare the decrease per restriction with the full model’s residual mean square instead.
- Forgetting \(q\). A block of several added coefficients is not treated as one restriction unless it truly represents one independent restriction. Use the number of independent restrictions in the hypothesis.
- Comparing models fitted to different observations. The residual-sum-of-squares reduction has the intended interpretation only when both models use the same response and sample.
- Calling a nonsignificant result proof of no effect. Failure to reject means the observed statistic is not sufficiently unusual under the null at the selected \(\alpha\); it does not establish that the null is true.
- Reporting a p-value without its context. State which added terms were tested, the degrees of freedom, and the conclusion in the outcome’s context. Statistical evidence of improved fit is not, by itself, evidence of practical importance or causation.
Check Your Understanding
Answer each question using the partial F-test logic and the definitions in this tutorial.
- A full model reduces residual sum of squares by 24 when six independent restrictions are released. What is the reduction per restriction?
- Why is the denominator of the partial F-statistic based on the full model’s residual sum of squares rather than the reduced model’s?
- Two models have \(RSS_R=150\), \(RSS_F=120\), \(q=3\), and \(df_F=24\). Calculate the partial F-statistic.
- What does a p-value of 0.08 imply about a test conducted at \(\alpha=0.05\), and what does it not imply?
- Name two conditions needed for a valid usual partial F-test, apart from the assumptions on the errors.