Estimation of the Lyapunov exponent spectrum

This function estimates the Lyapunov exponent spectrum through the QR decomposition method considering the argument set selected by the user.

lyapunov.spec(x, blocking = c("FULL", "NOVER", "EQS", "BOOT"), B = 100,
  doplot = TRUE)

a matrix or data frame containing the partial derivatives of jacobian.


a character denoting if you consider the full sample FULL, the non-overlapping sample NOVER, the equally spaced sample EQS or the bootstrap sample BOOT (Default FULL).


a non-negative integer denoting the number of bootstrap iterations (Default 100).


a logical value denoting if you want to draw a plot TRUE or not FALSE. If it is TRUE shows as many graphs as embedding dimensions have been considered. Each of them represents the evolution of the Lyapunov exponent spectrum values for the whole period considering the blocking method selected by the user (Default TRUE).


A list containing the Lyapunov exponent spectrum considering the parameter set selected by the user. The dataset saved by each blocking method are the estimated Lyapunov exponent value, the standard error, the z-test value and the p-value. If the user chooses the non-overlapping sample, the equally spaced sample or the bootstrap sample the mean and median values of the Lyapunov exponent are showed. Also some details about the embedding dimension, the sample size, the block length and the block number are recorded.


Ellner, S., Gallant, A., McCaffrey, D., Nychka, D. 1991 Convergence rates and data requirements for jacobian-based estimates of lyapunov exponents from data. Physics Letters A 153(6):357-363.

McCaffrey, D.F., Ellner, S., Gallant, A.R., Nychka, D.W. 1992 Estimating the lyapunov exponent of a chaotic system with nonparametric regression. Journal of the American Statistical Association 87(419):682-695.

Nychka, D., Ellner, S., Gallant, A.R., McCaffrey, D. 1992 Finding chaos in noisy systems. Journal of the Royal Statistical Society 54(2):399-426.

Whang, Y.J., Linton, O. 1999 The asymptotic distribution of nonparametric estimates of the lyapunov exponent for stochastic time series. Journal of Econometrics 91(1):1-42.

Shintani, M., Linton, O. 2004 Nonparametric neural network estimation of Lyapunov exponents and a direct test for chaos. Journal of Econometrics 120(1):1-33.

See Also

lyapunov.max, lyapunov

  • lyapunov.spec
## We show below an example considering time series from the
## logistic equation. We have estimated the Lyapunov exponent
## spectrum considering the bootstrap sample for an embedding
## dimension (m=4). First of all, we need to estimates the
## partial derivatives of the jacobian.
ts<-data$`Logistic 100`$time.serie
# }
Documentation reproduced from package DChaos, version 0.1-1, License: GPL (>= 2)

Community examples

Looks like there are no examples yet.