clear all set more off pause on set obs 1 expand 1000 gen id=_n gen b_human=. gen se_human=. gen t_human=. gen p_human=. gen b_non=. gen se_non=. gen t_non=. gen p_non=. save temp, replace forvalues k = 12345 12346 to 13345 { clear all set seed `k' set obs 1 expand 41 gen firmid=_n expand 15 bysort firmid: gen year=_n gen human=uniform()