关于结果输出

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

连老师，
我的模型是这样的:
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

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

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：结果输出 forvalues forvalue Cholesky preserve

forvalues i=1/N {
should be
forvalues i=1/`N '{

此外，b 和 se 矩阵在你的命令中都只能记录最后一次的结果，他们并未动起来。

初级视频打不开，可以联系人大的邓老师，他们负责处理技术问题。

arlionn 发表于 2011-9-24 20:08
forvalues i=1/N {
should be
forvalues i=1/`N '{

谢谢，能动了，只是se还有些问题