R version 3.0.3 (2014-03-06) -- "Warm Puppy" Copyright (C) 2014 The R Foundation for Statistical Computing 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('GTM' + ,110.39 + ,6.80 + ,52.90 + ,3.48 + ,21.38 + ,67.64 + ,26.64 + ,32.86 + ,27.47 + ,0.70 + ,12.70 + ,1.12 + ,2.11 + ,4.08 + ,49.87 + ,7.88 + ,2.84 + ,0.38 + ,5.01 + ,1.66 + ,67.64 + ,1.31 + ,0.99 + ,0.57 + ,77.34 + ,'PRO' + ,23.49 + ,6.35 + ,10.16 + ,1.90 + ,15.24 + ,58.42 + ,2.54 + ,2.54 + ,6.98 + ,1.27 + ,3.17 + ,1.90 + ,0.63 + ,0.63 + ,38.10 + ,3.17 + ,5.08 + ,1.90 + ,8.89 + ,1.27 + ,58.42 + ,0.63 + ,0.63 + ,0.63 + ,44.45 + ,'SAC' + ,66.63 + ,7.58 + ,39.79 + ,2.21 + ,15.79 + ,21.16 + ,13.89 + ,18.63 + ,22.74 + ,0.00 + ,9.16 + ,0.63 + ,0.95 + ,0.63 + ,10.10 + ,1.58 + ,1.26 + ,0.00 + ,6.32 + ,1.89 + ,21.16 + ,4.11 + ,0.32 + ,0.63 + ,36.95 + ,'CHM' + ,34.49 + ,3.47 + ,25.74 + ,2.64 + ,7.43 + ,14.36 + ,5.78 + ,4.79 + ,4.29 + ,0.33 + ,5.78 + ,0.00 + ,0.33 + ,0.17 + ,9.41 + ,0.83 + ,1.16 + ,0.33 + ,2.15 + ,0.50 + ,14.36 + ,0.83 + ,0.50 + ,0.17 + ,35.81 + ,'ESC' + ,74.73 + ,5.41 + ,200.57 + ,7.40 + ,26.33 + ,98.22 + ,18.08 + ,20.36 + ,2.14 + ,0.14 + ,4.13 + ,0.71 + ,1.42 + ,1.57 + ,70.75 + ,8.11 + ,5.98 + ,0.28 + ,11.25 + ,1.85 + ,98.22 + ,3.27 + ,1.00 + ,1.57 + ,73.88 + ,'SRO' + ,33.63 + ,3.48 + ,24.93 + ,2.32 + ,20.29 + ,82.34 + ,11.31 + ,12.18 + ,4.35 + ,0.87 + ,5.22 + ,0.87 + ,1.16 + ,0.00 + ,56.83 + ,4.64 + ,6.67 + ,1.74 + ,8.99 + ,3.48 + ,82.34 + ,0.58 + ,0.87 + ,0.00 + ,67.84 + ,'SOL' + ,2.55 + ,0.70 + ,1.16 + ,0.23 + ,0.93 + ,8.59 + ,3.25 + ,3.25 + ,1.86 + ,0.70 + ,2.55 + ,0.70 + ,0.23 + ,0.00 + ,3.72 + ,0.00 + ,0.70 + ,0.23 + ,3.48 + ,0.46 + ,8.59 + ,1.39 + ,0.23 + ,0.00 + ,18.35 + ,'TOT' + ,14.60 + ,1.06 + ,1.27 + ,0.00 + ,0.63 + ,7.83 + ,2.75 + ,2.96 + ,0.21 + ,0.00 + ,2.54 + ,0.42 + ,0.00 + ,0.00 + ,0.85 + ,0.42 + ,0.63 + ,0.21 + ,4.23 + ,1.48 + ,7.83 + ,0.85 + ,0.63 + ,0.42 + ,15.87 + ,'QUT' + ,42.78 + ,1.89 + ,30.92 + ,3.03 + ,9.47 + ,26.88 + ,18.43 + ,21.33 + ,3.66 + ,1.64 + ,4.04 + ,1.64 + ,0.63 + ,1.14 + ,17.54 + ,1.77 + ,1.77 + ,0.38 + ,4.29 + ,1.14 + ,26.88 + ,0.25 + ,0.76 + ,0.63 + ,34.83 + ,'SUC' + ,32.74 + ,1.45 + ,43.51 + ,1.86 + ,9.95 + ,30.87 + ,18.23 + ,20.93 + ,4.35 + ,1.24 + ,5.39 + ,1.04 + ,1.04 + ,0.41 + ,17.20 + ,2.28 + ,2.69 + ,0.21 + ,7.25 + ,1.24 + ,30.87 + ,0.41 + ,0.21 + ,0.41 + ,36.68 + ,'RET' + ,15.84 + ,1.32 + ,52.14 + ,1.98 + ,5.94 + ,35.97 + ,12.87 + ,16.50 + ,1.32 + ,1.65 + ,2.64 + ,1.98 + ,1.32 + ,0.66 + ,23.10 + ,1.32 + ,1.65 + ,0.33 + ,8.25 + ,1.32 + ,35.97 + ,1.32 + ,0.66 + ,0.66 + ,37.62 + ,'SMA' + ,8.02 + ,0.78 + ,7.14 + ,0.78 + ,4.70 + ,15.56 + ,4.21 + ,4.11 + ,1.66 + ,0.49 + ,1.66 + ,0.10 + ,0.10 + ,0.00 + ,10.67 + ,0.29 + ,1.37 + ,0.10 + ,2.74 + ,0.39 + ,15.56 + ,1.27 + ,0.10 + ,0.10 + ,12.72 + ,'HUE' + ,8.69 + ,0.70 + ,31.64 + ,1.30 + ,4.35 + ,6.43 + ,4.78 + ,4.87 + ,1.65 + ,0.52 + ,1.13 + ,0.09 + ,0.17 + ,0.00 + ,3.82 + ,0.26 + ,0.70 + ,0.09 + ,1.48 + ,0.09 + ,6.43 + ,1.83 + ,0.26 + ,0.70 + ,13.12 + ,'QUI' + ,4.50 + ,0.84 + ,12.66 + ,0.94 + ,2.09 + ,6.38 + ,2.93 + ,3.14 + ,0.73 + ,0.10 + ,2.41 + ,0.42 + ,0.21 + ,0.00 + ,2.30 + ,0.31 + ,0.73 + ,0.21 + ,2.41 + ,0.42 + ,6.38 + ,0.52 + ,0.42 + ,0.21 + ,14.44 + ,'BVP' + ,6.33 + ,1.86 + ,13.40 + ,2.23 + ,6.33 + ,18.62 + ,12.29 + ,14.89 + ,2.61 + ,0.37 + ,4.84 + ,2.61 + ,0.00 + ,0.00 + ,7.45 + ,1.49 + ,0.74 + ,0.37 + ,7.82 + ,0.74 + ,18.62 + ,0.37 + ,1.12 + ,0.74 + ,23.09 + ,'AVP' + ,8.48 + ,0.27 + ,14.56 + ,0.63 + ,6.97 + ,13.04 + ,8.13 + ,9.73 + ,1.43 + ,0.63 + ,1.07 + ,0.27 + ,0.18 + ,0.00 + ,7.41 + ,0.45 + ,1.61 + ,0.18 + ,2.86 + ,0.54 + ,13.04 + ,0.71 + ,0.18 + ,0.27 + ,14.47 + ,'PET' + ,3.82 + ,1.11 + ,38.03 + ,1.75 + ,11.30 + ,59.36 + ,5.89 + ,3.66 + ,1.27 + ,0.48 + ,1.91 + ,0.00 + ,0.48 + ,0.00 + ,41.54 + ,3.66 + ,7.48 + ,0.80 + ,4.93 + ,0.95 + ,59.36 + ,2.23 + ,0.64 + ,0.16 + ,25.30 + ,'IZA' + ,12.66 + ,1.95 + ,80.82 + ,3.65 + ,20.69 + ,79.85 + ,5.60 + ,5.84 + ,2.43 + ,0.73 + ,2.19 + ,0.49 + ,0.00 + ,0.00 + ,59.16 + ,6.82 + ,5.60 + ,1.46 + ,5.84 + ,0.97 + ,79.85 + ,0.49 + ,0.73 + ,0.73 + ,49.18 + ,'ZAC' + ,17.17 + ,0.90 + ,29.82 + ,2.26 + ,42.46 + ,92.61 + ,6.32 + ,4.97 + ,6.32 + ,0.00 + ,4.52 + ,0.00 + ,0.00 + ,0.00 + ,66.41 + ,7.68 + ,7.23 + ,0.45 + ,8.13 + ,2.71 + ,92.61 + ,1.81 + ,0.00 + ,1.81 + ,67.76 + ,'CHQ' + ,20.38 + ,2.45 + ,40.76 + ,0.82 + ,29.89 + ,96.20 + ,8.15 + ,7.61 + ,3.26 + ,0.27 + ,2.45 + ,0.27 + ,0.00 + ,0.00 + ,58.70 + ,6.25 + ,14.95 + ,1.90 + ,11.68 + ,2.72 + ,96.20 + ,2.45 + ,0.54 + ,0.82 + ,58.15 + ,'JAL' + ,1.27 + ,1.59 + ,2.87 + ,1.27 + ,12.42 + ,54.44 + ,4.78 + ,6.05 + ,1.59 + ,0.32 + ,1.59 + ,1.91 + ,0.32 + ,0.00 + ,33.11 + ,4.78 + ,7.96 + ,1.91 + ,5.09 + ,1.59 + ,54.44 + ,1.27 + ,0.00 + ,0.96 + ,35.02 + ,'JUT' + ,12.90 + ,1.15 + ,6.22 + ,1.15 + ,6.22 + ,61.26 + ,7.83 + ,8.52 + ,3.45 + ,0.23 + ,4.15 + ,0.69 + ,0.23 + ,0.00 + ,42.83 + ,7.37 + ,4.38 + ,1.38 + ,4.84 + ,0.46 + ,61.26 + ,4.38 + ,0.69 + ,0.23 + ,46.29) + ,dim=c(26 + ,22) + ,dimnames=list(c('DEP' + ,'Robo_vehiculos' + ,'Capt_Vehículos' + ,'Motocicletas' + ,'Capt_Motocicletas' + ,'Armas_robadas' + ,'Homicidios' + ,'Peatones' + ,'Capt_Peatones' + ,'Residencias' + ,'Capt_Residencias' + ,'Comercios' + ,'Capt_Comercios' + ,'Buses' + ,'Capt_Buses' + ,'HArmaFuego_M' + ,'HArmaFuego_F' + ,'HArmaBlanca_M' + ,'HArmaBlanca_F' + ,'HAsfixiaSec' + ,'MasfixiaSec' + ,'TOTHOM_' + ,'MA_DETERMINAR_M' + ,'MA_DETERMINAR_F' + ,'MA_DETERMINAR_IND' + ,'M_TOTAL_conIND') + ,1:22)) > y <- array(NA,dim=c(26,22),dimnames=list(c('DEP','Robo_vehiculos','Capt_Vehículos','Motocicletas','Capt_Motocicletas','Armas_robadas','Homicidios','Peatones','Capt_Peatones','Residencias','Capt_Residencias','Comercios','Capt_Comercios','Buses','Capt_Buses','HArmaFuego_M','HArmaFuego_F','HArmaBlanca_M','HArmaBlanca_F','HAsfixiaSec','MasfixiaSec','TOTHOM_','MA_DETERMINAR_M','MA_DETERMINAR_F','MA_DETERMINAR_IND','M_TOTAL_conIND'),1:22)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } There were 22 warnings (use warnings() to see them) > par1 = '3' > par1 <- '3' > #'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 Robo_vehiculos Capt_Veh\355culos Motocicletas Capt_Motocicletas GTM 110.39 6.80 52.90 3.48 PRO 23.49 6.35 10.16 1.90 SAC 66.63 7.58 39.79 2.21 CHM 34.49 3.47 25.74 2.64 ESC 74.73 5.41 200.57 7.40 SRO 33.63 3.48 24.93 2.32 SOL 2.55 0.70 1.16 0.23 TOT 14.60 1.06 1.27 0.00 QUT 42.78 1.89 30.92 3.03 SUC 32.74 1.45 43.51 1.86 RET 15.84 1.32 52.14 1.98 SMA 8.02 0.78 7.14 0.78 HUE 8.69 0.70 31.64 1.30 QUI 4.50 0.84 12.66 0.94 BVP 6.33 1.86 13.40 2.23 AVP 8.48 0.27 14.56 0.63 PET 3.82 1.11 38.03 1.75 IZA 12.66 1.95 80.82 3.65 ZAC 17.17 0.90 29.82 2.26 CHQ 20.38 2.45 40.76 0.82 JAL 1.27 1.59 2.87 1.27 JUT 12.90 1.15 6.22 1.15 Armas_robadas Homicidios Peatones Capt_Peatones Residencias GTM 21.38 67.64 26.64 32.86 27.47 PRO 15.24 58.42 2.54 2.54 6.98 SAC 15.79 21.16 13.89 18.63 22.74 CHM 7.43 14.36 5.78 4.79 4.29 ESC 26.33 98.22 18.08 20.36 2.14 SRO 20.29 82.34 11.31 12.18 4.35 SOL 0.93 8.59 3.25 3.25 1.86 TOT 0.63 7.83 2.75 2.96 0.21 QUT 9.47 26.88 18.43 21.33 3.66 SUC 9.95 30.87 18.23 20.93 4.35 RET 5.94 35.97 12.87 16.50 1.32 SMA 4.70 15.56 4.21 4.11 1.66 HUE 4.35 6.43 4.78 4.87 1.65 QUI 2.09 6.38 2.93 3.14 0.73 BVP 6.33 18.62 12.29 14.89 2.61 AVP 6.97 13.04 8.13 9.73 1.43 PET 11.30 59.36 5.89 3.66 1.27 IZA 20.69 79.85 5.60 5.84 2.43 ZAC 42.46 92.61 6.32 4.97 6.32 CHQ 29.89 96.20 8.15 7.61 3.26 JAL 12.42 54.44 4.78 6.05 1.59 JUT 6.22 61.26 7.83 8.52 3.45 Capt_Residencias Comercios Capt_Comercios Buses Capt_Buses HArmaFuego_M GTM 0.70 12.70 1.12 2.11 4.08 49.87 PRO 1.27 3.17 1.90 0.63 0.63 38.10 SAC 0.00 9.16 0.63 0.95 0.63 10.10 CHM 0.33 5.78 0.00 0.33 0.17 9.41 ESC 0.14 4.13 0.71 1.42 1.57 70.75 SRO 0.87 5.22 0.87 1.16 0.00 56.83 SOL 0.70 2.55 0.70 0.23 0.00 3.72 TOT 0.00 2.54 0.42 0.00 0.00 0.85 QUT 1.64 4.04 1.64 0.63 1.14 17.54 SUC 1.24 5.39 1.04 1.04 0.41 17.20 RET 1.65 2.64 1.98 1.32 0.66 23.10 SMA 0.49 1.66 0.10 0.10 0.00 10.67 HUE 0.52 1.13 0.09 0.17 0.00 3.82 QUI 0.10 2.41 0.42 0.21 0.00 2.30 BVP 0.37 4.84 2.61 0.00 0.00 7.45 AVP 0.63 1.07 0.27 0.18 0.00 7.41 PET 0.48 1.91 0.00 0.48 0.00 41.54 IZA 0.73 2.19 0.49 0.00 0.00 59.16 ZAC 0.00 4.52 0.00 0.00 0.00 66.41 CHQ 0.27 2.45 0.27 0.00 0.00 58.70 JAL 0.32 1.59 1.91 0.32 0.00 33.11 JUT 0.23 4.15 0.69 0.23 0.00 42.83 HArmaFuego_F HArmaBlanca_M HArmaBlanca_F HAsfixiaSec MasfixiaSec TOTHOM_ GTM 7.88 2.84 0.38 5.01 1.66 67.64 PRO 3.17 5.08 1.90 8.89 1.27 58.42 SAC 1.58 1.26 0.00 6.32 1.89 21.16 CHM 0.83 1.16 0.33 2.15 0.50 14.36 ESC 8.11 5.98 0.28 11.25 1.85 98.22 SRO 4.64 6.67 1.74 8.99 3.48 82.34 SOL 0.00 0.70 0.23 3.48 0.46 8.59 TOT 0.42 0.63 0.21 4.23 1.48 7.83 QUT 1.77 1.77 0.38 4.29 1.14 26.88 SUC 2.28 2.69 0.21 7.25 1.24 30.87 RET 1.32 1.65 0.33 8.25 1.32 35.97 SMA 0.29 1.37 0.10 2.74 0.39 15.56 HUE 0.26 0.70 0.09 1.48 0.09 6.43 QUI 0.31 0.73 0.21 2.41 0.42 6.38 BVP 1.49 0.74 0.37 7.82 0.74 18.62 AVP 0.45 1.61 0.18 2.86 0.54 13.04 PET 3.66 7.48 0.80 4.93 0.95 59.36 IZA 6.82 5.60 1.46 5.84 0.97 79.85 ZAC 7.68 7.23 0.45 8.13 2.71 92.61 CHQ 6.25 14.95 1.90 11.68 2.72 96.20 JAL 4.78 7.96 1.91 5.09 1.59 54.44 JUT 7.37 4.38 1.38 4.84 0.46 61.26 MA_DETERMINAR_M MA_DETERMINAR_F MA_DETERMINAR_IND M_TOTAL_conIND GTM 1.31 0.99 0.57 77.34 PRO 0.63 0.63 0.63 44.45 SAC 4.11 0.32 0.63 36.95 CHM 0.83 0.50 0.17 35.81 ESC 3.27 1.00 1.57 73.88 SRO 0.58 0.87 0.00 67.84 SOL 1.39 0.23 0.00 18.35 TOT 0.85 0.63 0.42 15.87 QUT 0.25 0.76 0.63 34.83 SUC 0.41 0.21 0.41 36.68 RET 1.32 0.66 0.66 37.62 SMA 1.27 0.10 0.10 12.72 HUE 1.83 0.26 0.70 13.12 QUI 0.52 0.42 0.21 14.44 BVP 0.37 1.12 0.74 23.09 AVP 0.71 0.18 0.27 14.47 PET 2.23 0.64 0.16 25.30 IZA 0.49 0.73 0.73 49.18 ZAC 1.81 0.00 1.81 67.76 CHQ 2.45 0.54 0.82 58.15 JAL 1.27 0.00 0.96 35.02 JUT 4.38 0.69 0.23 46.29 > fit <- principal(y, nfactors=par1, rotate='varimax') Loading required package: GPArotation The determinant of the smoothed correlation was zero. This means the objective function is not defined for the null model either. The Chi square is thus based upon observed correlations. In factor.scores, the correlation matrix is singular, an approximation is used Warning messages: 1: In log(det(m.inv.r)) : NaNs produced 2: In cor.smooth(r) : Matrix was not positive definite, smoothing was done 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 RC1 RC2 RC3 h2 u2 Robo_vehiculos 0.18 0.96 -0.08 0.95 0.047 Capt_Veh\355culos 0.25 0.74 -0.05 0.62 0.379 Motocicletas 0.48 0.47 -0.07 0.45 0.547 Capt_Motocicletas 0.46 0.60 0.06 0.58 0.415 Armas_robadas 0.86 0.23 -0.24 0.84 0.158 Homicidios 0.98 0.14 -0.03 0.98 0.025 Peatones 0.13 0.87 0.25 0.84 0.157 Capt_Peatones 0.07 0.89 0.27 0.86 0.136 Residencias 0.02 0.81 -0.23 0.71 0.285 Capt_Residencias -0.09 0.12 0.83 0.71 0.294 Comercios 0.04 0.88 -0.10 0.79 0.206 Capt_Comercios 0.02 0.19 0.78 0.64 0.357 Buses 0.14 0.84 0.26 0.80 0.199 Capt_Buses 0.09 0.89 0.11 0.82 0.184 HArmaFuego_M 0.95 0.17 -0.05 0.94 0.061 HArmaFuego_F 0.86 0.28 -0.15 0.85 0.152 HArmaBlanca_M 0.88 -0.22 -0.07 0.83 0.170 HArmaBlanca_F 0.68 -0.34 0.27 0.65 0.349 HAsfixiaSec 0.80 0.17 0.24 0.73 0.274 MasfixiaSec 0.73 0.21 0.01 0.58 0.421 TOTHOM_ 0.98 0.14 -0.03 0.98 0.025 MA_DETERMINAR_M 0.26 0.20 -0.62 0.49 0.512 MA_DETERMINAR_F 0.24 0.45 0.43 0.44 0.559 MA_DETERMINAR_IND 0.59 0.17 -0.18 0.41 0.593 M_TOTAL_conIND 0.81 0.52 0.00 0.93 0.073 RC1 RC2 RC3 SS loadings 8.45 7.58 2.40 Proportion Var 0.34 0.30 0.10 Cumulative Var 0.34 0.64 0.74 Proportion Explained 0.46 0.41 0.13 Cumulative Proportion 0.46 0.87 1.00 Test of the hypothesis that 3 components are sufficient. The degrees of freedom for the null model are 300 and the objective function was 155.94 The degrees of freedom for the model are 228 and the objective function was NaN The total number of observations was 22 with MLE Chi Square = NaN with prob < NaN Fit based upon off diagonal values = 0.97> fs <- factor.scores(y,fit) Warning messages: 1: In cor.smooth(r) : Matrix was not positive definite, smoothing was done 2: In invMatSqrt(t(L) %*% inv.r %*% L) : complex eigen values detected by invMatSqrt, results are suspect > fs $scores RC1 RC2 RC3 GTM 0.03545892 3.201586979 -0.004593578 PRO 0.62581570 -0.317190439 1.354360486 SAC -0.72278865 1.665972211 -1.812609552 CHM -0.90315369 0.146661518 -0.700559039 ESC 1.64514951 1.485854016 -0.564263656 SRO 1.21166760 0.007016348 1.087729828 SOL -1.10051331 -0.632259225 -0.071787878 TOT -0.91882123 -0.612739048 -0.424374531 QUT -0.51353188 0.617515375 1.655359165 SUC -0.42940514 0.488323077 0.897313632 RET -0.19770060 0.201042088 1.819534397 SMA -1.00483438 -0.641535670 -0.613935701 HUE -1.07888429 -0.489920834 -0.762677069 QUI -1.10420400 -0.595334611 -0.407219741 BVP -0.52713241 -0.021936882 1.343741642 AVP -0.95298753 -0.545011547 -0.161153966 PET 0.28886895 -0.674438663 -0.396917098 IZA 1.06880348 -0.595008350 0.262215290 ZAC 1.62368778 -0.470706214 -1.691591790 CHQ 1.98078516 -0.859739328 -0.407410746 JAL 0.64075586 -0.995934527 0.322628333 JUT 0.33296413 -0.362216272 -0.723788430 $weights RC1 RC2 RC3 Robo_vehiculos -0.011883649 0.162433515 -0.0605377093 Capt_Veh\355culos -0.010280172 0.093576501 -0.0466317021 Motocicletas 0.040318591 0.052336416 -0.0347835803 Capt_Motocicletas 0.030851589 0.061596979 0.0139676239 Armas_robadas 0.102576950 0.007257080 -0.0893250371 Homicidios 0.123625439 -0.027176570 0.0105229188 Peatones -0.022276643 0.112311301 0.0749967411 Capt_Peatones -0.031236781 0.116857135 0.0792597262 Residencias -0.043461438 0.135186357 -0.1276693086 Capt_Residencias 0.004481060 -0.011382636 0.3480631478 Comercios -0.043084419 0.134494190 -0.0787287777 Capt_Comercios 0.019513520 0.004985244 0.3295716487 Buses -0.013676915 0.115841056 0.0837931877 Capt_Buses -0.037667203 0.115599106 0.0102799206 HArmaFuego_M 0.123425797 -0.014749898 -0.0020014736 HArmaFuego_F 0.098242697 0.004802566 -0.0518688805 HArmaBlanca_M 0.130372271 -0.071309761 0.0008810816 HArmaBlanca_F 0.116890216 -0.096994639 0.1453094869 HAsfixiaSec 0.105791953 -0.023521508 0.1169035164 MasfixiaSec 0.085965012 -0.010669865 0.0146687811 TOTHOM_ 0.127890103 -0.025397932 0.0072209447 MA_DETERMINAR_M 0.007941942 0.041391951 -0.2661598235 MA_DETERMINAR_F 0.022981186 0.041607532 0.1722054400 MA_DETERMINAR_IND 0.066400287 0.003077523 -0.0674669785 M_TOTAL_conIND 0.082065645 0.035687047 -0.0010904374 $r.scores RC1 RC2 RC3 RC1 1.000000000 1.341838e-03 -1.601366e-03 RC2 0.001341838 1.000000e+00 6.506875e-05 RC3 -0.001601366 6.506875e-05 1.000000e+00 $R2 RC1 RC2 RC3 1.0079214 0.9996823 0.9997887 > postscript(file="/var/wessaorg/rcomp/tmp/1ohdq1396455957.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=16.666666666667,height=16.666666666667) > fa.diagram(fit) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/2kgfj1396455957.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=16.666666666667,height=16.666666666667) > 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/3rwkj1396455957.tab") > > try(system("convert tmp/1ohdq1396455957.ps tmp/1ohdq1396455957.png",intern=TRUE)) character(0) > try(system("convert tmp/2kgfj1396455957.ps tmp/2kgfj1396455957.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.932 0.809 3.745