关于matrix命令

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

clear
scalar beta0 = 1.5
scalar beta1 = 5.6
scalar sigma = 2
forvalues i = 1(1)50 {
    set obs 100
    gen x`i' =_n
    gen u`i' = invnorm(uniform())*sigma
    gen y`i'= beta0 + beta1*x`i'+u`i'      
    quietly reg y`i' x`i'
    matrix b`i' = e(b)
    matrix list b`i'
    matrix rmse`i'= e(rmse)
    matrix list rmse`i'
}
matrix b0 = (0,0)
matrix rmse0 = (0)
forvalues j = 1(1)50 {
   matrix b0 = b0\b`j'
   matrix b`j' = b0
   matrix rmse0 = rmse0\rmse`j'
   matrix rmse`j' = rmse0
   matrix list b`j'
   matrix list rmse`j'
}
matrix list b50
matrix list rmse50
matrix U = J(1,51,1)
matrix list U

以上，虽然看了输出结果和help，但是还是不怎么能理解matrix的意义。请教

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：matrix Mat forvalues forvalue Quietly matrix

这个是在做模拟啦。每次都保留估计出来的系数和rmse。最后b50和rmse50就是50次模拟里得到的系数和rmse值。

夏目贵志 发表于 2015-10-16 23:20
这个是在做模拟啦。每次都保留估计出来的系数和rmse。最后b50和rmse50就是50次模拟里得到的系数和rmse值。

查help的时候有点看不懂
matrix 表示什么？
还有e的function 看help也不大理解。。

夏目贵志 发表于 2015-10-16 23:20
这个是在做模拟啦。每次都保留估计出来的系数和rmse。最后b50和rmse50就是50次模拟里得到的系数和rmse值。

b和rmse的初始值是什么啊 看起来都没有设置过啊

陀花轮1995 发表于 2015-10-17 20:01
b和rmse的初始值是什么啊 看起来都没有设置过啊

不是在这儿吗
matrix b0 = (0,0)
matrix rmse0 = (0)

夏目贵志 发表于 2015-10-17 22:44
不是在这儿吗
matrix b0 = (0,0)
matrix rmse0 = (0)

第一个循环里面就用到了e(b)嘛……前面没有设置呀………………

陀花轮1995 发表于 2015-10-18 18:15
第一个循环里面就用到了e(b)嘛……前面没有设置呀………………

e(b)是Stata自动生成的。不用专门设置。可以直接用的。

夏目贵志 发表于 2015-10-18 22:16
e(b)是Stata自动生成的。不用专门设置。可以直接用的。

e(rmse)也是自动生成的吗？是不是只有在模拟的时候才会生成？
这里的b是专门留给模拟出来的betahat的吗？
那这些beta被stata存储在哪里，我能直接调用这列betahat的数吗
另外，`'（b`i'左右两边的符号）表示什么？
        \ （b0\b`J'）表示什么？（记得Java里\表示整除，这里说不通。。）

陀花轮1995 发表于 2015-10-18 23:54
e(rmse)也是自动生成的吗？是不是只有在模拟的时候才会生成？
这里的b是专门留给模拟出来的betahat的吗？ ...

这个说来话长了。看一下相关的帮助文件吧。
help return
help macro
另外，pdf帮助文件里说的比较详细。尽可能的看一下吧。