文章详情

短信预约-IT技能 免费直播动态提醒

请输入下面的图形验证码

提交验证

短信预约提醒成功

斯皮尔曼(spearman)相关系数python代码实现

2023-09-02 15:38

关注

简介

斯皮尔曼等级相关系数(简称等级相关系数,或称秩相关系数,英语:Spearman's rank correlation coefficient或Spearman's ρ)。一般用或者表示。它是衡量两个变量的相关性的无母数指标。它利用单调函数评价两个统计变量的相关性。若数据中没有重复值,且当两变量完全单调相关时,斯皮尔曼相关系数为+1或−1,而且位于-1到1之间如图所示。

更常用的一般为这个公式,但是比较麻烦。一般我们直接调用scipy.stats.spearman()直接调用。

备注:当所有的等级数值都为整数时,可以通过以下简单的公式计算等级相关系数。

斯皮尔曼(等级)系数主要是针对X,Y两个变量求相关性

适用范围

  1. 用于当数据不满足下列条件任意一个(线性关系,连续数据,正态分布)的时候

  1. 当数据为定序数据的时候

定序数据为反应登记的数据,代表了某种逻辑顺序,而且属于品质数据,(甲,乙,丙)(优,良,差)等类型都可以称为定序数据。

补充spearman系数数据排序

备注:如果变量数据相同的话,我们采取依次顺序排序,相同数值取位置的平均值

具体代码

import numpy as npimport pandas as pdimport scipy.stats as statsdef calculate_spearman_correlation(X, Y):    return stats.spearmanr(X, Y)[0]def calculate_spearman_correlation_p(X, Y):    return stats.spearmanr(X, Y)[1]if __name__=='__main__':    x = [1, 2, 3, 4, 5, 6, 7, 8, 9]    y = [2, 1, 2, 4.5, 7, 6.5, 6, 9, 9.5]    print(calculate_spearman_correlation_p(x, y))    print(calculate_spearman_correlation(x, y))

spearman检验

小样本检验(n<30)

在一般情况下,我们通过直接适用查表的方式进行验证,即为spearman系数要大于表中对应的临界值。即为我们认为我们的相关系数大于表中的临界值,我们认为相关系数是有显著性差异的,即为有相关性而且相关性不为0

大样本检验(n>30)

在一般情况下,我们通过构造统计量进行假设检验,假设检验如图,而且在统计量,如果统计量符合正态分布,具体公式如图所示。

H0:
H1:

(构造的统计量符合正态分布)

计算p值,如果p值大于0.05,即为证明有显著差异。即为有相关性,相关系数不为0。

求P值代码

def calculate_spearman_correlation_p(X, Y):    return stats.spearmanr(X, Y)[1]

来源地址:https://blog.csdn.net/xiaziqiqi/article/details/128890810

阅读原文内容投诉

免责声明:

① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的,并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据,供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。

② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341

软考中级精品资料免费领

  • 历年真题答案解析
  • 备考技巧名师总结
  • 高频考点精准押题
  • 2024年上半年信息系统项目管理师第二批次真题及答案解析(完整版)

    难度     813人已做
    查看
  • 【考后总结】2024年5月26日信息系统项目管理师第2批次考情分析

    难度     354人已做
    查看
  • 【考后总结】2024年5月25日信息系统项目管理师第1批次考情分析

    难度     318人已做
    查看
  • 2024年上半年软考高项第一、二批次真题考点汇总(完整版)

    难度     435人已做
    查看
  • 2024年上半年系统架构设计师考试综合知识真题

    难度     224人已做
    查看

相关文章

发现更多好内容

猜你喜欢

AI推送时光机
位置:首页-资讯-后端开发
咦!没有更多了?去看看其它编程学习网 内容吧
首页课程
资料下载
问答资讯