[R-sig-eco] Wavelet analysis using simulated distribution

[Sean Downey]

Dear all,

I am using the biwavelet package to run time series analysis, but I am 
trying to calculate a 'robust' version of the 95% CI using simulated 
distribution which I generate under special (non-linear) conditions. 
Stated differently, the standard 'white' and 'red' noise spectra 
normally used to calculate CIs in wavelet analysis are not appropriate 
for my data, but I have a separate simulation process that generates a 
distribution that is an appropriate null model.

So, I have a vector of observed (M) and simulated data (Si). I can 
easily perform the wavelet transform on M and extract the power matrix 
(Mp); similarly, it's also easy enough to iterate through Si and 
generate a distribution of power matrices for Mp. But is this the right 
approach? My question is, what is the best approach for calculating and 
plotting the wavelet CI using my simulated distribution, rather than the 
red or white noise spectrum?

Any help will be greatly appreciated.

-Sean




-- 
Sean S. Downey
Assistant Professor
Department of Anthropology
University of Maryland
1111 Woods Hall
College Park, MD 20742
USA

Office: 301-405-1427
Email: sdowney2 at umd.edu <mailto:sdowney2 at umd.edu>