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 ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.