连老师，
我的模型是这样的:
y=a+b1*x1+b2*x2+u
我对N个样本分别做回归，我想以矩阵形式输出系数b1， b2 和校正R square及其标准误（se）
下边的程序错在哪里呢？

qui tsset
bysort id: gen t = _n
order id t
qui tab id
local N = r(r)
mat Result = J(`N', 1, .)
forvalues i=1/N {
  preserve
    qui keep if id==`i'
    qui tsset t
    reg y x1 x2
    mat b = e(b)
    mat se = vecdiag(cholesky(diag(vecdiag(e(V)))))  
    mat Result[`i', 1] = (b,se,e(r2_a))
   restore  
}
  mat colnames Result = b1 b2 e(r2_a)
  mat list Result

另外：我打不开初级视频；高级的可以打开。
非常感谢！