用Python编程EP检验

Question

请问，如何用Python编程这段程序？这是统计里面的EP检验

Answer 1

计算部分应该如下：
import math


l = [141, 148, 132, 138, 154, 142, 150, 146, 155, 158,
     150, 140, 147, 148, 144, 150, 149, 145, 149, 158,
     143, 141, 144, 144, 126, 140, 144, 142, 141, 140,
     145, 135, 147, 146, 141, 136, 140, 146, 142, 137,
     148, 154, 137, 139, 143, 140, 131, 143, 141, 149,
     148, 135, 148, 152, 143, 144, 141, 143, 147, 146,
     150, 132, 142, 142, 143, 153, 149, 146, 149, 138,
     142, 149, 142, 137, 134, 144, 146, 147, 140, 142,
     140, 137, 152, 145]

n = len(l)
x_ = sum(l) / n
m2 = (1 / n) * sum((l[i] -x_) **2 for i in range(n))
A = sum(math.exp((-(l[i] -x_) **2) /(4 * m2)) for i in range(n))

B = 0
for k in range(1, n):
    for j in range(k):
        B += math.exp(-((l[j] - l[k]) **2) / (2 *m2))

T_EP = 1 + n / 3 ** 0.5 + 2 * B / n - 2 ** 0.5 * A
print(f"n={n}\nx_={x_}\nm2={m2}\nA={A}\nB={B}\nT_EP={T_EP}")