2015-04-23 259 views
0

我有三個數據框,我不知道這很令人興奮的數據,但這不是我想要解決的問題。在R圖中缺少x軸標籤

> 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軸上繪製,與11月 - 二月未能出現。我相信這是一個相當簡單的解決方法,但我似乎無法弄清楚。任何援助將不勝感激!謝謝!

+0

這是因爲你先暗算 「Ascidcv」,不具有所有月份。 – 2015-04-23 12:13:42

回答

3

你可以做你的情節分爲兩步:

第一,你的陰謀沒有軸

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) 

enter image description here

+1

@DavidArenburg,感謝編輯,使它更清晰! :-) – Cath