[R] Newbie Time Series Questions

Robert Schick robert.schick at noaa.gov
Thu Oct 17 23:04:10 CEST 2002


Thanks - That did the trick. Curious why the following doesn't quite
work:
klam.ts <- ts(data=klam.flow, start=c(1960,10), frequency=12)

The above doesn't make the vector.

Rob

Achim Zeileis wrote:
> 
> Robert Schick wrote:
> >
> > I have a data set of monthly river flows from 1960-2000, which are
> > similar in structure to the nottem data:
> >
> > > klam.flow
> >       Oct  Nov  Dec  Jan  Feb   Mar  Apr  May  Jun  Jul  Aug  Sep
> > 1961 1461 1716 2524 1773 1906  2005 1756 1575 1387  983 1094 1382
> > 1962 1907 2253 1985 1907 1769  1676 2634 1386  929  766  968 1309
> > ...
> >
> > I tried plotting with
> >
> > > ts.plot(klam.flow)
> >
> > Which quickly led me to realize the data aren't a "ts" object, as it
> > plotted 40 stacked lines instead of one long series. I've tried
> > converting to this using ts or as.ts, e.g.
> >
> > > klam.ts <- ts(data=klam.flow, start=1960,end=2000, frequency=12)
> 
> If I understand your data correctly they start in October 1961? Then you
> should do something like
> 
> R> klam.flow2 <- ts(as.vector(t(klam.flow)), start = c(1961, 10), freq =
> 12)
> R> plot(klam.flow2)
> 
> This turns your time series to one long vector first and then adds the
> "ts" properties.
> Z
> 
> > But this seems to explode the number of time series, e.g.
> >
> > > klam.ts
> >           Oct  Nov  Dec  Jan  Feb   Mar  Apr  May  Jun  Jul  Aug  Sep
> > Jan 1960 1461 1716 2524 1773 1906  2005 1756 1575 1387  983 1094 1382
> > Feb 1960 1907 2253 1985 1907 1769  1676 2634 1386  929  766  968 1309
> > Mar 1960 2511 2852 3661 2103 2189  2548 3841 2937  857  743 1058 1574
> > ...
> >
> > Any advice on how to properly create the ts object so it looks and plots
> > a la Figure 13.15 in VR?
> >
> > I'm using R 1.60 on a Windows 2000 box.
> >
> > Thanks
> > --
> > Rob Schick
> > Research Associate
> > NOAA Fisheries
> > Santa Cruz Lab
> > 110 Shaffer Road
> > Santa Cruz, CA 95060
> > Phone: 831.420.3960
> > -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> > 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
> > _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

-- 
Rob Schick
Research Associate
NOAA Fisheries
Santa Cruz Lab
110 Shaffer Road
Santa Cruz, CA 95060
Phone: 831.420.3960
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list