// Muster-Loesung Hausaufgabe 3
// kohler@wzb.eu

// Retrival
// --------

mkdat                                                       ///
    ap6801 bp9301 cp9601 dp9801 ep89 fp108 gp109 hp10901    ///
    ip10901 jp10901 kp10401 lp10401 mp11001 np11701 op12301 ///
    pp13501 qp14301 rp13501 sp13501 tp14201 up14501 vp154   ///
   using ../data/, netto(1) keep(gebjahr sex) clear         ///
   files(p) waves(a b c d e f g h i j k l m n o p q r s t u v )

holrein                                                         ///
 i1110284 i1110285 i1110286 i1110287 i1110288 i1110289 i1110290 ///
 i1110291 i1110292 i1110293 i1110294 i1110295 i1110296 i1110297 ///
 i1110298 i1110299 i1110200 i1110201 i1110202 i1110203 i1110204 ///
 i1110205 using ../data/ ///
 ,  files(pequiv) waves(a b c d e f g h i j k l m n o p q r s t u v )

// Reshape to long
// ---------------

soepren ap6801 bp9301 cp9601 dp9801 ep89 fp108 gp109 hp10901    ///
    ip10901 jp10901 kp10401 lp10401 mp11001 np11701 op12301 ///
    pp13501 qp14301 rp13501 sp13501 tp14201 up14501 vp154, newstub(lsat) waves(1984/2005)

soepren ///
 i1110284 i1110285 i1110286 i1110287 i1110288 i1110289 i1110290 ///
 i1110291 i1110292 i1110293 i1110294 i1110295 i1110296 i1110297 ///
 i1110298 i1110299 i1110200 i1110201 i1110202 i1110203 i1110204 ///
 i1110205, newstub(income) waves(1984/2005)

keep hhnr persnr gebjahr sex lsat* income*

reshape long lsat income, i(persnr) j(wave)

// Clean up
// --------

replace lsat = . if lsat < 0
replace income = . if income == 0

// Individual Dummy Variable Regression
// ------------------------------------

areg lsat income, absorb(persnr)


// De-meaning
// ----------

egen lsatmean = mean(lsat), by(persnr)
egen incomemean = mean(income), by(persnr)
gen lsatdm = lsat - lsatmean
gen incomedm = income - incomemean
reg lsatdm incomedm


exit