R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list('Case1' + ,2 + ,2 + ,2 + ,4 + ,2 + ,5 + ,4 + ,5 + ,2 + ,2 + ,4 + ,5 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,'Case2' + ,2 + ,2 + ,2 + ,2 + ,1 + ,4 + ,2 + ,5 + ,2 + ,2 + ,4 + ,2 + ,0 + ,1 + ,1 + ,1 + ,1 + ,0 + ,0 + ,0 + ,1 + ,1 + ,1 + ,1 + ,0 + ,'Case3' + ,2 + ,2 + ,2 + ,5 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,5 + ,5 + ,1 + ,1 + ,1 + ,0 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,'Case4' + ,1 + ,1 + ,1 + ,4 + ,1 + ,4 + ,4 + ,4 + ,2 + ,2 + ,4 + ,4 + ,0 + ,0 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,0 + ,1 + ,1 + ,0 + ,'Case5' + ,1 + ,1 + ,1 + ,2 + ,1 + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,4 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,'Case6' + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,2 + ,1 + ,1 + ,2 + ,4 + ,0 + ,1 + ,1 + ,1 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,1 + ,1 + ,1 + ,'Case7' + ,3 + ,4 + ,2 + ,NA + ,2 + ,4 + ,NA + ,4 + ,4 + ,4 + ,4 + ,NA + ,NA + ,NA + ,NA + ,0 + ,1 + ,1 + ,0 + ,0 + ,1 + ,1 + ,NA + ,1 + ,NA + ,'Case8' + ,2 + ,2 + ,2 + ,2 + ,NA + ,4 + ,4 + ,4 + ,NA + ,2 + ,4 + ,NA + ,NA + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,1 + ,1 + ,0 + ,1 + ,1 + ,'Case9' + ,1 + ,2 + ,1 + ,2 + ,1 + ,2 + ,2 + ,4 + ,2 + ,4 + ,4 + ,4 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,1 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,'Case10' + ,4 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,NA + ,2 + ,4 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,'Case12' + ,1 + ,1 + ,1 + ,4 + ,1 + ,2 + ,2 + ,1 + ,1 + ,1 + ,4 + ,5 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,NA + ,0 + ,0 + ,1 + ,0 + ,0 + ,'Case13' + ,2 + ,4 + ,4 + ,5 + ,2 + ,4 + ,4 + ,2 + ,4 + ,2 + ,4 + ,4 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,1 + ,1 + ,0 + ,1 + ,0 + ,0 + ,0 + ,'Case14' + ,2 + ,1 + ,1 + ,1 + ,1 + ,2 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,'Case15' + ,4 + ,2 + ,2 + ,2 + ,2 + ,4 + ,2 + ,2 + ,4 + ,4 + ,4 + ,4 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,'Case16' + ,2 + ,2 + ,4 + ,2 + ,2 + ,4 + ,NA + ,4 + ,4 + ,4 + ,4 + ,NA + ,1 + ,1 + ,1 + ,0 + ,1 + ,1 + ,1 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,'Case17' + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,4 + ,4 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,'Case18' + ,2 + ,2 + ,1 + ,2 + ,1 + ,1 + ,2 + ,1 + ,1 + ,1 + ,2 + ,1 + ,1 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,1 + ,0 + ,1 + ,0 + ,1 + ,'Case19' + ,1 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,1 + ,1 + ,1 + ,0 + ,0 + ,1 + ,1 + ,1 + ,1 + ,0 + ,1 + ,1 + ,1 + ,'Case20' + ,1 + ,2 + ,2 + ,2 + ,1 + ,2 + ,2 + ,4 + ,2 + ,2 + ,4 + ,4 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,'Case21' + ,2 + ,4 + ,1 + ,1 + ,1 + ,2 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,0 + ,0 + ,1 + ,1 + ,'Case22' + ,2 + ,2 + ,4 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,4 + ,0 + ,1 + ,1 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,'Case23' + ,2 + ,2 + ,2 + ,4 + ,2 + ,4 + ,NA + ,2 + ,2 + ,4 + ,4 + ,4 + ,1 + ,1 + ,1 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,1 + ,0 + ,'Case24' + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,2 + ,2 + ,2 + ,2 + ,4 + ,2 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,0 + ,0 + ,0 + ,1 + ,1 + ,1 + ,'Case25' + ,2 + ,2 + ,1 + ,2 + ,2 + ,4 + ,2 + ,1 + ,1 + ,2 + ,4 + ,4 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,'Case26' + ,2 + ,2 + ,2 + ,4 + ,4 + ,4 + ,4 + ,2 + ,2 + ,2 + ,4 + ,4 + ,1 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,1 + ,0 + ,'Case27' + ,2 + ,2 + ,2 + ,4 + ,4 + ,4 + ,4 + ,2 + ,2 + ,2 + ,4 + ,4 + ,1 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,NA + ,1 + ,0 + ,'Case28' + ,4 + ,2 + ,2 + ,5 + ,2 + ,5 + ,4 + ,4 + ,4 + ,5 + ,4 + ,4 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,'Case29' + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,1 + ,0 + ,'Case30' + ,2 + ,2 + ,1 + ,2 + ,1 + ,2 + ,2 + ,2 + ,4 + ,2 + ,4 + ,4 + ,0 + ,0 + ,1 + ,0 + ,1 + ,0 + ,1 + ,0 + ,1 + ,0 + ,1 + ,0 + ,0 + ,'Case31' + ,4 + ,5 + ,2 + ,4 + ,2 + ,4 + ,2 + ,4 + ,4 + ,2 + ,4 + ,4 + ,1 + ,1 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,1 + ,0 + ,'Case32' + ,4 + ,4 + ,4 + ,4 + ,2 + ,4 + ,2 + ,4 + ,4 + ,4 + ,4 + ,2 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,1 + ,0 + ,'Case33' + ,4 + ,2 + ,2 + ,2 + ,2 + ,4 + ,4 + ,4 + ,4 + ,4 + ,5 + ,4 + ,0 + ,0 + ,0 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,'Case34' + ,2 + ,2 + ,2 + ,3 + ,2 + ,4 + ,2 + ,4 + ,3 + ,3 + ,4 + ,5 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,1 + ,1 + ,1 + ,1 + ,1 + ,0 + ,'Case35' + ,2 + ,NA + ,2 + ,NA + ,2 + ,4 + ,4 + ,4 + ,2 + ,4 + ,5 + ,4 + ,1 + ,1 + ,1 + ,0 + ,1 + ,0 + ,1 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,'Case36' + ,4 + ,2 + ,4 + ,NA + ,4 + ,2 + ,5 + ,4 + ,4 + ,4 + ,5 + ,4 + ,0 + ,1 + ,1 + ,0 + ,1 + ,1 + ,NA + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,'Case37' + ,1 + ,4 + ,2 + ,4 + ,2 + ,4 + ,5 + ,2 + ,2 + ,2 + ,4 + ,4 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,1 + ,0 + ,0 + ,0 + ,1 + ,0 + ,'Case38' + ,4 + ,2 + ,4 + ,2 + ,2 + ,4 + ,2 + ,5 + ,4 + ,4 + ,5 + ,5 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,'Case39' + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,4 + ,1 + ,1 + ,4 + ,4 + ,0 + ,1 + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,NA + ,'Case40' + ,2 + ,2 + ,4 + ,NA + ,5 + ,NA + ,4 + ,5 + ,2 + ,2 + ,4 + ,4 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,NA + ,NA + ,1 + ,0 + ,'Case41' + ,2 + ,4 + ,2 + ,2 + ,NA + ,2 + ,2 + ,4 + ,4 + ,4 + ,4 + ,4 + ,0 + ,NA + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,0 + ,NA + ,0 + ,1 + ,0 + ,'Case42' + ,2 + ,4 + ,4 + ,2 + ,4 + ,4 + ,4 + ,2 + ,4 + ,4 + ,4 + ,4 + ,0 + ,0 + ,0 + ,0 + ,1 + ,1 + ,0 + ,0 + ,1 + ,0 + ,0 + ,1 + ,0 + ,'Case43' + ,4 + ,4 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,1 + ,1 + ,1 + ,1 + ,0 + ,1 + ,1 + ,1 + ,1 + ,0 + ,1 + ,1 + ,1 + ,'Case44' + ,2 + ,2 + ,4 + ,2 + ,2 + ,4 + ,2 + ,4 + ,2 + ,2 + ,4 + ,2 + ,0 + ,0 + ,1 + ,1 + ,0 + ,1 + ,1 + ,0 + ,1 + ,1 + ,1 + ,0 + ,0 + ,'Case45' + ,4 + ,4 + ,2 + ,4 + ,1 + ,5 + ,4 + ,5 + ,2 + ,4 + ,5 + ,5 + ,0 + ,0 + ,1 + ,0 + ,1 + ,1 + ,0 + ,0 + ,1 + ,1 + ,0 + ,1 + ,0) + ,dim=c(26 + ,44) + ,dimnames=list(c('label' + ,'1a' + ,'1b' + ,'1c' + ,'1d' + ,'1e' + ,'1f' + ,'1g' + ,'1h' + ,'1i' + ,'1j' + ,'1k' + ,'1l' + ,'2' + ,'3' + ,'4a' + ,'6' + ,'7' + ,'8' + ,'9' + ,'10' + ,'11' + ,'12' + ,'13' + ,'14' + ,'15') + ,1:44)) > y <- array(NA,dim=c(26,44),dimnames=list(c('label','1a','1b','1c','1d','1e','1f','1g','1h','1i','1j','1k','1l','2','3','4a','6','7','8','9','10','11','12','13','14','15'),1:44)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } There were 44 warnings (use warnings() to see them) > par1 = '4' > par1 <- '4' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa P., 2012, Factor Analysis (v1.0.2) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_factor_analysis.wasp/ > #Source of accompanying publication: > # > library(psych) > par1 <- as.numeric(par1) > x <- t(x) > nrows <- length(x[,1]) > ncols <- length(x[1,]) > y <- array(as.double(x[1:nrows,2:ncols]),dim=c(nrows,ncols-1)) > colnames(y) <- colnames(x)[2:ncols] > rownames(y) <- x[,1] > y 1a 1b 1c 1d 1e 1f 1g 1h 1i 1j 1k 1l 2 3 4a 6 7 8 9 10 11 12 13 14 Case1 2 2 2 4 2 5 4 5 2 2 4 5 0 0 0 0 0 1 1 0 1 0 1 1 Case2 2 2 2 2 1 4 2 5 2 2 4 2 0 1 1 1 1 0 0 0 1 1 1 1 Case3 2 2 2 5 2 2 2 2 2 2 5 5 1 1 1 0 1 1 1 1 1 1 1 1 Case4 1 1 1 4 1 4 4 4 2 2 4 4 0 0 1 1 1 1 1 1 1 0 1 1 Case5 1 1 1 2 1 2 2 2 2 2 4 4 0 1 0 0 0 0 0 0 1 0 0 1 Case6 1 1 1 2 2 4 2 2 1 1 2 4 0 1 1 1 0 1 0 0 0 1 1 1 Case7 3 4 2 NA 2 4 NA 4 4 4 4 NA NA NA NA 0 1 1 0 0 1 1 NA 1 Case8 2 2 2 2 NA 4 4 4 NA 2 4 NA NA 0 0 1 0 1 1 0 1 1 0 1 Case9 1 2 1 2 1 2 2 4 2 4 4 4 0 0 1 0 0 1 1 0 1 0 1 1 Case10 4 2 2 2 2 2 2 2 4 NA 2 4 0 0 0 0 1 1 0 0 0 0 0 1 Case12 1 1 1 4 1 2 2 1 1 1 4 5 0 0 1 0 1 1 0 NA 0 0 1 0 Case13 2 4 4 5 2 4 4 2 4 2 4 4 0 0 0 0 1 1 1 1 0 1 0 0 Case14 2 1 1 1 1 2 1 1 1 1 2 2 0 0 1 0 1 1 0 0 0 0 0 1 Case15 4 2 2 2 2 4 2 2 4 4 4 4 0 0 1 0 1 1 0 0 0 0 0 1 Case16 2 2 4 2 2 4 NA 4 4 4 4 NA 1 1 1 0 1 1 1 0 1 0 0 1 Case17 2 2 2 2 2 2 2 2 2 4 4 4 0 0 1 0 1 1 0 0 0 0 1 0 Case18 2 2 1 2 1 1 2 1 1 1 2 1 1 0 1 0 1 1 0 0 1 0 1 0 Case19 1 1 1 1 1 2 2 2 2 2 2 2 1 1 1 0 0 1 1 1 1 0 1 1 Case20 1 2 2 2 1 2 2 4 2 2 4 4 1 0 1 1 0 0 0 0 1 0 0 1 Case21 2 4 1 1 1 2 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 0 0 1 Case22 2 2 4 2 2 2 2 2 2 2 4 4 0 1 1 0 1 0 0 1 0 0 0 1 Case23 2 2 2 4 2 4 NA 2 2 4 4 4 1 1 1 0 1 0 0 0 0 1 1 1 Case24 2 2 2 2 2 4 2 2 2 2 4 2 1 1 1 1 1 1 1 0 0 0 1 1 Case25 2 2 1 2 2 4 2 1 1 2 4 4 0 0 1 0 0 0 0 0 0 0 0 1 Case26 2 2 2 4 4 4 4 2 2 2 4 4 1 0 1 0 1 1 0 0 0 1 0 1 Case27 2 2 2 4 4 4 4 2 2 2 4 4 1 0 1 0 1 1 0 0 0 0 NA 1 Case28 4 2 2 5 2 5 4 4 4 5 4 4 0 0 1 0 1 1 0 0 0 0 1 1 Case29 2 2 2 2 2 2 2 2 2 2 2 4 0 0 0 0 1 0 0 0 0 1 1 1 Case30 2 2 1 2 1 2 2 2 4 2 4 4 0 0 1 0 1 0 1 0 1 0 1 0 Case31 4 5 2 4 2 4 2 4 4 2 4 4 1 1 1 0 1 1 0 0 0 1 0 1 Case32 4 4 4 4 2 4 2 4 4 4 4 2 0 0 0 0 1 1 0 0 0 1 0 1 Case33 4 2 2 2 2 4 4 4 4 4 5 4 0 0 0 0 0 1 0 0 1 0 0 0 Case34 2 2 2 3 2 4 2 4 3 3 4 5 0 0 1 0 1 1 0 1 1 1 1 1 Case35 2 NA 2 NA 2 4 4 4 2 4 5 4 1 1 1 0 1 0 1 0 1 0 0 1 Case36 4 2 4 NA 4 2 5 4 4 4 5 4 0 1 1 0 1 1 NA 0 1 0 0 1 Case37 1 4 2 4 2 4 5 2 2 2 4 4 0 0 1 0 1 1 0 1 0 0 0 1 Case38 4 2 4 2 2 4 2 5 4 4 5 5 0 0 1 0 1 1 0 0 0 0 1 1 Case39 1 1 1 1 1 4 1 4 1 1 4 4 0 1 NA NA NA NA NA NA NA NA NA NA Case40 2 2 4 NA 5 NA 4 5 2 2 4 4 0 0 1 0 1 1 0 0 0 NA NA 1 Case41 2 4 2 2 NA 2 2 4 4 4 4 4 0 NA 0 0 1 1 0 0 0 NA 0 1 Case42 2 4 4 2 4 4 4 2 4 4 4 4 0 0 0 0 1 1 0 0 1 0 0 1 Case43 4 4 2 2 2 2 2 2 2 2 2 2 1 1 1 1 0 1 1 1 1 0 1 1 Case44 2 2 4 2 2 4 2 4 2 2 4 2 0 0 1 1 0 1 1 0 1 1 1 0 Case45 4 4 2 4 1 5 4 5 2 4 5 5 0 0 1 0 1 1 0 0 1 1 0 1 15 Case1 0 Case2 0 Case3 1 Case4 0 Case5 0 Case6 1 Case7 NA Case8 1 Case9 0 Case10 0 Case12 0 Case13 0 Case14 0 Case15 0 Case16 0 Case17 0 Case18 1 Case19 1 Case20 0 Case21 1 Case22 0 Case23 0 Case24 1 Case25 0 Case26 0 Case27 0 Case28 0 Case29 0 Case30 0 Case31 0 Case32 0 Case33 1 Case34 0 Case35 0 Case36 0 Case37 0 Case38 0 Case39 NA Case40 0 Case41 0 Case42 0 Case43 1 Case44 0 Case45 0 > fit <- principal(y, nfactors=par1, rotate='varimax') Loading required package: GPArotation Warning messages: 1: In log(det(m.inv.r)) : NaNs produced 2: In log(det(r)) : NaNs produced 3: In cor.smooth(r) : Matrix was not positive definite, smoothing was done > fit Principal Components Analysis Call: principal(r = y, nfactors = par1, rotate = "varimax") Standardized loadings (pattern matrix) based upon correlation matrix PC1 PC3 PC4 PC2 h2 u2 1a 0.09 0.75 -0.05 -0.09 0.587 0.41 1b 0.12 0.67 0.30 -0.14 0.568 0.43 1c 0.36 0.59 -0.02 -0.10 0.483 0.52 1d 0.75 -0.01 0.18 -0.39 0.745 0.26 1e 0.39 0.37 0.08 -0.39 0.445 0.55 1f 0.70 0.16 -0.03 0.07 0.524 0.48 1g 0.71 0.22 -0.01 -0.06 0.552 0.45 1h 0.60 0.27 -0.25 0.45 0.698 0.30 1i 0.25 0.78 -0.25 -0.02 0.728 0.27 1j 0.39 0.58 -0.35 0.04 0.610 0.39 1k 0.75 0.12 -0.27 0.07 0.664 0.34 1l 0.60 -0.13 -0.41 -0.23 0.600 0.40 2 -0.09 -0.02 0.73 0.05 0.537 0.46 3 -0.13 -0.05 0.53 0.21 0.344 0.66 4a 0.09 -0.43 0.36 -0.12 0.341 0.66 6 -0.14 -0.19 0.39 0.54 0.502 0.50 7 0.16 0.20 0.13 -0.65 0.511 0.49 8 0.13 0.33 0.27 -0.03 0.197 0.80 9 0.09 -0.12 0.43 0.60 0.566 0.43 10 0.02 -0.20 0.57 -0.02 0.367 0.63 11 0.07 -0.01 0.11 0.78 0.630 0.37 12 0.32 0.08 0.20 -0.04 0.153 0.85 13 0.05 -0.53 0.09 0.19 0.326 0.67 14 0.15 0.13 0.21 -0.03 0.084 0.92 15 -0.36 -0.03 0.59 0.39 0.625 0.37 PC1 PC3 PC4 PC2 SS loadings 3.72 3.34 2.76 2.57 Proportion Var 0.15 0.13 0.11 0.10 Cumulative Var 0.15 0.28 0.39 0.50 Proportion Explained 0.30 0.27 0.22 0.21 Cumulative Proportion 0.30 0.57 0.79 1.00 Test of the hypothesis that 4 components are sufficient. The degrees of freedom for the null model are 300 and the objective function was NaN The degrees of freedom for the model are 206 and the objective function was NaN The number of observations was 44 with Chi Square = NaN with prob < NaN Fit based upon off diagonal values = 0.85> fs <- factor.scores(y,fit) Warning message: In cor.smooth(r) : Matrix was not positive definite, smoothing was done > fs $scores PC1 PC3 PC4 PC2 Case1 1.83217230 -1.55582394 -0.66673196 0.81120201 Case2 0.48885383 -0.82456447 0.18070391 0.76780374 Case3 0.43094663 -0.58404892 0.60493577 1.10419101 Case4 1.62480657 -2.06415114 -0.37942310 1.21475384 Case5 0.13948929 -1.01149398 -0.56100977 -0.14625146 Case6 -0.12399882 -1.00940680 0.26036421 0.61286569 Case7 NA NA NA NA Case8 NA NA NA NA Case9 0.84766150 -1.35652290 -0.36423646 0.56338205 Case10 NA NA NA NA Case12 NA NA NA NA Case13 0.11435848 0.77113709 -0.01242767 -0.25655423 Case14 -1.54666527 0.29295112 -0.17831748 -1.09646791 Case15 -0.47705897 0.94672673 -0.38552713 -1.03146541 Case16 NA NA NA NA Case17 -0.06122791 -0.35528973 -0.92929057 -0.90777993 Case18 -1.86668081 0.58069966 0.94693290 -0.14266346 Case19 -1.13723903 -0.05054861 1.57498519 1.16062392 Case20 0.33436153 -0.98450086 -0.03538725 0.51588996 Case21 -2.42018929 1.53175734 2.52898539 0.93281951 Case22 -0.49650800 0.16629570 -0.30440433 -0.65114064 Case23 NA NA NA NA Case24 -0.79973692 0.41999754 1.34259907 0.89461943 Case25 0.20538183 -0.98019889 -1.02717923 -0.63599159 Case26 0.29277587 0.01896929 -0.56062156 -0.47353617 Case27 NA NA NA NA Case28 0.94883834 -0.03295913 -0.85446524 -0.43749211 Case29 -0.53922687 -0.07812115 -0.70015983 -1.00880622 Case30 0.27370626 -0.94372507 -0.59950903 -0.14698768 Case31 -0.69604428 1.66777580 0.64062229 -0.30484292 Case32 -0.81774535 2.18853369 0.41671512 -0.81137623 Case33 0.22601624 0.59115416 0.17200973 0.07671067 Case34 1.00560012 -0.86780003 -0.61989328 0.17176782 Case35 NA NA NA NA Case36 NA NA NA NA Case37 0.36950768 -0.19325994 -0.60111043 -0.39586617 Case38 0.64275352 0.19151155 -0.73802751 -0.55202629 Case39 NA NA NA NA Case40 NA NA NA NA Case41 NA NA NA NA Case42 -0.32940514 1.52087968 0.04562566 -0.70818674 Case43 -1.93615376 1.66912986 2.69696844 1.31029136 Case44 0.55446484 -0.63412796 0.48664050 1.06222502 Case45 1.24258306 -0.13548228 -0.50913826 0.12962391 $weights PC1 PC3 PC4 PC2 1a -0.14888441 0.274348823 0.1146772158 -0.05710295 1b -0.16603403 0.290606608 0.1683039582 -0.01380861 1c -0.03553950 0.156435328 0.0613623262 -0.01188237 1d 0.16001900 -0.103820455 -0.0839410645 0.02972408 1e -0.01791859 0.102494916 -0.0002395364 -0.04947338 1f 0.17794813 -0.098099031 -0.0388141010 0.08630641 1g 0.15556824 -0.071564370 -0.0471107944 0.05547540 1h 0.18391575 -0.091711820 -0.0028643250 0.12648532 1i -0.07641790 0.214228217 0.0683475563 -0.04731828 1j 0.02709298 0.088593472 0.0032789303 -0.01682767 1k 0.23523631 -0.168777547 -0.1122253769 0.06288864 1l 0.23325997 -0.226836872 -0.2096926333 -0.02211253 2 -0.10046070 0.111452823 0.1671231453 0.08771966 3 -0.06838351 0.071015510 0.1478081836 0.09886668 4a 0.06983387 -0.124533308 -0.0295360053 0.05135966 6 0.00665065 -0.024183070 0.1354969498 0.16085059 7 -0.07557802 0.111385951 -0.0343049432 -0.11965016 8 -0.07170182 0.140902556 0.1046796659 0.02293429 9 0.05952160 -0.047150332 0.1400643660 0.20167427 10 -0.01691042 0.001041801 0.0795347939 0.07751100 11 0.08427303 -0.061013462 0.1176154089 0.19543093 12 0.04915378 -0.007783671 0.0193939371 0.04907442 13 0.14335299 -0.222415443 -0.0695154028 0.08347607 14 -0.01106686 0.048392210 0.0500374752 0.03007975 15 -0.13121724 0.128566781 0.2174016337 0.12061628 $r.scores PC1 PC3 PC4 PC2 PC1 1.0000000 -0.6954926 -0.6559380 0.1640154 PC3 -0.6954926 1.0000000 0.5227935 -0.3059552 PC4 -0.6559380 0.5227935 1.0000000 0.5790302 PC2 0.1640154 -0.3059552 0.5790302 1.0000000 $R2 PC1 PC3 PC4 PC2 0.8008181 0.9472730 0.6487814 0.6035366 > postscript(file="/var/wessaorg/rcomp/tmp/1kjcj1352478640.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > fa.diagram(fit) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/2qg1w1352478640.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(fs$scores,pch=20) > text(fs$scores,labels=rownames(y),pos=3) > dev.off() null device 1 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Rotated Factor Loadings',par1+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variables',1,TRUE) > for (i in 1:par1) { + a<-table.element(a,paste('Factor',i,sep=''),1,TRUE) + } > a<-table.row.end(a) > for (j in 1:length(fit$loadings[,1])) { + a<-table.row.start(a) + a<-table.element(a,rownames(fit$loadings)[j],header=TRUE) + for (i in 1:par1) { + a<-table.element(a,round(fit$loadings[j,i],3)) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/3zpr11352478640.tab") > > try(system("convert tmp/1kjcj1352478640.ps tmp/1kjcj1352478640.png",intern=TRUE)) character(0) > try(system("convert tmp/2qg1w1352478640.ps tmp/2qg1w1352478640.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.847 0.490 3.414