

library(ade4)
# aravo$env 含定量變數和分類變數的環境測量資料
data(aravo)
head(aravo$env)
#環境資料的 PCA,結果中保留前 4 個排序軸的座標便於觀測資料
pca <- dudi.hillsmith(aravo$env, scannf = FALSE, nf = 4)
#PCA 概要
pca
summary(pca)
#展示前兩軸的排序圖,以及特徵值柱形圖
scatter(pca, choices = 1:2)
#或者
biplot(pca, choices = 1:2)
#提取結果,例如
#names(pca)
pca$l1 #標準化後的樣方得分(排序座標)
pca$c1 #標準化後的環境變數得分(排序座標)
pca$eig #各 PCA 軸特徵值
pca$eig/sum(pca$eig) #各 PCA 軸的貢獻度

參考資料
Hill, M. O., and A. J. E. Smith. 1976. Principal component analysis of taxonomic data with multi-state discrete characters. Taxon, 25, 249-255.



