[Stata]报表呈现（1）：相关分析

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

【stata代码模板】相关分析_correlate_pwcorr_pcorr_spearman_kendal 转自www.jiayounet.com 感谢数说工作室

（1）      Correlate计算pearson相关系数

pearson相关系数是最常用的相关分析依据，要求变量服从正态分布，代码为：

——————————————模板——————————————

1. correlate 变量1 变量2 … if var=value,means covariance

复制代码

——————————————模板——————————————

此代码用来计算变量1、变量2…之间的相关系数矩阵。

If用来筛选满足条件的数据，可缺省。

Means要求显示一半描述统计量（均值、标准差、最小最大值）。

Covariance要求显示协方差矩阵而不是相关系数矩阵。

（2）      pwcorr推断总体相关系数

只用correlate计算样本的相关系数，会受到抽样波动的影响，样本相关不能说明总体相关。Pwcorr不仅可以计算相关系数，还可对相关系数显著性进行检验，原假设是总体相关系数为0，即不相关。代码为：

————————————模板————————————

1. pwcorr 变量1 变量2 … if var=value,sig star(#)

复制代码

————————————模板————————————

Sig显示相关系数显著性的p值。

Star(#)是当相关系数显著性检验的p值小于设定的#%，则在相关系数旁标记星号。

比如，计算语文成绩与数学成绩的相关性，显著性水平为5%：

Pwcorr Chinese maths,sig star(5)

（3）      pcorr计算偏相关系数

   偏相关分析的任务就是在研究两变量之间的相关关系时控制可能对其影响的其他变量。

————————————模板——————————————

1. pcorr 变量1 其他变量1 其他变量2 … if var=value

复制代码

————————————模板——————————————

上面代码是考察其他变量与变量1的偏相关系数。

（4）      等级相关

   当变量联合正态分布不能满足，或者变量为定序与定类变量时，不宜使用pearson，而应使用spearman等级相关或kendal等级相关。

————————————模板————————————

1. spearman/kendal 变量1 变量2 … if var=value,star(#)

复制代码

————————————模板————————————

语句开头用spearman或kendal都可以，分别代表两种相关系数的检验。

Star(#)是当相关系数显著性检验的p值小于设定的#%，则在相关系数旁标记星号。

---

扫码加我 拉你入群

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

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

分享0 收藏12 回帖

关键词：Stata tata 相关分析 covariance correlate 正态分布 pearson 标准差 最大值 统计

呵呵，师弟有心，这个系列贴做得好。谢谢分享。

niuniuyiwan 发表于 2015-7-9 13:09
【stata代码模板】相关分析_correlate_pwcorr_pcorr_spearman_kendal转自www.jiayounet.com感谢数说工作室（ ...

怎么感觉手机上看这些代码很费劲，是软件的问题还是手机的问题哦

tkt718 发表于 2015-9-5 00:35
怎么感觉手机上看这些代码很费劲，是软件的问题还是手机的问题哦

用电脑来浏览比较好些，祝学习愉快。

1. //以下介绍将上述命令如pwcor执行时的结果输出到word/excel中方法：
   
2. //方法一：以auto.dta为例计算皮尔逊系数
   
3. 
   
4. clear all
   
5. sysuse auto
   
6. local vbles price rep78 weight length mpg
   
7. local vbles: list sort vbles
   
8. capture postutil clear
   
9. tempfile corrs
   
10. postfile handle str32 var1 str32 var2 float r_ using `corrs'
    
11. local nvars: word count `vbles'
    
12. forvalues i = 1/`nvars' {
    
13.      local v1: word `i' of `vbles'
    
14.      forvalues j = `=`i'+1'/`nvars' {
    
15.           local v2: word `j' of `vbles'
    
16.           corr `v1' `v2'
    
17.           post handle  ("`v1'") ("`v2'") (`r(rho)')
    
18.      }
    
19. }
    
20. postclose handle
    
21. 
    
22. //将长型数据转为宽型
    
23. use `corrs', clear
    
24. reshape wide r_, i(var1) j(var2) string
    
25. rename r_* *  
    
26. order _all, alphabetic
    
27. order var1, first
    
28. sort var1
    
29. browse //查看
    
30. 
    
31. *By Clyde Schechter 04 Aug 2014, 19:37

复制代码

1. //方法二：logout
   
2. 
   
3. ssc install logout
   
4. 
   
5. clear all
   
6. sysuse auto
   
7. local mylist price length weight rep78
   
8. logout, save ("d:\StataMP\ado\personal\data\auto.xml") excel replace:pwcorr `mylist', star(.05)
   
9. !=start auto.xml
   
10. //其中save()中内容为路径+文件名
    
11. //!=start 为windows下打开文件的方法(也可在窗口中点击dir)

复制代码

1. ssc install polychoric
   
2. sysuse auto
   
3. polychoric mpg length, pw
   
4. matrix a=matrix(r(R))
   
5. putexcel A7=matrix(a,names) using results.xls ,replace

复制代码

方法四：polychoric

1. //安装
   
2. ssc install matpwcorr
   
3. 
   
4. sysuse auto,clear
   
5. 
   
6. //创设矩阵
   
7. 
   
8. matpwcorr price mpg rep78 headroom trunk length weight
   
9. 
   
10. //列出矩阵
    
11. 
    
12. mat list corr

复制代码

方法五：matpwcorr

1. //方法6 ：makematrix
   
2. ssc install makematrix
   
3. use hsb2.dta
   
4. makematrix, from(r(rho)): corr read-socst
   
5. //指定单独列变量
   
6. makematrix, from(r(rho)) cols(socst) format(%9.2f): corr read-socst

复制代码

数据.rar (2.37 KB) 本附件包括：

• hsb2.dta

2015-9-14 17:36:32 上传

1. //使用plotmatrix
   
2. 
   
3. ssc install plotmatrix
   
4. 
   
5. sysuse auto  
   
6. 
   
7. matrix accum R = price mpg trunk weight length turn , nocons dev
   
8. matrix R = corr(R)
   
9. 
   
10. plotmatrix, m(R) s(-1(0.25)1) c(red) freq legend(off)

复制代码