我有4组值:y1,y2,y3,y4和一组x.y值具有不同的范围,我需要将它们绘制为单独的曲线,在y轴上具有单独的值集.
简单来说,我需要3个具有不同值(比例)的y轴,以便在同一个图上绘图.
我有几年的时间序列,我需要在一个图表中绘制.最大的系列平均值为340,最小值为245,最大值为900.最小系列的平均值为7,最小值为-28,最大值为31.其余系列的值均为6到700.该系列经历了多年的常规年度和季节性模式,直到突然出现一个月的温度升高,随后死亡人数大大增加.
我无法提供任何实际数据,但我已经模拟了以下数据并尝试了下面的代码,该代码基于此处的示例代码http://www.r-bloggers.com/multiple-y-axis-in-ar-情节/.但情节并没有产生我想要的东西.我有以下问题
非常感谢
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)
Run Code Online (Sandbox Code Playgroud)