如何得到pearson的p值矩阵？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如何得到pearson的p值矩阵？

stata似乎在-corr-里提供相关系数矩阵，而没有p值矩阵；而-pwcorr-结果输出窗口显示了p值，但没有p值矩阵。

因为计算相关系数的变量很多，从-pwcorr-结果输出窗口中一一粘贴p值，工作量很大。

请问如何得到p值矩阵？或者如何由相关系数计算p值？或者-pwcorr-结果输出窗口怎么设置，可以整个贴到excel或word里？（直接copy,paste，同一行的系数就是一个单元格，很郁闷！）

谢谢啦！

扫码加我 拉你入群

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

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

分享0 收藏11 回帖

关键词：pearson ear ARS RSO son 矩阵 pearson

*设数据库中各变量都没有缺失值

   
mata
X=st_data(.,.)
C=correlation(X)
for (i=2; i<=rows(C); i++) {
for (j=1; j<=i-1; j++) {
C[j,i]=2*ttail(rows(X)-2,abs(C[i,j]/sqrt((1-C[i,j]^2)/(rows(X)-2))))
}
}
C
end

*以上命令组可显示，左下角是相关系数，右上角是p值，对角线是全是1

[此贴子已经被作者于2009-5-17 10:00:21编辑过]

eblog  金币 +5  金钱 +50  魅力 +15  经验 +50  奖励 2009-5-17 7:26:04

非常感谢！

mata方面我不太懂，但我run一下是可以的。

可为什么p值与-pwcorr,sig－的输出结果不一样呢？

再次谢谢啦！

. pwcorr, sig

             |     code     year       yy
-------------+---------------------------
        code |   1.0000
             |
             |
        year |  -0.1565   1.0000
             |   0.0089
             |
          yy |  -0.0181  -0.0705   1.0000
             |   0.7633   0.2411

mata的结果：

 Cov
                  1              2              3
    +----------------------------------------------+
  1 |             1    .9955258129    .6183583394  |
  2 |  -.1565163265              1     .879447964  |
  3 |  -.0181433598   -.0705359956              1  |
    +----------------------------------------------+

p值不太一样。尤其像code和year，差距很大。

以下是引用jianlamhua在2009-5-16 22:56:00的发言：可为什么p值与-pwcorr,sig－的输出结果不一样呢？

修正了一处。加上了绝对值号。不过，应该还有问题。主要是关于rho的分布的。

可否把你的数据（三个变量）贴一下？

不能上传dta格式的附件吗？我把它转化成excel格式的了。

326210.xls (29 KB)

2009-5-17 00:03:00 上传

[此贴子已经被作者于2009-5-17 0:03:02编辑过]

以下是引用jianlamhua在2009-5-16 23:55:00的发言：不能上传dta格式的附件吗？我把它转化成excel格式的了。

又修正了一下。现在可以了。

（1）概率应乘2（双尾概率）；（2）统计量应取绝对值。

哇赛，sungmoo真是太猛了！而且还very very热心，呵呵。谢谢！

by the way，陈峰书上的相关系数p值是按你原来那样算的，而stata中pearson相关系数是按现在这样算的（命令组run的结果与pearson结果一致）。

哪一种正确的呢？还是有各自的适用前提？我印象中变量服从正态分布即可。。。

我知道了。陈峰书上是tprob,而这里是ttail.实际上是一样的。

thank you a lot again!

我是用matlab的，用corrcoef得到了相关矩阵和P值，相关矩阵的计算公式是知道了，R=C(i,j)/sart(C(i,i)C(j,j))，但是，不知道此处的P值怎么计算出来的，我不会其他的统计软件了，matlab也是才学的，有P值的数学公式就行了，恳请上面2位大虾赐教，我的邮箱是915449026@qq.com，如蒙两位大虾不介意辛苦，可以也回复我的邮箱吗？