对于未分组的变量，是否有一个高效的Python实现用于Somers d？

import pandas as pd
from scipy.stats import rankdata

def somers_d(x, y):
    n = len(x)
    x_ranks = rankdata(x)
    y_ranks = rankdata(y)
    d = sum((x_ranks[i] - y_ranks[i])**2 for i in range(n))
    return 1 - (6 * d) / (n * (n**2 - 1))

# 示例数据
x = [1, 2, 3, 4, 5]
y = [5, 4, 3, 2, 1]

# 计算Somers d
result = somers_d(x, y)
print("Somers d:", result)