statsby _b, by(year ind) saving(123.dta, replace): reg D_ta D_a D_r D_p, noco
> ns
(running regress on estimation sample)

      command:  regress D_ta D_a D_r D_p, nocons
           by:  year ind

Statsby groups
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
xxxxx.......................................

.
.
. merge m:1 year ind using 123.dta

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                             2,327  (_merge==3)
    -----------------------------------------

. gen yhat = D_a*_b_D_a + D_r*_b_D_r + D_p*_b_D_p+_b_cons
_b_cons not found
r(111);