<< Chapter < Page Chapter >> Page >
This module contains a brief introduction to the econometric problem of sample selectivity bias using Stata.

Sample selection bias

Introduction

These notes discuss how to handle one of the more common problems that arise in economic analyses—sample selection bias. Essentially, sample selection bias can arise whenever some potential observations cannot be observed. For instance, the students enrolled in an intermediate microeconomics course are not a random sample of all undergraduates. Students self-select when they enroll in any class or choose a major. While we do not know all of the reasons for this self-selection, we suspect that students choosing to take advanced economics courses have more quantitative skills than students choosing courses in the humanities. Since we do not observe the grades that students who did not enroll in the intermediate microeconomics class would have made had they enrolled, we can never observe the grades that they would have made. Under certain circumstances the omission of potential members of a sample will cause ordinary least squares (OLS) to give biased estimates of the parameters of a model.

In the 1970s James Heckman developed techniques that will correct the bias introduced by sample selection bias. Since then, most econometric computer programs include a command that automatically used Heckman’s method. However, blind use of these commands can lead to errors that would be avoided by a better understanding of his correction technique. This module is intended to provide this understanding.

In the first section I discuss the sources of sample selection bias by examining the basic economic model used to understand the problem. In the second section I present the estimation strategy first developed by Heckman. In the third section I discuss how to estimate the Heckman model in Stata . In the final section I examine an extended example of the technique. An exercise is included at the end of the discussion.

The model

Assume that there is an unobserved latent variable, y i , MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyEamaaDaaaleaacaWGPbaabaGaey4fIOcaaOGaaiilaaaa@39B6@ and an unobserved latent index, d i , MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizamaaDaaaleaacaWGPbaabaGaey4fIOcaaOGaaiilaaaa@39A1@ such that:

y i = x i β + ε i where i = 1 , , N ; MathType@MTEF@5@5@+=feaagyart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyEamaaDaaaleaacaWGPbaabaGaey4fIOcaaOGaeyypa0JabCiEayaafaWaaSbaaSqaaiaadMgaaeqaaOGaaCOSdiabgUcaRiabew7aLnaaBaaaleaacaWGPbaabeaakiaabccacaqG3bGaaeiAaiaabwgacaqGYbGaaeyzaiaabccacaWGPbGaeyypa0JaaGymaiaacYcacqWIMaYscaGGSaGaamOtaiaacUdaaaa@4DDA@
d i = z i γ + ν i where i = 1 , , N ; MathType@MTEF@5@5@+=feaagyart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizamaaDaaaleaacaWGPbaabaGaey4fIOcaaOGaeyypa0JabCOEayaafaWaaSbaaSqaaiaadMgaaeqaaOGaaC4SdiabgUcaRiabe27aUnaaBaaaleaacaWGPbaabeaakiaabccacaqG3bGaaeiAaiaabwgacaqGYbGaaeyzaiaabccacaWGPbGaeyypa0JaaGymaiaacYcacqWIMaYscaGGSaGaamOtaiaacUdaaaa@4DD9@
d i = { 1  if  d i > 0 0  if  d i 0   ;  and MathType@MTEF@5@5@+=feaagyart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizamaaBaaaleaacaWGPbaabeaakiabg2da9maaceaaeaqabeaacaaIXaGaaeiiaiaabMgacaqGMbGaaeiiaiaadsgadaqhaaWcbaGaamyAaaqaaiabgEHiQaaakiabg6da+iaaicdaaeaacaaIWaGaaeiiaiaabMgacaqGMbGaaeiiaiaadsgadaqhaaWcbaGaamyAaaqaaiabgEHiQaaakiabgsMiJkaaicdacaqGGaaaaiaawUhaaiaacUdacaqGGaGaaeyyaiaab6gacaqGKbaaaa@50BE@
y i = y i d i . MathType@MTEF@5@5@+=feaagyart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyEamaaBaaaleaacaWGPbaabeaakiabg2da9iaadMhadaqhaaWcbaGaamyAaaqaaiabgEHiQaaakiaadsgadaWgaaWcbaGaamyAaaqabaGccaGGUaaaaa@3EEC@

The matrix notation above means (1) that

  1. x i β = [ 1 x 1 i x K i ] [ β 0 β 1 β K ] = [ 1 x 1 i x K i ] [ β 0 β 1 β K ] = β 0 + β 1 x 1 i + + β K x K i MathType@MTEF@5@5@+=feaagyart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabCiEayaafaWaaSbaaSqaaiaadMgaaeqaaOGaaCOSdiabg2da9maadmaabaqbaeqabqqaaaaabaGaaGymaaqaaiaadIhadaWgaaWcbaGaaGymaiaadMgaaeqaaaGcbaGaeSO7I0eabaGaamiEamaaBaaaleaacaWGlbGaamyAaaqabaaaaaGccaGLBbGaayzxaaWaaWbaaSqabeaakiadacUHYaIOaaWaamWaaeaafaqabeabbaaaaeaacqaHYoGydaWgaaWcbaGaaGimaaqabaaakeaacqaHYoGydaWgaaWcbaGaaGymaaqabaaakeaacqWIUlstaeaacqaHYoGydaWgaaWcbaGaam4saaqabaaaaaGccaGLBbGaayzxaaGaeyypa0ZaamWaaeaafaqabeqaeaaaaeaacaaIXaaabaGaamiEamaaBaaaleaacaaIXaGaamyAaaqabaaakeaacqWIVlctaeaacaWG4bWaaSbaaSqaaiaadUeacaWGPbaabeaaaaaakiaawUfacaGLDbaadaWadaqaauaabeqaeeaaaaqaaiabek7aInaaBaaaleaacaaIWaaabeaaaOqaaiabek7aInaaBaaaleaacaaIXaaabeaaaOqaaiabl6Uinbqaaiabek7aInaaBaaaleaacaWGlbaabeaaaaaakiaawUfacaGLDbaacqGH9aqpcqaHYoGydaWgaaWcbaGaaGimaaqabaGccqGHRaWkcqaHYoGydaWgaaWcbaGaaGymaaqabaGccaWG4bWaaSbaaSqaaiaaigdacaWGPbaabeaakiabgUcaRiabl+UimjabgUcaRiabek7aInaaBaaaleaacaWGlbaabeaakiaadIhadaWgaaWcbaGaam4saiaadMgaaeqaaaaa@7E1C@
  2. z i γ = [ 1 z 1 i z L i ] [ γ 0 γ 1 γ L ] = γ 0 + j = 1 L γ j z j i . MathType@MTEF@5@5@+=feaagyart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabCOEayaafaWaaSbaaSqaaiaadMgaaeqaaOGaaC4Sdiabg2da9maadmaabaqbaeqabeabaaaabaGaaGymaaqaaiaadQhadaWgaaWcbaGaaGymaiaadMgaaeqaaaGcbaGaeS47IWeabaGaamOEamaaBaaaleaacaWGmbGaamyAaaqabaaaaaGccaGLBbGaayzxaaWaamWaaeaafaqabeabbaaaaeaacqaHZoWzdaWgaaWcbaGaaGimaaqabaaakeaacqaHZoWzdaWgaaWcbaGaaGymaaqabaaakeaacqWIUlstaeaacqaHZoWzdaWgaaWcbaGaamitaaqabaaaaaGccaGLBbGaayzxaaGaeyypa0Jaeq4SdC2aaSbaaSqaaiaaicdaaeqaaOGaey4kaSYaaabCaeaacqaHZoWzdaWgaaWcbaGaamOAaaqabaGccaWG6bWaaSbaaSqaaiaadQgacaWGPbaabeaaaeaacaWGQbGaeyypa0JaaGymaaqaaiaadYeaa0GaeyyeIuoakiaac6caaaa@617E@

Substituting (1), (2) and (3) into (4) gives:

y i = { x i β + ε i  if  z i γ + ν i > 0 0            if  z i γ + ν i 0   . MathType@MTEF@5@5@+=feaagyart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyEamaaBaaaleaacaWGPbaabeaakiabg2da9maaceaaeaqabeaaceWH4bGbauaadaWgaaWcbaGaamyAaaqabaGccaWHYoGaey4kaSIaeqyTdu2aaSbaaSqaaiaadMgaaeqaaOGaaeiiaiaabMgacaqGMbGaaeiiaiqahQhagaqbamaaBaaaleaacaWGPbaabeaakiaaho7acqGHRaWkcqaH9oGBdaWgaaWcbaGaamyAaaqabaGccqGH+aGpcaaIWaaabaGaaGimaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeyAaiaabAgacaqGGaGabCOEayaafaWaaSbaaSqaaiaadMgaaeqaaOGaaC4SdiabgUcaRiabe27aUnaaBaaaleaacaWGPbaabeaakiabgsMiJkaaicdacaqGGaaaaiaawUhaaiaac6caaaa@628C@

Note that N is the total sample size and n is the number of observations for which d i = 1. MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizamaaBaaaleaacaWGPbaabeaakiabg2da9iaaigdacaGGUaaaaa@3A74@

Since y i MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyEamaaDaaaleaacaWGPbaabaGaey4fIOcaaaaa@38FC@ is not observed for ( N n ) , MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWaaeaacaWGobGaeyOeI0IaamOBaaGaayjkaiaawMcaaiaacYcaaaa@3AE0@ the question becomes why are these observations missing. A concrete example of such a model is a model of female wage determination. Equation (1) would model the wage rate earned by women in the labor force and Equation (2) would model the decision by a female to enter the labor force. In this case, y i MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyEamaaBaaaleaacaWGPbaabeaaaaa@380C@ , the wage rate woman i receives, is a function of the variables in x i ; MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCiEamaaBaaaleaacaWGPbaabeaakiaacUdaaaa@38D8@ however, women not in the labor force are not included in the sample. If these missing observations are drawn randomly from the population, there is no need for concern. Selectivity bias arises if the ( N n ) MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWaaeaacaWGobGaeyOeI0IaamOBaaGaayjkaiaawMcaaaaa@3A30@ omitted observations have unobserved characteristics that affect the likelihood that d i = 1 MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizamaaBaaaleaacaWGPbaabeaakiabg2da9iaaigdaaaa@39C2@ and are correlated with the wage the woman would receive had she entered the labor force. For instance, a mentally unstable female is likely to earn relatively low wages and might be more unlikely to enter the labor force. In this case, the error terms, ε i MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTdu2aaSbaaSqaaiaadMgaaeqaaaaa@38B5@ and ν i MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyVd42aaSbaaSqaaiaadMgaaeqaaaaa@38C6@ would be independent and identically distributed N ( 0 , ) , MathType@MTEF@5@5@+=feaagyart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOtamaabmaabaGaaGimaiaacYcacqGHris5aiaawIcacaGLPaaacaGGSaaaaa@3C0E@ where

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Econometrics for honors students. OpenStax CNX. Jul 20, 2010 Download for free at http://cnx.org/content/col11208/1.2
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Econometrics for honors students' conversation and receive update notifications?

Ask