Jolliffe I. Principal Component Analysis (2ed., Springer, 2002)(518s)

302 12. PCA for Time Series and Other Non-Independent Dataencountered in this area, with observations corresponding to times andvariables to spatial position, they are not necessarily restricted to suchdata.Time series are usually measured at discrete points in time, but sometimesthe series are curves. The analysis of such data is known as functionaldata analysis (functional PCA is the subject of Section 12.3). The final sectionof the chapter collects together a number of largely unconnected ideasand references concerning PCA in the context of time series and othernon-independent data.12.2 PCA-Related Techniques for (Spatio-)Temporal Atmospheric Science DataIt was noted in Section 4.3 that, for a common type of data in atmosphericscience, the use of PCA, more often referred to as empirical orthogonalfunction (EOF) analysis, is widespread. The data concerned consist of measurementsof some variable, for example, sea level pressure, temperature,...,atp spatial locations (usually points on a grid) at n different times. Themeasurements at different spatial locations are treated as variables and thetime points play the rôle of observations. An example of this type was givenin Section 4.3. It is clear that, unless the observations are well-separatedin time, there is likely to be correlation between measurements at adjacenttime points, so that we have non-independence between observations. Severaltechniques have been developed for use in atmospheric science that takeaccount of correlation in both time and space, and these will be describedin this section. First, however, we start with the simpler situation wherethere is a single time series. Here we can use a principal component-liketechnique, called singular spectrum analysis (SSA), to analyse the autocorrelationin the series. SSA is described in Section 12.2.1, as is its extensionto several time series, multichannel singular spectrum analysis (MSSA).Suppose that a set of p series follows a multivariate first-order autoregressivemodel in which the values of the series at time t are linearly relatedto the values at time (t − 1), except for a multivariate white noise term.An estimate of the matrix defining the linear relationship can be subjectedto an eigenanalysis, giving insight into the structure of the series. Such ananalysis is known as principal oscillation pattern (POP) analysis, and isdiscussed in Section 12.2.2.One idea underlying POP analysis is that there may be patterns in themaps comprising our data set, which travel in space as time progresses, andthat POP analysis can help to find such patterns. Complex (Hilbert) empiricalorthogonal functions (EOFs), which are described in Section 12.2.3,are designed to achieve the same objective. Detection of detailed oscillatorybehaviour is also the aim of multitaper frequency-domain singularvalue decomposition, which is the subject of Section 12.2.4.