aa <- (structure(list(X.0.85 = c(-1.02, -1.17, -1.29, -1.39, -1.46, -1.5, -1.52, -1.5, -1.46, -1.39, -1.3, -1.19, -1.07, -0.93, -0.79, -0.65, -0.5, -0.36, -0.22, -0.08, 0.05, 0.18, 0.3, 0.41, 0.52, 0.62, 0.72, 0.81, 0.89, 0.98, 1.05, 1.13, 1.19, 1.25, 1.29, 1.31, 1.31, 1.29, 1.24, 1.16, 1.06, 0.93, 0.77, 0.58, 0.38, 0.16, -0.07, -0.31, -0.89, -1.05, -1.19, -1.31, -1.41, -1.47, -1.51, -1.51, -1.49, -1.44, -1.37, -1.28, -1.17, -1.04, -0.91, -0.76, -0.62, -0.47, -0.33, -0.19, -0.06, 0.08, 0.2, 0.32, 0.43, 0.54, 0.64, 0.73, 0.82, 0.91, 0.99, 1.07, 1.14, 1.2, 1.25, 1.29, 1.31, 1.31, 1.28, 1.23, 1.15, 1.04, 0.9, 0.73, 0.55, 0.34, 0.11, -0.12, -0.92, -1.08, -1.22, -1.33, -1.42, -1.48, -1.51, -1.51, -1.48, -1.43, -1.35, -1.26, -1.14, -1.01, -0.88, -0.74, -0.59, -0.45, -0.3, -0.16, -0.03, 0.1, 0.22, 0.34, 0.45, 0.56, 0.66, 0.75, 0.84, 0.93, 1.01, 1.08, 1.15, 1.21, 1.26, 1.3, 1.31, 1.3, 1.27, 1.21, 1.13, 1.01, 0.87, 0.7, 0.5, 0.29, 0.07, -0.17, -0.95, -1.11, -1.24, -1.35, -1.44, -1.49, -1.51, -1.51, -1.48, -1.42, -1.33, -1.23, -1.12, -0.99, -0.85, -0.71, -0.56, -0.42, -0.27, -0.14, 0, 0.13, 0.25, 0.36, 0.47, 0.58, 0.68, 0.77, 0.86, 0.94, 1.02, 1.1, 1.17, 1.23, 1.27, 1.3, 1.31, 1.3, 1.26, 1.2, 1.11, 0.98, 0.83, 0.66, 0.46, 0.25, 0.02, -0.22, -0.99, -1.14, -1.27, -1.37, -1.45, -1.5, -1.52, -1.5, -1.47, -1.4, -1.32, -1.21, -1.09, -0.96, -0.82, -0.68, -0.53, -0.39, -0.25, -0.11, 0.02, 0.15, 0.27, 0.39, 0.5, 0.6, 0.7, 0.79, 0.87, 0.96, 1.04, 1.11, 1.18, 1.24, 1.28, 1.31, 1.31, 1.29, 1.25, 1.18, 1.08, 0.96, 0.8, 0.62, 0.42, 0.2, -0.03, -0.27)), .Names = "X.0.85", class = "data.frame", row.names = c(NA, -240L))) library(wmtsa) library(fields) aats <- ts (aa, deltat =1/30, start = 0.0) aa.cwt <- wavCWT(aats) aa.tree <- wavCWTTree (aa.cwt) aa.peak <- wavCWTPeaks(aa.tree) sorry I didn't realize that only the mexican hat wavelet could be used for the peak function. This should work fine. and Use the tim.colors argument in the plot() function. I hope this helps.... This reason I suggested the morlet wavelet is that it is compactly suppoted
On Sat, Dec 6, 2008 at 1:39 AM, <mauede at alice.it> wrote:
aa.cwt <- wavCWT(aats,wavelet="morlet")
aa.cwt
Continuous Wavelet Transform of aats ------------------------------------ Wavelet : Morlet Wavelet frequency shift : 5 (rad/sec) Length of series : 241 Sampling interval : 0.03333333 Number of scales : 65 Range of scales : 0.0333333333333333 to 8.03333333333333
aa.cwt.coef <- as.matrix (aa.cwt) plot (aa.cwt, series=TRUE) plot (aa.cwt, type="persp")
aa.tree <- wavCWTTree (aa.cwt) aa.peak <- wavCWTPeaks (aa.tree)
Error in wavCWTPeaks(aa.tree) : Only CWT developed using the Mexican hat (gaussian2) filter are supported
Alice Messenger ;-) chatti anche con gli amici di Windows Live Messenger e tutti i telefonini TIM! Vai su http://maileservizi.alice.it/alice_messenger/index.html?pmk=footer
Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis