Step 10: Model Specification and Multicollinearity
- A specification error is a wrong choice of variables, functional form, measurement or sample. Some cause bias (omitted variables, bad controls, wrong form, error in \( x \), selection on \( y \)); others only cost precision.
- Multicollinearity (high but imperfect correlation among regressors) causes no bias. It inflates variances by \( \text{VIF}_j = 1/(1-R_j^2) \). Never “fix” it by dropping a relevant variable.
- A bad control is an outcome of the variable you care about; holding it fixed removes part of the effect you want.
- RESET adds \( \hat y^2 \) and \( \hat y^3 \) to test functional form. Classical measurement error in \( x \) biases the slope towards zero by the reliability ratio.
- Exam point: find influential points with leverage, studentised residuals and Cook’s distance, but correct only genuine errors. Never delete data because it is inconvenient.
- The main types of specification error and which ones bias your estimates.
- How to recognise multicollinearity, compute VIFs, and what to do (and not do) about it.
- How to choose controls, spot bad controls and judge a proxy variable.
- How to run Ramsey’s RESET, the \( J \) test and the AIC and BIC.
- What measurement error, missing data and selected samples do to OLS, and how to handle influential observations.
What is a specification error?
Every regression rests on choices: which variables, in what form, measured how, on which observations. Each can be wrong, and none shows up as an error message. The software reports neat coefficients and small standard errors whatever you feed it. This step is about the checks that stand between a regression and a believable result.
| Specification error | Consequence for OLS | How to detect or avoid it |
|---|---|---|
| Omitted relevant variable | Bias if correlated with an included \( x \) | Theory; sign the bias; proxies |
| Irrelevant variable included | No bias; larger variance | \( t \) and \( F \) tests; theory |
| Bad control (an outcome of \( x \)) | Bias: removes part of the effect | Ask what causes what |
| High multicollinearity | No bias; large variances | VIFs; joint tests |
| Wrong functional form | Bias in the partial effects | RESET; plots; \( J \) test |
| Measurement error in \( x \) | Bias towards zero | Better data; instruments |
| Measurement error in \( y \) | No bias; larger variance | Better data |
| Sample selected on \( y \) | Bias | Know how the sample was drawn |
| Influential observations | Results driven by a few points | Leverage; Cook’s distance |
Multicollinearity
Step 6 met perfect collinearity, an exact linear relationship such as exper = age − educ − 6, which breaks MLR.3. Multicollinearity is different: one regressor is highly, but not perfectly, correlated with a combination of the others. It breaks no assumption. OLS is still unbiased and BLUE. It makes the affected coefficients imprecise, because little independent variation is left to estimate each separate effect.
The standard error is inflated by \( \sqrt{\text{VIF}_j} \) compared with a regressor uncorrelated with the others. Inflation grows slowly: to multiply a standard error by 10 needs VIF = 100, that is \( R_j^2 = 0.99 \).
Chapter 10 of Econometrics Made Simple regresses house price on sqft, lot, bedrooms, bathrooms, age, dist_cbd and pool (500 sales). The other regressors explain \( R_j^2 = 0.771 \) of the variation in sqft. Compute its VIF, tolerance and standard-error inflation.
The four size measures overlap heavily. Age, distance and pool, with VIFs of about 1.01 to 1.02, are almost unrelated to size, so their coefficients are as precise as the sample allows.
The symptoms
Multicollinearity has three classic symptoms, all from one cause: the data cannot tell the correlated variables apart. (1) Large standard errors on the correlated variables, often with odd signs. (2) Joint but not individual significance. (3) Unstable coefficients across samples. It is above all a small-sample problem, because \( \text{SST}_j \) grows with \( n \).
The book shows this with the first 30 sales only. Floor area has \( p = 0.072 \), bedrooms \( p = 0.920 \) (with a standard error ten times its coefficient) and bathrooms \( p = 0.118 \). A careless reader would say size does not matter. Yet the \( F \) test that all four size variables are zero gives \( F = 20.21 \), \( p = 2.9 \times 10^{-7} \). Across neighbourhoods, the sqft coefficient stays between 0.107 and 0.135, while the bedrooms coefficient swings from −5.809 to 13.616.
What to do, and what not to do
First ask whether there is a problem at all. If the coefficient you care about is precise enough, collinearity among the other regressors does not matter. When it does matter: get more data; test the collinear group jointly (Step 7); combine the variables into an index; ask a question the data can answer, such as the value of a bundle of extra space and an extra bedroom; or centre polynomial terms to reduce the collinearity between \( x \) and \( x^2 \).
In the full house regression the bedrooms coefficient is 4.680 (se 5.422), not significant. A student drops sqft “to fix the multicollinearity”. The bedrooms coefficient becomes 40.877 (se 5.194), with \( t \) near 8. What went wrong?
This is the precision problem swapped for a bias problem. For the difference between the two kinds of collinearity, see Difference between Perfect and Imperfect Multicollinearity.
Which controls? Bad controls and proxies
Omitting a relevant variable biases the coefficient on \( x_1 \) by \( \beta_2\delta_1 \). Including an irrelevant one costs only precision: adding parents’ education to the wage equation raised the standard error of the education coefficient from 0.00786 to 0.00847. So when a variable is plausibly relevant and correlated with the regressor of interest, it usually belongs in. There is one important exception.
Controlling for occupation when estimating the return to education compares two accountants with different schooling and leaves out the main channel: better jobs. The same applies to the gender pay gap. The book simulates 20,000 workers in which being female lowers the log wage by 0.10 directly and a high-paying job raises it by 0.40. Half of men (50.6%) but a quarter of women (24.9%) hold such jobs. Without the job variable the gap is −0.2043: about −0.10 directly plus \( 0.40 \times (0.249 – 0.506) = -0.103 \) through jobs. Adding the job dummy gives −0.1032, the direct effect only. Both are correct; they answer different questions.
Proxy variables
A proxy is an observed variable included in place of an unobserved one, such as an IQ score for ability. It works only if, once the proxy is held fixed, the rest of the omitted factor is unrelated to the regressor. In practice a proxy usually reduces the bias rather than removing it. In the book’s simulation with a true return of 0.08, omitting ability gives 0.1065, adding a noisy test score gives 0.0857, and adding true ability gives 0.0773. Parents’ education, by contrast, moves the coefficient the wrong way (0.1048 to 0.1078), because in these data it shifts schooling but is unrelated to ability. Judge a proxy by whether it plausibly captures the omitted factor, not by whether it is available. Step 13 offers a stronger tool, instrumental variables.
Testing functional form
If the true relationship is curved and you fit a straight line, the omitted curvature sits in the error and is correlated with \( x \). When you have a specific alternative, test it directly (add \( x^2 \) and check its \( t \)). Ramsey’s RESET is a general check: if the model is right, no nonlinear function of the regressors should help.
The book runs RESET on wage and house models. Interpret: wage in levels \( F = 11.639 \) (\( p = 1.01 \times 10^{-5} \)); log wage \( F = 0.320 \) (\( p = 0.726 \)); log wage without exper² \( F = 0.783 \) (\( p = 0.458 \)).
RESET is a useful alarm with limited power against specific omissions. It never tells you what is wrong, and a pass is not proof the model is right.
For non-nested models, such as log price on sqft versus log price on log(sqft), use the Davidson–MacKinnon \( J \) test: add one model’s fitted values to the other and test them. Adding the log-size model’s fitted values to the level-size model gives \( t = 7.21 \) (reject the level model); the reverse gives \( t = -0.49 \) (keep the log model). Information criteria agree. For OLS, \( \text{AIC} = n\ln(\text{SSR}/n) + 2(k+1) \) and \( \text{BIC} = n\ln(\text{SSR}/n) + (k+1)\ln n \), up to a constant; smaller is better, and both must use the same dependent variable and observations. The log-size model has the lower AIC (−274.27 against −224.44).
Measurement error and selected samples
Error in \( y \) that is unrelated to the regressors only adds noise: OLS stays unbiased, with larger standard errors. Error in \( x \) is worse. If we observe \( x = x^* + e \) with classical error \( e \),
so the slope is biased towards zero (attenuation).
The variance of education in wages.csv is 3.4796. The book adds random reporting error with standard deviation 1 year. Predict the attenuated slope, starting from 0.1026.
The standard error does not warn you: it is smaller (0.0073 against 0.0081). Adding noise to log wage instead leaves the slope near 0.1026 (0.1090 in one draw) but raises the standard error to 0.0104.
Missing data matter according to why they are missing. Dropping 30% of workers at random leaves the education coefficient at 0.1060 (from 0.1048) with a larger standard error. Selecting on \( x \) (only workers with 12+ years of schooling) gives 0.1201, within sampling error. Selecting on \( y \) (only wages below $25) gives 0.0544, almost half: selection bias, with a small standard error (0.0067) that hides it.
Outliers and influential observations
OLS minimises squared residuals, so a few points can have a large say. Leverage \( h_i \) measures how far an observation’s \( x \) values are from the centre (flag \( h_i > 2(k+1)/n \)). The studentised residual scales the residual by a standard error estimated without that point (values beyond ±3 are rare). Cook’s distance combines both:
Influence needs both leverage and a large residual. In the book’s price regression on sqft, age, distance and pool, 41 sales have high leverage, 5 have studentised residuals beyond ±3, and 30 exceed \( 4/n = 0.008 \), but the largest Cook’s distance is only 0.082 (house 363). The biggest house has the highest leverage, 0.067, but sits near the line (Cook’s \( D = 0.005 \)). Dropping the five most influential sales moves the pool premium from 27.253 to 22.949 thousand dollars, a little over half a standard error.
Now change one digit: type house 1’s floor area as 27,050 instead of 2,705. The sqft coefficient falls from 0.1388 to 0.0352, \( R^2 \) from 0.676 to 0.286, and that sale’s Cook’s distance is 174.5. Diagnostics earn their keep when the data contain a mistake.
A specification checklist
Practise: 6 questions
Try each one on paper first, then check the worked answer.
EasyQ1. Regressing \( x_j \) on the other regressors gives \( R_j^2 = 0.90 \). Find the VIF, the tolerance, and the factor by which the standard error of \( \hat\beta_j \) is inflated. Show worked answer
Tolerance \( = 0.10 \).
Standard-error inflation: \( \sqrt{10} = 3.16 \).
Whether this matters depends on whether \( \hat\beta_j \) is still precise enough to answer the question.
EasyQ2. A study estimates the effect of smoking during pregnancy on birthweight. For each candidate control, say whether to include it: (a) the mother’s age; (b) the length of the pregnancy (gestation); (c) the mother’s years of education. Show worked answer
(b) Exclude: smoking shortens pregnancies, which lowers birthweight. Gestation is an outcome of smoking, a bad control; holding it fixed removes part of the effect.
(c) Include: determined before the pregnancy and likely related to both.
MediumQ3. The true slope is 0.5. The true regressor has variance 4 and is measured with classical error of variance 1. (a) What does OLS estimate in large samples? (b) What if the error variance doubles? Show worked answer
(b) \( 4/(4+2) = 0.667 \), so the estimate is \( 0.5 \times 0.667 = 0.333 \).
Noisier measurement means more attenuation towards zero.
MediumQ4. A model with \( k = 4 \) and \( n = 250 \) has \( R^2 = 0.360 \). Adding \( \hat y^2 \) and \( \hat y^3 \) raises it to 0.372. Compute the RESET \( F \) statistic and test at 5% (critical value of \( F_{2,243} \): 3.03). Show worked answer
\( F = \dfrac{(0.372-0.360)/2}{(1-0.372)/243} = \dfrac{0.006}{0.002584} = 2.32 \).
\( 2.32 < 3.03 \): no evidence of misspecification at 5%. That is not proof the form is right; RESET has limited power.
MediumQ5. In a regression with \( k = 3 \) and \( n = 100 \), observation 57 has leverage 0.10 and standardised residual 2.5. Compute its Cook’s distance and say what you would do. Show worked answer
This is above \( 4/n = 0.04 \), so worth a look, but well below 1.
Check the record for errors. If it is genuine, keep it, ask whether the model is missing something, and report results with and without it.
ExamQ6. Two models for the same \( y \) are fitted to \( n = 400 \) observations. Model A has \( k = 3 \) and SSR = 210; model B adds two regressors (\( k = 5 \)) and has SSR = 204. Using \( \text{AIC} = n\ln(\text{SSR}/n) + 2(k+1) \) and \( \text{BIC} = n\ln(\text{SSR}/n) + (k+1)\ln n \): (a) compute both AICs [2]; (b) compute both BICs [2]; (c) test the two extra regressors with an \( F \) test (5% critical value of \( F_{2,394} \): 3.02) [2]; (d) which model would you choose? [1] Show worked answer
(b) \( \ln 400 = 5.991 \). A: \( -257.74 + 4(5.991) = -233.78 \). B: \( -269.34 + 6(5.991) = -233.39 \). [M1 A1] A has the (slightly) lower BIC, because BIC penalises parameters more heavily.
(c) \( F = \dfrac{(210-204)/2}{204/394} = \dfrac{3}{0.5178} = 5.79 > 3.02 \): the two regressors are jointly significant. [M1 A1]
(d) B: the \( F \) test and AIC favour it, and the BIC difference is tiny (0.39), which is not decisive. [B1]
Total: 7 marks.
Check yourself
1. What does high (but imperfect) multicollinearity do to OLS?
2. Classical measurement error in an explanatory variable
3. Which sample restriction biases OLS?
Your score: 0 out of 3
The short version
- Omitted variables, bad controls, wrong form, error in \( x \) and selection on \( y \) cause bias; irrelevant variables, multicollinearity and error in \( y \) cost only precision.
- \( \text{VIF}_j = 1/(1-R_j^2) \); sqft has VIF 4.367 in the house regression. Symptoms are large standard errors and joint but not individual significance. Dropping sqft makes bedrooms jump from 4.680 to 40.877: bias, not a cure.
- Bad controls are outcomes of the variable of interest; controlling for job type halves the simulated gender gap from −0.204 to −0.103. A proxy helps only if it captures the omitted factor.
- RESET rejects the level-wage model (\( F = 11.64 \)) but not the log-wage model (\( F = 0.32 \)). Compare non-nested models with the \( J \) test and, for the same \( y \), AIC and BIC.
- Classical error in \( x \) attenuates the slope by the reliability ratio (0.777 in the wage simulation). Selection on \( y \) causes bias.
- Use leverage, studentised residuals and Cook’s distance; correct genuine errors, keep genuine unusual data, and report results both ways.
- Chapter 10 of Econometrics Made Simple covers all of this in full, with code and a full set of exercises with answers.
Studying more than one subject? The Complete Book Bundle includes every Economics, Statistics, Mathematics and Econometrics book we publish.
Further reading on EconTutorials
Related in-depth articles from our blog, and the neighbouring steps:
- Difference between Perfect and Imperfect Multicollinearity
- Omitted Variable Bias: How to Predict the Direction Before You Run the Regression
- Instrumental Variables and Two-Stage Least Squares: Solving the Endogeneity Problem
- Heteroskedasticity: What It Is, Why It Matters, and How to Fix It
- Previous step: Step 9: Dummy Variables and the Chow Test
- Next step: Step 11: Heteroskedasticity
