similar to: Wavelet Bootstrap Size Simulation

Given a time series of length N, I am trying to figure out its maximum DWT decomposition level. With reference to "dwt" function of R package "wavelets", running the provided example I get the following and wonder how the maximum decomposition level (which is not an integer number) is calculated. I tried myself the formula in the on-line documentation but could not get the

My goal is to store the DWT coefficients from a number of time series in such a way to be able to extract subsets satisfying giben conditions. I have defined the following 11 matrices whose x-axis represent time intervals and y-axis represent the number of time series: MS <- 11 d1.mat <- matrix (data=list(), nrow=TotNumCycles, ncol=2^(MS-1)) d2.mat <- matrix (data=list(),

This code consistenly segfaults for me. Can someone please take a look and tell me if the problem is due to something I am doing or is there a problems with the dwt (idwt) functions in the waveslim library. Thanks tom library(waveslim)

I take the liberty of pasting the on-line documentation for package "wavelets", functions "dwt" and "wt.filter", pertinent to the wavelet/filter choice. While I apologize for my verbose messages with huge attachments, I encourage whoever loads packages to CRAN to take some minutes to reread the attached documentation. Usage dwt(X, filter="la8", n.levels,

hi everyone . I am trying to use some packages but there are some functions have the same name in different package. for example dwt function both in packages wavelets and waveslim dwt(X, filter="la8", n.levels, boundary="periodic") How can I avoid mixing them up ? -- TANG Jie Email: totangjie@gmail.com Tel: 0086-2154896104 Shanghai Typhoon Institute,China

Can someone please help me understand the examples given in the dwt function of the waveslim library library(waveslim) ?dwt I'm having a problem understanding these lines. I assume they are required because of how the algorythm deals with the signal bounderies. Am I correct in thinking that for wavelets of scale 1 and 2, the coefiecints need to be shifted by 2 places, for scale 3 and 4 shift

Hi I have the following loop and would like to vectorize it. Any ideas if it is possible? Thanks, Rainer Tha Loop: for (i in 2:Result$NoSims) { ppp <- runifpoint(Result$NoPlants) K <- Kest(ppp) Result$LSim[i,] <- sqrt(K$iso / pi) - K$r CM <- (Result$LSim[i,] * Result$LSim[i,]) / abs(K$r[2] - K$r[1]) Result$SigCM[i] <- sum(CM, na.rm=TRUE) print(i) flush.console() }

Hi I have a loop which is doing time consuming calculations and I would like to be able to have some feedback on where it is in it's calculations. I tried to simply show the counter variable in the loop, but id doesn't work as all display seems to be delayed until the loop is completed. Is there any way of displaying the progress of a loop? Rainer The loop: for (i in

Dear all, Using R function "dwt", it seems that I cannot specify the wavelet transformation like Haar. What's the R code for wavelet decomposition which allows me to specify Haar wavelet transformation? Of course, if it can include "db2", that is even better. In general, I want an R function like matlab code "dwt". Thanks in advance! Zhen Li

Hello, I forgot in my last reply that my DWT implementation can be speed up, for example I'm doing for now: for (;_X1<_E_;_X1++,_RES+=2) //dilatation { _RES[0]=_X1[0]<<3; _RES[1]=(_X1[1]+_X1[0])<<2; } then for (;_X2<_E_;_X2++,_RES+=2) //details { _RES[0]-=(_X2[1]+_X2[0])<<1; _RES[1]+=6*_X2[1]-_X2[2]-_X2[0];

Hello, I have upgraded from Mandrake Linux 8.0 to Mandrake 8.1 and try to reinstall my favourite R... Everything is OK for the base software but I have trouble to get some packages installed. Specifically, for some packages, a "collect2" binary seems to be necessary during the compilation/linking process of the library, and LD complains about not finding it : Installing source

Using ordered probit model, I get errors from dwt and bptest. dwt: Error in durbinWatsonTest.default(...) : requires vector of residuals bptest: Error in storage.mode(y) <- "double" : invalid to change the storage mode of a factor I imagine I have to restate as an individual probit model for each category, but is there an easier way? thanks, bp [[alternative HTML version

Dear WizaRds, unfortunately, I have been unable to replace the '99' and '999' entries in library(UsingR) attach(babies) as definitions for missing values NA, because sometimes the 99 entry is indeed a correct value. Usually, or so I thought, NAs can easily replace a, say, 999 entry via mymat[mymat==999] <- "yodl" in a matrix or data frame. Alas, the babies'

After playing with the vorbis code for a while and doing tons of hacks and analysis on it, I've found it to perform very poorly with impulse signals. The MDCT seems to cause lots of spreading, and it seems to result in much worse impulse performance then mp3. What is the current plan on handling this? Will a smart quantizer be able to avoid it? I've been looking at various ways of

Dear R Users, For running wavelet functions using dwt( ), modwt( ), and mra( ), a wavelet filter algorithm is applied. For all these functions, default is "la8" and other possibility is "haar". In related documents, another possibilities like as symlet and coiflet ... are not cited. Besides "la8" and "haar", which wavelet filters can be used? Thank

Hello I do wavelet transform by using this code: dec=dwt(ld, filter='d8', n.levels=lev, boundary="reflaction"); dec consists of the decomposition coefficients and other How can I change the coefficients the decomposition manualy? -- View this message in context: http://r.789695.n4.nabble.com/wavlet-transform-tp2313427p2313427.html Sent from the R help mailing list archive at

Hi, I am interested in a applying wavelets as a smoothing tool for my (1-dimensional) data. I looked into wavetresh and waveslim packages but could not quite figure out an obvious way to do this after running dwt or wt functions. Would someone be able to point me in the right direction on how to denoise univariate data using one of wavelet packages available in R? Thank you very much Jane ps

Hello, I searched through the mailing list archives here, at securityfocus.com, at HP and google and have come up dry. Sorry in advance if my search was not complete enough or I missed something but I sure as heck didn't see anything. ENV: Slackware 10 w/ ssh (SSH-2.0-OpenSSH_3.8.1p1) connecting to an alpha with OpenVMS 7.3-2 w/ sshd (Remote protocol version 2.0, remote software version

> results <- mra(x$w1mcp, filter = "d4", n.levels = 3, boundary = "periodic", method = "dwt") > write.csv(results, "c:/mydata.csv") Error in as.data.frame.default(x[[i]], optional = TRUE) : cannot coerce class 'structure("mra", package = "wavelets")' into a data.frame I am not able to access data stored in results.

Hot Pennies.net Special Situation Alert Watch This Stock Trade!!! STOCK DETAILS Symbol: IPFS Exchange: OTCBB Shares Outstanding: 13,379,000 Current Price: $1.10 Industry: Wireless Manufacturing and Technology COMPANY DETAILS Integrated Performance Systems Inc 17300 N. Dallas Parkway, Ste. 2040 Dallas, Texas 75248 Phone: 972/381-1212 Web: