方差函数方程的计算

时间:2017-06-20 13:47:40

标签: statistics

我在此代码中有错误,因为我想计算(x1)和(x2)列表中值之间的差异。任何建议?!

def my_var(L):
    s = 0
    t = 0
    u = 0
    for i in range(0, len(L)):
        s += L[i]
        t = s/len(L)
        u += ((L[i]-t)*(L[i]-t))
    return u / len(L)
x1 = [1, 3, 4, -3, 8]
x2 = [1, -4, 7, 2]
v1 = my_var(x1)
v2 = my_var(x2)
print(v1)
print(v2)

3 个答案:

答案 0 :(得分:1)

根据我学习概率和统计数据的方式,你做错了很多事。您需要计算平均值(平均值),然后将每个值减去平均值减去平方。然后最后取该分子并除以样本大小(n-1)之前除以1。

def my_var(L):
    mean = float(sum(L) / Len(L))
    numerator = 0
    for i in range(0, len(L)):
        numerator += (L[i]-mean)**2
    return numerator / (len(L) - 1)
x1 = [1, 3, 4, -3, 8]
x2 = [1, -4, 7, 2]
v1 = my_var(x1)
v2 = my_var(x2)
print(v1)
print(v2)

不使用sum:

def my_var(L):
    my_sum = 0
    mean = 0
    numerator = 0
    for i in range(0, len(L)):
        my_sum += L[i]
    mean = float(my_sum / len(L))
    for i in range(0, len(L)):
        numerator += (L[i]-mean)**2
    return numerator / (len(L) - 1)
x1 = [1, 3, 4, -3, 8]
x2 = [1, -4, 7, 2]
v1 = my_var(x1)
v2 = my_var(x2)
print(v1)
print(v2)

答案 1 :(得分:0)

尝试numpy。

import numpy as np
x1 = [1, 3, 4, -3, 8]
x2 = [1, -4, 7, 2]
v1 = np.var(x1)
v2 = np.var(x2)

答案 2 :(得分:0)

谢谢@billy_ferguson。我修改了你的代码,它的工作原理。请问我,我仍然是一个业余爱好者但是你可以替换浮点和求和函数并使用更简单的算术运算符作为len(L)和+ =在这一行中的意思=浮点数(sum(L)/ len(L))

def my_var(L):
    mean = 0
    numerator = 0
    for i in range(0, len(L)):
        mean = float(sum(L) / len(L))
        numerator += (L[i]-mean)**2
    return numerator / len(L)
x1 = [1, 3, 4, -3, 8]
x2 = [1, -4, 7, 2]
v1 = my_var(x1)
v2 = my_var(x2)
print(v1)
print(v2)