1 つのグラフにプロットする必要がある数年の時系列があります。最大の系列の平均は 340、最小は 245、最大は 900 です。最小の系列の平均は 7、最小は -28、最大は 31 です。残りの系列の値は 6 ～ 700 の範囲です。このシリーズは、何年にもわたって定期的な年次および季節的なパターンに従いますが、突然、1か月間気温が上昇し、その後、通常よりもはるかに多くの死亡者が発生しました.
実際のデータを提供することはできませんが、次のデータをシミュレートし、ここにあるサンプル コードに基づいた以下のコードを試しましたhttp://www.r-bloggers.com/multiple-y-axis-in-ar-プロット/ . しかし、プロットは私が望んでいたものを生み出しませんでした. 次の質問があります
- プロットでは、シリーズのいずれかを明確に描写することは難しく、重要な事実は詳細に隠されています。このデータをより適切に提示するにはどうすればよいですか?
- Y 軸の長さが異なります。どうすれば同じ長さの軸を持つことができますか? このコードを改善し、より良いプロットを提示する方法についてのアイデアや提案に感謝します。極端な気象現象の期間を反映する極端な値をシミュレートできないため、シミュレートしたデータは私のデータを反映していません。
どうもありがとう
temp<- rnorm(365, 5, 10)
mort<- rnorm(365, 300, 45)
poll<- rpois(365, lambda=76)
date<-seq(as.Date('2011-01-01'),as.Date('2011-12-31'),by = 1)
df<-data.frame(date,mort,poll,temp)
windows(600,600)
par(mar=c(5, 12, 4, 4) + 0.1)
with(df, {
plot(date, mort, axes=F, ylim=c(170,max(mort)), xlab="", ylab="",type="l",col="black", main="")
points(date,mort,pch=20,col="black")
axis(2, ylim=c(170,max(mort)),col="black",lwd=2)
mtext(2,text="Mortality",line=2)
})
par(new=T)
plot(date, poll, axes=F, ylim=c(45,max(poll)), xlab="", ylab="",
type="l",col="red",lty=2, main="",lwd=1)
axis(2, ylim=c(45,max(poll)),lwd=1,line=3.5)
points(date, poll,pch=20)
mtext(2,text="PM10",line=5.5)
par(new=T)
plot(date, temp, axes=F, ylim=c(-28,max(temp)), xlab="", ylab="",
type="l",lty=3,col="brown", main="",lwd=1)
axis(2, ylim=c(-28,max(temp)),lwd=1,line=7)
points(date, temp,pch=20)
mtext(2,text="Temperature",line=9)
axis(1,pretty(range(date),10))
mtext("date",side=1,col="black",line=2)
