search

tgouhier / biwavelet

R package to conduct univariate and bivariate wavelet analyses
43 stars 17 forks source link

Significance level has no effect with cross-wavelet #85

Open JAFERRARI opened 3 years ago

JAFERRARI commented 3 years ago

Dear Sirs. The significance level works with wavelet coherence, but has no effect with cross-wavelet.

Thank you

See the example:

library(biwavelet)

t1 <- cbind(1:200, rnorm(200)) t2 <- cbind(1:200, rnorm(200))

cross Wavelet sig.level HAS NO EFFECT

biwv.t1t2 <- xwt(t1, t2,sig.level = 0.95, sig.test = 0) par(oma = c(0, 0, 0, 1), mar = c(5, 4, 4, 5) + 0.1) plot(biwv.t1t2, plot.cb = TRUE, plot.phase = TRUE)

biwv.t1t2 <- xwt(t1, t2,sig.level = 0.2, sig.test = 0) par(oma = c(0, 0, 0, 1), mar = c(5, 4, 4, 5) + 0.1) plot(biwv.t1t2, plot.cb = TRUE, plot.phase = TRUE)

wavelet coherence sig.level WORKS

biwv.t1t2 <- wtc(t1, t2,sig.level = 0.95, sig.test = 0,nrands = 100) par(oma = c(0, 0, 0, 1), mar = c(5, 4, 4, 5) + 0.1) plot(biwv.t1t2, plot.cb = TRUE, plot.phase = TRUE)

biwv.t1t2 <- wtc(t1, t2,sig.level = 0.2, sig.test = 0,nrands = 100) par(oma = c(0, 0, 0, 1), mar = c(5, 4, 4, 5) + 0.1) plot(biwv.t1t2, plot.cb = TRUE, plot.phase = TRUE)