similar to: How to interpret the phase spectrum?

I am having a hard time interpreting the phase and coherency plots. x is two timeseries that occur at the same time i.e. a b 1 11.2 12.3 16 11.3 12.4 31 11.4 12.5 46 11.5 12.6 ...etc even though my example is does not show this they are oscillating at more or less the same frequency just shifted by t=x (imagine two sine waves offset with the 2nd sine

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

I am reading some documentation about Cross Spectrum Analysis as a technique to compare spectra. My understanding is that it estimates the correlation strength between quasi-periodic structures embedded in two signals. I believe it may be useful for my signals analysis. I was referred to the R functions that implement this type of analysis. I tried all the examples which generated a series of

Dear useRs, I have a curve which is a mixture of Gaussian curves (for example UV emission or absorption spectrum). Do you have any suggestions how to implement searching for optimal set of Gaussian peaks to fit the curve? I know that it is very complex problem, but maybe it is a possibility to do it? First supposement is to use a nls() with very large functions, and compare AIC value, but it is

Dear all, I am searching the period of a time series usering R. Is there some lag window functions in R? Could you give me some books about spectral analysis usering R? best wishes, Wang [[alternative HTML version deleted]]

Dear R users I have two questions about estimating the spectral power of a time series: 1) I came across a funny thing with the following code: data(co2) par(mfrow=c(2,1)) co2.sp1<-spectrum(co2,detrend=T,demean=T,span=3) co2.sp2<-spectrum(co2[1:468],detrend=T,demean=T,span=3) The first plot displays the frequencies ranging from 0 to 6 whearas the second plot displays the same curve but

Dear All, How to do Hodrick-Prescott Filter in R? Thanks! [[alternative HTML version deleted]]

Dear all, What is the usage of the ".Call()"? What is the meaning of the follows: S=.Call("RS_fractal_spectral_density_function_direct",x,as.vector(taper.),as.logical(center),as.logical(recenter),TRUE,as.integer(npad),COPY=rep(FALSE,6),CLASSES=c(rep("matrix",2),rep("logical",3),"integer"),PACKAGE="ifultools") And is a function

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi I have 10.000 simulations for a sensitivity analysis. I have done a few sensitivity analysis for different response variables already, but now, as most of the simulations (if not all) show some cyclic behaviour, see how the independent input parameter influence the frequency of the cyclic changes and "how cyclic" they actually are. So

Hi, I'm wondering why the spectrum() phase of quadrature couple isn't purely +/-pi. But mostly, I'm looking for a recommended way to take a 2-D spectrum and convert it into a single complex array. Kindly consider: # 10 Hz sine wave 10 seconds long sampled at 50 Hz deltaT = 1/50 t = seq(0, 10, deltaT) w = 2 * pi * 10 x = ts( sin( w * t ), deltat = deltaT ) y = ts( sin(

Dear R users, I'm looking for a test of bimodality in order to make some decisions about how to procede with an analysis algorithm. I have not come across any such tests in my readings and discussions apart from the Rao which appears to be applicable to cyclic data. The data I'm interested in characterizing as uni- or bimodal are frequency x amplitude spectra of consonant speech sounds,

Hello. Does anyone know how to generate a power spectrum from the STFT function in package e1071? The output for this function supplies the Fourier coefficients, but I do not know how to relate these back to the power spectrum. Alternatively, if anyone is aware of other packages/functions that perform short time Fourier transforms, that would also be helpful. Thanks. -- View this message in

Hi, everyone, I tried to use "spectrum()" or "spec.pgram()" to get a periodogram of a time series but they didn't work. Even the examples given in the help file didn't work (all with the same error message, below). And the 'ts'ibrary was loaded with "library(ts)" or "library("ts"). I also tried library(tseries) but got the same problem.

hi R users I have a large 1D dataset and some of them is NA value . I found I cound get the spectrum by such a command. ua=c£¨10£¬30 £¬40£¬50£¬NA£© spectrum(ua) and I could not use na.rm just like mean or sd function How could I get the spectrum of ua ? thank you . -- TANG Jie [[alternative HTML version deleted]]

TB --- 2008-08-27 11:26:00 - tinderbox 2.3 running on freebsd-stable.sentex.ca TB --- 2008-08-27 11:26:00 - starting RELENG_7 tinderbox run for amd64/amd64 TB --- 2008-08-27 11:26:00 - cleaning the object tree TB --- 2008-08-27 11:26:30 - cvsupping the source tree TB --- 2008-08-27 11:26:30 - /usr/bin/csup -r 3 -g -L 1 -h localhost -s /tinderbox/RELENG_7/amd64/amd64/supfile TB --- 2008-08-27

Hi everybody! I am searching for an R implementation of multi-scale singular spectrum analysis (SSA), which was introduced by Yiou etal (2000) (Data-adaptive wavelets and multi-scale singular-spectrum analysis) and further described in Ghil etal (2002). For SSA alone I found an package recently published on r-forge: http://r-forge.r-project.org/projects/simsalabim/ I am happy about any

Dear list, I apologies first for my English, hope you will understand well my question. I am working on 1/2 hour piezometric data, time unit is second. They present daily oscillation when using the spectrum() function. What I am really interested in, is to find the amplitude corresponding to this oscillation. I work with a college using Matlab, and although we apply the same methodology, our

Hi, I have a plot that is essentially the same as that in Figure 5.6 of "Lattice - Multivariate Data Visualization with R". The key difference is that I would like to add a title to top of the grey color spectrum legend, but have thus far been unsuccessful. I've tried a variety of options, but to no avail. I'm new to R, and admittedly I don't understand the

Hi everybody! I am searching for an R implementation of multi-scale SSA (singular spectrum analysis), which was introduced by Yiou etal (2000) (Data-adaptive wavelets and *multi*-*scale *singular-spectrum analysis<http://linkinghub.elsevier.com/retrieve/pii/S0167278900000452>) and further described in Ghil etal (2002). For SSA alone I found an package recently published on r-forge:

Reading Monty's message, I better re-post this one :-). In another thread the following was said: > It might help if you also mentioned that the sound quality of Ogg is > significantly better than MP3 (assuming that you use regular MP3, not > MP3Pro), especially at 32kbps. Speaking of which, would there be anything to stop Vorbis, maybe not v1 but in the future, to use similar