Rでエラーバーが付いた棒グラフを描いてみる．
##棒グラフ エラーバー付き##
#文字サイズ拡大（デフォルトは12）
> par(ps=16)
#データセットはこんな感じ
> dat1
con1 con2 con3
1 9.12 12.68 8.85
2 9.55 11.22 11.36
3 8.35 15.06 12.38
4 13.96 18.71 16.65
5 16.53 31.17 21.74
6 9.64 9.26 7.02
7 10.46 28.45 12.62
8 8.96 13.78 15.76
9 19.29 21.07 26.91
10 4.74 16.14 8.00
11 22.70 31.77 29.02
12 16.77 14.46 22.13
13 16.81 8.94 21.36
14 11.23 17.26 9.01
#groupの定義付けをする
> group<-c("task1","task2","task3")
#まず，平均値を出す
> dat1means<-mean(dat1)
> dat1means
con1 con2 con3
12.72214 17.85500 15.91500
> mdat1means<-t(matrix(dat1means,3,1))
> colnames(mdat1means)=group
> rownames(mdat1means)="average"
> mdat1means
con1 con2 con3
average 12.72214 17.855 15.915
#標準偏差（SD）も出す（今回はエラーバーをSDにしただけ，別にSDじゃなくてもよい）
> dat1sd<-sd(dat1)
> mdat1sd<-t(matrix(dat1sd,3,1))
> colnames(mdat1sd)=group
> rownames(mdat1sd)="sd"
> mdat1sd
con1 con2 con3
sd 5.028285 7.634092 7.226356
#作図する
> dat1plot<-barplot(mdat1means,beside=TRUE,space=c(1,1),ylim=c(0,30),ylab="dat1 [m/sec]",col="grey",cex.lab=1.5,cex.axis=1.5)
> box(bty="L")
エラーバーを追加する
> arrows(dat1plot,mdat1means,dat1plot,mdat1means+mdat1sd,lwd=1.0,angle=90,length=0.1,code=3)
これでおしましい．
そのうち，図を添付します．
コメント