liuywustb 发表于 2019-5-21 22:09 
su e
然后利用r命令
g cancha=r(sd)
谢谢解答,那请问前面的分行业分年度的残差要怎么提取呢?
我以为这个命令是可以跑的,但实际操作上只是显示最后一次回归的残差,并且因为分组后最后一组样本量不足,还无法得到残差。
上网查了命令:
gen resid_stddev = .
forvalues year in 1980/2010 {
reg Y X1 X2 X3 if year == `year'
capture drop resid
predict resid if e(sample), r
su resid
replace resid_stddev = r(sd) if year == `year'
}
table year, contents(mean resid_stddev)
在实际操作中也出现了invalid syntax的情况。
谢谢你了!