sungmoo 发表于 2012-8-19 08:36 
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
谢谢~!