[R] output from STL

[Brian J. Cuffel, Ph.D.]

The following works for me, where 'outpatient' is a time series object

outpatstl<-stl(outpatient,"periodic")
seasonal<-outpatstl$time.series[,1]
outnoseason<-outpatient-seasonal
plot(outpatient,type='p')
lines(outnoseason,type='l')


-----Original Message-----
From:	Antonio Rodríguez Verdugo [SMTP:rod.chav at hsoft.es]
Sent:	Friday, May 11, 2001 8:14 AM
To:	r-help at stat.math.ethz.ch
Subject:	[R] output from STL

Hi All,

How do I can create a new vector, i.e. 'seasonal' or 'trend' from the resultant seasonal or trend component of the Time.Series object produced by STL, and how I could superpose in the same graphic, i.e. original data and trend or seasonality?
Thanks in advance!
Antonio
Antonio Rodríguez Verdugo
CICEM Agua del Pino
Huelva
Oceanography and Coastal Resources,
PhD Program, 
University of Huelva
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list-Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._