Jan.Verbesselt at csiro.au
2007-Jul-04 07:58 UTC
[R] how to plot a monthplot from a ts object where all individual years are shown (e.g. as lines) and can be compared with a "average or median " year?
Dear R help, I'm working with regular 8-daily time-series from 2000 up till now and would like to be able to compare years with each other. E.g. by creating a monthplot via the result of the stl() method it looks ok> but I was wondering whether there exist other methods to plot the different years as lines on top of each other such that years can be compared with each other (temporal change detection)? This is the type of time-series we are working with: forest <- ts(data, frequency=46, start=c(2000,8), end=c(2006,46)) I tried it as follows but this is not very clear. year0 <- window(forest, start=2000,end=c(2000,46)) year1 <- window(forest, start=2001,end=c(2001,46)) year2 <- window(forest, start=2002,end=c(2002,46)) year3 <- window(forest, start=2003,end=c(2003,46)) year4 <- window(forest, start=2004,end=c(2004,46)) year5 <- window(forest, start=2005,end=c(2005,46)) year6 <- window(forest, start=2006,end=c(2006,46)) plot(1:46,years[,1], col=2) lines(1:46,years[,2], col=3) lines(1:46,years[,3], col=4) lines(1:46,years[,4], col=5) lines(1:46,years[,5], col=6) lines(1:46,years[,6], col=7) lines(1:46,years[,7], col=8) Are there other options to be able to compare years with each other in order to detect change (e.g., per month)? Thanks a lot, J R 2.5, Win XP [[alternative HTML version deleted]]