legend関数の使い方でググったのでメモ。 par(mfrow=c(2, 2)) dat <- cbind(A=sample(10,3),B=sample(10,3),C=sample(10,3)) #1枚目 matplot(t(dat),col=1:3,lwd=1:3,pch=20,type="b",xlim=c(0,3)+.5,xaxt="n") axis(side=1,at=1:3,labels=colnames(dat)) le…