我想绘制时间序列的盒图(例如http://www.r-graph-gallery.com/146-boxplot-for-time-series/),并且可以让其他示例工作,禁止使用最后一个示例。我有六年(2011年至2016年)的平均值,并有2014年和2015年的数据(尽管数量很少),但由于某些原因,2014年和2015年的数据没有显示箱子。Boxplot不会绘制所有数据
我的输入数据有三列:年,月和居住索引(0到1之间的一个值)。有多人(本例中为37人),每人每年每月都有平均居住指数(包括2014年和2015年)。
例如:
year month RI
2015 1 NA
2015 2 NA
2015 3 NA
2015 4 NA
2015 5 NA
2015 6 NA
2015 7 0.387096774
2015 8 0.580645161
2015 9 0.3
2015 10 0.225806452
2015 11 0.3
2015 12 0.161290323
2016 1 0.096774194
2016 2 0.103448276
2016 3 0.161290323
2016 4 0.366666667
2016 5 0.258064516
2016 6 0.266666667
2016 7 0.387096774
2016 8 0.129032258
2016 9 0.133333333
2016 10 0.032258065
2016 11 0.133333333
2016 12 0.129032258
其被重复为每个单独的鱼。
我的代码:
#make boxplot
boxplot(RI$RI~RI$month+RI$year,
xaxt="n",xlab="",col=my_colours,pch=20,cex=0.3,ylab="Residency Index (RI)", ylim=c(0,1))
abline(v=seq(0,12*6,12)+0.5,col="grey")
axis(1,labels=unique(RI$year),at=seq(6,12*6,12))
平均趋势线工程按其他例子。
a=aggregate(RI$RI,by=list(RI$month,RI$year),mean, na.rm=TRUE)
lines(a[,3],type="l",col="red",lwd=2)
在这个问题上的任何帮助,将不胜感激。
难道你不能发布数据,或它的一个子集? –
我刚刚添加了一小部分(非常小)的数据。希望你遵循它。 – Rishard