如果matpe=(1,0.7,0.5\0.7,1,0.6\0.5,0.6,1)

matsp=(1,0.9,0.8\0.9,1,0.4\0.8,0.4,1)

那么希望得到的

matc=(1,0.9,0.8\0.7,1,0.4\0.5,0.6,1)

也就是对角线以下部分是pearson,以上部分是spearman，而矩阵阶数不变。

老登上不，用别人的帐号回复的，请见谅

*有一个比较笨的方法：

  
loc n=colsof(matsp)
forv j=1/`n'{
loc l=`j'+1
forv k=`l'/`n'{
mat matpe[`j',`k']=matsp[`j',`k']
}
}

[此贴子已经被作者于2009-5-14 4:58:34编辑过]

eblog  金币 +5  金钱 +20  魅力 +20  经验 +20  奖励 2009-5-14 7:33:22

谢谢啦。

我试试看。

我一般的做法是使用excel处理。先将输出结果用文本向导导入excel，如果是上三角的话，就再复制一下，用选择性粘贴中的转置选项将其转为上三角即可，但关于显著性水平的画星星，恐怕还得手工来过。

也学这个，谢谢啦~~~~~~~~~~~~~~~~~~~~~~~~~~~

加星号的问题，请参见：
http://www.pinggu.org/bbs/b121i32761p14.html

加星号的问题，请参见：
http://www.pinggu.org/bbs/b121i32761p14.html




Example：

. sysuse nlsw88, clear
(NLSW, 1988 extract)

. pwcorr_a grade wage hours ttl_exp  tenure age
             |    grade     wage    hours  ttl_exp   tenure      age
-------------+------------------------------------------------------
       grade |   1.000   
        wage |   0.325***  1.000   
       hours |   0.087***  0.159***  1.000   
     ttl_exp |   0.198***  0.265***  0.230***  1.000   
      tenure |   0.123***  0.178***  0.161***  0.577***  1.000   
         age |  -0.033    -0.036*   -0.028     0.124***  0.074***  1.000

以下程序主体用9楼sungmoo的代码，然后用estout程序包更加方便地获得格式化的输出（包括小数点位数、tab、rtf、excel等）：

// check if -estout- is installed
capt which estout
if _rc ssc inst estout
// get data
sysuse auto.dta, clear
// collect correlation matrices
corr trunk weight length
matrix matpe = r(C)
spearman trunk weight length
matrix matsp = r(Rho)
matrix both = matpe
// combine the matrices
local I = rowsof(matpe)
forval i = 1/`I' {
local j1 = `i' + 1
forval j = `j1'/`I' {
  matrix both[`i',`j'] = matsp[`i',`j']
}
}
// display and save the result in specified format, - help esttab -
mat l both
esttab matrix(both, fmt(3)) using example, nomtitle rtf

个人认为estout程序包的功能非常强大，值得大家好好看看。

想要加星号和格式化的输出，试试（加红部分可以因需调整，其它部分最好不要改动）：
// check if -estout- is installed
capt which estout
if _rc ssc inst estout

// clear stored estimates
eststo clear

// get data
sysuse auto, clear     // your .dta file here

// main: collect, display and save well-formated correlation matrices
local vlist trunk weight length    // your varlist here
local upper
local lower `vlist'
foreach v of local vlist {
    estpost correlate `v' `lower'
    foreach m in b rho p count {
        matrix `m' = e(`m')
    }
    if "`upper'"!="" {        // compute -spearman- with -corr-
                preserve
                tempvar `v'
                egen ``v'' = rank(`v')
                replace `v' = ``v''
                foreach var of local upper {
                        tempvar `var'
                        egen ``var'' = rank(`var')
                        replace `var' = ``var''
                }
                estpost correlate `v' `upper'
        foreach m in b rho p count {
            matrix `m' = e(`m'), `m'
        }
                restore
    }
    ereturn post b
    foreach m in rho p count {
        quietly estadd matrix `m' = `m'
    }
    eststo `v'
    local lower: list lower - v
    local upper `upper' `v'
}
esttab, nonumbers mtitles noobs not    // display
esttab using example, replace rtf nonumbers mtitles noobs not   // for use in Word

對不起......請問一下，我分別跑出pearson與spearman後，

分別執行 est stroe m1 與 est store m2，再執行

mat a=(0,0\1,0)
mat c=a#m1+a'#m2

stata出現m1 not found的訊息，請問我儲存有問題嗎？謝謝！

est sto不是用来生成矩阵的吧？

大大您好：

執行以下指令，出現錯誤，麻煩請指正，謝謝您！

pwcorr x1 x2 x3
mat pe=r(Rho)
spearman x1 x2 x3
mat sp=r(Rho)
mat a=(0,0\1,0)
mat b=(0,1\0,0)
mat c=a#pe+b#sp

此時出現「conformability error」，請問哪裡出問題了呢？thanks！

很好，受教了！非常感谢！

https://bbs.pinggu.org/thread-1089842-1-1.html
并帖

*还可以借助mata（设考察变量组x*的相关系数阵：左下角是cor x*的结果，右上角是spearman x*的结果；( )内是p值；***/0.01，**/0.05，*/0.1）：

spearman x*
n mata
x=st_data(.,"x*")
c=correlation(x)
n=rows(c)
b=strofreal(lowertriangle(c)+uppertriangle(st_matrix("r(Rho)")), "%9.3f")
p=st_matrix("r(P)")
for (i=2; i<=n; i++) {
for (j=1; j<=i-1; j++) {
p[i,j]=2*ttail(rows(x)-2,abs(c[i,j]/sqrt((1-c[i,j]^2)/(rows(x)-2))))
b[i,j]=b[i,j]+(p[i,j]<0.01?"***":(p[i,j]<0.05?"**":(p[i,j]<0.1?"*":"")))
b[j,i]=b[j,i]+(p[j,i]<0.01?"***":(p[j,i]<0.05?"**":(p[j,i]<0.1?"*":"")))
}
}
c=editvalue(b+J(n,n,"(")+strofreal(p,"%9.3f")+J(n,n,")"), "2.000(0.000)", "1")
c
end

命令倒是简便，但是如何保留小数点后三位和加星号问题？

往帖中不是给出相关链接了吗？https://bbs.pinggu.org/thread-32761-1-1.html

请教请教，我不懂这句话：
“把peason系数矩阵（存为matpe）和spearman系数矩阵(存为matsp)合并."  这里的matpe是将peason 系数矩阵存为matpe.mata文件吗？？？
我用sungmoo的命令依葫芦画瓢作，但系统不认matpe？？？

我知道，上面15和16楼voodoo 给出的命令可以保留小数点和加星号并且把pearson 和spearman 放到一起，但是那个命令很长的撒，我是想问你这个矩阵格式的命令能不能实现这个功能

21# sungmoo

谢谢提供这个命令，很棒。还想请教几个问题：
1.这个命令输出的结果小数位数很多，怎样才能修改一下变成只输出三位小数？
2.correlate 和 pwcorr这两个命令输出的结果不太一样，请问他们都是Pearson的命令吗？有何区别？在您的这个设计中用的是 correlate吗？（下三角的）
3.此外这个命令还能再加上p值吗？
4.这个可以用比如logout 这样的命令输出到excel吗？

非常感谢。

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

相当于用correlate。

大家为什么不愿意自己思考一下扩展别人的命令呢？

strofreal(real matrix R, string matrix format)

肯定可以。这个留给大家思考吧（比如，先想好把p值放在什么位置，p值保留几位小数）。

非常非常感谢。