无法手动计算PCA的分数，与princomp产生的分数相匹配

Question

使用此脚本：

d <- iris[1:4]
pca = princomp(d, cor=TRUE)

d2 <- scale(d, scale = FALSE)

m1 <- matrix(as.numeric(d[1,1:4]),1,4)
m2 <- matrix(pca$loadings,4,4)

mf1 <- m1 %*% m2
mf2 <- m2 %*% t(m1)

我无法获得矩阵乘法结果（mf1或mf2）以匹配运行pca$scores的结果。我的理解是，加载的值乘以得分，但我遗漏了一些东西（大概很简单）。

是否存在princomp适用的默认标准化？