许多paper里经常能看到Mean±SD(SE)这样的表达方式,或者在图表里用SD或者SE来表示error bar,用SD的居多,但是也有不少用SE的。初学者很容易混淆SD(standard deviation)和SE(standard error)。
SD
SD我们都很清楚,是表达数据的离散程度,然后实际应用中很多数据具有近似正态分布的概率分布,有了SD,我们就可以大致估计数据的范围,譬如经典的"68-95-99.7法则",即约 68% 数值分布在距离平均值有 1 个标准差之内的范围,约 95% 数值分布在距离平均值有 2 个标准差之内的范围,以及约 99.7% 数值分布在距离平均值有 3 个标准差之内的范围。如下图:
SE
SE是什么呢,一般来说,自然界里很难获得总体数据,我们只能用样本(无论是各种实验还是社会调查抽样)去近似估计总体,这样问题就来了,估计的准不准(平均值)?
我们可以理论上这样做,既然不能获得总体,我们可以尽可能多(无限)的从标准差为σ的总体数据里抽取大小为 n 的样本,每个样本各有一个平均值,所有样本平均值的标准差就可以用"68-95-99.7法则"评估准不准了(这就是所谓的置信区间),样本平均值的标准差可以被证明如下公式表达:
但由于通常σ为未知,此时可以用研究中取得样本的标准差 (S) 来估计:
这就是SE的来源,即样本平均值的SD,我们用MATLAB编程演示下。
MATLAB演示
以掷硬币为例,掷100次,统计正面(1)的次数,共统计1000次作为一个sample,然后我们这样采1000个sample(程序不怕累的)。
先贴结果,代码放后面。结果如下,我们能看到样本平均值的SD基本等于样本的SE!
完整代码如下:
%1000个样本(sample)
sample_mean = []
for m=1:1000
%掷硬币100次,统计正面的次数,共统计1000次作为一个sample
sample = []
for i = 1:1000
box = randi([0,1],100,1);
sample = [sample length(box(box==1))];
end
sample_mean = [sample_mean mean(sample)];
end
%sample的标准差和标准误
SD_sample = std(sample)
SE_sample = SD_sample/sqrt(1000)
%sample平均值的标准差
SD_sample_mean = std(sample_mean)
结论
1、标准差(SD)更能反应离散程度。
paper里需要Mean±SD这个信息,就是便于读者进行判断数据的离散性,e.g.,一般我们把偏离平均值2或3个SD的值作为outlier(i.e., 异常值)。
2、标准误则比较适合用于评估精确性或准确性的问题。
paper里根据需要也可以提供Mean±SE这个信息,就是便于读者进行判断数据的不确定性,e.g.,95%置信区间是用的Mean ± 2*SE。
无论用哪种表达方式,一定要注意标明,特别是error bar,好的paper都会说明这是什么的。