R version 2.11.0 (2010-04-22) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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(0.65321251377534 + ,0 + ,0.81954393554187 + ,1.6232492903979 + ,0.47712125471966 + ,0 + ,0.47712125471966 + ,0.79934054945358 + ,0.30102999566398 + ,1.83884909073726 + ,3.40602894496362 + ,3.66304097489397 + ,2.79518458968242 + ,0.47712125471966 + ,0.69897000433602 + ,0.60205999132796 + ,0.32221929473392 + ,0.25527250510331 + ,1.43136376415899 + ,1.02325245963371 + ,2.25406445291434 + ,2.25527250510331 + ,0.60205999132796 + ,0.60205999132796 + ,0.60205999132796 + ,0.95904139232109 + ,-0.15490195998574 + ,1.27875360095283 + ,-1.69897000433602 + ,-0.52287874528034 + ,1.54406804435028 + ,0 + ,0 + ,0 + ,1.19865708695442 + ,0.5910646070265 + ,1.48287358360875 + ,2.20411998265592 + ,2.22788670461367 + ,2.59328606702046 + ,0.60205999132796 + ,0.69897000433602 + ,0.60205999132796 + ,0.7160033436348 + ,0 + ,1.44715803134222 + ,0.51851393987789 + ,1.40823996531185 + ,1.79934054945358 + ,0 + ,0.30102999566398 + ,0 + ,1.03742649794062 + ,0.55630250076729 + ,1.69897000433602 + ,1.71733758272386 + ,2.64345267648619 + ,2.36172783601759 + ,0 + ,0 + ,0 + ,0.91907809237607 + ,0.14612803567824 + ,0.84509804001426 + ,-0.36653154442041 + ,0.80617997398389 + ,2.04921802267018 + ,0.69897000433602 + ,0.60205999132796 + ,0.60205999132796 + ,1.04139268515823 + ,0.17609125905568 + ,1.47712125471966 + ,2.66745295288995 + ,2.62634036737504 + ,2.44870631990508 + ,0.69897000433602 + ,0.69897000433602 + ,0.69897000433602 + ,0.50514997831991 + ,-0.15490195998574 + ,0.54406804435028 + ,-1.09691001300806 + ,0.079181246047625 + ,1.6232492903979 + ,0 + ,0 + ,0 + ,0.79934054945358 + ,0.32221929473392 + ,0.77815125038364 + ,-0.10237290870956 + ,0.54406804435028 + ,1.6232492903979 + ,0.30102999566398 + ,0.30102999566398 + ,0.30102999566398 + ,0.81954393554187 + ,0.61278385671974 + ,1.01703333929878 + ,-0.69897000433602 + ,0.69897000433602 + ,2.07918124604762 + ,0.30102999566398 + ,0.30102999566398 + ,0.30102999566398 + ,0.97772360528885 + ,0.079181246047625 + ,1.30102999566398 + ,1.44185217577329 + ,2.06069784035361 + ,2.17026171539496 + ,0.69897000433602 + ,0.69897000433602 + ,0.69897000433602 + ,0.51851393987789 + ,-0.30102999566398 + ,0.5910646070265 + ,-0.92081875395238 + ,0 + ,1.20411998265592 + ,0.47712125471966 + ,0 + ,0.30102999566398 + ,1.04139268515823 + ,0.53147891704226 + ,1.61278385671974 + ,1.92941892571429 + ,2.51188336097887 + ,2.49136169383427 + ,0 + ,0.47712125471966 + ,0 + ,0.67209785793572 + ,0.17609125905568 + ,0.95424250943932 + ,-1 + ,0.60205999132796 + ,1.44715803134222 + ,0.69897000433602 + ,0 + ,0.47712125471966 + ,1.01703333929878 + ,0.53147891704226 + ,0.88081359228079 + ,0.01703333929878 + ,0.74036268949424 + ,1.83250891270624 + ,0.69897000433602 + ,0.47712125471966 + ,0.60205999132796 + ,0.86923171973098 + ,-0.096910013008056 + ,1.66275783168157 + ,2.71683772329952 + ,2.81624129999178 + ,2.52633927738984 + ,0.69897000433602 + ,0.69897000433602 + ,0.69897000433602 + ,0.32221929473392 + ,-0.096910013008056 + ,1.38021124171161 + ,-2 + ,-0.60205999132796 + ,1.69897000433602 + ,0 + ,0 + ,0 + ,1.25285303097989 + ,0.30102999566398 + ,2 + ,1.79239168949825 + ,3.12057393120585 + ,2.42651126136458 + ,0 + ,0 + ,0 + ,0.78532983501077 + ,0.27875360095283 + ,0.50514997831991 + ,-1.69897000433602 + ,-0.39794000867204 + ,1.27875360095283 + ,0.60205999132796 + ,0 + ,0.47712125471966 + ,1.07554696139253 + ,0.11394335230684 + ,0.69897000433602 + ,0.23044892137827 + ,0.79934054945358 + ,1.07918124604762 + ,0.30102999566398 + ,0 + ,0 + ,1.13987908640124 + ,0.7481880270062 + ,0.81291335664286 + ,0.54406804435028 + ,1.03342375548695 + ,2.07918124604762 + ,0.30102999566398 + ,0 + ,0 + ,1.15533603746506 + ,0.49136169383427 + ,1.07918124604762 + ,-0.31875876262441 + ,1.19033169817029 + ,2.14612803567824 + ,0.30102999566398 + ,0.30102999566398 + ,0.30102999566398 + ,1.18184358794477 + ,0.25527250510331 + ,1.30535136944662 + ,1 + ,2.06069784035361 + ,2.23044892137827 + ,0.60205999132796 + ,0.60205999132796 + ,0.60205999132796 + ,1 + ,-0.045757490560675 + ,1.11394335230684 + ,0.20951501454263 + ,1.05690485133647 + ,1.23044892137827 + ,0.30102999566398 + ,0 + ,0.30102999566398 + ,1.07554696139253 + ,0.25527250510331 + ,1.43136376415899 + ,2.28330122870355 + ,2.25527250510331 + ,2.06069784035361 + ,0.60205999132796 + ,0.60205999132796 + ,0.60205999132796 + ,0.81291335664286 + ,0.27875360095283 + ,1.25527250510331 + ,0.39794000867204 + ,1.08278537031645 + ,1.49136169383427 + ,0.69897000433602 + ,0.69897000433602 + ,0.69897000433602 + ,0.8750612633917 + ,-0.045757490560675 + ,0.67209785793572 + ,-0.55284196865778 + ,0.27875360095283 + ,1.32221929473392 + ,0.47712125471966 + ,0 + ,0.47712125471966 + ,1.02530586526477 + ,0.41497334797082 + ,0.99122607569249 + ,0.62736585659273 + ,1.70243053644553 + ,1.7160033436348 + ,0 + ,0 + ,0 + ,0.86923171973098 + ,0.38021124171161 + ,1.46239799789896 + ,0.83250891270624 + ,2.25285303097989 + ,2.2148438480477 + ,0.30102999566398 + ,0.47712125471966 + ,0.30102999566398 + ,0.92427928606188 + ,0.079181246047625 + ,0.84509804001426 + ,-0.1249387366083 + ,1.0899051114394 + ,2.35218251811136 + ,0.30102999566398 + ,0.30102999566398 + ,0.30102999566398 + ,0.75587485567249 + ,-0.045757490560675 + ,0.77815125038364 + ,0.55630250076729 + ,1.32221929473392 + ,2.35218251811136 + ,0.47712125471966 + ,0.30102999566398 + ,0.47712125471966 + ,0.69019608002851 + ,-0.30102999566398 + ,1.30102999566398 + ,1.74429298312268 + ,2.24303804868629 + ,2.17897694729317 + ,0.69897000433602 + ,0.69897000433602 + ,0.69897000433602 + ,0.50514997831991 + ,-0.22184874961636 + ,0.65321251377534 + ,-0.045757490560675 + ,0.41497334797082 + ,1.77815125038364 + ,0.30102999566398 + ,0 + ,0.30102999566398 + ,1.04139268515823 + ,0.36172783601759 + ,0.8750612633917 + ,0.30102999566398 + ,1.0899051114394 + ,2.30102999566398 + ,0.47712125471966 + ,0 + ,0.47712125471966 + ,0.69019608002851 + ,-0.30102999566398 + ,0.36172783601759 + ,-1 + ,0.39794000867204 + ,1.66275783168157 + ,0.47712125471966 + ,0.30102999566398 + ,0.30102999566398 + ,1.12057393120585 + ,0.41497334797082 + ,1.38021124171161 + ,0.6222140229663 + ,1.76342799356294 + ,2.32221929473392 + ,0.60205999132796 + ,0.47712125471966 + ,0.60205999132796 + ,0.98677173426624 + ,-0.22184874961636 + ,0.47712125471966 + ,0.54406804435028 + ,0.5910646070265 + ,1.14612803567824 + ,0.30102999566398 + ,0 + ,0 + ,1.10720996964787 + ,0.81954393554187) + ,dim=c(9 + ,39) + ,dimnames=list(c('logL' + ,'logWb' + ,'logWbr' + ,'logtg' + ,'logP' + ,'logS' + ,'logD' + ,'logSWS' + ,'logPS') + ,1:39)) > y <- array(NA,dim=c(9,39),dimnames=list(c('logL','logWb','logWbr','logtg','logP','logS','logD','logSWS','logPS'),1:39)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '1 1 0' > par12 = '1 0 0' > par11 = '0 0 0' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '8 9' > par4 = '1 2 3 4' > par3 = '5 6 7' > par2 = 'A A A' > par1 = 'ecol const sleep' > library(plspm) Loading required package: amap Loading required package: diagram Loading required package: shape > library(diagram) > y <- as.data.frame(t(y)) > is.data.frame(y) [1] TRUE > head(y) logL logWb logWbr logtg logP logS logD 1 0.6532125 0.0000000 0.8195439 1.623249 0.4771213 0.00000 0.4771213 2 1.8388491 3.4060289 3.6630410 2.795185 0.4771213 0.69897 0.6020600 3 1.4313638 1.0232525 2.2540645 2.255273 0.6020600 0.60206 0.6020600 4 1.2787536 -1.6989700 -0.5228787 1.544068 0.0000000 0.00000 0.0000000 5 1.4828736 2.2041200 2.2278867 2.593286 0.6020600 0.69897 0.6020600 6 1.4471580 0.5185139 1.4082400 1.799341 0.0000000 0.30103 0.0000000 logSWS logPS 1 0.7993405 0.3010300 2 0.3222193 0.2552725 3 0.9590414 -0.1549020 4 1.1986571 0.5910646 5 0.7160033 0.0000000 6 1.0374265 0.5563025 > trim <- function(char) { + return(sub('s+$', '', sub('^s+', '', char))) + } > (latnames <- strsplit(par1,' ')[[1]]) [1] "ecol" "const" "sleep" > (n <- length(latnames)) [1] 3 > (L1 <- as.numeric(strsplit(par3,' ')[[1]])) [1] 5 6 7 > (L2 <- as.numeric(strsplit(par4,' ')[[1]])) [1] 1 2 3 4 > (L3 <- as.numeric(strsplit(par5,' ')[[1]])) [1] 8 9 > (L4 <- as.numeric(strsplit(par6,' ')[[1]])) numeric(0) > (L5 <- as.numeric(strsplit(par7,' ')[[1]])) numeric(0) > (L6 <- as.numeric(strsplit(par8,' ')[[1]])) numeric(0) > (L7 <- as.numeric(strsplit(par9,' ')[[1]])) numeric(0) > (L8 <- as.numeric(strsplit(par10,' ')[[1]])) numeric(0) > (S1 <- as.numeric(strsplit(par11,' ')[[1]])) [1] 0 0 0 > (S2 <- as.numeric(strsplit(par12,' ')[[1]])) [1] 1 0 0 > (S3 <- as.numeric(strsplit(par13,' ')[[1]])) [1] 1 1 0 > (S4 <- as.numeric(strsplit(par14,' ')[[1]])) numeric(0) > (S5 <- as.numeric(strsplit(par15,' ')[[1]])) numeric(0) > (S6 <- as.numeric(strsplit(par16,' ')[[1]])) numeric(0) > (S7 <- as.numeric(strsplit(par17,' ')[[1]])) numeric(0) > (S8 <- as.numeric(strsplit(par18,' ')[[1]])) numeric(0) > if (n==1) sat.mat <- rbind(S1) > if (n==2) sat.mat <- rbind(S1,S2) > if (n==3) sat.mat <- rbind(S1,S2,S3) > if (n==4) sat.mat <- rbind(S1,S2,S3,S4) > if (n==5) sat.mat <- rbind(S1,S2,S3,S4,S5) > if (n==6) sat.mat <- rbind(S1,S2,S3,S4,S5,S6) > if (n==7) sat.mat <- rbind(S1,S2,S3,S4,S5,S6,S7) > if (n==8) sat.mat <- rbind(S1,S2,S3,S4,S5,S6,S7,S8) > sat.mat [,1] [,2] [,3] S1 0 0 0 S2 1 0 0 S3 1 1 0 > if (n==1) sat.sets <- list(L1) > if (n==2) sat.sets <- list(L1,L2) > if (n==3) sat.sets <- list(L1,L2,L3) > if (n==4) sat.sets <- list(L1,L2,L3,L4) > if (n==5) sat.sets <- list(L1,L2,L3,L4,L5) > if (n==6) sat.sets <- list(L1,L2,L3,L4,L5,L6) > if (n==7) sat.sets <- list(L1,L2,L3,L4,L5,L6,L7) > if (n==8) sat.sets <- list(L1,L2,L3,L4,L5,L6,L7,L8) > sat.sets [[1]] [1] 5 6 7 [[2]] [1] 1 2 3 4 [[3]] [1] 8 9 > (sat.mod <- strsplit(par2,' ')[[1]]) [1] "A" "A" "A" > res <- plspm(x=y, sat.mat, sat.sets, sat.mod, scheme='centroid', scaled=TRUE, boot.val=TRUE) > (r <- summary(res)) PARTIAL LEAST SQUARES PATH MODELING (PLS-PM) ---------------------------------------------------------- MODEL SPECIFICATION 1 Number of Cases 39 2 Latent Variables 3 3 Manifest Variables 9 4 Scale of Data Standardized Data 5 Weighting Scheme centroid 6 Tolerance Crit 1e-05 7 Max Num Iters 100 8 Convergence Iters 3 9 Paths by PLS-R FALSE 10 Bootstrapping TRUE 11 Bootstrap samples 100 ---------------------------------------------------------- BLOCKS DEFINITION Block Type NMVs Mode 1 S1 Exogenous 3 Reflective 2 S2 Endogenous 4 Reflective 3 S3 Endogenous 2 Reflective ---------------------------------------------------------- BLOCKS UNIDIMENSIONALITY Type.measure MVs C.alpha DG.rho eig.1st eig.2nd S1 Reflective 3 0.896 0.936 2.49 0.449 S2 Reflective 4 0.927 0.948 3.29 0.355 S3 Reflective 2 0.716 0.876 1.56 0.442 ---------------------------------------------------------- OUTER MODEL weights std.loads communal redundan S1 logP 0.227 0.868 0.753 0.000 logS 0.505 0.901 0.813 0.000 logD 0.371 0.940 0.883 0.000 S2 logL 0.193 0.838 0.703 0.156 logWb 0.301 0.927 0.859 0.191 logWbr 0.292 0.970 0.941 0.209 logtg 0.313 0.881 0.776 0.173 S3 logSWS 0.588 0.892 0.796 0.517 logPS 0.545 0.873 0.762 0.495 ---------------------------------------------------------- CORRELATIONS BETWEEN MVs AND LVs S1 S2 S3 S1 logP 0.868 0.129 -0.464 logS 0.901 0.649 -0.670 logD 0.940 0.311 -0.659 S2 logL 0.274 0.838 -0.448 logWb 0.488 0.927 -0.640 logWbr 0.451 0.970 -0.645 logtg 0.448 0.881 -0.726 S3 logSWS -0.555 -0.711 0.892 logPS -0.663 -0.509 0.873 ---------------------------------------------------------- INNER MODEL $S2 concept value 1 R2 0.2226 2 Intercept 0.0000 3 path_S1 0.4718 $S3 concept value 1 R2 0.6495 2 Intercept 0.0000 3 path_S1 -0.4628 4 path_S2 -0.4767 ---------------------------------------------------------- CORRELATIONS BETWEEN LVs S1 S2 S3 S1 1.000 0.472 -0.688 S2 0.472 1.000 -0.695 S3 -0.688 -0.695 1.000 ---------------------------------------------------------- SUMMARY INNER MODEL LV.Type Measure MVs R.square Av.Commu Av.Redun AVE S1 Exogen Rflct 3 0.000 0.816 0.000 0.816 S2 Endogen Rflct 4 0.223 0.820 0.182 0.819 S3 Endogen Rflct 2 0.649 0.779 0.506 0.779 ---------------------------------------------------------- GOODNESS-OF-FIT GoF value 1 Absolute 0.5941 2 Relative 0.7695 3 Outer.mod 0.9965 4 Inner.mod 0.7722 ---------------------------------------------------------- TOTAL EFFECTS relationships dir.effects ind.effects tot.effects 1 S1->S2 0.472 0.000 0.472 2 S1->S3 -0.463 -0.225 -0.688 3 S2->S3 -0.477 0.000 -0.477 --------------------------------------------------------- BOOTSTRAP VALIDATION weights Original Mean.Boot Std.Error perc.05 perc.95 logP 0.227 0.220 0.0666 0.0877 0.290 logS 0.505 0.512 0.0888 0.4019 0.686 logD 0.371 0.366 0.0235 0.3275 0.403 logL 0.193 0.181 0.0464 0.0914 0.238 logWb 0.301 0.305 0.0304 0.2566 0.354 logWbr 0.292 0.292 0.0259 0.2477 0.335 logtg 0.313 0.320 0.0363 0.2749 0.383 logSWS 0.588 0.584 0.0587 0.5002 0.681 logPS 0.545 0.550 0.0501 0.4736 0.634 loadings Original Mean.Boot Std.Error perc.05 perc.95 logP 0.868 0.860 0.0686 0.757 0.934 logS 0.901 0.908 0.0276 0.867 0.955 logD 0.940 0.930 0.0422 0.850 0.976 logL 0.838 0.821 0.0910 0.644 0.929 logWb 0.927 0.929 0.0268 0.877 0.962 logWbr 0.970 0.970 0.0096 0.956 0.983 logtg 0.881 0.881 0.0278 0.838 0.919 logSWS 0.892 0.890 0.0364 0.826 0.933 logPS 0.873 0.872 0.0508 0.776 0.927 paths Original Mean.Boot Std.Error perc.05 perc.95 S1->S2 0.472 0.484 0.1080 0.273 0.648 S1->S3 -0.463 -0.458 0.0860 -0.596 -0.305 S2->S3 -0.477 -0.479 0.0737 -0.584 -0.327 rsq Original Mean.Boot Std.Error perc.05 perc.95 S2 0.223 0.245 0.1012 0.0746 0.42 S3 0.650 0.656 0.0827 0.5163 0.78 total.efs Original Mean.Boot Std.Error perc.05 perc.95 S1->S2 0.472 0.484 0.1080 0.273 0.648 S1->S3 -0.688 -0.687 0.0760 -0.801 -0.555 S2->S3 -0.477 -0.479 0.0737 -0.584 -0.327 > myr <- res$path.coefs > myind <- 1 > for (j in 1:(length(sat.mat[1,])-1)) { + for (i in 1:length(sat.mat[,1])) { + if (sat.mat[i,j] == 1) { + if ((res$boot$path[myind,'perc.05'] < 0) && (res$boot$path[myind,'perc.95'] > 0)) { + myr[i,j] = 0 + } + myind = myind + 1 + } + } + } > postscript(file="/var/www/html/rcomp/tmp/1oi101292493392.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plotmat(round(myr,4), pos = NULL, curve = 0, name = latnames, + lwd = 1, box.lwd = 1, cex.txt = 1, box.type = 'circle', + box.prop = 0.5, box.cex = 1, arr.type = 'triangle', + arr.pos = 0.5, shadow.size = 0.01, prefix = '', arr.lcol = 'blue', + arr.col = 'blue', arr.width = 0.2, main = c('Inner Model', + 'Path Coefficients')) > dev.off() null device 1 > myr <- res$path.coefs > myind <- 1 > myi <- 1 > for (j in 1:(length(sat.mat[1,])-1)) { + for (i in 1:length(sat.mat[,1])) { + if (i > j) { + myr[i,j] = res$boot$total.efs[myi,'Original'] + myi = myi + 1 + if ((res$boot$total.efs[myind,'perc.05'] < 0) && (res$boot$total.efs[myind,'perc.95'] > 0)) { + myr[i,j] = 0 + } + myind = myind + 1 + } + } + } > postscript(file="/var/www/html/rcomp/tmp/29iz61292493392.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plotmat(round(myr,4), pos = NULL, curve = 0, name = latnames, + lwd = 1, box.lwd = 1, cex.txt = 1, box.type = 'circle', + box.prop = 0.5, box.cex = 1, arr.type = 'triangle', + arr.pos = 0.5, shadow.size = 0.01, prefix = '', arr.lcol = 'blue', + arr.col = 'blue', arr.width = 0.2, main = c('Inner Model', + 'Total Effects')) > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'PARTIAL LEAST SQUARES PATH MODELING (PLS-PM)',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'MODEL SPECIFICATION',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Number of Cases',header=TRUE) > a<-table.element(a,r$xxx$obs) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Latent Variables',header=TRUE) > a<-table.element(a,n) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Manifest Variables',header=TRUE) > a<-table.element(a,length(y[1,])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Scaled?',header=TRUE) > a<-table.element(a,r$xxx$scaled) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Weighting Scheme',header=TRUE) > a<-table.element(a,r$xx$scheme) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Bootstrapping?',header=TRUE) > a<-table.element(a,r$xx$boot.val) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Bootstrap samples',header=TRUE) > a<-table.element(a,r$xx$br) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/3yjwh1292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BLOCKS DEFINITION',4,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > a<-table.element(a,'Type',header=TRUE) > a<-table.element(a,'NMVs',header=TRUE) > a<-table.element(a,'Mode',header=TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],header=TRUE) + a<-table.element(a,r$input$Type[i]) + a<-table.element(a,r$unidim$MVs[i]) + a<-table.element(a,r$unidim$Type.measure[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/49sd21292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BLOCKS UNIDIMENSIONALITY',7,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > a<-table.element(a,'Type.measure',header=TRUE) > a<-table.element(a,'MVs',header=TRUE) > a<-table.element(a,'eig.1st',header=TRUE) > a<-table.element(a,'eig.2nd',header=TRUE) > a<-table.element(a,'C.alpha',header=TRUE) > a<-table.element(a,'DG.rho',header=TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],header=TRUE) + a<-table.element(a,r$unidim$Type.measure[i]) + a<-table.element(a,r$unidim$MVs[i]) + a<-table.element(a,r$unidim$eig.1st[i]) + a<-table.element(a,r$unidim$eig.2nd[i]) + a<-table.element(a,r$unidim$C.alpha[i]) + a<-table.element(a,r$unidim$DG.rho[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/5utcq1292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'OUTER MODEL',5,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > a<-table.element(a,'weights',header=TRUE) > a<-table.element(a,'std.loads',header=TRUE) > a<-table.element(a,'communal',header=TRUE) > a<-table.element(a,'redundan',header=TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],5,header=TRUE) + a<-table.row.end(a) + for (j in 1:length(r$outer.mod[[i]][,1])) { + a<-table.row.start(a) + a<-table.element(a,rownames(r$outer.mod[[i]])[j],header=T) + a<-table.element(a,r$outer.mod[[i]][j,1]) + a<-table.element(a,r$outer.mod[[i]][j,2]) + a<-table.element(a,r$outer.mod[[i]][j,3]) + a<-table.element(a,r$outer.mod[[i]][j,4]) + a<-table.row.end(a) + } + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/6gtaw1292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'CORRELATIONS BETWEEN MVs AND LVs',n+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > for (iii in 1:n) { + a<-table.element(a,latnames[iii],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],n+1,header=TRUE) + a<-table.row.end(a) + for (j in 1:length(r$outer.cor[[i]][,1])) { + a<-table.row.start(a) + a<-table.element(a,rownames(r$outer.cor[[i]])[j],header=T) + for (iii in 1:n) { + a<-table.element(a,r$outer.cor[[i]][j,iii]) + } + a<-table.row.end(a) + } + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/7jc9k1292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'INNER MODEL',3,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > a<-table.element(a,'Concept',header=TRUE) > a<-table.element(a,'Value',header=TRUE) > a<-table.row.end(a) > for (i in 1:(length(labels(r$inner.mod)))) { + a<-table.row.start(a) + print (paste('i=',i,sep='')) + a<-table.element(a,labels(r$inner.mod)[i],3,header=TRUE) + a<-table.row.end(a) + for (j in 1:length(r$inner.mod[[i]][,1])) { + print (paste('j=',j,sep='')) + a<-table.row.start(a) + a<-table.element(a,rownames(r$inner.mod[[i]])[j],header=T) + a<-table.element(a,r$inner.mod[[i]][j,1],header=T) + a<-table.element(a,r$inner.mod[[i]][j,2]) + a<-table.row.end(a) + } + } [1] "i=1" [1] "j=1" [1] "j=2" [1] "j=3" [1] "i=2" [1] "j=1" [1] "j=2" [1] "j=3" [1] "j=4" > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/8nc7q1292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'CORRELATIONS BETWEEN LVs',n+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (iii in 1:n) { + a<-table.element(a,latnames[iii],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],header=T) + for (j in 1:n) { + a<-table.element(a,r$latent.cor[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/98v6e1292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'SUMMARY INNER MODEL',8,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > a<-table.element(a,'LV.Type',header=TRUE) > a<-table.element(a,'Measure',header=TRUE) > a<-table.element(a,'MVs',header=TRUE) > a<-table.element(a,'R.square',header=TRUE) > a<-table.element(a,'Av.Commu',header=TRUE) > a<-table.element(a,'Av.Redun',header=TRUE) > a<-table.element(a,'AVE',header=TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],header=T) + a<-table.element(a,r$inner.sum[i,1]) + a<-table.element(a,r$inner.sum[i,2]) + a<-table.element(a,r$inner.sum[i,3]) + a<-table.element(a,r$inner.sum[i,4]) + a<-table.element(a,r$inner.sum[i,5]) + a<-table.element(a,r$inner.sum[i,6]) + a<-table.element(a,r$inner.sum[i,7]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/10td4j1292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'GOODNESS-OF-FIT',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'GoF',header=TRUE) > a<-table.element(a,'Value',header=TRUE) > a<-table.row.end(a) > for (i in 1:4) { + a<-table.row.start(a) + a<-table.element(a,r$gof[i,1],header=T) + a<-table.element(a,r$gof[i,2]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/11mn3m1292493392.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'TOTAL EFFECTS',4,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'relationships',header=TRUE) > a<-table.element(a,'dir.effect',header=TRUE) > a<-table.element(a,'ind.effect',header=TRUE) > a<-table.element(a,'tot.effect',header=TRUE) > a<-table.row.end(a) > for (i in 1:length(r$effects[,1])) { + a<-table.row.start(a) + a<-table.element(a,r$effects[i,1],header=T) + a<-table.element(a,r$effects[i,2]) + a<-table.element(a,r$effects[i,3]) + a<-table.element(a,r$effects[i,4]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12752a1292493392.tab") > dum <- r$boot$weights > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - WEIGHTS',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13t61y1292493392.tab") > dum <- r$boot$loadings > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - LOADINGS',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14mfij1292493392.tab") > dum <- r$boot$paths > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - PATHS',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/157yzp1292493392.tab") > dum <- r$boot$rsq > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - RSQ',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/16i7ys1292493392.tab") > dum <- r$boot$total.efs > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - TOTAL EFFECTS',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/17l7eg1292493392.tab") > try(system("convert tmp/1oi101292493392.ps tmp/1oi101292493392.png",intern=TRUE)) character(0) > try(system("convert tmp/29iz61292493392.ps tmp/29iz61292493392.png",intern=TRUE)) character(0) > > #-SERVER-193.190.124.10:1001 > > > proc.time() user system elapsed 5.972 0.236 6.352