R でエラーバー付きの折れ線グラフを作成する．

Rで折れ線グラフを作ってみる．

##散布図もしくは折れ線図　エラーバー付き##

#文字サイズ拡大（デフォルトは12）

> par(ps=16)

#データセットはこんな感じ

> dat

tsk1 tsk2 tsk3

1 876.90 979.8 868.3

2 902.20 882.9 941.9

3 843.10 863.9 932.3

4 971.10 955.9 988.5

5 873.40 1022.3 930.6

6 756.80 809.9 766.2

7 779.22 857.5 882.8

8 948.10 999.8 953.0

9 1135.40 1085.7 1186.3

10 1082.60 1034.6 1102.2

11 898.20 994.1 922.7

12 822.50 868.0 911.3

13 950.70 876.5 950.6

14 901.80 963.2 1021.8

> tsk

[1] tsk1 tsk1 tsk1 tks1 tsk1 tsk1 tsk1 tsk1 tsk1

[10] tsk1 tsk1 tsk1 tsk1 tsk1 tsk2 tsk2 tsk2 tsk2

[19] tsk2 tsk2 tsk2 tsk2 tsk2 tsk2 tsk2 tsk2 tsk2

[28] tsk2 tsk3 tsk3 tsk3 tsk3 tsk3 tsk3 tsk3 tsk3

[37] tsk3 tsk3 tsk3 tsk3 tsk3 tsk3

Levels: tsk3 tsk3 tsk3

>

> tapplydat<-tapply(dat,tsk,mean)

> tapplydat

tsk1 tsk2 tsk3

910.1443 942.4357 954.1786

> tapplydatsd<-tapply(dat,tsk,sd)

> tapplydatsd

tsk1 tsk2 tsk3

104.82401 82.15649 101.29027

#作図

> datplot<-plot(1:3,tapplydat,type="b",xlim=c(0.75,3.25),cex=1.5,pch=1,xlab="",ylab="dat [mm]",ylim=c(800,1100), xaxt="n",byt="l")

#エラーバー追加

> arrows(1:3,tapplydat-tapplydatsd,1:3,tapplydat+tapplydatsd,lwd=1.0,angle=90,length=0.1,code=3)

#X軸を追加

> axis(1, 1:3, c("task1", "task2", "task3"))