主对角线上面的部分是Pearson相关系数
主对角线下面的部分是Spearman相关系数
当然，二者互换是一样的喽

请教高手，spearman系数可以一次全部算得吗，还是需要两两去算

https://bbs.pinggu.org/thread-455856-3-1.html

非常感谢~！

不懂mata。。。对你的这个程序不太理解，只能机械套用了。运行中stata报错，调试了半天也没太明白。。。。还得麻烦你帮我看一下例子如下：
sysuse auto, clear   
local vlist trunk weight length  //我把变量列表都放在这里了

spearman `vlist'
n mata
x=st_data(.,"`vlist'")
c=correlation(`vlist')    //这里的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))))           //这里的“x”该放什么？看不懂，惭愧。。。。。
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



谢谢~！