[英]R plotting PCA scores with ggbiplot
我正在尝试使用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
是否应该基于忽略了任何带有NA的行的GCelem副本重新计算GCelem.rxtp?
您可以尝试:
GCelem.rxtp <- GCelem[complete.cases(GCelem_trans.log), 9]
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.