similar to: Harmonic Regression in R

Hi: I have encountered problems with imodwt and universal.thresh.modwt and cannot find any reference in R Search. Hope someone can give me some ideas: Starting with modwt.la8 <- modwt(xdata, "la8", n.level=6) <-- this seems to work fine (1) ydata <- imodwt(modwt.la8) will always give ydata as numeric(0) (no values) instead of being a time series data with

Hello list; I am working with a dataset containing animal locations over time for a large number of individuals. I would like to compute a center of activity by finding the minimum harmonic mean on a grid overlaying the points of an animal (a "standard" to express center of activity in animal ecology). I have searched the archives at Jonathan Barons server but have not been able

An idea would be like for WMA 9 lower bitrates (32-42-48Kbps) to use a 'crystallizer'; which is basically a harmonics booster focussed at transposing sharp tones some octaves higher. Eg: A file has been recorded @ 20khz computer (or 10khz real life) to preserve space. While playing back the file sounds a bit mushy, almost as if someone was speaking through a cardboard wall. The higher

I' need to smooth a serie univariate with harmonic regression, but i do not know if it exists somepeople know in R this type of regression, not in S-plus. Thank's _________________________________________________________

I am looking for a package to perform harmonic analysis with the goal of estimating the period of the dominant high frequency component in some mono-channel signals. I guess there are presumably a number of CRAN packages allowing for such analysis. However, my search with keywords was not successfull. It brought up a lot of Fourier miscellanea but nothing specifically geared for my needs. I would

Dear All, I have some time series data where X=month and Y=nutrient concentration (I can have several concentration data for one month). Is there a way to fit for it an Harmonic Function. Is there a package, script,etc which I can use? Thx -- View this message in context: http://r.789695.n4.nabble.com/Fitting-a-Harmonic-Function-to-Time-Series-Data-tp3901266p3901266.html Sent from the R help

Dear R Users, I have a CO2 time series. I want to fit this series seasonal cycle and trend with fourth harmonic function, and then compute residuals. I am doing something like: file<-read.csv("co2data.csv") names(file) attach(file) fit<-lm(co2~1+time+I(time^2)+sin(2*pi*time)+cos(2*pi*time)+sin(4*pi*time)+cos(4*pi*time)+

Hello, I will try to explain the problem, sorry if it will be a little long... I'm using R to analyze results of cyclic mechanical testing, like this: - apply quasi-sinusoidal load - measure quasi-sinusoidal vertical and horizontal deformations (quasi-sinusoidal load means that load "should be" sinusoidal, but testing machine puts in some noise...) I enclose a sample of data at

Hello, I try to synthesize harmonic functions with a subset of frequencies determined by fast fourier transform. I wrote two different functions, a looped version and a matrix multiplication version. As an example the looped version takes 5 sec and the matrix-version takes 3-4 sec (R 1.4, Athlon 1.2 GHz, 256 MB, Win ME), but the latter needs huge amount of memory due to the matrices. So, none of

Hi, Sorry for the very short explanation about the problem of intersection. I have a wave function monitored from the heart beat in a particular interval of times.?Apart?fom that, there is?a line with positive slope (e.g: y = x+2) which lies across the wave and intersect on a number of points. My problem is i have no exact equation for?such a?complex harmonic wave produced by the heart

Hi Has anybody done THD or THD-N measurements with the CELT Codec (bext would be on various bit rates) If someone could share results for Mono at 64kBit and Stereo at 128 and 196kBit it would be great. thank you very much Jochen -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.xiph.org/pipermail/opus/attachments/20100817/b1da17cc/attachment-0002.htm

Hi, I need to derive a time series that represents the first derivative of an original time series. The function coefDeriv in the cyclones package seemed to be the ticket, but I'm not sure if I am interpreting the output of the function correctly...or even using the function correctly. This is a snipbit of what I've been trying: -------- library(cyclones) ## read in my 1-column of

Airon, I used R2.4.0 on a Windows XP (SP2) (not Chinese :-)) and it still works: > data = read.csv("u:/airon.csv") > xdata = data$Adj..Close > modwt.la8 = modwt(xdata, "la8", n.level=6) > summary(modwt.la8) Length Class Mode d1 1467 -none- numeric d2 1467 -none- numeric d3 1467 -none- numeric d4 1467 -none- numeric d5 1467 -none- numeric d6 1467

I upgraded from Fedora Core 4 to Fedora Core 5 and I find a lot of previously installed packages won't run because shared libraries or other system things have changed "out from under" the installed R libraries. I do not know for sure if the R version now from Fedora-Extras (2.2.1) is exactly the same one I was using in FC4. I see problems in many packages. Example, Hmisc: unable

Airon, I don't think you have to find an English computer 'cause the following must work in your Chinese one :-) Let me explain. First of all, change your lines to xdata <- ckhdat$Adj..Close[1:1447] #names(ckhdwt.la8) <- c("w1", "w2", "w3", "w4", "w5","w6", "v6") note the # sign, i.e., DO NOT change the names

Is there some reason why > (1:2)+(1:3) [1] 2 4 4 Warning message: In (1:2) + (1:3) : longer object length is not a multiple of shorter object length can't be made into an error? I realise it was there in S-PLUS, but since it produces a warning there can't be many examples on CRAN or Bioconductor using it, and I can't think of any situation where it would be used deliberately.

Skipped content of type multipart/alternative-------------- next part -----= --------- A non-text attachment was scrubbed... Name: not available Type: image/gif Size: 18867 bytes Desc: not available Url : http://lists.linux-foundation.org/pipermail/virtualization/attachment= s/20070408/f3caf0e4/attachment.gif

Hi, I was wondering if it is possible to harmonize the ouput of symbols() and legend() both from the graphics package. Let us take this example: x<-runif(10) y<-runif(10) z<-runif(10) leg<-round(seq(min(z),max(z),l=4),2) # 4 values rounded up to 2 decimals for the legend symbols(x,y,circles=z,inches=0.2) legend("topright",legend=leg,pch=1,pt.cex=leg/max(leg)*2) #

Skipped content of type multipart/alternative-------------- next part -----= --------- A non-text attachment was scrubbed... Name: not available Type: image/gif Size: 18867 bytes Desc: not available Url : http://lists.linux-foundation.org/pipermail/virtualization/attachment= s/20070408/f3caf0e4/attachment.gif

I just finished testing and integrating Speex into our test project. Speex works excellent and sounds great for the bitrate. I have noticed, however, on cheaper microphones that Speex will produce loud static noise at high harmonics/amplitude in the sound. Is there some setting to prevent this static? I don't have this problem with various other manufacturer voice coders that I am using