[R] Newbie Time Series Questions
Jason Turner
jasont at indigoindustrial.co.nz
Thu Oct 17 10:41:25 CEST 2002
On Thu, Oct 17, 2002 at 12:39:08PM -0700, 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
> ...
>
> 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.
looks like a data frame, then. You've got to R it's not a multivariate
time series; R thinks it's one series for Oct, one for Nov, etc.
So, convert to matrix, and you'll need to:
1) transpose it, so the *column* vectors are sequential in time (by default,
when R unrolls a matrix to a vector, it stacks the columns, not the rows).
2) convert the matrix to a vector, so R knows it's univariate.
3) *then* convert to time series. whew. ;-)
klam.ts <- ts(as.vector(t(klam.flow)),frequency=12,start=c(1961,10))
klam.ts
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1961 1461 1716 2524
1962 1773 1906 2005 1756 1575 1387 983 1094 1382 1907 2253 1985
1963 1907 1769 1676 2634 1386 929 766 968 1309
If that's wrong - i.e. each row in klam.flow is one year and the columns
are just in a funny order - you'll have to juggle the columns before step
1 above. Like so:
> klam.ts <- ts(as.vector(t(klam.flow[,c(4:12,1:3)])),frequency=12,start=c(1961,1))
> klam.ts
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1961 1773 1906 2005 1756 1575 1387 983 1094 1382 1461 1716 2524
1962 1907 1769 1676 2634 1386 929 766 968 1309 1907 2253 1985
Cheers
Jason
--
Indigo Industrial Controls Ltd.
64-21-343-545
jasont at indigoindustrial.co.nz
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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