Modes of variation

In statistics, modes of variation[1] are a continuously indexed set of vectors or functions that are centered at a mean and are used to depict the variation in a population or sample. Typically, variation patterns in the data can be decomposed in descending order of eigenvalues with the directions represented by the corresponding eigenvectors or eigenfunctions. Modes of variation provide a visualization of this decomposition and an efficient description of variation around the mean. Both in principal component analysis (PCA) and in functional principal component analysis (FPCA), modes of variation play an important role in visualizing and describing the variation in the data contributed by each eigencomponent.[2] In real-world applications, the eigencomponents and associated modes of variation aid to interpret complex data, especially in exploratory data analysis (EDA).

## Contents

* 1 Formulation
* 1.1 Modes of variation in PCA
* 1.2 Modes of variation in FPCA
* 2 Estimation
* 2.1 Modes of variation in PCA
* 2.2 Modes of variation in FPCA
* 3 Applications
* 4 References

## Formulation

Modes of variation are a natural extension of PCA and FPCA.

### Modes of variation in PCA

If a random vector 𝐗=(X1,X2,⋯,Xp)T has the mean vector μp, and the covariance matrix Σp×p with eigenvalues λ1≥λ2≥⋯≥λp≥0 and corresponding orthonormal eigenvectors 𝐞1,𝐞2,⋯,𝐞p, by eigendecomposition of a real symmetric matrix, the covariance matrix Σ can be decomposed as

Σ=𝐐Λ𝐐T,

where 𝐐 is an orthogonal matrix whose columns are the eigenvectors of Σ, and Λ is a diagonal matrix whose entries are the eigenvalues of Σ. By the Karhunen–Loève expansion for random vectors, one can express the centered random vector in the eigenbasis

𝐗−μ=∑k=1pξk𝐞k,

where ξk=𝐞kT(𝐗−μ) is the principal component[3] associated with the k-th eigenvector 𝐞k, with the properties

E⁡(ξk)=0,Var⁡(ξk)=λk, and E⁡(ξkξl)=0 for l≠k.

Then the k-th mode of variation of 𝐗 is the set of vectors, indexed by α,

𝐦k,α=μ±αλk𝐞k,α∈[−A,A],

where A is typically selected as 2 or 3.

### Modes of variation in FPCA

For a square-integrable random function X(t),t∈𝒯⊂Rp, where typically p=1 and 𝒯 is an interval, denote the mean function by μ(t)=E⁡(X(t)), and the covariance function by

G(s,t)=Cov⁡(X(s),X(t))=∑k=1∞λkφk(s)φk(t),

where λ1≥λ2≥⋯≥0 are the eigenvalues and {φ1,φ2,⋯} are the orthonormal eigenfunctions of the linear Hilbert–Schmidt operator

G:L2(𝒯)->L2(𝒯),G(f)=∫𝒯G(s,t)f(s)ds.

By the Karhunen–Loève theorem, one can express the centered function in the eigenbasis,

X(t)−μ(t)=∑k=1∞ξkφk(t),

where

ξk=∫𝒯(X(t)−μ(t))φk(t)dt

is the k-th principal component with the properties

E⁡(ξk)=0,Var⁡(ξk)=λk, and E⁡(ξkξl)=0 for l≠k.

Then the k-th mode of variation of X(t) is the set of functions, indexed by α,

mk,α(t)=μ(t)±αλkφk(t), t∈𝒯, α∈[−A,A]

that are viewed simultaneously over the range of α, usually for A=2 or 3.[2]

## Estimation

The formulation above is derived from properties of the population. Estimation is needed in real-world applications. The key idea is to estimate mean and covariance.

### Modes of variation in PCA

Suppose the data 𝐱1,𝐱2,⋯,𝐱n represent n independent drawings from some p-dimensional population 𝐗 with mean vector μ and covariance matrix Σ. These data yield the sample mean vector 𝐱‾, and the sample covariance matrix 𝐒 with eigenvalue-eigenvector pairs (λ^1,𝐞^1),(λ^2,𝐞^2),⋯,(λ^p,𝐞^p). Then the k-th mode of variation of 𝐗 can be estimated by

𝐦^k,α=𝐱‾±αλ^k𝐞^k,α∈[−A,A].

### Modes of variation in FPCA

Consider n realizations X1(t),X2(t),⋯,Xn(t) of a square-integrable random function X(t),t∈𝒯 with the mean function μ(t)=E⁡(X(t)) and the covariance function G(s,t)=Cov⁡(X(s),X(t)). Functional principal component analysis provides methods for the estimation of μ(t) and G(s,t) in detail, often involving point wise estimate and interpolation. Substituting estimates for the unknown quantities, the k-th mode of variation of X(t) can be estimated by

m^k,α(t)=μ^(t)±αλ^kφ^k(t),t∈𝒯,α∈[−A,A].

## Applications

The first and second modes of variation of female mortality data from 41 countries in 2003[4] The first and second modes of variation of male mortality data from 41 countries in 2003[4]

Modes of variation are useful to visualize and describe the variation patterns in the data sorted by the eigenvalues. In real-world applications, modes of variation associated with eigencomponents allow to interpret complex data, such as the evolution of function traits[5] and other infinite-dimensional data.[6] To illustrate how modes of variation work in practice, two examples are shown in the graphs to the right, which display the first two modes of variation. The solid curve represents the sample mean function. The dashed, dot-dashed, and dotted curves correspond to modes of variation with α=±1,±2, and ±3, respectively.

The first graph displays the first two modes of variation of female mortality data from 41 countries in 2003.[4] The object of interest is log hazard function between ages 0 and 100 years. The first mode of variation suggests that the variation of female mortality is smaller for ages around 0 or 100, and larger for ages around 25. An appropriate and intuitive interpretation is that mortality around 25 is driven by accidental death, while around 0 or 100, mortality is related to congenital disease or natural death.

Compared to female mortality data, modes of variation of male mortality data shows higher mortality after around age 20, possibly related to the fact that life expectancy for women is higher than that for men.

## References

1. ↑ Castro, P. E.; Lawton, W. H.; Sylvestre, E. A. (November 1986). "Principal Modes of Variation for Processes with Continuous Sample Curves". Technometrics 28 (4): 329. doi:10.2307/1268982. ISSN 0040-1706.
2. ↑ 2.0 2.1 Wang, Jane-Ling; Chiou, Jeng-Min; Müller, Hans-Georg (June 2016). "Functional Data Analysis". Annual Review of Statistics and Its Application 3 (1): 257–295. doi:10.1146/annurev-statistics-041715-033624. ISSN 2326-8298. https://zenodo.org/record/895750.
3. ↑ Kleffe, Jürgen (January 1973). "Principal components of random variables with values in a seperable HILBERT space". Mathematische Operationsforschung und Statistik 4 (5): 391–406. doi:10.1080/02331887308801137. ISSN 0047-6277.
4. ↑ 4.0 4.1 4.2 "Human Mortality Database". https://www.mortality.org/.
5. ↑ Kirkpatrick, Mark; Heckman, Nancy (August 1989). "A quantitative genetic model for growth, shape, reaction norms, and other infinite-dimensional characters". Journal of Mathematical Biology 27 (4): 429–450. doi:10.1007/bf00290638. ISSN 0303-6812. PMID 2769086.
6. ↑ Jones, M. C.; Rice, John A. (May 1992). "Displaying the Important Features of Large Collections of Similar Curves". The American Statistician 46 (2): 140–145. doi:10.1080/00031305.1992.10475870. ISSN 0003-1305.

0.00

(0 votes)

Original source: https://en.wikipedia.org/wiki/Modes of variation. Read more

|

Retrieved from "https://handwiki.org/wiki/index.php?title=Modes_of_variation&oldid=4150257"