R의 princomp () 객체에 대한 summary ()와 loading ()의 차이점은 무엇입니까?

예제 코드 :

(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

추신 : 요약 (pc.cr)으로 만든 테이블에 어떻게 액세스 할 수 있습니까 ?? (나는 str에서 그것을 찾을 수없는 것 같습니다.)

첫 번째 출력은 정확하고 가장 유용한 출력입니다. loadings()객체를 호출 하면 SS가 항상 1과 같은 요약이 반환되므로 % 분산은 SS 하중을 변수 수로 나눈 값입니다. 요인 분석을 사용할 때만 (예 :와 같이 factanal) 의미가 있습니다. 나는 princompSVD 기반 대안 ( prcomp)을 사용하지 않으며 훨씬 더 강력한 FactoMineR 또는 ade4 패키지를 선호합니다 !

두 번째 질문에 대해 summary()함수는 각 구성 요소 ( pc.cr$sdev귀하의 경우) 에 대해 SD를 반환하고 나머지 테이블은 나중에 계산됩니다 ( print또는 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").