I'm trying to perform a principal component analysis using more than one dataset of compositional data. I have two sets of data:
data(expenditure)
data(wellbeing)
Usually to calculate pca on a single composition I use the compositions
package:
library(compositions)
#conversion into acomp object
cdata <- acomp(expenditure)
#principal component analysis
pca <- princomp(cdata)
pca
But I can't figure out how to use also thewellbeing
compositional data to perform a pca along with expenditure
data.
Any idea or suggestions would be appreciated, thanks!