use https://www3.nd.edu/~rwilliam/statafiles/bollenbrand, clear
set matsize 7500
set more off
timer clear


* Replicate Bollen & Brand 2010 Social Forces Paper, Most of Table 3, 
* some Table 5. The bigger models can take hours so be patient if you
* try to run them all!!! If you have Mplus it can run these models
* far more quickly.

* Random Effects Model 1
timer on 1
xtdpdml lnwg hchild, ylag(0) re fiml tfix errorinv /// 
	 ti(Bollen & Brand Social Forces 2010 Random Effects Model 1 Table 3)
est store re1
timer off 1
timer list 1


* Random Effects Model 2
timer on 2
xtdpdml lnwg hchild marr div, ylag(0) re fiml tfix errorinv /// 
	ti(Bollen & Brand Social Forces 2010 Random Effects Model 2 Table 3)
est store re2
timer off 2
timer list 2

* Random Effects Model 3
timer on 3
xtdpdml lnwg hchild marr div eduatt cursc snrpt snrft exppt expft break , ylag(0) re fiml tfix errorinv /// 
	ti(Bollen & Brand Social Forces 2010 Random Effects Model 3 Table 3)
est store re3
timer off 3
timer list 3
	
* Fixed Effects Model 1
timer on 4
xtdpdml lnwg hchild, ylag(0) fiml tfix errorinv /// 
	ti(Bollen & Brand Social Forces 2010 Fixed Effects Model 1 Table 3)
est store fe1
timer off 4
timer list 4

* Fixed Effects Model 2
timer on 5
xtdpdml lnwg hchild marr div, ylag(0) fiml tfix errorinv /// 
	ti(Bollen & Brand Social Forces 2010 Fixed Effects Model 2 Table 3)
est store fe2
timer off 5
timer list 5

* Fixed Effects Model 3
timer on 6
xtdpdml lnwg hchild marr div eduatt cursc snrpt snrft exppt expft break , ylag(0) fiml tfix errorinv /// 
	ti(Bollen & Brand Social Forces 2010 Fixed Effects Model 3 Table 3)
est store fe3
timer off 6
timer list 6

*********************************************

* Some lagged IV models from Table 5.

* Model 5b2, Lagged dv, Tables 5 & 6
timer on 7
xtdpdml lnwg hchild marr div eduatt cursc snrpt snrft exppt expft break , alphafree errorinv fiml tfix /// 
	inv(black hisp) ti(Bollen & Brand Social Forces 2010 Table 6 last model)
est store m5b2
timer off 7
timer list 7