STATA编程回归问题，在线等！

我要做N个子样本回归，最后把生成的变量合起来再做回归。DO-FILE如下：
gen lnz=.
forvalues i=1(1)31{
reg lnyperl lnkperl if id==i
gen aa=lnyperl-_b[lnkperl]*lnkperl if id==i
replace lnz=aa if id==i
drop aa
}
xtreg return lnz llr lnlperk,fe


结果STATA给我跑了31个全样本回归。。完全不知道怎么回事，在线求助大家！高额悬赏！谢谢！

看帮助里面的例子啊
明显你哪里少写了 引号。
* Generate 100 uniform random variables named x1, x2, ..., x100.

forvalues i = 1(1)100 {
                       generate x`i' = runiform()
                                }

https://bbs.pinggu.org/thread-1422107-1-1.html

呵呵。蓝色正解！