大家看文献,或者自己做文章的时候应该都接触过误差线。误差线是通常用于统计或数据科学,用来显示潜在的误差或相对于系列中每个数据的不确定程度。误差线可以用标准差或标准误差,一般用标准差(standard deviation)。
样本标准差的计算公式是,n是样本的个数
标准误的计算公式是
下面来举个实际的例子,比如说我们手上有4个男性的身高和4个女性的身高。如下表所示
我们先来看无代码版画误差线,需要用到Excel
1)首先我们来计算均值,在Excel里面计算均值使用average函数
2)计算标准差,在Excel里面计算标准差用stdev函数
3)计算标准误,根据前面的公式计算标准误只需要将标准差除以根号n就可以了,sqrt函数是计算平方根的函数。
3)画柱形图
添加误差线
指定误差线的值,这里以标准差为例,正错误值和负错误值都选择E4:F4
这样就能得到下面这张图了,如果只想上方有误差线,在方向的地方选择正偏差就可以了。也可以根据标准误来画误差线。留给大家自己练习吧!
下来我们来用R代码画误差线
1)读入数据
这个数据存放在height.txt文件中
代码语言:javascript复制a=read.table("height.txt",header=T,sep="t")
2)计算均值,方差和标准误
代码语言:javascript复制#计算均值
ave=apply(a,2,mean)
#计算标准差
std=apply(a,2,sd)
#计算标准误
se=apply(a,2,function(x){sd(x)/sqrt(length(x))})
3)画图
代码语言:javascript复制#画柱形图
x<-barplot(ave,xlim=c(0,1),width=0.1,space=0.8,ylim=c(0,200),col="lightblue")
#画误差线
arrows(x0=x,y0=ave std,
x1=x,y1=ave-std, #如果只画上误差线,这里改成y1=ave
angle=90,
code=3,
length=0.04,
lwd=0.4)
也可以根据标准误来画误差线。留给大家自己练习吧!