2015-04-23 208 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