肯德尔秩相关系数(Kendall‘s Tau)排名

肯德尔秩相关系数(Kendall's Tau)是一种用于衡量两个排列之间相似性的统计指标,它考虑了元素之间的顺序关系而不考虑具体数值。该系数被广泛用于排序、排名和比较不同实验结果的相关性等领域。

具体而言,肯德尔秩相关系数衡量了两个排列之间逆序对(倒置对)的数量,从而反映了它们的相似程度。逆序对是指在一个排列中,两个元素的顺序与另一个排列中它们的顺序相反。肯德尔秩相关系数的计算步骤如下:

  1. 给定两个排列,比较它们的每一对元素,计算其中一个排列中元素在另一个排列中的顺序关系(逆序或非逆序)。

  2. 计算逆序对的数量。

  3. 根据公式计算肯德尔秩相关系数:

    τ = Number of concordant pairs − Number of discordant pairs Total number of pairs \tau = \frac{\text{Number of concordant pairs} - \text{Number of discordant pairs}}{\text{Total number of pairs}} τ=Total number of pairsNumber of concordant pairs−Number of discordant pairs

    其中,"concordant pairs" 是在两个排列中顺序关系相同的元素对数,而 "discordant pairs" 是顺序关系不同的元素对数。

肯德尔秩相关系数的取值范围在 -1 到 1 之间,具体解释如下:

  • τ = 1 \tau = 1 τ=1:表示两个排列完全一致,没有逆序对。
  • τ = − 1 \tau = -1 τ=−1:表示两个排列完全相反,每个元素都是一个逆序对。
  • τ = 0 \tau = 0 τ=0:表示两个排列之间没有明显的相关性,逆序对和非逆序对的数量大致相等。

肯德尔秩相关系数具有直观性和易于理解的特点,适用于各种排序和排名场景中,但在某些情况下,它可能对一些特殊关系的敏感性较强。

代码

python 复制代码
def kendall_tau(arr1, arr2):
    if len(arr1) != len(arr2):
        raise ValueError("Input arrays must have the same length.")
    
    n = len(arr1)
    concordant_pairs = 0
    discordant_pairs = 0
    
    for i in range(n - 1):
        for j in range(i + 1, n):
            # Compare the order of elements in both arrays
            if (arr1[i] < arr1[j] and arr2[i] < arr2[j]) or (arr1[i] > arr1[j] and arr2[i] > arr2[j]):
                concordant_pairs += 1
            elif (arr1[i] < arr1[j] and arr2[i] > arr2[j]) or (arr1[i] > arr1[j] and arr2[i] < arr2[j]):
                discordant_pairs += 1
    
    tau = (concordant_pairs - discordant_pairs) / (0.5 * n * (n - 1))
    return tau

# Example arrays
arr1 = [1, 2, 3, 4]
arr2 = [2, 1, 4, 3]

# Calculate Kendall's Tau
tau = kendall_tau(arr1, arr2)
print("Kendall's Tau:", tau)
相关推荐
BK贝壳4 小时前
亿级用户签到位图选型复盘:用 bool-hybrid-array 扛冷启动打卡,也把它和 Roaring/Redis 的边界测清楚了
数据结构·python·性能优化·numpy
2601_966949656 小时前
如何根据研究的起止时间截取分钟级行情?Python 实战与 QuantDash 数据处理指南
开发语言·python·pandas·量化交易·1024程序员节·股票数据·quantdash
在世修行10 小时前
智学文枢:页面切换与跨页数据预填
python·桌面应用·signal·1024程序员节·pyside·跨页跳转·状态预填
2601_9628857210 小时前
如何用 Python 把 K 线重采样成任意周期?
开发语言·python
打工仔折腾 AI10 小时前
在 Ubuntu ARM64 部署 Q1 V3.0 MuJoCo 模型:加载、单步推理与执行器响应解析
linux·运维·人工智能·python·深度学习·ubuntu·性能优化
信誓旦旦的程序猿11 小时前
【Python量化实战宝典 #21】北向资金跟随手册:聪明钱择时信号(Python 可跑)
python·1024程序员节·股票数据api·股票数据api接口·股票api数据接口·股票量化数据api·股票量化数据接口
严谨的麻辣烫11 小时前
静态 ISP 代理请求超时怎么排查?Python 与 curl 实战
python·http·requests·1024程序员节·网络排障
KeepRobots13 小时前
11-01 温哥华不回拨:用 60 行 Python 查你的工单系统会不会差一小时(KeepRobots 工程笔记)
python
xiao77mwq13 小时前
你以为 OCR 慢?实测 81.3 s 里全是故障代价:如何把「定格态」消掉,让耗时从 6.5 s 到 13 倍差距显出原形
python
深圳市益普科技有限公司13 小时前
半导体良率管理YMS:良率与设备参数的相关性,先拉哪几组数据
人工智能·python·算法