获取R中SVD的列名

Question

我在R中对一个名为data的数据帧执行SVD

svd1 <- svd(scale(data))

我使用以下命令绘制结果：

plot(svd1$d^2/sum(svd1$d^2),xlab="Column",ylab="Percent of variance explained",pch=19)

我发现在图中，列号并不与数据帧中的列号相对应(无论我在SVD中使用哪个列子集，第一列总是显示最高方差)。

我的问题是，如何在图中获得列名(或“真实”索引)？

Answer 1

根据R文档(.../library/base/html/svd.html)：

d is a vector containing the singular values of x, of length min(n, p).

这些奇异值是通过SVD算法从整个输入矩阵中计算出来的，因此无法根据列名来标记奇异值。