similar to: Help with wavCWTPeaks

aa <- (structure(list(X.0.85 = c(-1.02, -1.17, -1.29, -1.39, -1.46, -1.5, -1.52, -1.5, -1.46, -1.39, -1.3, -1.19, -1.07, -0.93, -0.79, -0.65, -0.5, -0.36, -0.22, -0.08, 0.05, 0.18, 0.3, 0.41, 0.52, 0.62, 0.72, 0.81, 0.89, 0.98, 1.05, 1.13, 1.19, 1.25, 1.29, 1.31, 1.31, 1.29, 1.24, 1.16, 1.06, 0.93, 0.77, 0.58, 0.38, 0.16, -0.07, -0.31, -0.89, -1.05, -1.19, -1.31, -1.41, -1.47, -1.51, -1.51,

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,

Here is another occurrence of wmTSA internal error. My time series is a short breathing cycle (2425-Cyle_9.txt). Since wmtsa 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

I have attached the signal that causes the error message in this email subject. Only columns 1 and 3 have to be considered. It is the work trajectory of a molecule migrating between two equilibrium conformations. The curve 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

Dear All, I would like to ask a couple of questions on a LME model. I tested 4 selection lines at 4 food concentrations against a standard competitor stock. I had 3 replicate cages per selection line. In each cage I have 10 vials. I counted the number of wild type flies and competitor stock emerging in each vial. My main question is: is there any difference between selection lines? I did fit

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 =

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

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.3091308821262, -1.42473112931594, 0.234125312118165, -0.307047554490597,

Hi folks, I'm attempting to use the EMD package to analyze some neuroimaging data (timeseries with 64 channels sampled across 1 million time points within each of 20 people). I found that processing a single channel of data using EMD::emd() took about 8 hours. Exploration using Rprof() suggested that most of the compute time was spent in EMD::extrema(). Looking at the code for EMD:extrema(),

This is an example from chapter 11 of the 6th edition of Devore's engineering statistics text. It happens to be a balanced data set in two factors but the calculations will also work for unbalanced data. I create a factor called 'cell' from the text representation of the Variety level and the Density level using '/' as the separator character. The coefficients for the linear

I have a vector with about 100.000 values representing a quite regular function (sinusoid like). I would like to find all local maxima of this function (should be about 4000). Is there a native routine for R? Thanks in advance Eryk. -- _|_ \|/ \|/ Eryk Witold Wolski tel :0049-(0)30-8413-1543 w w ?v? 'v? \'v'/ MPI Moleculare Genetik fax :0049-(0)30-8413-1139 |

I have a data frame with two columns, a factor and a numeric. I want to create data frame with the factor, its frequency and the median of the numeric column > head(motifList) events score 1 aeijm -0.25000000 2 begjm -0.25000000 3 afgjm -0.25000000 4 afhjm -0.25000000 5 aeijm -0.25000000 6 aehjm 0.08333333 To get the frequency table of events: > motifTable <-

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

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

Hi! I have written an R-package (http://tocsy.agnld.uni-potsdam.de/wavelets/index.html) in R 2.4.1 that requires the ts() function. Users using R 2.5.0 now have a problem installing this package. I checked the package using R 2.5.0: _______________________________________________________ * Installing *source* package 'sowas' ... ** libs gcc -std=gnu99

Hello, # Full disclosure. I am not sure if my problem is a bug(s) in the code, or a fundamental misunderstanding on my part about what I am trying to do with these statistics. I am not familiar with maximum likelihood tests. # I currently have two vectors Aequipecten<-c(0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0,

Cc'd back to the list... always use reply-all. You say res should have been LSS but LSS is a scalar so the for loop will only run once. What does a successful output look like for a sample input? How do you (we) know when success has been achieved? In fact, what is the formula you want to implement? If your code below is the definition of your "formula" then we are left with no

Hi all, as I am trying to move slowly from just "working" to "good" code, I'd like to ask if there's a smarter way than using a for-loop in tasks like the example below. I need to obtain the extrema of the cumulated sum of a detrended time series. The following code is currently used, please have a look at the comments for my questions and remarks: system.time({ X

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!