如何在R中使用2个不同的y轴和X轴上的DATE进行绘图？

Question

基本上问题与此问题相同：How can I plot with 2 different y-axes? 除了我在x轴上有一个日期和漂亮（range（））对我使用日期不起作用。

（抱歉，我没有足够的“声誉”来评论上述问题以询问更多信息，所以我需要开始一个新问题）。

我写了以下函数：

  PvsRef<-function(id,feature){
  ymax= max(data[[feature]], na.rm=TRUE)
  plotdata = data[data$ID==id,]

  plot(plotdata$date,plotdata$ref, type="l", col="red", axes=FALSE, ylim=c(0,48), main=paste0("Reference vs.",feature), sub=as.character(id),  xlab="", ylab="")
  axis(2, ylim=c(0,48),col="red",las=1)
  mtext("Reference",side=2,col="red",line=2.5)
  box()

  par(new=TRUE)

  plot(plotdata$date,plotdata[[feature]], type="b", col="blue", axes=FALSE, ylim=c(0,ymax),  xlab="", ylab="")
  mtext(feature,side=4,col="blue",line=4) 
  axis(4, ylim=c(0,ymax), col="blue",col.axis="blue",las=1)

  axis(1,pretty(range(plotdata$date),10))
  mtext("Date",side=1,col="black",line=2.5)  


  axis(1,pretty(range(plotdata$date),10))
  mtext("Date",side=1,col="black",line=2.5)  


  }

然而，x轴显示奇怪的数字而不是日期。

示例数据：

data = data.frame(ID=c(1,1,1,1,1,1), date=as.Date(c("2000-01-01","2000-02-01","2000-03-01","2000-04-01","2000-05-01","2000-06-01")), ref=c(30,23,43,12,34,43), other=c(120,140,230,250,340,440))

然后我想使用

运行该功能

PvsRef(1,"other")

但没有正确的x轴。

编辑：

如果我使用

  axis.Date(1,pretty(range(plotdata$date),10))
  mtext("Date",side=1,col="black",line=2.5)

正确的日期出现了，但是漂亮的要求没有10个蜱。

Answer 1

对于axis.Date，您需要在设置刻度/标签位置时指定参数名称at。如果您未能执行此操作，axis.Date会将第二个参数视为x，这会导致选择合适的标签网格（请参阅?axis.Date）。< / p>

将该行更改为，例如：

axis.Date(1, at=pretty(range(plotdata$date), 10), format='%d %b', las=2, 
          cex.axis=0.8)`

产生你正在寻找的10个滴答。但是，间隔不一致。也许最好指定一个日期序列，例如：

seq(min(plotdata$date), max(plotdata$date), length.out=10)

所以你的功能可以调整为：

f <- function(id, feature) {
  opar <- par()
  on.exit(suppressWarnings(par(opar)))
  par(mar=c(6, 4, 4, 5))
  ymax <- max(data[[feature]], na.rm=TRUE)
  plotdata <- data[data$ID==id,]

  plot(plotdata$date, plotdata$ref, type="l", col="red", axes=FALSE, 
       ylim=c(0, 48), main=paste0("Reference vs.", feature), 
       sub=as.character(id), xlab="", ylab="")
  mtext("Reference", side=2, col="red", line=2.5)
  axis(2, ylim=c(0, 48), col='red', col.axis="red", las=1)

  par(new=TRUE)  
  plot(plotdata$date, plotdata[[feature]], type="b", col="blue", axes=FALSE, 
       ylim=c(0,ymax),  xlab="", ylab="")
  mtext(feature, side=4, col="blue", line=3) 
  axis(4, ylim=c(0, ymax), col="blue", col.axis="blue", las=1)

  axis.Date(1, at=seq(min(plotdata$date), max(plotdata$date), length.out=10), 
            format='%d %b', las=2, cex.axis=0.8)
  box(lwd=2)
}  

f(1, 'other')