PCA得分与Varimax旋转的PCA得分

Question

我已经在R中使用prcomp对75-76个指标变量数据库和7232个公司（包括NA）进行了PCA。在应用函数之前，我将数据居中，但没有重新缩放它们，因为它们都是指标变量。 （我的推理正确吗？）

此后，我按照 amoeba here的指示varimax旋转了2个或3个第一主成分的载荷。

由于我居中但未缩放数据，因此将代码更改为：

Varimax_results <- varimax(rawLoadings,normalize = FALSE)
invLoadings     <- t(pracma::pinv(VarimaxLoadings))
scores          <- scale(DatosPCA, scale = FALSE) %*% invLoadings

现在，我试图弄清楚为什么“ prcomp”给出的分数与使用上面的代码获得的分数为何不相同。

我可能缺少一些理论背景，因此，如果有人可以告诉我分数是否应该相同以及在那种情况下我在代码中做错了什么，我将不胜感激。如果不应该相同，应该使用哪个？

非常感谢您！