在主成分分析中获取 "PC1" insted 变量名

Question

我有一些数据如下所示：

head(data)
   net1re net2re net3re net4re net5re net6re
24      3      2      1      2      3      3
33      1      1      1      1      1      2
30      3      1      1      1      1      3
22      2      1      1      1      1      1
31      3      2      1      1      1      2
1       2      1      1      1      1      2

我运行主成分分析如下：

library(psych)
fit <- principal(data[,1:6], rotate="varimax")
data$friendship=fit$scores

这将创建我可以在控制台上调用的变量 "friendship"：

> colnames(data)
[1] "net1re"     "net2re"     "net3re"     "net4re"     "net5re"    
[6] "net6re"     "friendship"

但是当我想查看我的数据时，我得到的不是变量名 "PC1":

> head(data)
   net1re net2re net3re net4re net5re net6re         PC1
24      3      2      1      2      3      3  1.29231531
33      1      1      1      1      1      2 -0.68448111
30      3      1      1      1      1      3  0.02783916
22      2      1      1      1      1      1 -0.67371031
31      3      2      1      1      1      2  0.10251282
1       2      1      1      1      1      2 -0.44345075

这变成了一个大麻烦，因为我需要用不同的变量重复它，所有的结果都会得到 "PC1"。

为什么会这样，我如何分配变量名而不是 "PC1"。

谢谢

Answer 1

出现这种不寻常的效果是因为fit$scores是一个矩阵：

str(data)
#'data.frame':  6 obs. of  7 variables:
# $ net1re    : int  3 1 3 2 3 2
# $ net2re    : int  2 1 1 1 2 1
# $ net3re    : int  1 1 1 1 1 1
# $ net4re    : int  2 1 1 1 1 1
# $ net5re    : int  3 1 1 1 1 1
# $ net6re    : int  3 2 3 1 2 2
# $ friendship: num [1:6, 1] 1.1664 -1.261 0.0946 -0.5832 1.1664 ...
#  ..- attr(*, "dimnames")=List of 2
#  .. ..$ : chr  "24" "33" "30" "22" ...
#  .. ..$ : chr "PC1"

要获得想要的结果，您可以使用

data$friendship=as.vector(fit$scores)

或

data$friendship=fit$scores[,1]

无论哪种情况，输出都将是：

data
#   net1re net2re net3re net4re net5re net6re  friendship
#24      3      2      1      2      3      3  1.16635312
#33      1      1      1      1      1      2 -1.26098965
#30      3      1      1      1      1      3  0.09463653

str(data)
#'data.frame':  6 obs. of  7 variables:
# $ net1re    : int  3 1 3 2 3 2
# $ net2re    : int  2 1 1 1 2 1
# $ net3re    : int  1 1 1 1 1 1
# $ net4re    : int  2 1 1 1 1 1
# $ net5re    : int  3 1 1 1 1 1
# $ net6re    : int  3 2 3 1 2 2
# $ friendship: num  1.1664 -1.261 0.0946 -0.5832 1.1664 ...

在主成分分析中获取 "PC1" insted 变量名

getting "PC1" insted of variable name in principal component analysis

r

pca

psych