similar to: Lyapunov exponent code for time series

Hello, Does anybody know if there is somewhere in R a function to calculate the Lyapunov exponent in a time series? Thanks, Philippe Grosjean .......................................................<??}))><.... ) ) ) ) ) ( ( ( ( ( Prof. Philippe Grosjean \ ___ ) \/ECO\ ( Numerical Ecology of Aquatic Systems /\___/ ) Mons-Hainaut University, Pentagone / ___ /( 8, Av. du

Dear all, I need to solve the following Lyapunov Matrix equation: C=ACA' + B, with A and B given square symmetric matrices. Does anyone knows of a package that can solve the lyapunov matrix equation in R? Or even a C/Fortran implementation? I did not find one on netlib. Thank you.

Has anyone writen R code for estimating Hurst exponent with R/S method or other methods? or any other source of R code available? Many thanks Catherine Wang

Hello, I am trying to collect several global measures or statistics for time-series as well as packages of R that can compute them. I have found several of them in papers and books, but the literature is so big i am sure i am missing several of them. skewness kurtosis min max mean SD trend seasonality periodicity chaos (Lyapunov Exponent) / Largest Lyapunov Exponent (i think is the same

Dear Sir, I’am a frensh student and i’am a new user of the R software. After using the command (x<-read.delim(“clipboard”) to read a spreadsheet of Excel, I want to run the bds test and calculate the Lyapunov exponent. I have charged the R software by the packages tseries and tseriesChaos. when i run bds.test(x,m=2) Unfortunately the R software displays “error in as.vector(x,mode= “double”) :

R-users -- I would be interested in tools in R to fit the following gam model: logit(p) = a + f(x1) + f(x2) + f(x1,x2), where f(x1,x2) defines a surface. I have looked into the mgcv library, but it seems only to fit models of the form: logit(p) = a + f(x1) + f(x2) Any ideas? Cheers, Dan =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dan Powers Associate Professor,

Dear All, I am a novice when it comes to time-series analysis and at the moment I am actually interested in calculating the Hurst exponent of a time series. This question has already been asked quite some time ago http://bit.ly/98dZsi and I trust some progress has been made ever since. I was able to find some functions in the packages http://cran.r-project.org/web/packages/Rwave/index.html

> From: llvmdev-bounces at cs.uiuc.edu [mailto:llvmdev-bounces at cs.uiuc.edu] > On Behalf Of Mikhail Lyapunov > Subject: [LLVMdev] Multiple modules JITting > We're looking for a way to use LLVM to JIT many modules, assuming that > the full list of modules and their content are not available at a time > when some of jitted pieces are already in use. We do this in our

Hi, I wonder if someone could help. I needed to transfer (copy) a workspace file that had been generated in linux (R 2.11) to windows running the same version of R 2.11 (but of course windows binary). Usually, there is no problem in doing this and all objects work as expected. I am often doing this to be able to produce wmf or emf graphic files that I need. This time I had some spectra that I

Dear all, I have uploaded to CRAN a new package: tseriesChaos. This is an early version (0.1) with basic tools for the explorative analysis of nonlinear time series motivated by chaos theory. Until now, the package is largely inspired by the TISEAN project (by Rainer Hegger, Holger Kantz and Thomas Schreiber: http://www.mpipks-dresden.mpg.de/~tisean/ ). This version includes: - Method of false

Dear all, I have uploaded to CRAN a new package: tseriesChaos. This is an early version (0.1) with basic tools for the explorative analysis of nonlinear time series motivated by chaos theory. Until now, the package is largely inspired by the TISEAN project (by Rainer Hegger, Holger Kantz and Thomas Schreiber: http://www.mpipks-dresden.mpg.de/~tisean/ ). This version includes: - Method of false

Hi, I?m trying to use the STLperArea function in the ndvits package. I can run the sample data (?SLPSAs_full?) without any problem. However, when I come to run the function on my own data, I get the following message. Waiting to confirm page change... Error in .checkPath(path) : no TISEAN executables found in that directory. Please set a proper TISEAN executables path using

Hi R-Users I need to estimate Lyapunov exponent of my time series. After reading description of all functions available I still don't know how to determine time delay. My time series length is 4200. Is it possible to determine time delay with other function's output or I can choose any random value? Thanks for your help! -- View this message in context:

hi all, I have a simple question that does power spectral analysis related to capacity dimension, information dimension, lyapunov exponent, hurst exponent. If yes then please show me the way. I am newbie in the world of chaos. Sayonara With Smile & With Warm Regards :-) G a u r a v Y a d a v Senior Executive Officer, Economic Research & Surveillance Department, Clearing

For a Monte Carlo study I need to extract from an lme model the estimated standard deviation of a random effect and store it in a vector. If I do a print() or summary() on the model, the number I need is displayed in the Console [it's the 0.1590195 in the output below] >print(fit) >Linear mixed-effects model fit by maximum likelihood > Data: datag2 > Log-likelihood:

Hi, Attempting to use LLVM in jitting mode for AMD64, we met a problem. When the jitted routine needs a big stack frame (> 1 page), the system attempts to call __chkstk to probe the stack. This attempt results in assertion in RuntimeDyldELF::resolveX86_64Relocation(), case ELF::R_X86_64_PC32, because the RealOffset does not fit in 32 bits. Same happens with __alloca (when

Hi, I spotted an inconsistency between the output of funfits under R, more specifically "nnreg", and the one obtained under S-plus. Notice the difference in the "Root MSE", "GVC" and "GVC cost=2". As a check mark, I am also including the data summaries, confirming that I was in fact using the same data set. Here's my output from both systems:

Hello, i have a little problem with R and i hope you can help me. I want to use splines to estimate a function but i want to force the interpolation to be monotone. Is this possible with R ? Thank you, Rémi. --------------------------------- [[alternative HTML version deleted]]

Hi, I met troubles with jitting x86 codes when using Intrinsic::x86_sse_ldmxcsr. The target code must execute some SSE2 instruction with DAZ/FTZ modes enabled and others with DAZ/FTZ disabled. I'm trying to get this by emitting LDMXCSR instructions with proper flag words. It appeared however that execution engine sometimes reorders these instructions with computational ones (say with

Hello, I would like to plot two superposed surfaces (original data, and smoothed models). Basically, I would like to create a usual persp plot with box, and then superpose a second surface (without box and axes) on the top plane of the box defined by the first graph. I don't know if it is clear :-( Is it possible to draw such a second surface in an existing persp plot in R? Another (related)