r - PCA 分数与 Varimax 旋转 PCA 分数

Question

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

之后，我按照amoeba herevarimax的说明旋转了 2 或 3 个第一主成分的负载。

由于我已经居中，但没有重新调整我的数据，我将代码更改为：

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

现在我想弄清楚为什么“prcomp”给出的分数和使用上面代码获得的分数不一样。

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

非常感谢！