[R] squared coherency and cross-spectrum

[Ling Jin]

Hi All,

I have two time series, each has length 354. I tried to calculate the 
coherency^2 between them, but the value I got is always 1. On a website, 
it says: " Note that if the ensemble averaging were to be omitted, the 
coherency (squared) would be 1, independent of the data". Does any of 
you know how to specify properly in R in order to get more useful 
coherency? The examples in the help do give coherencies that are not 1s, 
but I did not notice any special specification.

Next question is on co-spectrum. When I supply "spectrum" function with 
multiple time series, it only gives me spectrum (smoothed periodogram) 
of individual time series. Is there any way I can get the 
cross-spectrum? I believe R has calculated it, but I could not find in 
the returned values.

Attached is the smoothed periodogram of the two time series.

Thanks a lot!

Ling

-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.spectrum.png
Type: image/png
Size: 4596 bytes
Desc: not available
Url : https://stat.ethz.ch/pipermail/r-help/attachments/20051201/1beca4ef/test.spectrum.png