误差棒的添加首先要计算每个分组下的误差值,然后添加。
还是以小编近三周的生活消费数据为例子(炫穷开始)
数据中,money变量为一天的总消费,Lunch为午餐消费
先编写个代码计算每一周的总消费和午餐消费的平均值与误差值,这里需要说明,误差值为计算的标准误,而不是标准差,两者具体的差异可以去网上查阅一下。
mystats <- function(x,no.omit=FALSE){
if(no.omit)
x <- x[!is.na(x)]
mean <- mean(x)
n <- length(x)
se <- sd(x/sqrt(n))
return(c(n=n,mean=mean,se=se))
}
计算代码写好后,可以计算数据的参数
aggregion <- summaryBy(money~week_n,data=df2,FUN=mystats)
结果如下:
然后根据此数据绘制:
ggplot