Làm cách nào để xuất bảng xếp hạng yếu tố được sắp xếp?

Question

nét Vấn đề: xuất khẩu một phiên bản sắp xếp của một đối tượng của lớp "tải trọng"

Sau khi chạy một phân tích yếu tố với fa -function của psych gói, tôi nhận được một bảng tải trọng yếu tố đó trông giống như một trong những đưa ra ở đây:

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   

Sau khi chạy chức năng print trên bảng này (được lưu trữ trong biến f.loadings), tôi nhận được một bảng sắp xếp 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 tuy nhiên trả lại bản sao "ẩn" của đối tượng, vì vậy tôi không thể xuất kết quả này theo định dạng được yêu cầu. Tuy nhiên, tôi muốn tìm cách xuất phiên bản .csv của bảng này.

Tôi không thể tìm cách xác định tham số write.csv để thực hiện sắp xếp chính xác đối tượng của lớp "đang tải". Việc gán kết quả của hàm in cũng không giải quyết được điều này, vì nó chỉ trả về phiên bản chưa được phân loại. Do đó x <- print(f.loadings, digits=2, cutoff=.3, sort=TRUE) và sau đó gọi cho biến mới x, vẫn trả về phiên bản chưa phân loại của bảng.

Chức năng nào sẽ phù hợp để phân loại đối tượng "nạp" và trả lại đối tượng này một cách rõ ràng? Nói cách khác, làm thế nào tôi có thể xuất một bảng được sắp xếp như vậy?

Mã để tạo bảng:

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

Dưới đây là cách bạn có thể hack phương thức in mặc định cho loadings. Tôi đã gán câu lệnh in cho newx và xuất khẩu. Khi bạn gán kết quả printLoadings cho một biến, bạn có đối tượng được sắp xếp theo ý của mình. Kết quả bây giờ là một bảng các ký tự. Tôi sẽ để nó cho bạn như một bài tập để chuyển nó sang dạng số.

> 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

Khi một hàm trả về kết quả của nó với invisible() nó chỉ đơn giản có nghĩa là kết quả không được in. Tuy nhiên, bạn vẫn có thể gán kết quả cho một biến và thao tác nó thành bất kỳ đối tượng nào khác.

Vì vậy:

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 

Tương tự, str(x) chỉ ra kết quả là một ma trận:

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" 

Bây giờ bạn có thể sử dụng write.csv(x) export các kết quả.

Answer 3

Thử chức năng fa.sort thay vì sắp xếp chức năng in.

Nghĩa là, nếu f < - fa (Thurstone, 3) bạn có thể in (f, sắp xếp = TRUE) mà không làm bản in nhưng không trả lại đối tượng được sắp xếp, hoặc fs < - fa. sắp xếp (f) trả về đối tượng đã sắp xếp mà bạn có thể xuất dưới dạng csv.

Bill