1 2 3 4 5 6 7 8 9 library(waveletComp) my.w = analyze.wavelet(mydata, "Index", loess.span = 0, dt = 0.25, dj = 1/250, lowerPeriod = 0.25, upperPeriod = 32, make.pval = TRUE, n.sim = 30) wt.image(my.w, n.levels = 250, legend.params = list(lab = "wavelet power levels") )