各位老师，我想按行业和年份分组回归求残差，应该是使用循环语句，但是我编的程序不对
gen resid = .
forvalues indus=1/21 {
forvalues year=2003/2012 {
        reg var1 var2 if year == `year' & indus == `indus'
        capture drop resid
        predict resid if e(sample), r
        su resid
        replace resid= r(sd) if year == `year' & indus == `indus'
}
}

请问应该怎么修改？