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