data(coffee)
p1 <- pcaCoDa(coffee[,-1])
p1
biplot(p1)
## with labels for the scores:
data(arcticLake)
rownames(arcticLake) <- paste(sample(letters[1:26], nrow(arcticLake), replace=TRUE),
1:nrow(arcticLake), sep="")
pc <- pcaCoDa(arcticLake, method="classical")
biplot(pc, xlabs=rownames(arcticLake))
Run the code above in your browser using DataLab