R中princomp（）对象的summary（）和loadings（）之间有什么区别？

示例代码：

(pc.cr <- princomp(USArrests))  
summary(pc.cr)
loadings(pc.cr)  ## note that blank entries are small but not zero

我从每个人那里得到了不同的输出，而且我不确定我理解有什么不同。

这是输出：

> summary(pc.cr)
Importance of components:
                           Comp.1      Comp.2      Comp.3       Comp.4
Standard deviation     82.8908472 14.06956001 6.424204055 2.4578367034
Proportion of Variance  0.9655342  0.02781734 0.005799535 0.0008489079
Cumulative Proportion   0.9655342  0.99335156 0.999151092 1.0000000000


> loadings(pc.cr)  ## note that blank entries are small but not zero

...

               Comp.1 Comp.2 Comp.3 Comp.4
SS loadings      1.00   1.00   1.00   1.00
Proportion Var   0.25   0.25   0.25   0.25
Cumulative Var   0.25   0.50   0.75   1.00

PS：如何访问由summary（pc.cr）创建的表？（我似乎在str中找不到它。）

第一个输出是正确且最有用的输出。调用loadings()您的对象只会返回一个汇总，其中SS始终等于1，因此％方差就是SS负载除以变量数量。仅当使用因素分析时（如中的factanal），它才有意义。我从不使用princomp它或基于SVD的替代方法（prcomp），而我更喜欢FactoMineR或ade4软件包，它们的功能更强大！

关于您的第二个问题，该summary()函数仅返回每个组件的SD（pc.cr$sdev在您的情况下），该表的其余部分似乎是在以后计算的（通过printor show方法，我没有对此进行详细调查）。

> getS3method("summary","princomp")
function (object, loadings = FALSE, cutoff = 0.1, ...)
{
    object$cutoff <- cutoff
    object$print.loadings <- loadings
    class(object) <- "summary.princomp"
    object
}
<environment: namespace:stats>

princomp()使用可以查看本身所做的事情getAnywhere("princomp.default")。