Python自由度

Question

我试图计算自由度，我使用了以下公式：

不幸的是，它给了我一个错误的答案。在

我的代码：

def stdev(X):
    m = mean(X)
    return math.sqrt(sum((x-m)**2 for x in X) / len(X))

def degreesOfFreedom(X, Y):
    s1 = (stdev(X)**2)
    s2 = (stdev(Y)**2)
    df = (s1 / len(X) + s2 / len(Y))**2 / ((s1 / len(X))**2 / (len(X) - 1) + (s2 / len(Y))**2 / len(Y) - 1)
    return(df)

print('Degrees of freedom for Student-t distribution: ' + str(degreesOfFreedom(A, B)))

其中A和B是包含数值的独立数组。 我没有获得有价值的值，而是看到了以下形式的结果：

Degrees of freedom for Student-t distribution: -0.00004981053615260864

Answer 1

 len(Y))**2 / len(Y) - 1)

因为缺少括号

^{pr2}$

根据博德马斯的说法，除法是首先发生的，因此是错误的