R中的主成分分析中的内积

时间:2018-10-29 11:48:10

标签: r pca

为此,我使用data(banknote)给出的R中的钞票数据,该数据显示200张瑞士钞票的尺寸。我的数据矩阵称为X,我已经通过pca.banknote<-prcomp(X)执行了PCA。

我试图证明X[i,]给出的每个观察值pca.banknote$rot[,3]与主成分负荷3之间的内积与pca.banknote$x[,3]给出的第三次PC分数相同。

我尝试过:

all.equal(as.matrix(X[,])%*%banknote.pca$rot[,3], as.matrix(banknote.pca$x[,3]), check.attributes=FALSE)

但这只是平均差为1,即它们不相等。

我是否需要将其中之一的格式更改为矢量/数据帧等才能正常工作?还是对问题出在哪里有任何想法?

任何反馈将不胜感激。谢谢。

0 个答案:

没有答案