为此,我使用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,即它们不相等。
我是否需要将其中之一的格式更改为矢量/数据帧等才能正常工作?还是对问题出在哪里有任何想法?
任何反馈将不胜感激。谢谢。