Statistical methods to correct for endogeneity problems
Control functions (also known as two-stage residual inclusion) are statistical methods to correct for endogeneity problems by modelling the endogeneity in the error term. The approach thereby differs in important ways from other models that try to account for the same econometric problem. Instrumental variables, for example, attempt to model the endogenous variable "X" as an often invertible model with respect to a relevant and exogenous instrument "Z". Panel analysis uses special data properties to difference out unobserved heterogeneity that is assumed to be fixed over time.
Control functions were introduced by Heckman and Robb although the principle can be traced back to earlier papers. A particular reason why they are popular is because they work for non-invertible models (such as discrete choice models) and allow for heterogeneous effects, where effects at the individual level can differ from effects at the aggregate. A well-known example of the control function approach is the Heckman correction.
Formal definition.
Assume we start from a standard endogenous variable setup with additive errors, where "X" is an endogenous variable, and "Z" is an exogenous variable that can serve as an instrument.
 ;templatestyles src="Numbered block/styles.css" / ;
 ;templatestyles src="Numbered block/styles.css" / ;
 ;templatestyles src="Numbered block/styles.css" / ;
 ;templatestyles src="Numbered block/styles.css" / ;
A popular instrumental variable approach is to use a two-step procedure and estimate equation (2) first and then use the estimates of this first step to estimate equation (1) in a second step. The control function, however, uses that this model implies
 ;templatestyles src="Numbered block/styles.css" / ;
The function "h"("V") is effectively the control function that models the endogeneity and where this econometric approach lends its name from.
In a Rubin causal model potential outcomes framework, where "Y"1 is the outcome variable of people for who the participation indicator "D" equals 1, the control function approach leads to the following model
 ;templatestyles src="Numbered block/styles.css" / ;
as long as the potential outcomes "Y"0 and "Y"1 are independent of "D" conditional on "X" and "Z".
Variance correction.
Since the second-stage regression includes generated regressors, its variance-covariance matrix needs to be adjusted.
Examples.
Endogeneity in Poisson regression.
Wooldridge and Terza provide a methodology to both deal with and test for endogeneity within the exponential regression framework, which the following discussion follows closely. While the example focuses on a Poisson regression model, it is possible to generalize to other exponential regression models, although this may come at the cost of additional assumptions (e.g. for binary response or censored data models).
Assume the following exponential regression model, where formula_1 is an unobserved term in the latent variable. We allow for correlation between formula_1 and formula_3 (implying formula_3 is possibly endogenous), but allow for no such correlation between formula_1 and formula_6. 
 formula_7
The variables formula_6 serve as instrumental variables for the potentially endogenous formula_3. One can assume a linear relationship between these two variables or alternatively project the endogenous variable formula_3 onto the instruments to get the following reduced form equation:
 ;templatestyles src="Numbered block/styles.css" / ;
The usual rank condition is needed to ensure identification. The endogeneity is then modeled in the following way, where formula_11 determines the severity of endogeneity and formula_12 is assumed to be independent of formula_13. 
 formula_14
Imposing these assumptions, assuming the models are correctly specified, and normalizing formula_15, we can rewrite the conditional mean as follows:
 ;templatestyles src="Numbered block/styles.css" / ;
If formula_12 were known at this point, it would be possible to estimate the relevant parameters by quasi-maximum likelihood estimation (QMLE). Following the two step procedure strategies, Wooldridge and Terza propose estimating equation (1) by ordinary least squares. The fitted residuals from this regression can then be plugged into the estimating equation (2) and QMLE methods will lead to consistent estimators of the parameters of interest. Significance tests on formula_17 can then be used to test for endogeneity within the model.
Extensions.
The original Heckit procedure makes distributional assumptions about the error terms, however, more flexible estimation approaches with weaker distributional assumptions have been established. Furthermore, Blundell and Powell show how the control function approach can be particularly helpful in models with nonadditive errors, such as discrete choice models. This latter approach, however, does implicitly make strong distributional and functional form assumptions.
References.
 ;templatestyles src="Reflist/styles.css" / ;