https://cran.r-project.org/web/packages/ggfortify/vignettes/plot_pca.html

> install.packages("ggfortify")

> library（“ggfortify")

{ggfortify}let{ggplot2}know how to interpret PCA objects. After loading{ggfortify}, you can useggplot2::autoplot function for stats::prcomp and stats::princomp objects.

讀入文件中的數(shù)據(jù)

> df<-read.csv("pom.csv",header=T)

> df

> autoplot(prcomp(df))?????? #出圖

圖形編輯

>help(autoplot.prcomp)

> autoplot(prcomp(df), data = df, colour = 'code')

> autoplot(prcomp(df), data = df, colour = 'code', label = TRUE, label.size = 3)

autoplot(prcomp(df), data = df, colour = 'code',

+? loadings = TRUE, loadings.colour = 'blue',

+ loadings.label = TRUE, loadings.label.size = 3)