similar to: Writing Complex Formulas

I am trying to create a set of wavelets in frequency space--namely Cauchy wavelets for an intensity analysis (von Tscharner, 2000). The wavelets are defined by the following formula: [(f/cf)^(cf*scale)]*[e^((-f/cf)+1)^(cf*scale)] where *f *is frequency of length *n*, *cf* is center frequency (defined below) and is an array of *j *columns and *n* row, and scale is a constant. cf = (1/scale)*(j

Full_Name: Jörg Polzehl Version: 1.5.1 OS: Windows 2000 Submission from: (NULL) (193.175.148.198) When calculating MLE's in a variance component model using constrained optimization, i.e. optim(...,method="L-BFGS-B",...) I observed an inproper behaviour in cases where the likelihood function was evalueted at the constraint. Parameters and value of the function at the constraint

Hi, I am trying to recreate the right graph on page 524 of Gelman's 2006 paper "Prior distributions for variance parameters in hierarchical models" in Bayesian Analysis, 3, 515-533. I am only interested, however, in recreating the portion of the graph for the overlain prior density for the half-Cauchy with scale 25 and not the posterior distribution. However, when I try:

Hi!! I'm looking for a package in R for time series analysis using wavelets( I'm interested on extractingt information from electrocardiograms with the aim of study the heart rate variability) so i want to know which package is the best for me. I've thinking about using some of those(you can recommend me others,xD): *wavelets: A package of funtions for computing wavelet filters,

I have question (curiosity) regarding returned values of R's qcauchy () function, for nonexceedance probability (F). It seems the ideal returned range of cauchy distribution should be [-Inf,Inf]. For F=0 > qcauchy(0) [1] -Inf but for F=1 > qcauchy(1) [1] 8.16562e+15 It seems to me that the proper return value should be Inf??? For default (location=0,scale=1) quantile function of

Hi I'm using R(2.9.0) and gstat package under Windows to plot sample variograms. When I want to fit them, I don't have the possibility to choose Cauchy functions. Under Mathematica for example, we have it and it seems to be the best one to fit the variograms I have I tested the others functions, they can fit but I would,if possible, use the Cauchy one. If anybody knows something about

Hi all, I am using the function COV.WT to estimate the estimators (location and scale) of a bivariate cauchy distribution. My doubt is about the option WT (weight), cause at the R-help shows that the weight is uniform according to the number of observations. But, checking the theory, for example, the mean is given by mean_estimator=mean(u(s)x)/mean(u(s)), where x=my data (bivariate)

Can anyone explain such different output: > stableFit(s,alpha = 1.75, beta = 0, gamma = 1, delta = 0, + type = c("q", "mle"), doplot = TRUE, trace = FALSE, title = NULL, + description = NULL) Title: Stable Parameter Estimation Call: .qStableFit(x = x, doplot = doplot, title = title, description = description) Model: Student-t Distribution Estimated

I'm new to the topic of wavelets. When I tried to use the mra function in the wavelets package, the data is not getting compressed. eg. if the original data has 500 values , the output data also has the same. However in MATLAB, depending on the level of decompositon, the data gets compressed. How do I implement this in R? -- View this message in context:

Hello, In my last post, I seemed to say that there were no directional wavelet transforms that would give improvement to the NHW codec, but are directional wavelets an active research field today? I have also seen some recent patents on advanced and improved SPIHT, and I am a little surprised that there is no "up-to-date" state-of-the-art wavelet codec.Rududu in April 2008 was the last

Following a post on vqf.com forum, I found out that some company called Adapted Wave Technology is ready to release its audio codec based on wavelets. They say that this wavelet-based codec is superior to MDCT codecs by 30-50%. One sentence caught my attention: "Adapted Wave Technologies holds exclusive patent rights for applying wavelet packets to audio applications." - sounds like a

waveslim 1.4 has recently been uploaded to CRAN and is fully compatible with Rv2.0. 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

waveslim 1.4 has recently been uploaded to CRAN and is fully compatible with Rv2.0. 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

I found this a few months ago, but for the life of me I can't remember what the function or package was, and I have had no luck finding it this week. I have found, again, the functions for working with distributions like Cauchy, F, normal, &c., and ks.test, but I have not found the functions for estimating the distribution parameters given a vector of values. What I need to do is

Todd Day wrote: > How many players (players that can pull an MP3 stream) are using > Vorbis? I have used it a couple of times for testing and was > impressed by the results. Does Winamp have Vorbis decode built-in? Winamp just put the finishing touches on their own vorbis plugin, and we expect it to ship in the next version. Sonique has already had Vorbis for some time (millions of

Is there a wavelet toolbox in R available that can be compared to the Matlab toolbox for wavelets. Does this toolbox offer the possibility to calculated approximations with different wavelets such as Daubechies 2, Daubechies 3, etc, Symmlet 3, Symmlets 5, etc, Coiflets and so on. I have tried to proceed in this way but I was only able to use one sort of wavelets with a pre-fixed order (i.e.

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

I am attempting to run a Poisson EWMA model using Patrick Brandt's source code. I get the following error when I run the code: Error in dimnames(x) <- dn : length of 'dimnames' [1] not equal to array extent Dimnames(x) looks like this: [[1]] NULL [[2]] [1] "mip" "div" "nom" "unity" "mood"

Hi user's Does anybody work with wavelets on R? Please I need some help.....!!!! Atte Marize Sim?es -- View this message in context: http://r.789695.n4.nabble.com/Package-wavelets-tp2526023p2526023.html Sent from the R help mailing list archive at Nabble.com.

hi; i'm a lurker & trying to make sense of the discussion... is there a glossary of terms available somewhere? mdct, wavelets, 4th dimensions & preechoes remind me of my time in the twilight zone;) thanks peter --- >8 ---- List archives: http://www.xiph.org/archives/ Ogg project homepage: http://www.xiph.org/ogg/ To unsubscribe from this list, send a message to