6

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

使用 psych-package的 -function 运行因子分析后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以对“加载”类的对象进行正确排序。分配 print 函数的结果也不能解决这个问题，因为它只返回未排序的版本。因此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")
4

3 回答 3

5

以下是如何破解loadings. 我已将打印语句分配给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
于 2012-08-24T19:59:39.317 回答
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)来导出结果。

于 2012-08-24T18:57:46.150 回答
0

尝试 fa.sort 函数，而不是在 print 函数中排序。

也就是说，如果 f <- fa(Thurstone,3) 您可以 print(f,sorted=TRUE) 不打印但不返回排序的对象，或者 fs <- fa.sort(f) 返回排序的对象，然后您可以将其导出为 csv。

账单

于 2013-06-09T18:10:05.917 回答