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)

由於某些原因，x軸上僅繪制了9月，10月和3月，而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)