Displaying 6 results from an estimated 6 matches for "wavcwttree".
2008 Dec 09
1
package "wmtsa": wavCWTPeaks error (PR#13381)
...n 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 (aa.cwt, type="minima")
aa.maxpeak <- wavCWTPeaks (aa.maxtree)
aa.minpeak <- wavCWTPeaks (aa.mintree)
Error in `row.names<-.data.frame`(`*tmp*`, value = c("1", "0")) :
invalid 'row.n...
2008 Dec 06
1
Morlet wavelet not supportd by wavCWTPeaks
..., 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...
2010 Apr 12
0
WMTSA "wavCWTPeaks" : Error in if (time.start < times.range[1]) time.start <- times.range[1]
...t:
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.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...
2008 Dec 12
0
Is there anyone in charge of package wmtsa ?
...<- 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 (Bad-Tree-Configuration) has only 4 branches. I believe this may cause the error I get when I try to g...
2008 Dec 05
2
Help with wavCWTPeaks
...332766 -0.266518818
[241] -0.314798550
I convert it into a time series and then I get the CWT coefficients. Then I build the tree
that exhibits only 3 branches (see attached plot)
aats <- ts (aa, deltat =1/30, start = 0.0)
aa.cwt <- wavCWT(aats)
aa.tree <- wavCWTTree (aa.cwt)
I can get the data for each of the 3 branches:
> aa.tree[[1]]
$itime
[1] 135 135 134 133 132 130 128 126 123 122 122 122 122 123 126
$iscale
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
$time
[1] 4.466667 4.466667 4.433333 4.400000 4.366667 4.300000 4.233333 4.166667
[9] 4....
2010 Mar 24
0
package "wmtsa": how to set up the Tolerance vector to detect time series peaks blindly
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 see attachments).
I am looking for some criteria, possibly based on the CWT results, to make this...