Question

我有一个带有列标签的数据框*。例如，在使用prcomp时，如何使用或显示那些 labels 列而不是 names 列？（下面的功能get_label来自程序包sjlabelled。）

> get_label(spss.wvs) %>% head(2)
           wvs_e069_01                wvs_e069_02 
"Confidence: churches" "Confidence: armed forces"

（通过prcomp，我得到的是：

> pca.wvs$rotation %>% head(2)
                PC1        PC2
wvs_e069_01 -0.08513771 0.45688379 
wvs_e069_02 -0.23062304 0.05508813

我想要的是：

> pca.wvs$rotation %>% head(2)
                PC1        PC2
"Confidence: churches" -0.08513771 0.45688379 
"Confidence: armed forces" -0.23062304 0.05508813

*数据集是按照以下方式从spss（dropbox link to .sav）通过包haven导入的：

library(haven)
haven.imp <- read_spss("qog_std.sav")
library(dplyr)
spss.wvs <- haven.imp %>% dplyr:: select(starts_with("wvs_e069_0"), starts_with("wvs_e069_1"), grep("wvs_e069_20", names(haven.imp)))

（我确信有更有效的方法来获取“ wvs_e0690_00和wvs_e0690_20之间的所有内容”，但是至少可以做到这一点。指针值得赞赏。如您所知，我是{{3}的初学者}。

str(spss.wvs[1])的输出是：

Classes ‘tbl_df’, ‘tbl’ and 'data.frame':   193 obs. of  1 variable:
 $ wvs_e069_01: num  NA NA NA 3.17 NA ...
  ..- attr(*, "label")= chr "Confidence: churches"
  ..- attr(*, "format.spss")= chr "F8.2"

Answer 1

啊！我可以的：

> library(sjlabelled)
> rownames(spss.wvs$rotation) <- get_label(wvs)

哪个给我

> spss.wvs$rotation %>% head(2)
                                 PC1        PC2
Confidence: churches     -0.08513771 0.45688379
Confidence: armed forces -0.23062304 0.05508813

使用prcomp

1 个答案: