2
$\begingroup$

I have a question about conducting a PCA between variables that are measured in different units. I understand the importance of using a correlation matrix versus a covariance matrix to minimize variance. The data I'm working with is not normally distributed and has not been transformed in other tests.

For example, there are three variables A, B, and C, and 20 observations, where 10 observations are measured using 1 set of units, and the other 10 observations are measured using another set of units*. The values between the units are quite different in in value and variance (expected). The data is not normal in either units and has not been transformed.

The measurements using the first set of units is 2 to 3 orders of magnitude higher than those measured using the other units (expected). I have conducted a PCA using a correlation matrix and interpreted results. However a non-statistician recommended I `standardize' the measured data, such that I'm using ratio or fractions for all the observations for each of the variables: Variable A/(sum of all 3 variables) and so on and so forth for Variables B and C.

However, a PCA using the contributing fraction of each variable is different from PCA using the measured value in different units leading to different eigenvalues and eigenvectors, thus leading to two different scientific interpretations.

Beyond using the appropriate association matrix, Is this "standardization" step valid and or necessary from a statistical perspective? Update: Should PCA be done on compositional data?

$\endgroup$
4
  • $\begingroup$ I didn't understand about the different sets of units - units of measurement or experimental units (two groups of individuals with A, B, & C measured on each individual)? $\endgroup$ Jul 15, 2013 at 9:53
  • $\begingroup$ Yes different units of measuring the variable a, b, and c in different samples. As an example, variable A is measured in units of mass, while variable B is measured in units of force. This leads to differences in absolute values and variances in the data set. $\endgroup$
    – Oleic
    Jul 16, 2013 at 0:36
  • $\begingroup$ That different variables are measured in different units is the situation that typically leads people to want to standardize them i.e. use the correlation matrix for PCA. (I'm still not sure quite what you meant by 20 observations being made of three variables, 10 using one set of units & 10 using another.) $\endgroup$ Jul 16, 2013 at 8:36
  • $\begingroup$ It's also hard to see what meaning you'd give to $\frac{A}{A+B+C}$ when $A$, $B$, & $C$ are each measured in different units. $\endgroup$ Jul 16, 2013 at 8:47

1 Answer 1

3
$\begingroup$

First, whether you use the covariance matrix, or the correlation matrix (equivalent to standardizing each variable before carrying out PCA on the covariance matrix), or transform the data in any other way before carrying out PCA, the results of the PCA apply to that transformed data. So you should not be surprised to see different eigensystems using different transformations; any interpretations you may make may of course be different, but are are not conflicting. If they seem to conflict you must be misinterpreting them.

Second, whether it's more meaningful to express each variable as a fraction of the sum of variables for each individual is for you to decide, before thinking about principal components. If it is more meaningful, PCA on the data thus transformed may not be what you want: any one variable is expressible in terms of the other two, which are still constrained not to exceed unity in total. A scatterplot would be an obvious method to look at three variables, using barycentric co-ordinates if you like. If you still need PCA for something, Aitchison (1983), Biometrika 70 (1) discusses the issues, & gives useful transformations to use for vectors of proportions, & you may be interested in the R packages compositions & robCompositions.

$\endgroup$
2
  • $\begingroup$ Okay thank you @Scortchi. I'm trying to take a closer look at the differences between the PCA on the data sets. I looked at the article link but think it is a beyond my knowledge level, could you explain a little bit more about the component of each vector summing to unity? I think this will help me understand why or why not proportions are appropriate way to standardize the data. Thank you. $\endgroup$
    – Oleic
    Jul 16, 2013 at 0:39
  • 1
    $\begingroup$ Just that if if you measure $A$ as say 50%, & $B$ as say 30%, then you know $C$ is 20% because the total must be 100%. So the last principal component will always have an eigenvalue of zero: you might as well just carry out PCA on A & B; & even then linear correlations are probably not so relevant to you. $\endgroup$ Jul 16, 2013 at 8:22

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.