*program to simulate OLS results program reg_simple, rclass gen dis=rnormal(0,1) gen y = 2 + 5*x + dis regress y x return scalar b2=_b[x] return scalar seb2=_se[x] return scalar b1=_b[_cons] return scalar seb1=_se[_cons] drop dis y end