Dear r-help list: In the following R session, I seem to be somehow breaking the spectrum() function, but I'm not sure how. Could somebody please point out my mistake? My apologies if it's something that should be obvious. The mysterious bit is how spectrum(my.ts) at first works, but then later fails with an error message that I find difficult to interpret. I'm running R 1.3.0 that I compiled myself on what was once Caldera Linux (2.2.14), but I dare say I've deviated somewhat from the original Caldera distribution (added GNOME and generally updated software). Thanks in advance. Regards, Mike ---- session transcript follows ---- R : Copyright 2001, The R Development Core Team Version 1.3.0 (2001-06-22) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or `licence()' for distribution details. R is a collaborative project with many contributors. Type `contributors()' for more information. Type `demo()' for some demos, `help()' for on-line help, or `help.start()' for a HTML browser interface to help. Type `q()' to quit R.> ls()character(0)> library(tseries)`tseries' version: 0.7-4 `tseries' is a package for time series analysis with emphasis on non-linear modelling. See `library (help=tseries)' for details.> my.ts<-ts(rnorm(200)) > spectrum(my.ts) # nice graphic display appears > plot(my.ts) # another nice graphic display appears > spectrum(my.ts) # As expected, yet another nice graphic displayappears> perf<-read.delim('vmout3.txt') > attach(perf) > ls(perf)[1] "X" "avm" "b" "cs" "cy" "fr" "fre" "idl" "in" "pi" "po" "r" [13] "re" "sr" "sy" "sys" "t" "usr" "wai"> idl.ts<-ts(idl)> spectrum(idl.ts) # No graphic appearsError in "[<-"(*tmp*, , i, value = c(w, rep(1, nr - 2 * m), rev(w)) * : number of items to replace is not a multiple of replacement length In addition: Warning messages: 1: longer object length is not a multiple of shorter object length in: pi * seq(1, 2 * m - 1, by = 2) 2: longer object length is not a multiple of shorter object length in: c(w, rep(1, nr - 2 * m), rev(w)) * x[, i]> spectrum(my.ts) # I would have thought this would have workedError in "[<-"(*tmp*, , i, value = c(w, rep(1, nr - 2 * m), rev(w)) * : number of items to replace is not a multiple of replacement length In addition: Warning messages: 1: longer object length is not a multiple of shorter object length in: pi * seq(1, 2 * m - 1, by = 2) 2: longer object length is not a multiple of shorter object length in: c(w, rep(1, nr - 2 * m), rev(w)) * x[, i]>-------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20010719/4bf2b5bc/attachment.html
MichaelNielsen@synergy.com.au
2001-Jul-19 01:28 UTC
[R] Re: Strange behaviour of spectrum()?
Speaking of spectrum, my face has now been red-shifted.
I guess it's pretty dumb having a variable named 'pi' in a data
frame that
I've attached. This is obviously newbie sloppiness, but it raises the
question about how R-guru's manage inadvertent clobbering of magic names.
Constructive suggestions would be warmly received.
Thanks.
Regards,
Mike
----- Forwarded by Michael Nielsen/Syd/Synergy on 19/07/2001 11:22 AM
-----
Michael Nielsen
19/07/2001 11:13 AM
To: r-help at stat.math.ethz.ch
cc:
Subject: Strange behaviour of spectrum()?
Dear r-help list:
In the following R session, I seem to be somehow breaking the spectrum()
function, but I'm not sure how. Could somebody please point out my
mistake? My apologies if it's something that should be obvious.
The mysterious bit is how spectrum(my.ts) at first works, but then later
fails with an error message that I find difficult to interpret.
I'm running R 1.3.0 that I compiled myself on what was once Caldera Linux
(2.2.14), but I dare say I've deviated somewhat from the original Caldera
distribution (added GNOME and generally updated software).
Thanks in advance.
Regards,
Mike
---- session transcript follows ----
R : Copyright 2001, The R Development Core Team
Version 1.3.0 (2001-06-22)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type `license()' or `licence()' for distribution details.
R is a collaborative project with many contributors.
Type `contributors()' for more information.
Type `demo()' for some demos, `help()' for on-line help, or
`help.start()' for a HTML browser interface to help.
Type `q()' to quit R.
> ls()
character(0)> library(tseries)
`tseries' version: 0.7-4
`tseries' is a package for time series analysis with emphasis
on non-linear modelling.
See `library (help=tseries)' for details.
> my.ts<-ts(rnorm(200))
> spectrum(my.ts) # nice graphic display appears
> plot(my.ts) # another nice graphic display appears
> spectrum(my.ts) # As expected, yet another nice graphic display
appears> perf<-read.delim('vmout3.txt')
> attach(perf)
> ls(perf)
[1] "X" "avm" "b" "cs"
"cy" "fr" "fre" "idl" "in"
"pi" "po" "r"
[13] "re" "sr" "sy" "sys"
"t" "usr" "wai"> idl.ts<-ts(idl)
> spectrum(idl.ts) # No graphic appears
Error in "[<-"(*tmp*, , i, value = c(w, rep(1, nr - 2 * m), rev(w))
* :
number of items to replace is not a multiple of replacement length
In addition: Warning messages:
1: longer object length
is not a multiple of shorter object length in: pi * seq(1, 2 * m -
1, by = 2)
2: longer object length
is not a multiple of shorter object length in: c(w, rep(1, nr - 2
* m), rev(w)) * x[, i]
> spectrum(my.ts) # I would have thought this would have worked
Error in "[<-"(*tmp*, , i, value = c(w, rep(1, nr - 2 * m), rev(w))
* :
number of items to replace is not a multiple of replacement length
In addition: Warning messages:
1: longer object length
is not a multiple of shorter object length in: pi * seq(1, 2 * m -
1, by = 2)
2: longer object length
is not a multiple of shorter object length in: c(w, rep(1, nr - 2
* m), rev(w)) * x[, i] >
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
https://stat.ethz.ch/pipermail/r-help/attachments/20010719/a6ede6d0/attachment.html