使用此脚本：
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适用的默认标准化？