/* This is the Presentation Do-file */
capture log close
log using xt_pres.log, replace

clear
set matsize 800
* set mem 10m
use abdata
tsset id year

forvalues yr = 1(1)9 {
	local yr2 = 1975+`yr'
	rename yr`yr' yr`yr2'
}

gen ln=l.n
gen l2n=l2.n

gen lw=l.w

gen lk=l.k
gen l2k=l2.k

gen lys=l.ys
gen l2ys=l2.ys














tsset id year



more
xtdes
more
xtsum n w k ys
more
regress d.n  d.w dl.w d.k dl.k dl2.k d.ys dl.ys dl2.ys 
predict double res_d ,res
more
xtreg n w lw k lk l2k ys lys l2ys, fe
more
predict double e, e
gen le = l.e
more
xtreg n w lw k lk l2k ys lys l2ys, re sa
more
xtreg n w lw k lk l2k ys lys l2ys, mle nolog
more

di "on to testing"
more
xtreg n w lw k lk l2k ys lys l2ys, re sa
xthausman

regress e le w lw k lk l2k ys lys l2ys
more
regress res_d l.res_d  d.w dl.w d.k dl.k dl2.k d.ys dl.ys dl2.ys 
more

di "xtregar section"
more
xtregar n w lw k lk l2k ys lys l2ys, fe lbi 
more
hausman, save
more
xtregar n w lw k lk l2k ys lys l2ys, re lbi 
more
hausman

more
di "xtivreg section"
more
xtivreg n l2.n l(0/1).w l(0/2).(k ys) yr1981-yr1984 (l.n = l3.n), i(id) fd
more

xtivreg n l(0/1).w l(0/2).(k ys) yr1981-yr1984 (l.n= l3.n), fe
more
xtivreg n l(0/1).w l(0/2).(k ys) yr1981-yr1984 (l.n= l3.n), ec2sls
more
xtivreg n l(0/1).w l(0/2).(k ys) yr1981-yr1984 (l.n= l3.n), re
more

di "panel data with lagged endogenous variables section"
more
xtivreg n l2.n l(0/1).w l(0/2).(k ys) yr1981-yr1984 (l.n = l3.n), i(id) fd
hausman, save
more

xtabond n l(0/1).w l(0/2).(k ys) yr1980-yr1984, lags(2)
hausman
more
xtabond n l(0/1).w l(0/2).(k ys) yr1980-yr1984, lags(2) robust
more
xtabond n l(0/1).w l(0/2).(k ys) yr1980-yr1984, lags(2) twostep
more
xtabond n l(0/1).ys yr1980-yr1984 , lags(2) twostep pre(w, lag(1,.)) /*
	*/ pre(k, lag(2,.)) 
more

di "THE END"