library(fmsb) MorA <- select.list(choices=c("ASFR","ASMFR"), preselect="ASFR", title="どちらのパタンを描画?") if (MorA=="ASFR") { matplot(Jfert[,1], Jfert[,2:14], type="b", col=1:13, pch=1:13, lty=1, main="日本人女性の年齢別出生率の推移", xlab="女性の年齢", ylab="ASFR") legend("right",col=1:13,pch=1:13,lty=rep(1,13),legend=colnames(Jfert[,2:14])) } else { matplot(Jfert[,1], Jfert[,15:27], type="b", col=1:13, pch=1:13, lty=1, main="日本人女性の年齢別有配偶出生率の推移", xlab="女性の年齢", ylab="ASMFR") legend("right",col=1:13,pch=1:13,lty=rep(1,13),legend=colnames(Jfert[,15:27])) }