Spearman and Pearson相关系数矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

右上Spearman左下Pearson相关系数矩阵

1. 1. 输出为csv格式
   
2. 2. 星号为* 0.1 ** 0.05 *** 0.01

复制代码

代码：

//修改部分

local        vlist "Y X1 X2"

local        expFile "Corr_Matrix.csv"

//不用修改部分

preserve

eststo clear

set more off

local upper

local lower `vlist'

expand 2, gen(version)

foreach v of local vlist {

    egen rank = rank(`v') if version == 1

    replace `v' = rank if version ==1

    drop rank

}

foreach v of local vlist {

   estpost correlate `v' `lower' if version == 0

   foreach m in b rho p count {

       matrix `m' = e(`m')

   }

   if "`upper'"!="" {

   estpost correlate `v' `upper' if version == 1

       foreach m in b rho p count {

           matrix `m' = e(`m'), `m'

       }

   }

   ereturn post b

   foreach m in rho p count {

       quietly estadd matrix `m' = `m'

   }

   eststo `v', title(`v')

   local lower: list lower - v

   local upper `upper' `v'

}

/* Export the matrix */

esttab using `expFile', nonumbers mtitles noobs not title("Correlations") star(* 0.1 ** 0.05 *** 0.01) replace

shellout using `expFile'

restore

参考来源：https://stackoverflow.com/questions/25211192/combined-pearson-spearman-rank-correlation-matrix-with-significance-stars-in-sta
https://www.titanwolf.org/Network/q/a9bdb93d-f4f6-45dd-9f89-41f9ef893f95/y
http://repec.org/bocode/e/estout/advanced.html#advanced908

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：spearman pearman pearson 相关系数矩阵 ARMA

格式有问题，此处再写一下，代码内容同上

1. //修改部分
   
2. 
   
3. local        vlist "Y X1 X2"
   
4. local        expFile "Corr_Matrix.csv"
   
5. 
   
6. //不用修改部分
   
7. preserve
   
8. 
   
9. eststo clear
   
10. set more off
    
11. 
    
12. local upper
    
13. local lower `vlist'
    
14. 
    
15. expand 2, gen(version)
    
16. 
    
17. foreach v of local vlist {
    
18.     egen rank = rank(`v') if version == 1
    
19.     replace `v' = rank if version ==1
    
20.     drop rank
    
21. }
    
22. 
    
23. foreach v of local vlist {
    
24. 
    
25.    estpost correlate `v' `lower' if version == 0
    
26.    foreach m in b rho p count {
    
27.        matrix `m' = e(`m')
    
28.    }
    
29. 
    
30.    if "`upper'"!="" {
    
31.    estpost correlate `v' `upper' if version == 1
    
32.        foreach m in b rho p count {
    
33.            matrix `m' = e(`m'), `m'
    
34.        }
    
35.    }
    
36.    ereturn post b
    
37.    foreach m in rho p count {
    
38.        quietly estadd matrix `m' = `m'
    
39.    }
    
40.    eststo `v', title(`v')
    
41.    local lower: list lower - v
    
42.    local upper `upper' `v'
    
43. }
    
44. 
    
45. /* Export the matrix */
    
46. esttab using `expFile', nonumbers mtitles noobs not title("Correlations") star(* 0.1 ** 0.05 *** 0.01) replace
    
47. shellout using `expFile'
    
48. 
    
49. restore

复制代码