[求助] 求协方差矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

不知道为什么结果提取不出来

corr x1 x2 x3

matrix a=r(C)

matirx list a

结果显示a是1*1的对称矩阵，元素为空，即‘.’

不知道哪错了？

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：协方差矩阵 协方差 matirx matrix 结果显示 矩阵 协方差

这个命令并不保留（save）方差阵的结果，所以引用（提取）不了。

corr x1 x2 x3,c

可以看到方差阵，但它没有被存成一个以后可以调用的矩阵。

没问题

clear all
set obs 10
set seed 12345
forvalues i=1/3 {
gen x`i'=uniform()
}
corr x1 x2 x3
mat list r(C)
mat a=r(C)
mat list a

结果：



symmetric a[3,3]
           x1         x2         x3
x1          1
x2  .22581843          1
x3  .22128659  -.0648734          1

以下是引用eblog在2009-6-2 17:01:00的发言：mat list r(C)

可能是版本不同。我的stata提供不了这样的r(C)。

我用的是10.1版

https://bbs.pinggu.org/thread-455771-1-1.html

这里好像与楼主谈到了同一问题。

嗯，所以要保持更新
更新会修补一些小bug，尤其是一些外部命令，变化更大

好像eblog版主提供的是相关系数矩阵吧，不是协方差阵？

语句本身没有问题，正版下可以出结果

同样版本的stata，

在同学的机上就可以，在我这不行

不知道是不是系统的原因，我的是vista，同学的是xp

以下是引用rayality在2009-6-3 0:49:00的发言：

好像eblog版主提供的是相关系数矩阵吧，不是协方差阵？

呵呵，是的，我只是对按楼主提供的代码进行了验证