Sự khác biệt giữa tóm tắt () và loadings () cho đối tượng Princeomp () trong R là gì?

Mã ví dụ:

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

Tôi đang nhận được kết quả đầu ra khác nhau từ mỗi loại, và tôi không chắc tôi hiểu sự khác biệt là gì.

Đây là đầu ra:

> 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: làm thế nào tôi có thể truy cập vào bảng được tạo bởi tóm tắt (pc.cr) ?? (Tôi dường như không thể tìm thấy nó trong str.)

Đầu ra đầu tiên là đầu ra chính xác và hữu ích nhất. Gọi loadings()đối tượng của bạn chỉ trả về một bản tóm tắt trong đó SS luôn bằng 1, do đó phương sai% chỉ là các tải SS chia cho số lượng biến. Nó chỉ có ý nghĩa khi sử dụng Phân tích nhân tố (như trong factanal). Tôi không bao giờ sử dụng princomphoặc nó thay thế SVD-based ( prcomp), và tôi thích FactoMineR hoặc ade4 gói mà là của xa mạnh mẽ hơn!

Về câu hỏi thứ hai của bạn, summary()hàm chỉ trả về SD cho từng thành phần ( pc.cr$sdevtrong trường hợp của bạn) và phần còn lại của bảng dường như được tính toán sau đó (thông qua phương thức printhoặc show, tôi không điều tra chi tiết điều này).

> 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>

Những gì princomp()bản thân nó có thể được xem bằng cách sử dụng getAnywhere("princomp.default").