search for: wavcwtpeaks

Full_Name: Maura Monville Version: 2.8 OS: Mac OS/X 10.5 Submission from: (NULL) (87.4.122.234) Here is the code that causes wavCWTPeaks error aats <- create.signalSeries(aa, pos=list(from=0.0, by=0.033)) aa.cwt <- wavCWT(aats) x11 (width=10,height=12) plot (aats,main=paste(insig," Cycle: ",j,sep="")) aa.maxtree <- wavCWTTree (aa.cwt, type="maxima") aa.mintree <- wavCWTTree...

...le/Alanine Dipeptide/Work_curves-Dir") work <- as.ts(read.table ("calc_work_332")[,c(1,3)]) plot(work[,1],work[,2]) work.wt <- wavCWT(work, n.scale=1000) scales <- attr(work.wt, "scale") work.tree <- wavCWTTree(work.wt,tolerance=0.5/sqrt(scales)) work.peak <- wavCWTPeaks(work.tree) > summary(work.tree) End Time Length Octaves Min Max Mean SD Var MAD 1 1043 79 11.78 -6.46e-09 1.96e+06 3.58e+05 6.24e+05 3.89e+11 0.102 2 NA 17 1.95 6.30e-01 1.26e+01 8.09e+00 3.13e+00 9.81e+00 3.249 > work.peak <- wavCWTPe...

...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, D...

I keep getting the following error when I look for minima in the series: > aa.peak <- wavCWTPeaks (aa.tree) Error in `row.names<-.data.frame`(`*tmp*`, value = c("1", "0")) : invalid 'row.names' length How can I work it around ? Thank you. Regards, Maura Alice Messenger ;-) chatti anche con gli amici di Windows Live Messenger e tutti i telefonini TIM! Vai su...

...nalSeries(aa3, pos=list(from=0.0, by=0.033)) #CONVERT AMPLITUDE INTO TIME-SERIES aa3.cwt <- wavCWT(aa3ts) #CWT aa3.maxtree <- wavCWTTree (aa3.cwt, type="maxima") # GENERATE CWT MAXIMA TREE aa3.maxpeak <- wavCWTPeaks (aa3.maxtree) #GET MAXIMUM PEAKS Notice that in this case the generated tree (Bad-Tree-Configuration) has only 4 branches. I believe this may cause the error I get when I try to get the maxime: > aa3.maxpeak <- wavCWTPeaks (aa3.maxtree) Error in `row.name...

I cannot understand the following error printed out when I try to get the extrema of my time series. I would appreciate some suggestion as I really cannot interpret the error. I might not be using a proper set of parameters in calling such functions. I am learning by doing ... > aa.peak <- wavCWTPeaks (aa.tree) Error in `row.names<-.data.frame`(`*tmp*`, value = c("1", "0")) : invalid 'row.names' length My data vector is as follows: > aa [1] -0.852625404 -0.886941142 -0.920699269 -0.953858240 -0.986377090 [6] -1.018215725 -1.049335086 -1.079697274 -1.109...

In my previous experience with wmtsa function wavCWTPeaks I realized the Tolerance vector passed to the multi-scale decomposition function "wavCWTTree" is crucial for detecting the time series extrema and prevent the R script from aborting. By trial & error I recently managed to figure out the Tolerance vector values in a single case (please...

There is a series of data contains time in fixed step and energy varying with time, how to test its periodicity?In R, it seems there is no direct tools since I have search the R manual with periodic and I have not found any related topic. Thanks a lot

...f the time series) so to get to your time units (which you set in the deltat or frequency argument when you create a timeseries with ts() ) 2^(value of the scale). I hope this helps Stephen 2008/9/21 yuankun shi <shiyuankun.debian@gmail.com>: > Thanks, I have succeeded to do this, first wavCWTPeaks to get every peaks' > coordinate, then calculated their horizontal distance, finally,bkde output > the distance's distribution, that's what I want. > On the contrary, picture of wavCWT seems hard to understand, I am not sure > what the y axis and the color mean. Could you do...