Heckman correction

The Heckman correction (the two-stage method, Heckman's lambda or the Heckit method^[1]) is any of a number of related statistical methods developed by James Heckman at the University of Chicago in 1976 to 1979 which allow the researcher to correct for selection bias.^[2] Selection bias problems are endemic to applied econometric problems,^[3]^[4] which make Heckman’s original technique, and subsequent refinements by both himself and others, indispensable to applied econometricians. Heckman received the Economics Nobel Prize in 2000 for this achievement.

Method

Statistical analyses based on non-randomly selected samples can lead to erroneous conclusions and poor policy. The Heckman correction, a two-step statistical approach, offers a means of correcting for non-randomly selected samples.

Heckman discussed bias from using nonrandom selected samples to estimate behavioral relationships as a specification error. He suggests a two-stage estimation method to correct the bias. The correction uses a control function idea and is easy to implement. Heckman’s correction involves a normality assumption, provides a test for sample selection bias and formula for bias corrected model.

Suppose that a researcher wants to estimate the determinants of wage offers, but has access to wage observations for only those who work. Since people who work are selected non-randomly from the population, estimating the determinants of wages from the subpopulation who work may introduce bias. The Heckman correction takes place in two stages.

In the first stage, the researcher formulates a model, based on economic theory, for the probability of working. The canonical specification for this relationship is a probit regression of the form

\operatorname {Prob}(D=1|Z)=\Phi (Z\gamma ),\,

where D indicates employment (D = 1 if the respondent is employed and D = 0 otherwise), Z is a vector of explanatory variables, $\gamma$ is a vector of unknown parameters, and Φ is the cumulative distribution function of the standard normal distribution. Estimation of the model yields results that can be used to predict this employment probability for each individual.

In the second stage, the researcher corrects for self-selection by incorporating a transformation of these predicted individual probabilities as an additional explanatory variable. The wage equation may be specified,

w^{*}=X\beta +u\,

where $w^{*}$ denotes an underlying wage offer, which is not observed if the respondent does not work. The conditional expectation of wages given the person works is then

E[w|X,D=1]=X\beta +E[u|X,D=1].\,

Under the assumption that the error terms are jointly normal, we have

E[w|X,D=1]=X\beta +\rho \sigma _{u}\lambda (Z\gamma ),\,

where ρ is the correlation between unobserved determinants of propensity to work $\varepsilon$ and unobserved determinants of wage offers u, σ_u is the standard deviation of $u$ , and $\lambda$ is the inverse Mills ratio evaluated at $Z\gamma$ . This equation demonstrates Heckman's insight that sample selection can be viewed as a form of omitted-variables bias, as conditional on both X and on $\lambda$ it is as if the sample is randomly selected. The wage equation can be estimated by replacing $\gamma$ with Probit estimates from the first stage, constructing the $\lambda$ term, and including it as an additional explanatory variable in linear regression estimation of the wage equation. Since $\sigma _{u}>0$ , the coefficient on $\lambda$ can only be zero if $\rho =0$ , so testing the null that the coefficient on $\lambda$ is zero is equivalent to testing for sample selectivity.

Heckman's achievements have generated a large number of empirical applications in economics as well as in other social sciences. The original method has subsequently been generalized, by Heckman and by others.^[5]

Disadvantages

The two-step estimator discussed above is a limited information maximum likelihood (LIML) estimator. In asymptotic theory and in finite samples as demonstrated by Monte Carlo simulations, the full information (FIML) estimator exhibits better statistical properties. However, the FIML estimator is more computationally difficult to implement.^[6]
The covariance matrix generated by OLS estimation of the second stage is inconsistent. Correct standard errors and other statistics can be generated from an asymptotic approximation or by resampling, such as through a bootstrap.
The canonical model assumes the errors are jointly normal. If that assumption fails, the estimator is generally inconsistent and can provide misleading inference in small samples.^[7] Semiparametric and other robust alternatives can be used in such cases.^[8]
The model obtains formal identification from the normality assumption when the same covariates appear in the selection equation and the equation of interest, but identification will be tenuous unless there are many observations in the tails where there is substantial nonlinearity in the inverse Mills ratio. Generally, an exclusion restriction is required to generate credible estimates: there must be at least one variable which appears with a non-zero coefficient in the selection equation but does not appear in the equation of interest, essentially an instrument. If no such variable is available, it may be difficult to correct for sampling selectivity.^[6]

Implementations in statistics packages

R: Heckman-type procedures are available as part of the sampleSelection package.^[9]^[10]
Stata: the command heckman provides the Heckman selection model.^[11]^[12]

References

↑ Heckit: 'Heck-' from Heckman and '-it' as in probit, tobit, and logit.
↑ Heckman, J. (1979). "Sample selection bias as a specification error". Econometrica. 47 (1): 153–61. doi:10.2307/1912352. JSTOR 1912352.
↑ Gronau, Reuben (1974). "Wage Comparisons—A Selectivity Bias". Journal of Political Economy. 82 (6): 1119–1143. doi:10.1086/260267. JSTOR 1830664.
↑ Lewis, H. Gregg (1974). "Comments on Selectivity Biases in Wage Comparisons". Journal of Political Economy. 82 (6): 1145–1155. doi:10.1086/260268. JSTOR 1830665.
↑ Lee, Lung-Fei (2001). "Self-selection". In Baltagi, B. A Companion to Theoretical Econometrics. Oxford: Blackwell. doi:10.1002/9780470996249.ch19.
1 2 Puhani, P. (2000). "The Heckman Correction for sample selection and its critique". Journal of Economic Surveys. 14 (1): 53–68. doi:10.1111/1467-6419.00104.
↑ Goldberger, A. (1983). "Abnormal Selection Bias". In Karlin, Samuel; Amemiya, Takeshi; Goodman, Leo. Studies in Econometrics, Time Series, and Multivariate Statistics. New York: Academic Press. pp. 67–84. ISBN 0-12-398750-4.
↑ Newey, Whitney; Powell, J.; Walker, James R. (1990). "Semiparametric Estimation of Selection Models: Some Empirical Results". American Economic Review. 80 (2): 324–28.
↑ Toomet, O.; Henningsen, A. (2008). "Sample Selection Models in R: Package sampleSelection". Journal of Statistical Software. 27 (7): 1–23. doi:10.18637/jss.v027.i07.
↑ "sampleSelection: Sample Selection Models". R Project.
↑ "heckman — Heckman selection model" (PDF). Stata Manual.
↑ Cameron, A. Colin; Trivedi, Pravin K. (2010). Microeconometrics Using Stata (Revised ed.). College Station: Stata Press. pp. 556–562. ISBN 978-1-59718-073-3.

External links

Encyclopædia Britannica article
Article on Heckman-MacFadden Nobel Prize.
Nobel prize official description.

Statistics

Descriptive statistics

Continuous data

Center	Mean arithmetic geometric harmonic Median Mode

Dispersion	Variance Standard deviation Coefficient of variation Percentile Range Interquartile range

Shape	Moments Skewness Kurtosis L-moments

Count data

Index of dispersion

Summary tables

Dependence

Graphics

Data collection

Study design	Population Statistic Effect size Statistical power Sample size determination Missing data

Survey methodology	Sampling Standard error stratified cluster Opinion poll Questionnaire

Controlled experiments	Design control optimal Controlled trial Randomized Random assignment Replication Blocking Interaction Factorial experiment

Uncontrolled studies	Observational study Natural experiment Quasi-experiment

Statistical inference

Statistical theory

Frequentist inference

Point estimation	Estimating equations Maximum likelihood Method of moments M-estimator Minimum distance Unbiased estimators Mean-unbiased minimum-variance Rao–Blackwellization Lehmann–Scheffé theorem Median unbiased Plug-in

Interval estimation	Confidence interval Pivot Likelihood interval Prediction interval Tolerance interval Resampling Bootstrap Jackknife

Testing hypotheses	1- & 2-tails Power Uniformly most powerful test Permutation test Randomization test Multiple comparisons

Parametric tests	Likelihood-ratio Wald Score

Specific tests

Z (normal) Student's t-test F

Goodness of fit	Chi-squared Kolmogorov–Smirnov Anderson–Darling Normality (Shapiro–Wilk) Likelihood-ratio test Model selection Cross validation AIC BIC

Rank statistics	Sign Sample median Signed rank (Wilcoxon) Hodges–Lehmann estimator Rank sum (Mann–Whitney) Nonparametric anova 1-way (Kruskal–Wallis) 2-way (Friedman) Ordered alternative (Jonckheere–Terpstra)

Bayesian inference

Correlation	Pearson product–moment Partial correlation Confounding variable Coefficient of determination

Regression analysis	Errors and residuals Regression model validation Mixed effects models Simultaneous equations models Multivariate adaptive regression splines (MARS)

Linear regression	Simple linear regression Ordinary least squares General linear model Bayesian regression

Non-standard predictors	Nonlinear regression Nonparametric Semiparametric Isotonic Robust Heteroscedasticity Homoscedasticity

Generalized linear model	Exponential families Logistic (Bernoulli) / Binomial / Poisson regressions

Partition of variance	Analysis of variance (ANOVA, anova) Analysis of covariance Multivariate ANOVA Degrees of freedom

Categorical / Multivariate / Time-series / Survival analysis

Categorical

Multivariate

Time-series

General	Decomposition Trend Stationarity Seasonal adjustment Exponential smoothing Cointegration Structural break Granger causality

Specific tests	Dickey–Fuller Johansen Q-statistic (Ljung–Box) Durbin–Watson Breusch–Godfrey

Time domain	Autocorrelation (ACF) partial (PACF) Cross-correlation (XCF) ARMA model ARIMA model (Box–Jenkins) Autoregressive conditional heteroskedasticity (ARCH) Vector autoregression (VAR)

Frequency domain	Spectral density estimation Fourier analysis Wavelet

Survival

Survival function	Kaplan–Meier estimator (product limit) Proportional hazards models Accelerated failure time (AFT) model First hitting time

Hazard function	Nelson–Aalen estimator

Test	Log-rank test

Applications

Biostatistics	Bioinformatics Clinical trials / studies Epidemiology Medical statistics

Engineering statistics	Chemometrics Methods engineering Probabilistic design Process / quality control Reliability System identification

Social statistics	Actuarial science Census Crime statistics Demography Econometrics National accounts Official statistics Population statistics Psychometrics

Spatial statistics	Cartography Environmental statistics Geographic information system Geostatistics Kriging

Category
Portal
Commons
WikiProject

Least squares and regression analysis

Computational statistics

Correlation and dependence

Regression analysis

Regression as a
statistical model

Linear regression	Simple linear regression Ordinary least squares Generalized least squares Weighted least squares General linear model

Predictor structure	Polynomial regression Growth curve (statistics) Segmented regression Local regression

Non-standard	Nonlinear regression Nonparametric Semiparametric Robust Quantile Isotonic

Non-normal errors	Generalized linear model Binomial Poisson Logistic

Decomposition of variance

Model exploration

Background

Design of experiments

Numerical approximation