请问,如何用Python编程这段程序?这是统计里面的EP检验
计算部分应该如下:
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}")