如何用一维数组计算方差？编程新手求均值等统计量实现指导

你的代码核心错误在于方差计算的最后一行多了无意义的- (1)，同时需要明确你要计算的是总体方差还是样本方差，两者的计算逻辑有细微区别：

错误分析 #

你当前代码的最终计算行：

variance_1 = list1_sum / (number_of_elements) - (1);

这里的- (1)完全不符合方差的数学定义，是导致结果错误的直接原因。

修正后的代码 #

1. 计算总体方差 #

如果你的数据代表完整的数据集（总体），方差公式为平方和除以元素总数：

// 循环累加平方和的逻辑是正确的，保留不变
for (element = first_element; element < number_of_elements; element++) {
  list1_sum = list1_sum + pow((list1_numbers[element] - list1_mean), 2);
}

// 修正后的总体方差计算
variance_1 = list1_sum / number_of_elements;

2. 计算样本方差 #

如果你的数据只是总体的一个样本，为了得到无偏估计，分母需要用元素总数减1：

// 循环部分保持不变
for (element = first_element; element < number_of_elements; element++) {
  list1_sum = list1_sum + pow((list1_numbers[element] - list1_mean), 2);
}

// 修正后的样本方差计算
variance_1 = list1_sum / (number_of_elements - 1);

额外优化提示 #

• 变量类型：建议将list1_sum声明为double类型，避免整数除法导致的精度丢失。
• 平方运算效率：计算平方时，直接用(list1_numbers[element] - list1_mean) * (list1_numbers[element] - list1_mean)替代pow函数，效率更高且减少浮点数函数调用的潜在问题。

内容的提问来源于stack exchange，提问作者Clownfish305