iaaft | R Documentation |
Performs the iterated amplitude adjusted Fourier transformation (IAAFT) upon a time-series x
.
iaaft(x, xdist=x, N=1, ...)
x |
a numeric vector containing the values of the time-series. |
xdist |
a numeric vector containing the values of the distribution to match. Default is set to |
N |
number of surrogates to create. Default is 1. |
... |
additional parameters to control for accuracy. |
The algorithm can be used to remove a time-series' x
non-linear correlation structure (multi-fractality) while keeping the degree of linear correlation (persistence). It can further transform \itemx to match any empirical distribution xdist
.
Returns a numeric vector of the surrogated time-series.
Code is based on Henning Rust (2006).
Wolfgang Schadner
T. Schreiber and A. Schmitz (2000), Surrogate time series, Physica D, 142, pp.346-382.
iaawt
n <- 1000 # random normal variable with positive auto-correlation: x <- sort(rnorm(n)) # random variable from t-distribution: y <- rt(n, 5) z <- iaaft(x, y)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.