Re: [R] rollmean and stl

From: Achim Zeileis <Achim.Zeileis_at_wu-wien.ac.at>
Date: Tue, 22 Jul 2008 23:10:11 +0200 (CEST)

On Tue, 22 Jul 2008, stephen sefick wrote:

> I need to investigate how rollmean and the trend returned from stl differ.
> I am trying to find out exactly what the trend part of stl is (I have just
> started coding in R and do not know fortran). I need to extract this
> because it will be used in further calculations, and it needs to be verified
> to make sure that I am using the right process. I would like to use this
> to remove the seasonal component of a time series to investigate the
> advective transport of oxygen and the local production.
>
> I tried
> d <- stl(data, s.window="periodic")
> d$trend
> and I get NULL

Look at the "Value" section of ?stl:

time.series: a multiple time series with columns 'seasonal', 'trend'

           and 'remainder'.

i.e., use

d$time.series[,"trend"]

> thanks
>
> Stephen
>
> --
> Let's not spend our time and resources thinking about things that are so
> little or so large that all they really do for us is puff us up and make us
> feel like gods. We are mammals, and have not exhausted the annoying little
> problems of being mammals.
>
> -K. Mullis
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help_at_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.
>
>



R-help_at_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. Received on Tue 22 Jul 2008 - 21:16:03 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 22 Jul 2008 - 21:32:51 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive