python 正态分布

Question

我有一个数字列表，其中包含这些数字的样本平均值和标准差。现在我正在尝试找出平均值+-SD、平均值+-2SD 和平均值+-3SD 中的数字。例如，在mean+-SD部分，我编写了这样的代码:

ND1 = [np.mean(l)+np.std(l,ddof=1)]    
ND2 = [np.mean(l)-np.std(l,ddof=1)]

m=sorted(l)

print(m)

ND68 = []

if ND2 > m and m< ND1:

    ND68.append(m<ND2 and m>ND1)
    print (ND68)

这是我的问题:1. 能否通过列表计算数量并排列。如果是这样，我哪一部分做错了。或者有一些包我可以用来解决这个问题。

Answer 1

这可能会有所帮助。我们将使用 numpy 来获取您正在寻找的值。在我的示例中，我创建一个正态分布数组，然后使用 bool 切片返回 +/- 1、2 或 3 个标准差之外的元素。

import numpy as np

# create a random normally distributed integer array
my_array = np.random.normal(loc=30, scale=10, size=100).astype(int)

# find the mean and standard dev
my_mean = my_array.mean()
my_std = my_array.std()

# find numbers outside of 1, 2, and 3 standard dev
# the portion inside the square brackets returns an
# array of True and False values.  Slicing my_array
# with the boolean array return only the values that
# are True
out_std_1 = my_array[np.abs(my_array-my_mean) > my_std]
out_std_2 = my_array[np.abs(my_array-my_mean) > 2*my_std]
out_std_3 = my_array[np.abs(my_array-my_mean) > 3*my_std]