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

   
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编辑过]

非常感谢！

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，差距很大。

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

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

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

326210.xls
大小:(29 KB)

 马上下载

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