103のコードを改変して、2群のデータを折れ線グラフで描いてみます。
#Data
x<-c(0,1,2,4)
y1<-c(1.1,1,0.9,1.2)#対象群
y2<-c(1,2,3,4)#実験群
#Plot
plot(x,y1,type="o",pch=21,bg="white",
xlim=c(0,4),ylim=c(0,5),axes=FALSE,ann=FALSE)
box()
par(new=T)
plot(x,y2,type="o",pch=21,bg="black",
xlim=c(0,4),ylim=c(0,5),axes=FALSE,ann=FALSE)
#Y-axis
axis(side=2,at=seq(0,5,1),las=1
mtext("Plasma Y (ng/ml)",side=2,line=2)
#X-axis
axis(side=1,at=seq(0,4,1))
mtext("Time (h)",side=1,line=2)
#Title
mtext("Title",side=3,line=1)
Rの新規文書ファイルにコピーして、すべての行を選択して、編集メニューから実行を選択します。
２つのデータをひとつのグラフにまとめて表示させるときは、二つ目のデータをplotする前に、par(new=T)と書きます。
par(new=T)とplot(...)の行を増やせば、3群のデータも表示できます。
#Data
x<-c(0,1,2,4)
y1<-c(1.1,1,0.9,1.2)
y2<-c(1,2,3,4)
y3<-c(1.2,3,4,5)
#Plot
plot(x,y1,type="o",pch=21,bg="white",
xlim=c(0,4),ylim=c(0,5),axes=FALSE,ann=FALSE)
box()
par(new=T)
plot(x,y2,type="o",pch=21,bg="black",
xlim=c(0,4),ylim=c(0,5),axes=FALSE,ann=FALSE)
par(new=T)
plot(x,y3,type="o",pch=21,bg="red",
xlim=c(0,4),ylim=c(0,5),axes=FALSE,ann=FALSE)
#Y-axis
axis(side=2,at=seq(0,5,1))
mtext("Plasma Y (ng/ml)",side=2,line=2)
#X-axis
axis(side=1,at=seq(0,4,1),las=1)
mtext("Time (h)",side=1,line=2)
#Title
mtext("Title",side=3,line=1)
3群のデータを表示できました。