search for: cwt

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

...a functions that extract extrema seem to expect longer series than I have, I tried the folowing two tricks: 1) I prolong the 1-cycle series on both ends through duplicating the first value (on the left end( and the last value on the right end as any ties as the 1-cycle series length. Then I apply CWT to such prolonged series and then 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...

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

Hello, Listers I'm trying to run blloean logit model with R. My code is: > library(boolean) > library(foreign) > pr <- read.dta ("prcore1.dta") > bp <- boolprep ("(a&b)|c", "cwt", a="O1", b="t", c="DM2 > + ah + md + con + n3 + rel + slo + pyrs > + sp1 + sp2 + spl3") > answer <- boolean (bp, link = "logit", method = "nlm") But I got an error message here as follows: Error in eval(expr, envir, enclos) : obj...

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

Dear, I am using "cwt "function from Rwave package to perform Morlet wavelet analysis. d1<-c(1.31673209591515, -0.171333455797346, -1.67618546079420, -0.931604651010362, -0.954614183224057, -1.19821794547758, 0.516096363353144, -0.0432811977638559, 0.737764943619919, 0.438046629673177, -0.208607167908743, -0...

...this will be an adequate example as I had the code already created/data uploaded). However, I have other plots that have fewer space-time strata (fewer tiles) and the effect is very noticeable. Code: ##You will need to load the following packages: latticeExtra, gridExtra and RCurl #Upload data cwtb_csv<-getURL("https://docs.google.com/spreadsheet/pub?key=0AjzYZNH9Dw9qd Fc5c1FXZ19uYXE3U1QwU1MxVkR4dGc&single=true&gid=0&output=csv",ssl.verifypeer = FALSE) cwtb<-read.csv(textConnection(cwtb_csv),header=T) gsi_csv<-getURL("https://docs.google.com/spreadsheet/pu...

...text html latex createar text html latex example createwgn text html latex example criticalvaluesWCO text html latex example criticalvaluesWSP text html latex example cwt.ts text html latex example nao text html latex nino3 text html latex plot.wt text html latex example plotwt text html latex...

...Besides ensuring usability with the most recent version of R, two additional "flavors" of wavelet methodology have been added to the package: (1) Hilbert wavelet pairs and (2) the dual-tree complex wavelet transform [only 1D and 2D ported from Matlab code by Selesnick]. The dual-tree CWT code has not been optimized for R (only a "direct" port from Matlab) so reasonable speed-ups are anticipated in the future. Comments/bug reports are welcome. References for (1) include: Selesnick, I.W. (2001). Hilbert transform pairs of wavelet bases. IEEE Signal Processing Letters...

I am looking for some explanations about the usage of the poorely documented R paclkage Rwave. Has anyone ever tried out its functions for Wavelet Analysis ? Thank you so much. Maura 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 [[alternative HTML version deleted]]

...Besides ensuring usability with the most recent version of R, two additional "flavors" of wavelet methodology have been added to the package: (1) Hilbert wavelet pairs and (2) the dual-tree complex wavelet transform [only 1D and 2D ported from Matlab code by Selesnick]. The dual-tree CWT code has not been optimized for R (only a "direct" port from Matlab) so reasonable speed-ups are anticipated in the future. Comments/bug reports are welcome. References for (1) include: Selesnick, I.W. (2001). Hilbert transform pairs of wavelet bases. IEEE Signal Processing Letters...

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

Ok, here's a weird one. I've attached a spa2000 to asterisk, and got the two phones to register as exten 706 and 707. I can call exten 708 (a cisco 7940) from 707 and everything works fine. I can call exten 708 from 706 and everything works fine. When I make a call to either 706 or 707 from any phone, the phone attached to the spa does not ring. However, if I pick up the appropriate

Hello, can anybody tell me how to produce a plot like the one in http://cran.r-project.org/web/packages/lme4/vignettes/Implementation.pdf on page 13, Figure 6? The data is stored in: library(nlme) data(Oats) Cheers -- View this message in context: http://r.789695.n4.nabble.com/Help-needed-in-reproducing-a-plot-tp4119603p4119603.html Sent from the R help mailing list archive at

Anyone have an example XML file for the PAP2T? Cheers, j

...it Cc: r-help@stat.math.ethz.ch Message-ID:     <c502a9e10905271556j1f96b9e1j86898964ae082e36@mail.gmail.com> Content-Type: text/plain; charset=UTF-8 Do you need time localization, or are you only interested in the period of the high frequency?  If you do need time localization why not use a CWT to look at the signal? Stephen On Wed, May 27, 2009 at 4:59 PM,  <mauede@alice.it> wrote: > Well, the time series I am? dealing with are non-linear and not-stationary. > Maura > > -----Messaggio originale----- > Da: r-help-bounces@r-project.org per conto di stephen sefick &gt...