R help - Jul 2004 - stl,package=stats

Greetings:
 
I'm using the time series decomposition routine "stl" from the
package "stats".
But how do I get the results into a vector to work with them?
example:
 
data(AirPassengers)
m<-stl(AirPassengers,"per")
print(m)
 
This lists the output but can't figure out how to extract the individual
series like seasonal, trend, irregular.
 
Thanks,
Bob

		
---------------------------------

	[[alternative HTML version deleted]]

bob mccall wrote:
> Greetings:
>  
> I'm using the time series decomposition routine "stl" from
the package "stats".
> But how do I get the results into a vector to work with them?
> example:
>  
> data(AirPassengers)
> m<-stl(AirPassengers,"per")
> print(m)
>  
> This lists the output but can't figure out how to extract the
individual series like seasonal, trend, irregular.
There are two ways:
Either use str(m) to look at the structure of the object or read the 
help page ?stl which tells you that

   ... 'stl returns an object of class "stl" with components
    time.series     a multiple time series with columns
                    seasonal, trend and remainder.' ...

I don't know whether there are any specific extractor methods, but to 
extract the decomposed seasonal values, you can simply write:

    m$time.series[,"seasonal"]

Uwe Ligges

> Thanks,
> Bob
> 
> 		
> ---------------------------------
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html

On Sun, 2004-07-18 at 18:38, bob mccall wrote:
> Greetings:
>  
> I'm using the time series decomposition routine "stl" from
the package "stats".
> But how do I get the results into a vector to work with them?
> example:
>  
> data(AirPassengers)
> m<-stl(AirPassengers,"per")
> print(m)
>  
> This lists the output but can't figure out how to extract the
individual series like seasonal, trend, irregular.
# Seasonal
as.vector(m$time.series[,1])
# Trend
as.vector(m$time.series[,2])
# Remainder
as.vector(m$time.series[,3])
>  
> Thanks,
> Bob
> 
> 		
> ---------------------------------
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html