Multivariate methods Principal component analysis PCA The idea

  • Slides: 13
Download presentation
Multivariate methods Principal component analysis (PCA). The idea we replace a number of intercorrelated

Multivariate methods Principal component analysis (PCA). The idea we replace a number of intercorrelated variables with a smaller number of variables, we will make the data and our life simpler. We will draw a new axis according to the maximal spread of the data cloud, Z 1 = a 1 X 1 + a 2 X 2 + a 3 X 3+. . Component score – coordinate on this new axis. The new variable is a linear combination of the old ones. We do not test anything, make preparations for further tests.

Calculating one principal component The (adult) age of a butterfly caught in the wild,

Calculating one principal component The (adult) age of a butterfly caught in the wild, - cannot be directly measured, Following variables correlate with it but do not measure it: - wear; - damage; - capture date. no wear damage date size eggslaid

All correlate with each ohter: wear damage date size We will „merge“ them to

All correlate with each ohter: wear damage date size We will „merge“ them to one principal component. Coordinate on the new axis is the estimated age, will be used in further analyses. size

PCA gives the following results: wear damage date size age = 0, 56*wear+0, 58*damage+0,

PCA gives the following results: wear damage date size age = 0, 56*wear+0, 58*damage+0, 59*date+0, 02*size for each observation: 1 2 3 4 5 -2. 48 -2. 74 -0. 64 -0. 28 0. 30, do not worry about the minus

no wear damage date size eggslaid „age“

no wear damage date size eggslaid „age“

eggs laid estimated age Use when the value of a variable cannot be measured

eggs laid estimated age Use when the value of a variable cannot be measured - for technical reasons; - in principle.

Principal component is the axis of data cloud, there can be several, orthogonal to

Principal component is the axis of data cloud, there can be several, orthogonal to each other! Two principal components: classifying plant communities. Objects: forest patches; measured variables: abundances of different plant species; variables being estimated: parameters of the forest patches. Let there be two – humidity and fertility, calculate two PC and respective PC scores so that: Z 1 = a 1, 1 X 1 + a 1, 2 X 2 + a 1, 3 X 3+. . Z 2 = a 2, 1 X 1 + a 2, 2 X 2 +. . …

Ordination plot: fertility humidity

Ordination plot: fertility humidity

fertility Biplot: nettle heather humidity

fertility Biplot: nettle heather humidity

Canonical correspondence analysis fertility CCA biplot humidity light

Canonical correspondence analysis fertility CCA biplot humidity light

Discriminant function analysis (DFA) - we know in advance in which classes objects belong;

Discriminant function analysis (DFA) - we know in advance in which classes objects belong; - which new variable - a linear combination of original variables best allows to discriminate between the species. Z 1 = a 1, 1 X 1 + a 1, 2 X 2 + a 1, 3 X 3+. . - each object has a a coordinate on this axis, and there is a critical value according to which the decision is made.