search for: wavcwt

Displaying 10 results from an estimated 10 matches for "wavcwt".

2008 May 16
0
wmtsa (wavCWT)
I can send data along, but it is too large for the list: I am using wmtsa and the function wavCWT I would like to use the range.scale to only evaluate a subset of frequencies. deltat of the time series that I am trying to analyze is 1 d = wavCWT(RM202.ts, range.scale = c(96, 2688)) and I get the error message Error in wavCWT(RM202.ts, scale.range = c(96, 2688) : Minimum scale must be great...
2008 May 28
1
Writing intermediate analysis to disk
Is there a way to write and analysis to disk and then reconstruct the whole thing back into an object. wavCWT() #wmtsa package I am running out of memory on my computer and I was wondering if there was a way to iterate through this process (as it is an iterative process anyway- it just stores the whole thing to memory). Or is there a way to set the scale that I want to look at so that wavCWT can use some...
2008 Dec 09
1
package "wmtsa": wavCWTPeaks error (PR#13381)
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 <- wavCWT...
2008 Dec 06
1
Morlet wavelet not supportd by wavCWTPeaks
....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 sugges...
2008 Sep 11
3
periodicity validation
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
2008 Dec 12
0
Is there anyone in charge of package wmtsa ?
...n obtain the maxima trees and the peaks as follows: aal <- rep (aa[1], length(aa)) aar <- rep (aa[length(aa)], length(aa)) aa3 <- c(aal,aa,aar) aa3ts <- create.signalSeries(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...
2008 Dec 05
2
Help with wavCWTPeaks
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 -...
2010 Apr 12
0
WMTSA "wavCWTPeaks" : Error in if (time.start < times.range[1]) time.start <- times.range[1]
...has 2 peaks, as shown in its plot. But I keep missing the 2nd one. Here is my short script: library(wmtsa) setwd("C:/Documents and Settings/Monville/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....
2008 Sep 09
1
help on wavelet
Hi, I have little experience using wavelet and I would like to know if it is possible,using R wavelet package, to have a plot of frequency versus time. thank you giov -- View this message in context: http://www.nabble.com/help-on-wavelet-tp19395583p19395583.html Sent from the R help mailing list archive at Nabble.com.
2008 Sep 22
1
R-help Digest, Vol 67, Issue 23
...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 y...