similar to: How-To extract data from a matrix of lists subject to a set of given constraints

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,

I am using some functions from package clusterSim to evaluate the best clusters layout. Here is the features vector I am using to cluater 12 signals: > alpha.vec [1] 0.8540039 0.8558350 0.8006592 0.8066406 0.8322754 0.8991699 0.8212891 [8] 0.8815918 0.9050293 0.9174194 0.8613281 0.8425293 In the following I pasted an excerpt of my program:

Hello, I am trying to call a FORTRAN subroutine from R. The Fortran code is @: http://lib.stat.cmu.edu/apstat/206 It performs a bivariate isotonic regression on a rectangular grid (m X n) matrix. I used the g77 compiler and successfully created a dll file and it also loads successfully from R. But somehow the programs fails to run properly. (I do get the correct result when I compile the

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

I have written a short script to estimate the size of a test of non-constant mean in an AR(1) time series. When I run the simulation on my PC (R version 2.7.1), I get the expected result: an empirical size much larger than the nominal size. On the Red Hat machine (R version 2.7.2) in my department, however, I get p-values > .45 for every simulated test. Below is my simulation code (uses

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)

Hi There: I was wondering if there is a way to read FORTRAN list data (similar to IDL's readf function). I often use FORTRAN for most of my number crunching, and use something like IDL to visualize and perform statistical analysis on that data. Since the each file is rather large (>100 Mb), formatting the output into columns or tables is impractical, hence the "list"

Dear all, I'm working on a design with rotating panels. Each site is sampled every X year. Have a look at the code below the get an idea of the design. In reality the number of sites will be (much) higher. However I'm not sure if I coded the design correctly in svydesign(). Can someone confirm that it is either correct or wrong? I've added the lmer() equivalent that I want to

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

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

hi there ! the current version of tinc crashes with a segfault. it seems, as if has to to with fact, that i don`t have configured any subnets.... if the client tries to start up subnet_update() ist called for myself. the for iteration does not get an entry (because the subnet avl-tree is empty). the later call of net2str causes a segfault, because it gets an empty subnet and tries to access

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'

Dear R Gurus, I am currently working on the two dataset ( A and B), they both have the same fields: ID , REGION, OFFICE, CSTART, CEND, NCYCLE, STATUS and CB. I want to merge the two data set by ID. The problem I have is that the in data A, the ID's are unique. However in the data set B, the ID's are not unique, thus some repeat themselves. How do I the merge or retrieve the

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