0

我正在尝试使用 ggbiplot 绘制 PCA 分数，但由于我的分数和我的分组不匹配，我不能。我认为不匹配源于我的原始对数转换数据中的 NA 值（我在计算 PC 时忽略了它）。有没有办法解决这个问题，以便我可以使用 ggbiplot 进行绘图？

GCelem_trans.pca<-prcomp(na.omit(GCelem_trans.log), center=TRUE, scale=TRUE)
GCelem.rxtp <- GCelem[, 9]


g <- ggbiplot(GCelem_trans.pca, obs.scale = 1, var.scale = 1, 
                  groups = GCelem.rxtp, ellipse = TRUE, 
                  circle = TRUE)
    Error in `$<-.data.frame`(`*tmp*`, "groups", value = c(2L, 5L, 5L, 2L,  : 
      replacement has 33337 rows, data has 30804

我是否应该根据 Gcelem 的副本重新计算 Gcelem.rxtp，该副本忽略任何带有 NA 的行？

4

1 回答 1

0

您可以尝试：

GCelem.rxtp <- GCelem[complete.cases(GCelem_trans.log), 9]
于 2016-03-09T00:10:10.877 回答