corsp——教你输出相关系数矩阵

2017-02-28 赵文泽 爬虫俱乐部 爬虫俱乐部
相关关系矩阵  

在实证研究中,变量间相关关系分析是必不可少的。

无论是皮尔森相关系数,还是斯皮尔曼相关系数,

爬虫君今天给你介绍corsp命令,帮你统统搞定!


使用corsp命令可以计算四种相关关系统计指标:皮尔森相关系数(r)、斯皮尔曼相关系数(rho)、肯德尔tau-a相关系数(tau_a)和肯德尔tau-b相关系数(tau_b)。命令每次运算最多同时输出两种相关系数结果,使用格式是:


corsp varlist [if] [in] [ , options ]


corsp命令的使用非常简单,安装完毕后,打开需要处理的数据文件:

sysuse auto,clear

corsp price weight length,lower(rho) upper(r)

上 图输出矩阵就是选择price、weight、length三个变量进行相关分析的结果,相关系数结果保留小数点后四位。其中,lower(rho)和 upper(r)选项对输出矩阵的排列方向进行了选择,lower(rho)指定了斯皮尔曼相关系数(rho)排列在输出矩阵对角线的左下方,同理 upper(r)指定了皮尔森相关系数排列在矩阵对角线右上方。使用lower和upper选项,你也可以指定只输出一种相关系数,以皮尔森相关系数为 例:

如果不使用lower和upper选项,在corsp命令的默认输出矩阵中,对角线右上方是斯皮尔曼相关系数,对角线左下方是皮尔森相关系数。你也可以使用switch选项调换两种相关系数的排列位置。

corsp price weight length

corsp price weight length,switch

你还可以使用pv或sig选项为相关系数添加P值,两个选项作用完全相同。

corsp price weight length,sig

corsp price weight length,pv

添加p值后,可以使用print(#)指定想要显示的相关系数显著水平。如果一个相关系数的p值显著水平超过了设定值,则该相关系数会在输出矩阵中变成空白。

corsp price weight length,sig print(0.00001)

Bonferroni选项和sidak选项分析计算了Bonferroni校正后的p值和sidak校正后的p值。矫正后,可以使用returnrp选项返回没有经过校正的相关系数矩阵。

以上就是corsp命令及全部选项的使用方法,如果你想把相关系数结果保存在word中,并且添加表示显著性的*号,请关注爬虫俱乐部的后续推文。


接下来是空气质量报告

全国空气质量如下

山西的朋友请注意

以上就是今天给大家分享的内容了,说得好就赏个铜板呗!有钱的捧个钱场,有人的捧个人场~,点赞打赏随您心意,么么哒~

应广大粉丝要求,爬虫俱乐部的推文公众号打赏功能可以开发票啦,累计打赏超过1000元我们即可给您开具发票,发票类别为“咨询费”。用心做事,只为做您更贴心的小爬虫。第一批发票已经寄到各位小主的手中,大家快来给小爬虫打赏呀~

文字编辑:梅洁瓷傲

技术总编:刘贝贝



往期推文推荐:

1.合并输出回归结果和其他检验结果——esttab和estadd

2.关于RTF你不知道的命令

3.关于RTF你不知道的命令之番外篇

4.免费事件研究,一片片从邮局寄来

5.免费的股价同步性,一片片从邮局寄来

6.Stata叫你回家听音乐了!

7.一言不合就用stata写邮件(Outlook/Foxmail)

8.玩转stata之调用浏览器

9.I have a Stata, I have a python

10.I have a Stata, I have a Python之二——pdf转word




关于我们

微信公众号“爬虫俱乐部”分享实用的stata命令,欢迎转载、打赏。爬虫俱乐部是由李春涛教授领导下的研究生及本科生组成的大数据分析和数据挖掘团队。

此外,欢迎大家踊跃投稿,介绍一些关于stata的数据处理和分析技巧。

投稿邮箱:xueyuan19920310@163.com

投稿要求:
1)必须原创,禁止抄袭;
2)必须准确,详细,有例子,有截图;
注意事项:
1)所有投稿都会经过本公众号运营团队成员的审核,审核通过才可录用,一经录用,会在该推文里为作者署名,并有赏金分成。
2)邮件请注明投稿,邮件名称为“投稿”+“推文名称”。
3)应广大读者要求,现开通有偿问答服务,如果大家遇到关于stata分析数据的问题,可以在公众号中提出,只需支付少量赏金,我们会在后期的推文里给予解答。

欢迎关注爬虫俱乐部

微信扫一扫
关注该公众号