R图中缺少x轴标签

Question

我有三个数据框，不是我知道的非常令人兴奋的数据，但这不是我想解决的问题。

  > Ascidcv
             Date   Average       SE
    3  2014-09-01 37.250000 6.326673
    15 2014-10-02  6.285714 2.738613
    > Ascidc1
              Date Average SE
    3   2014-10-15       0  0
    34  2014-11-12       0  0
    62  2014-12-11       0  0
    88  2015-02-11       0  0
    119 2015-03-09       0  0
    > Ascidc2
              Date Average SE
    18  2014-10-15       0  0
    48  2014-11-12       0  0
    75  2014-12-11       0  0
    103 2015-02-11       0  0
    135 2015-03-09       0  0

我使用这些数据框来制作情节：

plot(Ascidcv$Date, Ascidcv$Average, type='p', pch=4, col="red", xlab='Date', ylab='', main=expression(italic('Ascidiella sp.')), xlim=c(as.Date("2014-09-01"), as.Date("2015-03-09")), ylim=c(0,120))
points(Ascidc1$Date, Ascidc1$Average, type='p', pch=19, xlab='Date', main=expression(italic('Ascidiella sp.')), xlim=c(as.Date("2014-09-01"),as.Date("2014-12-11")), ylim=c(0,100))
points(Ascidc2$Date, Ascidc2$Average, type='p', pch=2, col="blue", xlab='Date', ylab='Average num ind.', main=expression(italic('Bugula sp.')), xlim=c(as.Date("2014-09-01"),as.Date("2014-12-11")), ylim=c(0,100))
mtext("Average % cover",side=2,line=3)

出于某种原因，仅在9月，10月和3月的月份上绘制了x轴，11月至2月未能出现。我确信这是一个非常简单的修复，但我似乎无法弄明白。任何帮助将不胜感激！谢谢！

Answer 1

你可以分两步完成你的情节：

首先，没有轴的情节

plot(Ascidcv$Date, Ascidcv$Average, type='p', pch=4, col="red", xlab='Date',
     ylab='', main=expression(italic('Ascidiella sp.')), xlim=c(as.Date("2014-09-01"), 
     as.Date("2015-03-09")), ylim=c(0,120),xaxt="n")

然后，在2014年9月至2015年3月的每个月的第一天添加带有刻度的轴 使用axis.Date和相应日期的向量

axis.Date(1, at=as.Date(c(paste0("2014-",c("09",10:12),"-01"),
          paste0("2015-0",1:3,"-01"))))

最后，添加其他点和文字

points(Ascidc1$Date, Ascidc1$Average, type='p', pch=19, xlab='Date', 
       main=expression(italic('Ascidiella sp.')), xlim=c(as.Date("2014-09-01"),
       as.Date("2014-12-11")), ylim=c(0,100))

points(Ascidc2$Date, Ascidc2$Average, type='p', pch=2, col="blue", xlab='Date', 
       ylab='Average num ind.', main=expression(italic('Bugula sp.')), 
       xlim=c(as.Date("2014-09-01"),as.Date("2014-12-11")), ylim=c(0,100))

mtext("Average % cover",side=2,line=3)