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