R、 如何可视化PCA

我面临着一堵墙，我有我的数据，我想展示PCA。数据如下所示：

              mean_HL  mean_LL  mean_FL  mean_ML
 AT1G01080.1  6659961  2481626  3334466  1917721
 AT1G01090.1 14305565 15664978 16938777 13772388
 AT1G01320.1  2417028  7043427 15705767  6236322
 AT1G01790.1  1689920  4840314  6193037  1870497
 AT1G02140.1  5878921  3766975  3536914  3674864
 AT1G02150.1  3271086  2408324  1741409  4520847
 > 

我是这样计算的

tbl_log <- log(MERGE[, 1:4])
Light_mean_PCA <- prcomp(tbl_log,
                       center = TRUE,
                       scale. = TRUE)

我是这样做的。但是我想要颜色取决于：原始表格的平均值，平均值，平均值，平均值

---

在实验中，如何在四种不同的光照条件下得到四种不同的颜色

使用示例数据，可以绘制6个点。您希望如何使用四种条件为它们着色？除了使用ggplot，你还可以尝试使用base R

plot（Light\u mean\u PCA$x[，1:2]）

当我起诉Jimbou提供的代码时，我得到了一个这样的图，我想让它像这样，但要根据条件进行着色。所以也许我提供了这么小的示例。但这只是我的数据中的一个头部，所以可以从原始表格colnames中注释颜色吗？使用您的示例数据，可以绘制6个点。您希望如何使用四种条件为它们着色？除了使用ggplot，你还可以尝试使用base R

plot（Light\u mean\u PCA$x[，1:2]）

当我起诉Jimbou提供的代码时，我得到了一个这样的图，我想让它像这样，但要根据条件进行着色。所以也许我提供了这么小的示例。但这只是我数据中的一个标题，所以有可能从原始表格colnames中注释颜色吗？

library(devtools)
install_github("ggbiplot", "vqv")

library(ggbiplot)
g <- ggbiplot(Light_mean_PCA, obs.scale = 1, var.scale = 1, 
          groups = row.names(Light_mean_PCA), ellipse = TRUE, 
          circle = TRUE)
g <- g + scale_color_discrete(name = "")
g <- g + theme(legend.direction = 'horizontal', 
           legend.position = 'top')
print(g)

library(ggfortify); library(ggplot2)
autoplot(prcomp(tbl_log), data = MERGE, colour = MERGE[,c(1,2,3,4)])