Jolliffe I. Principal Component Analysis (2ed., Springer, 2002)(518s)
Jolliffe I. Principal Component Analysis (2ed., Springer, 2002)(518s) Jolliffe I. Principal Component Analysis (2ed., Springer, 2002)(518s)
12.2. PCA and Atmospheric Time Series 311gether with the variance accounted for by each PC. Figure 12.5 displayssimilar plots for the real and imaginary parts of the first Hilbert EOF (labelledCEOFs on the plots). It can be seen that the real part of the firstHilbert EOF in Figure 12.5 looks similar to EOF1 in Figure 12.4. There arealso similarities between EOF3 and the imaginary part of the first HilbertEOF.Figures 12.6, 12.7 show plots of time series (scores), labelled temporalcoefficients in 12.7, for the first and third ordinary PCs, and the real andimaginary parts of the first HEOF, respectively. The similarity betweenthe first PC and the real part of the first HEOF is obvious. Both representthe same oscillatory behaviour in the series. The imaginary part of thefirst HEOF is also very similar, but lagged by π 2, whereas the scores on thethird EOF show a somewhat smoother and more regular oscillation. Cai andBaines (2001) note that the main oscillations visible in Figures 12.6, 12.7can be identified with well-known El Niño-Southern Oscillation (ENSO)events. They also discuss other physical interpretations of the results ofthe HEOF analysis relating them to other meteorological variables, andthey provide tests for statistical significance of HEOFs.The main advantage of HEOF analysis over ordinary PCA is its abilityto identify and reconstruct propagating waves, whereas PCA only findsstanding oscillations. For the first and second Hilbert EOFs, and for theirsum, this propagating behaviour is illustrated in Figure 12.8 by the movementof similar-valued coefficients from west to east as time progresses inthe vertical direction.12.2.4 Multitaper Frequency Domain-Singular ValueDecomposition (MTM SVD)In a lengthy paper, Mann and Park (1999) describe MTM-SVD, developedearlier by the same authors. It combines multitaper spectrum estimationmethods (MTM) with PCA using the singular value decomposition (SVD).The paper also gives a critique of several of the other techniques discussed inthe present section, together with frequency domain PCA which is coveredin Section 12.4.1. Mann and Park (1999) provide detailed examples, bothreal and artificial, in which MTM-SVD is implemented.Like MSSA, POP analysis, HEOF analysis and frequency domain PCA,MTM-SVD looks for oscillatory behaviour in space and time. It is closestin form to PCA of the spectral matrix F(λ) (Section 12.4.1), as it operatesin the frequency domain. However, in transferring from the time domainto the frequency domain MTM-SVD constructs a set of different taperedFourier transforms (hence the ‘multitaper’ in its name). The frequency domainmatrix is then subjected to a singular value decomposition, giving‘spatial’ or ‘spectral’ EOFs, and ‘principal modulations’ which are analogousto principal component scores. Mann and Park (1999) state that the
312 12. PCA for Time Series and Other Non-Independent DataFigure 12.4. The first four EOFs for Southern Hemisphere SST.Figure 12.5. Real and imaginary parts of the first Hilbert EOF for SouthernHemisphere SST.
- Page 292 and 293: 10.3. Sensitivity and Stability 261
- Page 294 and 295: 10.4. Robust Estimation of Principa
- Page 296 and 297: 10.4. Robust Estimation of Principa
- Page 298 and 299: 10.4. Robust Estimation of Principa
- Page 300 and 301: 11Rotation and Interpretation ofPri
- Page 302 and 303: 11.1. Rotation of Principal Compone
- Page 304 and 305: oot of the corresponding eigenvalue
- Page 306 and 307: 11.1. Rotation of Principal Compone
- Page 308 and 309: 11.1. Rotation of Principal Compone
- Page 310 and 311: 11.2. Alternatives to Rotation 279w
- Page 312 and 313: 11.2. Alternatives to Rotation 281F
- Page 314 and 315: 11.2. Alternatives to Rotation 283F
- Page 316 and 317: 11.2. Alternatives to Rotation 285T
- Page 318 and 319: 11.2. Alternatives to Rotation 287T
- Page 320 and 321: 11.2. Alternatives to Rotation 289A
- Page 322 and 323: 11.2. Alternatives to Rotation 291
- Page 324 and 325: 11.3. Simplified Approximations to
- Page 326 and 327: 11.3. Simplified Approximations to
- Page 328 and 329: 11.4. Physical Interpretation of Pr
- Page 330 and 331: 12Principal Component Analysis forT
- Page 332 and 333: 12.1. Introduction 301series is alm
- Page 334 and 335: 12.2. PCA and Atmospheric Time Seri
- Page 336 and 337: 12.2. PCA and Atmospheric Time Seri
- Page 338 and 339: and a typical row of the matrix is1
- Page 340 and 341: 12.2. PCA and Atmospheric Time Seri
- Page 344 and 345: 12.2. PCA and Atmospheric Time Seri
- Page 346 and 347: 12.2. PCA and Atmospheric Time Seri
- Page 348 and 349: 12.3. Functional PCA 317A key refer
- Page 350 and 351: 12.3. Functional PCA 319The sample
- Page 352 and 353: 12.3. Functional PCA 321speed (mete
- Page 354 and 355: 12.3. Functional PCA 323of the data
- Page 356 and 357: 12.3. Functional PCA 325subject to
- Page 358 and 359: 12.3. Functional PCA 327series than
- Page 360 and 361: 12.4. PCA and Non-Independent Data
- Page 362 and 363: 12.4. PCA and Non-Independent Data
- Page 364 and 365: 12.4. PCA and Non-Independent Data
- Page 366 and 367: 12.4. PCA and Non-Independent Data
- Page 368 and 369: 12.4. PCA and Non-Independent Data
- Page 370 and 371: 13.1. Principal Component Analysis
- Page 372 and 373: 13.1. Principal Component Analysis
- Page 374 and 375: 13.2. Analysis of Size and Shape 34
- Page 376 and 377: 13.2. Analysis of Size and Shape 34
- Page 378 and 379: 13.3. Principal Component Analysis
- Page 380 and 381: 13.3. Principal Component Analysis
- Page 382 and 383: 13.4. Principal Component Analysis
- Page 384 and 385: 13.4. Principal Component Analysis
- Page 386 and 387: 13.5. Common Principal Components 3
- Page 388 and 389: 13.5. Common Principal Components 3
- Page 390 and 391: 13.5. Common Principal Components 3
12.2. PCA and Atmospheric Time Series 311gether with the variance accounted for by each PC. Figure 12.5 displayssimilar plots for the real and imaginary parts of the first Hilbert EOF (labelledCEOFs on the plots). It can be seen that the real part of the firstHilbert EOF in Figure 12.5 looks similar to EOF1 in Figure 12.4. There arealso similarities between EOF3 and the imaginary part of the first HilbertEOF.Figures 12.6, 12.7 show plots of time series (scores), labelled temporalcoefficients in 12.7, for the first and third ordinary PCs, and the real andimaginary parts of the first HEOF, respectively. The similarity betweenthe first PC and the real part of the first HEOF is obvious. Both representthe same oscillatory behaviour in the series. The imaginary part of thefirst HEOF is also very similar, but lagged by π 2, whereas the scores on thethird EOF show a somewhat smoother and more regular oscillation. Cai andBaines (2001) note that the main oscillations visible in Figures 12.6, 12.7can be identified with well-known El Niño-Southern Oscillation (ENSO)events. They also discuss other physical interpretations of the results ofthe HEOF analysis relating them to other meteorological variables, andthey provide tests for statistical significance of HEOFs.The main advantage of HEOF analysis over ordinary PCA is its abilityto identify and reconstruct propagating waves, whereas PCA only findsstanding oscillations. For the first and second Hilbert EOFs, and for theirsum, this propagating behaviour is illustrated in Figure 12.8 by the movementof similar-valued coefficients from west to east as time progresses inthe vertical direction.12.2.4 Multitaper Frequency Domain-Singular ValueDecomposition (MTM SVD)In a lengthy paper, Mann and Park (1999) describe MTM-SVD, developedearlier by the same authors. It combines multitaper spectrum estimationmethods (MTM) with PCA using the singular value decomposition (SVD).The paper also gives a critique of several of the other techniques discussed inthe present section, together with frequency domain PCA which is coveredin Section 12.4.1. Mann and Park (1999) provide detailed examples, bothreal and artificial, in which MTM-SVD is implemented.Like MSSA, POP analysis, HEOF analysis and frequency domain PCA,MTM-SVD looks for oscillatory behaviour in space and time. It is closestin form to PCA of the spectral matrix F(λ) (Section 12.4.1), as it operatesin the frequency domain. However, in transferring from the time domainto the frequency domain MTM-SVD constructs a set of different taperedFourier transforms (hence the ‘multitaper’ in its name). The frequency domainmatrix is then subjected to a singular value decomposition, giving‘spatial’ or ‘spectral’ EOFs, and ‘principal modulations’ which are analogousto principal component scores. Mann and Park (1999) state that the