R 图中缺少 x 轴标签

Missing x axis labels in R plot

我有三个数据框,我知道不是很令人兴奋的数据,但这不是我要解决的问题。

  > 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)

出于某种原因,x 轴上仅绘制了九月、十月和三月,十一月至二月未出现。我确信这是一个非常简单的修复,但我似乎无法弄清楚。任何帮助将不胜感激!谢谢!

您可以分两步完成绘图:

首先,你的情节没有轴

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)