Извлечение значений (недоступно через «str») из объекта в R

Question

У меня проблема с извлечением значения из объекта PCA. Мне нужно «Пропорция разницы» значение от «Comp.1». В этом примере это значение 0.7056111.

Я могу видеть его в резюме, но не может извлечь его, потому что его нет в str(pca_test).

pca_test<-data.frame(var1=rnorm(20), var2=rnorm(20)) 
summary(princomp(pca_test)) 

Importance of components: 
          Comp.1 Comp.2 
Standard deviation  1.0200426 0.6588649 
Proportion of Variance 0.7056111 0.2943889 
Cumulative Proportion 0.7056111 1.0000000 


# value 0.7056111 is absent here 
> str(summary(princomp(pca_test))) 
List of 9 
$ sdev   : Named num [1:2] 1.02 0.659 
    ..- attr(*, "names")= chr [1:2] "Comp.1" "Comp.2" 
$ loadings  : loadings [1:2, 1:2] 0.289 0.957 -0.957 0.289 
    ..- attr(*, "dimnames")=List of 2 
    .. ..$ : chr [1:2] "var1" "var2" 
    .. ..$ : chr [1:2] "Comp.1" "Comp.2" 
$ center  : Named num [1:2] 0.1497 -0.0897 
    ..- attr(*, "names")= chr [1:2] "var1" "var2" 
$ scale   : Named num [1:2] 1 1 
    ..- attr(*, "names")= chr [1:2] "var1" "var2" 
$ n.obs   : int 20 
$ scores  : num [1:20, 1:2] -1.8965 1.4866 -1.5019 0.0841 0.4751 ... 
    ..- attr(*, "dimnames")=List of 2 
    .. ..$ : NULL 
    .. ..$ : chr [1:2] "Comp.1" "Comp.2" 
$ call   : language princomp(x = pca_test) 
$ cutoff  : num 0.1 
$ print.loadings: logi FALSE 
- attr(*, "class")= chr "summary.princomp" 

Answer 1

Но как agstudy найти его? Вот мой общий подход, чтобы узнать, как некоторый параметр кодируется

pca_test<-data.frame(var1=rnorm(20), var2=rnorm(20)) 
pc = princomp(pca_test) 
str(pc) # proportion not found 
sp = summary(pc) 
str(sp) # The proportion is still not there 
# It is of class summary.princomp 
# So it probably comes with the print 
# This is naughty behavioR, the print function should not compute, 
# but rather focus on display 
getAnywhere(print.summary.princomp) 
# There it is, at the top 
#vars <- x$sdev^2 
#vars <- vars/sum(vars) 
#cat("Importance of components:\n") 
#print(rbind(`Standard deviation` = x$sdev, `Proportion of Variance` = vars, 

Answer 2

Вы можете получить компоненты пропорция, как это:

object <- princomp(pca_test) 
vars <- object$sdev^2 
vars/sum(vars)