<br><font size=2 face="sans-serif">Dear r-help list:</font>
<br>
<br><font size=2 face="sans-serif">In the following R session, I seem to be somehow breaking the spectrum() function, but I'm not sure how. &nbsp;Could somebody please point out my mistake? &nbsp;My apologies if it's something that should be obvious.</font>
<br>
<br><font size=2 face="sans-serif">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.</font>
<br>
<br><font size=2 face="sans-serif">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).</font>
<br>
<br><font size=2 face="sans-serif">Thanks in advance.</font>
<br>
<br><font size=2 face="sans-serif">Regards,</font>
<br>
<br><font size=2 face="sans-serif">Mike</font>
<br>
<br><font size=2 face="sans-serif">---- session transcript follows ----</font>
<br>
<br><font size=2 face="sans-serif">R : Copyright 2001, The R Development Core Team<br>
Version 1.3.0 &nbsp;(2001-06-22)</font>
<br>
<br><font size=2 face="sans-serif">R is free software and comes with ABSOLUTELY NO WARRANTY.<br>
You are welcome to redistribute it under certain conditions.<br>
Type `license()' or `licence()' for distribution details.</font>
<br>
<br><font size=2 face="sans-serif">R is a collaborative project with many contributors.<br>
Type `contributors()' for more information.</font>
<br>
<br><font size=2 face="sans-serif">Type `demo()' for some demos, `help()' for on-line help, or<br>
`help.start()' for a HTML browser interface to help.<br>
Type `q()' to quit R.</font>
<br>
<br><font size=2 face="sans-serif">&gt; ls()<br>
character(0)<br>
&gt; library(tseries)</font>
<br>
<br><font size=2 face="sans-serif">&nbsp; &nbsp; &nbsp; &nbsp; `tseries' version: 0.7-4 <br>
 <br>
 &nbsp; &nbsp; &nbsp; &nbsp;`tseries' is a package for time series analysis with emphasis<br>
 &nbsp; &nbsp; &nbsp; &nbsp;on non-linear modelling.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;See `library (help=tseries)' for details.</font>
<br>
<br><font size=2 face="sans-serif">&gt; my.ts&lt;-ts(rnorm(200))<br>
&gt; spectrum(my.ts) &nbsp; &nbsp; &nbsp; # nice graphic display appears<br>
&gt; plot(my.ts) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; # another nice graphic display appears<br>
&gt; spectrum(my.ts) &nbsp; &nbsp; &nbsp; # As expected, yet another nice graphic display appears<br>
&gt; perf&lt;-read.delim('vmout3.txt')<br>
&gt; attach(perf)<br>
&gt; ls(perf)<br>
 [1] &quot;X&quot; &nbsp; &quot;avm&quot; &quot;b&quot; &nbsp; &quot;cs&quot; &nbsp;&quot;cy&quot; &nbsp;&quot;fr&quot; &nbsp;&quot;fre&quot; &quot;idl&quot; &quot;in&quot; &nbsp;&quot;pi&quot; &nbsp;&quot;po&quot; &nbsp;&quot;r&quot; &nbsp;<br>
[13] &quot;re&quot; &nbsp;&quot;sr&quot; &nbsp;&quot;sy&quot; &nbsp;&quot;sys&quot; &quot;t&quot; &nbsp; &quot;usr&quot; &quot;wai&quot;<br>
&gt; idl.ts&lt;-ts(idl)</font>
<br>
<br><font size=2 face="sans-serif">&gt; spectrum(idl.ts) &nbsp; &nbsp; &nbsp;# No graphic appears<br>
Error in &quot;[&lt;-&quot;(*tmp*, , i, value = c(w, rep(1, nr - 2 * m), rev(w)) * &nbsp;: <br>
 &nbsp; &nbsp; &nbsp; &nbsp;number of items to replace is not a multiple of replacement length<br>
In addition: Warning messages: <br>
1: longer object length<br>
 &nbsp; &nbsp; &nbsp; &nbsp;is not a multiple of shorter object length in: pi * seq(1, 2 * m - 1, by = 2) <br>
2: longer object length<br>
 &nbsp; &nbsp; &nbsp; &nbsp;is not a multiple of shorter object length in: c(w, rep(1, nr - 2 * m), rev(w)) * x[, i] </font>
<br>
<br>
<br><font size=2 face="sans-serif">&gt; spectrum(my.ts) &nbsp; &nbsp; &nbsp; # I would have thought this would have worked<br>
Error in &quot;[&lt;-&quot;(*tmp*, , i, value = c(w, rep(1, nr - 2 * m), rev(w)) * &nbsp;: <br>
 &nbsp; &nbsp; &nbsp; &nbsp;number of items to replace is not a multiple of replacement length<br>
In addition: Warning messages: <br>
1: longer object length<br>
 &nbsp; &nbsp; &nbsp; &nbsp;is not a multiple of shorter object length in: pi * seq(1, 2 * m - 1, by = 2) <br>
2: longer object length<br>
 &nbsp; &nbsp; &nbsp; &nbsp;is not a multiple of shorter object length in: c(w, rep(1, nr - 2 * m), rev(w)) * x[, i] <br>
&gt; </font>