next up previous contents
Next: Bayesian methods for data Up: Practical considerations Previous: Delay coordinates in practice   Contents


Prediction algorithms

Assume we are analysing the scalar time series $ x(1), \ldots, x(T)$. Using the time-delay embedding with delay $ d$ transforms the series to vectors of the form $ \mathbf{y}(t) = (x(t-d+1), \ldots, x(t))$. Prediction in these coordinates corresponds to predicting the next sample $ x(t+1)$ from the $ d$ previous ones, as all the other components of $ \mathbf{y}(t+1)$ are directly available in $ \mathbf{y}(t)$. Thus the problem reduces to finding a predictor of the form

$\displaystyle x(t+1) = f(\mathbf{y}(t)).$ (2.8)

Using a simple linear function for $ f$ leads to a linear auto-regressive (AR) model [20]. There are also several generalisations of the same model leading to other algorithms for the same purpose.

Farmer and Sidorowich [14] propose a locally linear predictor in which there are several linear predictors for different areas of the embedded state-space. The achieved performance is comparable with the global linear predictor for small embedding dimensionalities but as the embedding dimension grows, the local method clearly outperforms the global one. The problem with the locally linear predictor is that it is not continuous.

Casdagli [8] presents a review of several different predictors including a globally linear predictor, a locally linear predictor and a global nonlinear predictor using a radial basis function (RBF) network [21] as the nonlinearity. According to his results, the nonlinear methods are best for small data sets whereas the locally linear method of Farmer and Sidorowich gives the best results for large data sets.

For noisy data, the choice of coordinates can make a big difference on the success of the prediction algorithm. The embedding approach is by no means the only possible alternative, and as Casdagli et al. note in [9], there are often clearly better alternatives. There are many approaches in the field of neural computation where the choice is left to the algorithm. This corresponds to trying to blindly invert Equation (2.7). Such models are called (nonlinear) state-space models and they are studied in detail in Section 4.2.


next up previous contents
Next: Bayesian methods for data Up: Practical considerations Previous: Delay coordinates in practice   Contents
Antti Honkela 2001-05-30