楼主: xtcdw
12769 21

[其他] 请教:怎么求n个变量两两间的协方差? [推广有奖]

  • 1关注
  • 0粉丝

高中生

92%

还不是VIP/贵宾

-

威望
0
论坛币
52 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
490 点
帖子
39
精华
0
在线时间
37 小时
注册时间
2009-11-2
最后登录
2016-7-12

楼主
xtcdw 发表于 2012-4-7 19:17:48 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
有x1-xn共n个数据t期的数据,现在要求两两变量间的协方差,得到如下形式:
cov(x1,x2),cov(x1,x3),cov(x1,x4).......cov(x1,xn),cov(x2,x3),cov(x2,x4),......,cov(x2,xn),cov(x3,x4),cov(x3,x5)......cov(x3,xn),请高人帮忙,感激不尽!
二维码

扫码加我 拉你入群

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

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

关键词:协方差 感激不尽

沙发
sungmoo 发表于 2012-4-8 07:23:27
cor x1-xn,c
已有 1 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
Sunknownay + 20 + 8 + 1 + 1 + 1 热心帮助其他会员

总评分: 经验 + 20  论坛币 + 8  学术水平 + 1  热心指数 + 1  信用等级 + 1   查看全部评分

藤椅
xtcdw 发表于 2012-4-8 12:00:06
谢谢sungmoo 的回答,还有个问题就是,我现在的数据形式是
stkd        r
1        0.024445976
1        0.069237111
1        0.020535734
1        0.090597795
1        0.017095258
1        0.072745767
2        0.043986666
2        0.021858773
2        0.035048333
2        0.062093105
2        0.033950999
2        0.058640246
3        0.067889655
3        0.051035785
3        0.096212028
3        0.034610521
3        0.008738863
3        0.07686777
......
现在要得到两两间的协方差该怎么做?

板凳
sungmoo 发表于 2012-4-8 13:39:27
bys stkd: g i=_n
reshape wide r,i(i) j(stkd)
cor r*,c

*若要把方差阵生成为一个r.xls文件(在Stata默认目录中):
mat r=r(C)
clear
svmat r
outsheet using r.xls


已有 1 人评分学术水平 热心指数 信用等级 收起 理由
ywh19860616 + 1 + 1 + 1 精彩帖子

总评分: 学术水平 + 1  热心指数 + 1  信用等级 + 1   查看全部评分

报纸
xtcdw 发表于 2012-4-8 16:06:15
sungmoo 发表于 2012-4-8 13:39
bys stkd: g i=_n
reshape wide r,i(i) j(stkd)
cor r*,c
谢谢sungmoo版主的回复,可以了。如果我想将协方差的数据输出到excel表中该怎么做呢?

地板
sungmoo 发表于 2012-4-8 18:12:56
不太复杂的结果,直接用右键copy table,粘贴到excel中即可。

7
xtcdw 发表于 2012-4-8 18:23:46
sungmoo 发表于 2012-4-8 18:12
不太复杂的结果,直接用右键copy table,粘贴到excel中即可。
我试过,因为数据量太多,这样直接复制的结果不是太好

8
xtcdw 发表于 2012-4-8 18:34:10
我用logout命令也试了,输出的结果如下:
r600795        r600825        r600827        r600875        r600886        r600887        r600900       
r600795        0.000486                                               
r600825        0.000506        0.00244                                       
r600827        0.000304        0.000673        0.000390                               
r600875        0.000124        0.000351        0.000402        0.000658                       
r600886        0.000203        0.000767        0.000329        0.000325        0.000318               
r600887        -6.50e-05        -9.10e-05        -3.10e-05        1.00e-05        -2.40e-05        1.10e-05       
r600900        8.10e-05        0.000304        0.000181        0.000237        0.000164        -4.00e-06        9.90e-05
r600970        0.000605        2.50e-05        0.000143        -0.000178        -1.40e-05        -8.70e-05        -5.30e-05
r600971        0.000210        -0.000122        0.000176        0.000243        6.80e-05        -1.30e-05        6.90e-05
                                                       
r600970        r600971                                               
r600970        0.00106                                               
r600971        0.000270        0.000239                                       
上面还有很大一部分数据

9
xtcdw 发表于 2012-4-9 10:58:22
sungmoo 发表于 2012-4-8 13:39
bys stkd: g i=_n
reshape wide r,i(i) j(stkd)
cor r*,c
. mat a=r(C)
. svmat a
“number of observations will be reset to 1
Press any key to continue, or Break to abort
obs was 0, now 1“
输出的文件里面没有数据,请问是怎么回事

10
sungmoo 发表于 2012-4-9 11:06:12
有没有clear?

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群
GMT+8, 2026-1-3 06:44