主页 / coding / 问题 / 79591853
Accepted
Vicktor Krum
Asked: 2025-04-25 13:23:17 +0800 CST

Scipy 上的 Z 分数

  • 772

我需要使用 Scipy 找出属于 1 个特定点（即 X 的 1 个值）的 Zscore。

以下是手动代码：

data = [25, 37, 15, 36, 92, 28, 33, 40]
mean = sum(data)/len(data)
summation = 0
for i in range(0, len(data)):
  summation += (data[i]-mean)**2
std = ((1/len(data))*summation))**(1/2)
Z = (x-mean)/std

接下来，当我尝试使用 Scipy 执行相同操作时：

Z = stats.zscore(data)

我得到输出：

[-0.61219538 -0.05775428 -1.07422964 -0.10395771  2.48343411 -0.47358511
 -0.24256798  0.08085599]

可能是因为我只传递了一个参数，也就是数据本身。如何获取只有一个 X 值的 Z 分数？

python

2 个回答

  1. Best Answer

    这个解决方案可能适合你。这个解决方案可能比其他方案更好。

    from scipy import stats
import numpy as np

data = [25, 37, 15, 36, 92, 28, 33, 40]
x = 40 # Please assign a specific value

mean = np.mean(data)
std = np.std(data, ddof=0)
Z = (x - mean) / std
print(Z)

    您可能也喜欢这种方法（scipy.stats.zscore）：

    from scipy import stats
data = [25, 37, 15, 36, 92, 28, 33, 40]
x = 40
all_zscores = stats.zscore(data, ddof=0)
x_index = data.index(x)
Z = all_zscores[x_index]
print(Z)

    输出：

    0.08085599417810461
    • 2

  2. script.stats.zscore 将返回一个 array_like 对象（在本例中为 ndarray）。

    然后，您可以在源列表中搜索您感兴趣的值的索引，并使用该索引从 ndarray 中提取相关数据。

    下面实现了计算 zscores 的“手动”和 scipy 技术，并从一组精心挑选的值中选择值来获得它们的 zscores：

    from scipy.stats import zscore

data = [25, 37, 15, 36, 92, 28, 33, 40]
ld = len(data) # length of the list
mu = sum(data) / ld # mean

# note that we use ld rather than (ld-1) to align with 
# how scipy calculates zscores
sigma = (sum((n - mu) ** 2 for n in data) / ld) ** 0.5

# list of zscores
mz = [(n - mu) / sigma for n in data]

# ndarray of zscores
sz = zscore(data)

# demonstrate that both techniques produce the same results
for a, b in zip(mz, sz):
    sa = f"{a:>7.4f}"
    sb = f"{b:>7.4f}"
    print(sa, sb)
    assert sa == sb

# a selection of values to test
for x in (25, 15, 33, 99):
    try:
        i = data.index(x)
        print(x, sz[i])
    except ValueError:
        print(x, "is not in the source list")

    输出：

    -0.6122 -0.6122
-0.0578 -0.0578
-1.0742 -1.0742
-0.1040 -0.1040
 2.4834  2.4834
-0.4736 -0.4736
-0.2426 -0.2426
 0.0809  0.0809
25 -0.6121953844913635
15 -1.0742296369376756
33 -0.24256798253431383
99 is not in the source list
    • 1

相关问题