当手动计算浮点数列表的方差时,我使用“reduce”但发现它与我的预期略有偏差(给定numpy.var)。然后我使用列表理解重新计算它,并得到我期望的值。
sumSqrdReduce = reduce((lambda total, val: total+(val - mean)**2), lst)
sumSqrdComprehension = sum([(val-mean)**2 for val in lst])
演示此问题的示例列表:
lst = [0.53839998, 4.36650467, 3.64258786, 3.62987329, -0.33371547, 10.16436997, 3.11141481, 4.62991016, 0.72292498, -2.9477603, 4.0144724, 7.14428721, -3.05925725, 4.83175576, 5.55112354, 5.03295696, -2.40226829, 1.87662003, -1.02187228, 5.25553533, 1.54985611, 2.71460086, 0.83336707, -3.3935002, 3.88551682, -2.47155389, 1.76985117, 3.57110149, -5.17191153, 4.80879124, -0.97037815, 0.99500531, -0.22062183, 9.96261967, 3.31320864, 0.39606156, -2.71492665, 0.31085669, -1.82838686, 0.38113291, 2.7265862, 6.46300302, 3.11995554, 0.15073258, 12.03547416, 4.82310128, 2.43649615, 3.2195886, 2.84891094, 9.75191341]
使用上面的列表(平均值= 2.4008863134):
sumSqrdReduce = 671.241430925
sumSqrdComprehension = 674.171886287
我是否错误地执行了减少?或者这是一个通用的“浮点积累错误”,如果是这样,为什么这两种方法不会复制相同的浮点不准确?我希望每个方法与“真相”的任何差异都是一样的,希望不会有太大的不同。
答案 0 :(得分:5)
您确实正确执行了缩减。您将lst
的第一个元素作为初始total
,而不考虑其与平均值的平方偏差。
如果您真的想使用传统的函数编程函数,(val - mean)**2
的{{1}}部分在概念上更适合reduce
:
map
或者您可以指定初始累加器值0.0:
reduce(lambda x, y: x+y, map(lambda x: (x-mean)**2, lst))
由于您已将问题标记为NumPy,因此,如果出于某种原因您希望避免使用NumPy数组reduce((lambda total, val: total+(val - mean)**2), lst, 0.0)
,请执行以下操作:在arr
:
numpy.var