Singular spectrum analysis
In time series analysis, singular spectrum analysis is a nonparametric spectral estimation method. It combines elements of classical time series analysis, multivariate statistics, multivariate geometry, dynamical systems and signal processing. Its roots lie in the classical Karhunen –Loève spectral decomposition of time series and random fields and in the Mañé –Takens embedding theorem. SSA can be an aid in the decomposition of time series into a sum of components, each having a meaningful interpretation. The name "singular spectrum analysis" relates to the spectrum of eigenvalues in a singular value decomposition of a covariance matrix, and not directly to a frequency domain decomposition.
Brief history
The origins of SSA and, more generally, of subspace-based methods for signal processing, go back to the eighteenth century. A key development was the formulation of the spectral decomposition of the covariance operator of stochastic processes by Kari Karhunen and Michel Loève in the late 1940s.Broomhead and King and Fraedrich proposed to use SSA and multichannel SSA in the context of nonlinear dynamics for the purpose of reconstructing the attractor of a system from measured time series. These authors provided an extension and a more robust application of the idea of reconstructing dynamics from a single time series based on the embedding theorem. Several other authors had already applied simple versions of M-SSA to meteorological and ecological data sets.
Ghil, Vautard and their colleagues noticed the analogy between the trajectory matrix of Broomhead and King, on the one hand, and the Karhunen–Loeve decomposition, on the other. Thus, SSA can be used as a time-and-frequency domain method for time series analysis — independently from attractor reconstruction and including cases in which the latter may fail. The survey paper of Ghil et al. is the basis of the section of this article. A crucial result of the work of these authors is that SSA can robustly recover the "skeleton" of an attractor, including in the presence of noise. This skeleton is formed by the least unstable periodic orbits, which can be identified in the eigenvalue spectra of SSA and M-SSA. The identification and detailed description of these orbits can provide highly useful pointers to the underlying nonlinear dynamics.
The so-called ‘Caterpillar’ methodology is a version of SSA that was developed in the former Soviet Union, independently of the mainstream SSA work in the West. This methodology became known in the rest of the world more recently. ‘Caterpillar-SSA’ emphasizes the concept of separability, a concept that leads, for example, to specific recommendations concerning the choice of SSA parameters. This method is thoroughly described in of this article.
Methodology
In practice, SSA is a nonparametric spectral estimation method based on embedding a time series in a vector space of dimension. SSA proceeds by diagonalizing the lag-covariance matrix of to obtain spectral information on the time series, assumed to be stationary in the weak sense. The matrix can be estimated directly from the data as a Toeplitz matrix with constant diagonals, i.e., its entries depend only on the lag :An alternative way to compute, is by using the "trajectory matrix" that is formed by lag-shifted copies of, which are long; then
The eigenvectors of the lag-covariance matrix are called temporal empirical orthogonal functions. The eigenvalues of account for the partial variance in the
direction and the sum of the eigenvalues, i.e., the trace of
, gives the total variance of the original time series
. The name of the method derives from the singular values of
Decomposition and reconstruction
Projecting the time series onto each EOF yields the correspondingtemporal principal components :
An oscillatory mode is characterized by a pair of
nearly equal SSA eigenvalues and associated PCs that are in approximate phase quadrature. Such a pair can represent efficiently a nonlinear, anharmonic oscillation. This is due to the fact that a single pair of data-adaptive SSA eigenmodes often will capture better the basic periodicity of an oscillatory mode than methods with fixed basis functions, such as the sines and cosines used in the Fourier transform.
The window width determines the longest periodicity captured by SSA. Signal-to-noise separation can be obtained by merely inspecting the slope break in a "scree diagram" of eigenvalues or singular values vs.. The point at which this break occurs should not be confused with a "dimension" of the underlying deterministic dynamics.
A Monte-Carlo test can be applied to ascertain the statistical significance of the oscillatory pairs detected by SSA. The entire time series or parts of it that correspond to trends, oscillatory modes or noise can be reconstructed by using linear combinations of the PCs and EOFs, which provide the reconstructed components :
here is the set of EOFs on which the reconstruction is based. The values of the normalization factor, as well as of the lower and upper bound of summation and, differ between the central part of the time series and the vicinity of its endpoints.
Multivariate extension
Multi-channel SSA is a natural extension of SSA to an -channel time series of vectors or maps with data points. In the meteorological literature, extended EOF analysis is often assumed to be synonymous with M-SSA. The two methods are both extensions of classical principal component analysis but they differ in emphasis: EEOF analysis typically utilizes a number of spatial channels much greater than the number of temporal lags, thus limiting the temporal and spectral information. In M-SSA, on the other hand, one usually chooses. Often M-SSA is applied to a few leading PCs of the spatial data, with chosen large enough to extract detailed temporal and spectral information from the multivariate time series. However, Groth and Ghil have demonstrated possible negative effects of this variance compression on the detection rate of weak signals when the number of retained PCs becomes too small. This practice can further affect negatively the judicious reconstruction of the spatio-temporal patterns of such weak signals, and Groth et al. recommend retaining a maximum number of PCs, i.e.,.Groth and Ghil have demonstrated that a classical M-SSA analysis suffers from a degeneracy problem, namely the EOFs do not separate well between distinct oscillations when the corresponding eigenvalues are similar in size. This problem is a shortcoming of principal component analysis in general, not just of M-SSA in particular. In order to reduce mixture effects and to improve the physical interpretation, Groth and Ghil have proposed a subsequent VARIMAX rotation of the spatio-temporal EOFs of the M-SSA. To avoid a loss of spectral properties, they have introduced a slight modification of the common VARIMAX rotation that does take the spatio-temporal structure of ST-EOFs into account. Alternatively, a closed matrix formulation of the algorithm for the simultaneous rotation of the EOFs by iterative SVD decompositions has been proposed.
M-SSA has two forecasting approaches known as recurrent and vector. The discrepancies between these two approaches are attributable to the organization of the single trajectory matrix of each series into the block trajectory matrix in the multivariate case. Two trajectory matrices can be organized as either vertical or horizontal as was recently introduced in Hassani and Mahmoudvand, and it was shown that these constructions lead to better forecasts. Accordingly, we have four different forecasting algorithms that can be exploited in this version of MSSA.
Prediction
In this subsection, we focus on phenomena that exhibit a significant oscillatory component: repetition increases understanding and hence confidence in a prediction method that is closely connected with such understanding.Singular spectrum analysis and the maximum entropy method have been combined to predict a variety of phenomena in meteorology, oceanography and climate dynamics. First, the “noise” is filtered out by projecting the time series onto a subset of leading EOFs obtained by SSA; the selected subset should include statistically significant, oscillatory modes. Experience shows that this approach works best when the partial variance associated with the pairs of RCs that capture these modes is large.
The prefiltered RCs are then extrapolated by least-square fitting to an autoregressive model AR, whose coefficients give the MEM spectrum of the remaining “signal”. Finally, the extended RCs are used in the SSA reconstruction process to produce the forecast values. The reason why this approach – via SSA prefiltering, AR extrapolation of the RCs, and SSA reconstruction – works better than the customary AR-based prediction is explained by the fact that the individual RCs are narrow-band signals, unlike the original, noisy time series X. In fact, the optimal order p obtained for the individual RCs is considerably lower than the one given by the standard Akaike information criterion or similar ones.
Spatio-temporal gap filling
The gap-filling version of SSA can be used to analyze data sets that are unevenly sampled or contain missing data. For a univariate time series, the SSA gap filling procedure utilizes temporal correlations to fill in the missing points. For a multivariate data set, gap filling by M-SSA takes advantage of both spatial and temporal correlations. In either case: estimates of missing data points are produced iteratively, and are then used to compute a self-consistent lag-covariance matrix and its EOFs ; and cross-validation is used to optimize the window width and the number of leading SSA modes to fill the gaps with the iteratively estimated "signal," while the noise is discarded.SSA as a model-free tool
The areas where SSA can be applied are very broad: climatology, marine science, geophysics, engineering, image processing, medicine, econometrics among them. Hence different modifications of SSA have been proposed and different methodologies of SSA are used in practical applications such as trend extraction, periodicity detection, seasonal adjustment, smoothing, noise reduction .Basic SSA
SSA can be used as a model-free technique so that it can be applied to arbitrary time series including non-stationary time series. The basic aim of SSA is to decompose the time series into the sum of interpretable components such as trend, periodic components and noise with no a-priori assumptions about the parametric form of these components.Consider a real-valued time series of length. Let be some integer called the window length and.
Main algorithm of SSA
1st step: Embedding.Form the trajectory matrix of the series, which is the matrix
where are lagged vectors of size. The matrix is a Hankel matrix which means that has equal elements on the anti-diagonals.
2nd step: Singular Value Decomposition.
Perform the singular value decomposition of the trajectory matrix. Set and denote by the eigenvalues of taken in the decreasing order of magnitude and by the orthonormal system of the eigenvectors of the matrix corresponding to these eigenvalues.
Set and . In this notation, the SVD of the trajectory matrix can be written as
where
are matrices having rank 1; these are called elementary matrices. The collection will be called the th eigentriple of the SVD. Vectors are the left singular vectors of the matrix, numbers are the singular values and provide the singular spectrum of ; this gives the name to SSA. Vectors are called vectors of principal components.
3rd step: Eigentriple grouping.
Partition the set of indices into disjoint subsets.
Let. Then the resultant matrix corresponding to the group is defined as. The resultant matrices are computed for the groups and the grouped SVD expansion of can now be written as
4th step: Diagonal averaging.
Each matrix of the grouped decomposition is hankelized and then the obtained Hankel matrix is transformed into a new series of length using the one-to-one correspondence between Hankel matrices and time series.
Diagonal averaging applied to a resultant matrix produces a reconstructed series. In this way, the initial series is decomposed into a sum of reconstructed subseries:
This decomposition is the main result of the SSA algorithm. The decomposition is meaningful if each reconstructed
subseries could be classified as a part of either trend or some periodic component or noise.
Theory of SSA separability
The two main questions which the theory of SSA attempts to answer are: what time series components can be separated by SSA, and how to choose the window length and make proper grouping for extraction of a desirable component. Many theoretical results can be found in Golyandina et al..Trend, periodic components and noise are asymptotically separable as. In practice is fixed and one is interested in approximate separability between time series components. A number of indicators of approximate separability can be used, see Golyandina et al.. The window length determines the resolution of the method: larger values of provide more refined decomposition into elementary components and therefore better separability. The window length determines the longest periodicity captured by SSA. Trends can be extracted by grouping of eigentriples with slowly varying eigenvectors. A sinusoid with frequency smaller than 0.5 produces two approximately equal eigenvalues and two sine-wave eigenvectors with the same frequencies and -shifted phases.
Separation of two time series components can be considered as extraction of one component in the presence of perturbation by the other component. SSA perturbation theory is developed in Nekrutkin and Hassani et al..
Forecasting by SSA
If for some series the SVD step in Basic SSA gives, then this series is called time series of rank . The subspace spanned by the leading eigenvectors is called signal subspace. This subspace is used for estimating the signal parameters in signal processing, e.g. ESPRIT for high-resolution frequency estimation. Also, this subspace determines the linear homogeneous recurrence relation governing the series, which can be used for forecasting. Continuation of the series by the LRR is similar to forward linear prediction in signal processing.Let the series be governed by the minimal LRR. Let us choose, be the eigenvectors, which are provided by the SVD step of SSA. Then this series is governed by an LRR, where are expressed through , and can be continued by the same LRR.
This provides the basis for SSA recurrent and vector forecasting algorithms. In practice, the signal is corrupted by a perturbation, e.g., by noise, and its subspace is estimated by SSA approximately. Thus, SSA forecasting can be applied for forecasting of a time series component that is approximately governed by an LRR and is approximately separated from the residual.
Multivariate extension
Multi-channel, Multivariate SSA is a natural extension of SSA to for analyzing multivariate time series, where the size of different univariate series does not have to be the same. The trajectory matrix of multi-channel time series consists of linked trajectory matrices of separate times series. The rest of the algorithm is the same as in the univariate case. System of series can be forecasted analogously to SSA recurrent and vector algorithms. MSSA has many applications. It is especially popular in analyzing and forecasting economic and financial time series with short and long series length.Other multivariate extension is 2D-SSA that can be applied to two-dimensional data like digital images. The analogue of trajectory matrix is constructed by moving 2D windows of size.
MSSA and causality
A question that frequently arises in time series analysis is whether one economic variable canhelp in predicting another economic variable. One way to address this question was proposed by
Granger, in which he formalized the causality concept. A comprehensive causality test based on MSSA has recently introduced for causality measurement. The test is based on the forecasting accuracy and predictability of the direction of change of the MSSA algorithms.
MSSA and EMH
The MSSA forecasting results can be used in examining the efficient market hypothesis controversy.The EMH suggests that the information contained in the price series of an asset is reflected “instantly, fully, and perpetually” in the asset’s current price. Since the price series and the information contained in it are available to all market participants, no one can benefit by attempting to take advantage of the information contained in the price history of an asset by trading in the markets. This is evaluated using two series with different series length in a multivariate system in SSA analysis.
MSSA, SSA and business cycles
Business cycles plays a key role in macroeconomics, and are interest for a variety of players in the economy, including central banks, policy-makers, and financial intermediaries. MSSA-based methods for tracking business cycles have been recently introduced, and have been shown to allow for a reliable assessment of the cyclical position of the economy in real-time.MSSA, SSA and unit root
SSA's applicability to any kind of stationary or deterministically trending series has been extended to the case of a series with a stochastic trend, also known as a series with a unit root. In Hassani and Thomakos and Thomakos the basic theory on the properties and application of SSA in the case of series of a unit root is given, along with several examples. It is shown that SSA in such series produces a special kind of filter, whose form and spectral properties are derived, and that forecasting the single reconstructed component reduces to a moving average. SSA in unit roots thus provides an `optimizing' non-parametric framework for smoothing series with a unit root. This line of work is also extended to the case of two series, both of which have a unit root but are cointegrated. The application of SSA in this bivariate framework produces a smoothed series of the common root component.Gap-filling
The gap-filling versions of SSA can be used to analyze data sets that are unevenly sampled or contain missing data.Schoellhamer shows that the straightforward idea to formally calculate approximate inner products omitting unknown terms is workable for long stationary time series.
Golyandina and Osipov uses the idea of filling in missing entries in vectors taken from the given subspace. The recurrent and vector SSA forecasting can be considered as particular cases of filling in algorithms described in the paper.
Detection of structural changes
SSA can be effectively used as a non-parametric method of time series monitoring and change detection. To do that, SSA performs the subspace tracking in the following way. SSA is applied sequentially to the initial parts of the series, constructs the corresponding signal subspaces and checks the distances between these subspaces and the lagged vectors formed from the few most recent observations. If these distances become too large, a structural change is suspected to have occurred in the series.In this way, SSA could be used for change detection not only in trends but also in the variability of the series, in the mechanism that determines dependence between different series and even in the noise structure. The method have proved to be useful in different engineering problems.
Relation between SSA and other methods
SSA and Autoregression. Typical model for SSA is, where and is noise. The model of AR is. Despite these two models look similar they are very different. SSA considers AR as a noise component only. AR, which is red noise, is typical model of noise for Monte-Carlo SSA.SSA and spectral Fourier Analysis. In contrast with Fourier analysis with fixed basis of sine and cosine functions, SSA uses an adaptive basis generated by the time series itself. As a result, the underlying model in SSA is more general and SSA can extract amplitude-modulated sine wave components with frequencies different from. SSA-related methods like ESPRIT can estimate frequencies with higher resolution than spectral Fourier analysis.
SSA and Linear Recurrence Relations. Let the signal be modeled by a series, which satisfies a linear recurrence relation ; that is, a series that can be represented as sums of products of exponential, polynomial and sine wave functions. This includes the sum of dumped sinusoids model whose complex-valued form is. SSA-related methods allow estimation of frequencies and exponential factors . Coefficients can be estimated by the least squares method. Extension of the model, where are replaced by polynomials of, can be also considered within the SSA-related methods.
SSA and Signal Subspace methods. SSA can be considered as a subspace-based method, since it allows estimation of the signal subspace of dimension by.
SSA and State Space Models. The main model behind SSA is, where and is noise. Formally, this model belongs to the general class of state space models. The specifics of SSA is in the facts that parameter estimation is a problem of secondary importance in SSA and the data analysis procedures in SSA are nonlinear as they are based on the SVD of either trajectory or lag-covariance matrix.
SSA and Independent Component Analysis. SSA is used in blind source separation by ICA as a preprocessing step. On the other hand, ICA can be used as a replacement of the SVD step in the SSA algorithm for achieving better separability.
SSA and Regression. SSA is able to extract polynomial and exponential trends. However, unlike regression, SSA does not assume any parametric model which may give significant advantage when an exploratory data analysis is performed with no obvious model in hand.
SSA and Linear Filters. The reconstruction of the series by SSA can be considered as adaptive linear filtration. If the window length is small, then each eigenvector generates a linear filter of width for reconstruction of the middle of the series,. The filtration is non-causal. However, the so-called Last-point SSA can be used as a causal filter.
SSA and Density Estimation. Since SSA can be used as a method of data smoothing it can be used as a method of non-parametric density estimation.