Jan.Verbesselt at csiro.au
2008-Feb-26 04:02 UTC
[R] adjusting monthplot() towards a seasonal diagnostic plot for stl()
Hi all, I would like to adjust the monthplot() of an stl() so that for a time-series with freq=12 (months): a) the curve on the panel for the k-th month graphs the seasonal values minus their monthly mean values b) add to the fig. the values of the k-th month of the seasonal + remainder, also minus their monthly mean values which corresponds to the 'seasonal diagnostic plot as described by Cleveland et al. 1990 to select ns parameters of the stl() function. How can I use the R code for the monthplot() so that I can adjust the function to graph this? fit <- stl(log(co2), s.window = 20, t.window = 20) plot(fit) monthplot(fit, choice = "seasonal", cex.axis = 0.8, base=mean, type='l') Thanks, Jan