我的PANEL DATA有1万家公司的每个月的股票收益（monthly data), 我想针对每家公司的股票收益和市场收益做OLS，但是要每家公司股票收益的第1到第36个观察值与相应时间的市场收益做OLS，求出 市场收益的系数, 然后接着用这个公司的第2到第37个观察值与对应的市场收益做OLS，这样一直下去，直到这个公司最后36个观察值与对应的市场收益做OLS，然后换下一家公司做同样的操作，最后把所得到的所有 市场收益的系数存在一个新的变量里，要对应相应的时间点（第一个系数对应第36个观察值的时间点）。我自己写了个命令，别说最后存在对应时间点了，连估计系数都不能开始，求连教授帮我看一下 xret 代表公司的股票价格， xrm代表市场收益，公司代码已经变成从1到10000

xtset id time
gen a = .
forvalues i = 1/10000 {
     local  i = 1
     local  j = 36
     quietly reg xret xrm in `i++'/`j++' if id == `i'
     replace a = _b[xrm]
}