Question

问题定义：导出类“加载”对象的已排序版本

使用psych-package的fa函数运行因子分析后，我得到一个因子加载表，看起来像这里显示的那样：

Loadings:
         Factor1 Factor2 Factor3
TH_Q1     0.173   0.548   0.403 
TH_Q2     0.306   0.291   0.825 
TH_Q3     0.334   0.203   0.825 
TH_Q4     0.262   0.536   0.171 
TH_Q5     0.235   0.686         
TH_Q6     0.125   0.836         
TH_Q7     0.200   0.838         
TH_Q8_A1                        
TH_Q8_A2          0.155         
TH_Q9     0.644   0.133   0.171 
TH_Q10    0.608   0.208   0.157 
TH_Q11    0.569   0.161   0.306 
TH_Q12    0.722           0.127 
TH_Q13    0.661   0.311         
TH_Q14    0.562   0.407         
TH_Q15    0.675   0.422

在此表上运行函数print（存储在变量 f.loadings 中）之后，我得到一个已排序的表print(f.loadings, digits=2, cutoff=.3, sort=TRUE)：

Loadings:
         Factor1 Factor2 Factor3
TH_Q9     0.64                  
TH_Q10    0.61                  
TH_Q11    0.57            0.31  
TH_Q12    0.72                  
TH_Q13    0.66    0.31          
TH_Q14    0.56    0.41          
TH_Q15    0.68    0.42          
TH_Q1             0.55    0.40  
TH_Q4             0.54          
TH_Q5             0.69          
TH_Q6             0.84          
TH_Q7             0.84          
TH_Q2     0.31            0.82  
TH_Q3     0.33            0.83  
TH_Q8_A1                        
TH_Q8_A2

然而，

print会返回对象的“不可见”副本，因此我无法以请求的格式导出此结果。但是，我希望找到一种方法来导出此表的.csv版本。

我无法找到指定write.csv参数的方法来对“loading”类对象进行正确排序。分配打印功能的结果也不能解决这个问题，因为它只返回未排序的版本。因此x <- print(f.loadings, digits=2, cutoff=.3, sort=TRUE)并随后调用新变量x仍然返回表的未排序版本。

什么函数适合于对“加载”对象进行排序并明显返回此对象？换句话说，我如何导出这样的排序表？

生成表格的代码：

f.loadings <- structure(c(0.172693322885797, 0.306277415972136, 0.334012445825371, 
0.261822356615649, 0.234600824098634, 0.124541887813939, 0.200125976802047, 
0.0199775267669519, 0.0771905784767979, 0.643886342785064, 0.608004298828405, 
0.569498016145868, 0.722454442131503, 0.660683752725898, 0.561975379133291, 
0.675119271585253, 0.548184083921831, 0.291215413974386, 0.20334622551054, 
0.535545380240845, 0.685635981787823, 0.836401389336655, 0.837525597359627, 
0.0186113870539496, 0.154659865540958, 0.132908227837058, 0.20832344061795, 
0.160657979843522, 0.0933961709813049, 0.311465272208257, 0.406860675137862, 
0.421946817384512, 0.402664774610544, 0.824934582975472, 0.825220077707656, 
0.170809720550637, -0.0486225264368695, 0.0612401518170266, 0.052596915030506, 
-0.0463868732056794, 0.0208945338424677, 0.171412077700389, 0.156524506151013, 
0.306203004564158, 0.127377474768802, -0.0869197819037828, -0.0962274476959987, 
-0.0465278761105364), .Dim = c(16L, 3L), .Dimnames = list(c("TH_Q1",  "TH_Q2", "TH_Q3", "TH_Q4", "TH_Q5", "TH_Q6", "TH_Q7", "TH_Q8_A1",  "TH_Q8_A2", "TH_Q9", "TH_Q10", "TH_Q11", "TH_Q12", "TH_Q13",  "TH_Q14", "TH_Q15"), c("Factor1", "Factor2", "Factor3")), class = "loadings")

Answer 1

以下是如何破解loadings的默认打印方法。我已将print语句分配给newx并将其导出。将printLoadings结果分配给变量时，您可以使用已排序的对象。结果现在是一个字符表。我将把它留给你作为练习将其转换为数字。

> getS3method("print","loadings") #get the hidden method and modify it
printLoadings <- function (x, digits = 3, cutoff = 0.1, sort = FALSE, ...) 
{
   Lambda <- unclass(x)
   p <- nrow(Lambda)
   factors <- ncol(Lambda)
   if (sort) {
      mx <- max.col(abs(Lambda))
      ind <- cbind(1L:p, mx)
      mx[abs(Lambda[ind]) < 0.5] <- factors + 1
      Lambda <- Lambda[order(mx, 1L:p), ]
   }
   cat("\nLoadings:\n")
   fx <- format(round(Lambda, digits))
   names(fx) <- NULL
   nc <- nchar(fx[1L], type = "c")
   fx[abs(Lambda) < cutoff] <- paste(rep(" ", nc), collapse = "")
   newx <- print(fx, quote = FALSE, ...) # I assigned this to a variable
   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))
   }
   cat("\n")
   print(round(varex, digits))
   invisible(newx) #previously returned x
}

mmm <- printLoadings(f.loadings)
> str(mmm)
 chr [1:16, 1:3] " 0.173" " 0.306" " 0.334" " 0.262" " 0.235" ...
 - attr(*, "dimnames")=List of 2
  ..$ : chr [1:16] "TH_Q1" "TH_Q2" "TH_Q3" "TH_Q4" ...
  ..$ : chr [1:3] "Factor1" "Factor2" "Factor3"

> as.table(mmm)
         Factor1 Factor2 Factor3
TH_Q1     0.173   0.548   0.403 
TH_Q2     0.306   0.291   0.825 
TH_Q3     0.334   0.203   0.825 
TH_Q4     0.262   0.536   0.171

Answer 2

当函数使用invisible()返回其结果时，它只表示不打印结果。但是，您仍然可以将结果分配给变量并将其作为任何其他对象进行操作。

所以：

x <- print(f.loadings)
x
Loadings:
         Factor1 Factor2 Factor3
TH_Q1     0.173   0.548   0.403 
TH_Q2     0.306   0.291   0.825 
TH_Q3     0.334   0.203   0.825 
TH_Q4     0.262   0.536   0.171 
TH_Q5     0.235   0.686         
TH_Q6     0.125   0.836         
TH_Q7     0.200   0.838         
TH_Q8_A1                        
TH_Q8_A2          0.155         
TH_Q9     0.644   0.133   0.171 
TH_Q10    0.608   0.208   0.157 
TH_Q11    0.569   0.161   0.306 
TH_Q12    0.722           0.127 
TH_Q13    0.661   0.311         
TH_Q14    0.562   0.407         
TH_Q15    0.675   0.422         

               Factor1 Factor2 Factor3
SS loadings      3.259   3.145   1.747
Proportion Var   0.204   0.197   0.109
Cumulative Var   0.204   0.400   0.509

相似，str(x)表示结果是矩阵：

str(x)
 loadings [1:16, 1:3] 0.173 0.306 0.334 0.262 0.235 ...
 - attr(*, "dimnames")=List of 2
  ..$ : chr [1:16] "TH_Q1" "TH_Q2" "TH_Q3" "TH_Q4" ...
  ..$ : chr [1:3] "Factor1" "Factor2" "Factor3"

您现在可以使用write.csv(x)导出结果。

Answer 3

尝试使用fa.sort函数而不是在print函数中进行排序。

也就是说，如果 f＆lt; - fa（瑟斯顿，3）你可以打印（f，sorted = TRUE），它不打印但不返回已排序的对象，或者 fs＆lt; - fa.sort（f）返回已排序的对象，然后可以将其导出为csv。

比尔

如何导出已排序的因子加载表？

3 个答案: