Why Test Predictors Together?
When predictors are correlated, a regression model can estimate their individual coefficients imprecisely. Each coefficient describes the association with the outcome when the other predictors are held fixed; if the predictors tend to vary together, the data may contain limited information about changing one while holding the others constant. Yet the predictors may still have a detectable combined association with the outcome. A joint test addresses that different question by testing whether a set of coefficients is zero all at once.
A joint result and a collection of individual results are therefore not competing answers to the same question. A joint test can reject the hypothesis that all coefficients in a group are zero even though none of the individual coefficient tests rejects its own null hypothesis. The joint test does not, by itself, identify which coefficient is nonzero or establish a causal effect.
The Model and the Joint Hypothesis
Consider a linear regression with an outcome \(Y_i\), two predictors of interest \(X_{i1}\) and \(X_{i2}\), and any other predictors that are kept in both models:
The coefficients \(\beta_1\) and \(\beta_2\) describe the conditional associations of their respective predictors with the mean outcome, holding the other included predictors fixed. To test the two predictors jointly, the null hypothesis is \(H_0:\beta_1=\beta_2=0\). The alternative is that at least one of those two coefficients is not zero. This is a test of two restrictions, so \(q=2\). More generally, a joint test can impose \(q\) coefficient restrictions.
A partial F-test compares two nested models fitted to the same observations. The full model contains the predictors whose coefficients are being tested, along with the predictors to be retained. The reduced model omits the tested predictors but retains the other predictors. Let \(RSS_F\) and \(RSS_R\) be the residual sums of squares for the full and reduced models, respectively. Let \(df\) be the full model’s residual degrees of freedom. The statistic is:
The numerator measures the reduction in residual sum of squares per restriction when moving from the reduced model to the full model. The denominator is the full model’s residual mean square. Under the null hypothesis and the assumptions below, the statistic follows an \(F\) distribution with \(q\) numerator degrees of freedom and \(df\) denominator degrees of freedom.
Assumptions for the Partial F-Test
The test uses the usual linear regression model assumptions for exact small-sample inference: the conditional mean is correctly described by the specified linear model; errors are independent, normally distributed with mean zero and constant variance; and the predictors have no exact linear dependence in the full model. Equivalently, conditional on the observed predictors, the errors satisfy \(\varepsilon_i\sim N(0,\sigma^2)\) independently. The full and reduced models must be nested and must use the same observations. Correlation among predictors is permitted: it is near or exact redundancy—not correlation by itself—that affects what can be estimated.
Normality and constant variance matter for the exact \(F\) reference distribution stated here. Independence and a suitable mean model matter for the regression inference more generally. Before interpreting a test, consider whether the data collection and model make these assumptions reasonable; a calculation does not repair a poorly specified model.
Why the Statistic Has an F Distribution
The comparison works because adding predictors to a nested model cannot increase its residual sum of squares, so \(RSS_R-RSS_F\) is nonnegative. Under \(H_0\), the added predictors do not improve the conditional mean beyond random variation. With normally distributed errors and a full-rank design, the scaled reduction in residual sum of squares, \((RSS_R-RSS_F)/\sigma^2\), has a chi-square distribution with \(q\) degrees of freedom. The scaled residual sum of squares from the full model, \(RSS_F/\sigma^2\), has a chi-square distribution with \(df\) degrees of freedom. These quantities are independent for nested least-squares models under the stated assumptions. Dividing each chi-square quantity by its degrees of freedom gives the partial F-statistic:
The p-value is the probability, assuming \(H_0\) and the model assumptions, of obtaining an F-statistic at least as large as the observed one. It is not the probability that \(H_0\) is true. A small p-value is evidence against the joint null; it does not tell us which coefficient is responsible.
Worked Example: A Joint Test That Is Not Significant
Suppose a hypothetical analysis uses 30 observations and compares a full model with two predictors against a reduced model containing only an intercept. Assume the linear-model conditions are reasonable, the full model has no exact collinearity, and the predictors are not strongly correlated. The full model has residual sum of squares \(RSS_F=108\) and \(df=27\); the reduced model has \(RSS_R=117\). Test \(H_0:\beta_1=\beta_2=0\) at \(\alpha=0.05\).
There are \(q=2\) restrictions. The full model’s residual mean square is \(108/27=4\), and the reduction in residual sum of squares is \(117-108=9\). Therefore:
Under the null, this is compared with an \(F(2,27)\) distribution. For two numerator degrees of freedom, the upper-tail probability at \(F=1.125\) is \((27/(27+2(1.125)))^{27/2}=(12/13)^{13.5}\approx0.3394\). Since \(0.3394>0.05\), the test does not reject the joint null. These data do not provide sufficient evidence at the 0.05 level that the two predictors, considered together and conditional on the specified model, improve the mean model. This is not proof that both population coefficients are zero.
Joint Evidence Can Coexist With Imprecise Individual Tests
An individual coefficient test usually uses a \(t\)-statistic: the estimated coefficient divided by its standard error. When predictors are correlated, the coefficient estimates can have substantial standard errors and can be negatively correlated with each other. A joint test uses the estimated coefficients and their covariance together. If the estimates move in a coordinated direction, the joint evidence can be stronger than either separate \(t\)-test suggests.
The following example gives values that illustrate this pattern. The coefficient covariance matters: a pair of individual standard errors alone is not enough to calculate a joint test. In practice, the partial F-test from the nested-model comparison incorporates that covariance information.
Worked Example: Significant Joint Test, Inconclusive Individual Tests
In a hypothetical analysis of 50 observations, the full model includes two highly correlated predictors of interest and one adjustment predictor, plus an intercept. The reduced model retains the adjustment predictor but omits the two predictors being tested. Suppose the full model has \(RSS_F=180\) and \(df=46\), while the reduced model has \(RSS_R=240\). Assume both models use the same observations, the full design has no exact collinearity, and the error and mean-model assumptions are reasonable.
The joint null is \(H_0:\beta_1=\beta_2=0\), so \(q=2\). The numerator mean square is \((240-180)/2=30\). Keeping the unrounded denominator in the calculation gives:
The reference distribution is \(F(2,46)\). Its upper-tail p-value is approximately \(0.00134\), so the joint null is rejected at \(\alpha=0.05\). The full model provides evidence that at least one of the two tested coefficients is nonzero, conditional on the adjustment predictor and the model assumptions.
Suppose the coefficient estimates are both \(1.20\) outcome units per predictor unit, each with standard error \(0.80\). Each individual statistic is \(t=1.20/0.80=1.50\), with 46 degrees of freedom; each two-sided p-value is about \(0.140\). Neither individual test rejects at 0.05. This does not conflict with the joint result. For instance, if the estimated correlation between the two coefficient estimates is about \(-0.7065\), the joint statistic based on their covariance can be large even though each separate ratio is only 1.50. The individual tests ask whether each coefficient, by itself, is distinguishable from zero; the joint test asks whether the two coefficients are simultaneously zero.
The proper conclusion is not that both predictors are individually established as important. The joint test identifies evidence against the hypothesis that both coefficients are zero, but it does not identify which coefficient differs from zero or show that either predictor causes the outcome to change.
The Comparison Depends on What Is Retained
A partial F-test is conditional on the predictors included in both models. For example, testing two predictors while retaining an adjustment variable asks whether the two add explanatory information beyond that adjustment variable. Omitting the adjustment variable from both models would ask a different question and could produce a different result. The reduced model must represent the particular null hypothesis being tested; it is not simply any simpler model.
Worked Example: Testing Two Added Predictors Beyond a Baseline Model
A hypothetical data set has 52 observations. A baseline model contains an intercept and one predictor; a full model adds two predictors of interest. The baseline model’s residual sum of squares is \(365\), and the full model’s residual sum of squares is \(275\). Suppose the regression assumptions are reasonable, the full model has no exact collinearity, and the models use the same observations. The full model has four fitted coefficients including the intercept, so its residual degrees of freedom are \(52-4=48\).
The null hypothesis is that both added coefficients are zero, so \(q=2\). The full-model residual mean square is \(275/48\approx5.729\). The partial F-statistic is:
For \(F(2,48)\), the upper-tail probability is \((48/(48+2F))^{24}\), which here is \((48/63.71)^{24}\approx0.00055\). The p-value is below 0.05, so the data provide evidence that the two added coefficients are not both zero, conditional on the baseline predictor. This conclusion does not say that both added predictors have individually significant coefficients.
Common Mistake: Treating Separate Tests as the Joint Answer
- Wrong reasoning: “Neither individual coefficient has a p-value below 0.05, so the two predictors have no association with the outcome.”
- Correction: The individual tests address separate conditional coefficients. Test the joint null directly with the partial F-test if the question concerns whether the group of coefficients is zero.
- Wrong reasoning: “The joint test is significant, so both predictors are important.”
- Correction: Rejecting \(H_0:\beta_1=\beta_2=0\) means at least one coefficient is nonzero under the model; it does not establish that each is.
- Wrong reasoning: “The joint p-value gives the probability that the null hypothesis is true.”
- Correction: A p-value is calculated under the null and describes the tail probability of a statistic at least as extreme as the observed one. It is not a probability assigned to the null itself.
These distinctions matter when reporting results. State which coefficients were tested together, which predictors remained in both models, the test statistic and degrees of freedom, and the p-value. Then describe the evidence in context without turning an association into a causal claim. Correlated predictors may make individual coefficient estimates unstable or imprecise, even when a group-level test indicates that the predictors add information collectively.
Check Your Understanding
Answer each question using the distinction between a joint hypothesis and individual coefficient hypotheses.
- For \(H_0:\beta_2=\beta_3=0\), how many restrictions are tested, and what does rejection establish?
- A reduced model has \(RSS_R=150\), and a full model has \(RSS_F=120\), \(df=40\), and two added predictors. Write the partial F-statistic before calculating it.
- Why can both individual coefficient tests be inconclusive while the joint test rejects?
- When testing two coefficients, why must the reduced model retain the other predictors from the full model?
- What does a small joint-test p-value fail to establish about the individual predictors?