我有一个如下数据框:
x <- c(1, 2, 1, 2)
y <- c(1, 2, 3, 4)
z <- c(4, 3, 2, 1)
df <- data.frame(x, y, z)
我正在使用fa
软件包中的psych
函数进行因子分析:
fit <- fa(df, nfactors = 2)
fit$loadings
这将产生以下输出:
Loadings:
MR1 MR2
x 0.448
y 0.999
z -0.999
MR1 MR2
SS loadings 2.195 0.000
Proportion Var 0.732 0.000
Cumulative Var 0.732 0.732
我想将带有MR1和MR2的表保存为数据帧。有谁知道该怎么做?谢谢。
答案 0 :(得分:1)
x <- c(1, 2, 1, 2)
y <- c(1, 2, 3, 4)
z <- c(4, 3, 2, 1)
df <- data.frame(x, y, z)
fit <- psych::fa(df, nfactors = 2)
x <- fit$loadings
通过stats:::print.loadings
:
Lambda <- unclass(x)
p <- nrow(Lambda)
factors <- ncol(Lambda)
vx <- colSums(x^2)
varex <- rbind(`SS loadings` = vx)
if (is.null(attr(x, "covariance"))) {
varex <- rbind(varex, `Proportion Var` = vx/p)
if (factors > 1)
varex <- rbind(varex, `Cumulative Var` = cumsum(vx/p))
}
tibble::rownames_to_column(as.data.frame(varex), "x")
## x MR1 MR2
## 1 SS loadings 2.1954555 3.000000e-30
## 2 Proportion Var 0.7318185 1.000000e-30
## 3 Cumulative Var 0.7318185 7.318185e-01
然后,对于第一个表格:
cutoff <- 0.1 # (the default for the `print.loadings()` function)
Lambda <- unclass(x)
p <- nrow(Lambda)
fx <- setNames(Lambda, NULL)
fx[abs(Lambda) < cutoff] <- NA_real_
fx <- as.data.frame(fx)
rownames(fx) <- NULL
fx
## MR1 MR2
## 1 0.4476761 NA
## 2 0.9987596 NA
## 3 -0.9987596 NA