无法在 R 中绘制 PCA 数据。分数是否由给定的 object/name 定义以专门绘制它们？

Question

我已经使用机构传下来的代码完成了一个简单的 PCA 函数。它输出分数、载荷、特征值、%特征值、主成分的数量、列的平均值、标准偏差，最后是起始数据。在输出文件中，分数在显示分数之前用 [[1]] 标记。我正在尝试绘制这些分数，但我不确定如何从这一点获取数据。我假设它被分配给这个 [[1]] 或代码中定义这些分数的东西。这行代码如下所示：

 "#"perform pca on x
 x.svd <- svd(x);
 x.R <- x.svd$u %*% diag(x.svd$d);
 x.C <- t(x.svd$v);
 x.EV <- x.svd$d * x.svd$d
 x.EVpct <- x.EV/sum(x.EV);
 x.EV <- x.EV[1:sm];
 x.EVpct <- x.EVpct[1:sm];
 x.CumEVpct <- x.EVpct;

x.R 是执行分数的代码部分，但它也不适用于 plot 函数。希望有人明白我在努力问什么。非常感谢任何帮助。谢谢你的时间。

Answer 1

最简单的做法是：

pc <- prcomp(x)
plot(pc$x[, 1:2]

无法在 R 中绘制 PCA 数据。分数是否由给定的 object/name 定义以专门绘制它们？

Unable to plot PCA data in R. Are scores defined by a given object/name to plot them specifically?

plot

r

object

pca