R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-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(6 + ,6 + ,6 + ,5 + ,4 + ,4 + ,4 + ,4 + ,6 + ,2 + ,2 + ,5 + ,4 + ,2 + ,2 + ,3 + ,2 + ,2 + ,2 + ,2 + ,5 + ,5 + ,5 + ,4 + ,1 + ,1 + ,1 + ,1 + ,6 + ,4 + ,4 + ,5 + ,3 + ,5 + ,5 + ,4 + ,4 + ,4 + ,4 + ,3 + ,3 + ,1 + ,1 + ,5 + ,2 + ,2 + ,2 + ,4 + ,4 + ,4 + ,4 + ,3 + ,3 + ,2 + ,2 + ,1 + ,5 + ,6 + ,6 + ,6 + ,3 + ,3 + ,3 + ,2 + ,3 + ,2 + ,2 + ,2 + ,6 + ,6 + ,6 + ,6 + ,1 + ,2 + ,2 + ,1 + ,1 + ,6 + ,6 + ,4 + ,1 + ,1 + ,1 + ,2 + ,5 + ,4 + ,4 + ,5 + ,2 + ,1 + ,1 + ,2 + ,3 + ,4 + ,4 + ,3 + ,3 + ,2 + ,2 + ,4 + ,5 + ,5 + ,5 + ,1 + ,5 + ,5 + ,5 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,1 + ,1 + ,4 + ,3 + ,4 + ,4 + ,4 + ,5 + ,4 + ,4 + ,4 + ,3 + ,5 + ,5 + ,5 + ,5 + ,5 + ,5 + ,6 + ,3 + ,3 + ,3 + ,3 + ,2 + ,1 + ,1 + ,3 + ,4 + ,5 + ,5 + ,4 + ,3 + ,2 + ,2 + ,1 + ,3 + ,3 + ,3 + ,5 + ,6 + ,6 + ,6 + ,6 + ,6 + ,4 + ,4 + ,2 + ,3 + ,4 + ,4 + ,4 + ,3 + ,2 + ,2 + ,3 + ,4 + ,3 + ,3 + ,4 + ,4 + ,5 + ,5 + ,5 + ,2 + ,5 + ,5 + ,2 + ,2 + ,2 + ,2 + ,3 + ,3 + ,2 + ,2 + ,2 + ,2 + ,1 + ,1 + ,4 + ,5 + ,4 + ,4 + ,5 + ,5 + ,1 + ,1 + ,2 + ,6 + ,4 + ,4 + ,4 + ,3 + ,4 + ,4 + ,3 + ,3 + ,5 + ,5 + ,4 + ,2 + ,4 + ,4 + ,2 + ,5 + ,5 + ,5 + ,5 + ,2 + ,2 + ,2 + ,4 + ,4 + ,3 + ,3 + ,4 + ,4 + ,2 + ,2 + ,2 + ,1 + ,1 + ,1 + ,1 + ,6 + ,5 + ,5 + ,4 + ,2 + ,2 + ,2 + ,4 + ,3 + ,3 + ,3 + ,3 + ,4 + ,4 + ,4 + ,4 + ,3 + ,2 + ,2 + ,5 + ,4 + ,5 + ,5 + ,5 + ,3 + ,2 + ,2 + ,5 + ,6 + ,6 + ,6 + ,6 + ,4 + ,5 + ,5 + ,4 + ,5 + ,4 + ,4 + ,5 + ,3 + ,2 + ,2 + ,2 + ,4 + ,4 + ,4 + ,6 + ,3 + ,5 + ,5 + ,5 + ,3 + ,4 + ,4 + ,2 + ,6 + ,6 + ,6 + ,6 + ,6 + ,6 + ,6 + ,5 + ,2 + ,4 + ,4 + ,4 + ,2 + ,4 + ,4 + ,4 + ,6 + ,3 + ,3 + ,4 + ,5 + ,6 + ,6 + ,7 + ,1 + ,2 + ,2 + ,1 + ,2 + ,2 + ,2 + ,2 + ,5 + ,2 + ,2 + ,3 + ,3 + ,3 + ,3 + ,3 + ,3 + ,2 + ,2 + ,4 + ,4 + ,4 + ,4 + ,4 + ,6 + ,5 + ,5 + ,5 + ,4 + ,4 + ,4 + ,4 + ,2 + ,3 + ,3 + ,3 + ,4 + ,4 + ,4 + ,2 + ,4 + ,3 + ,3 + ,3 + ,2 + ,2 + ,2 + ,3 + ,4 + ,5 + ,5 + ,5 + ,3 + ,3 + ,3 + ,4 + ,6 + ,5 + ,5 + ,4 + ,4 + ,3 + ,3 + ,3 + ,4 + ,3 + ,3 + ,3 + ,3 + ,4 + ,4 + ,4 + ,3 + ,3 + ,3 + ,4 + ,3 + ,3 + ,3 + ,3 + ,4 + ,6 + ,6 + ,6 + ,2 + ,3 + ,3 + ,3 + ,2 + ,1 + ,1 + ,1 + ,3 + ,5 + ,5 + ,5 + ,3 + ,6 + ,6 + ,5 + ,4 + ,4 + ,4 + ,3 + ,4 + ,4 + ,4 + ,4 + ,2 + ,5 + ,5 + ,2 + ,5 + ,4 + ,4 + ,6 + ,5 + ,4 + ,4 + ,5 + ,4 + ,3 + ,3 + ,5 + ,4 + ,1 + ,1 + ,1 + ,3 + ,2 + ,2 + ,2 + ,4 + ,3 + ,3 + ,3 + ,4 + ,2 + ,2 + ,3 + ,5 + ,5 + ,5 + ,5 + ,4 + ,2 + ,2 + ,3 + ,5 + ,2 + ,2 + ,3 + ,3 + ,4 + ,4 + ,4 + ,2 + ,2 + ,2 + ,1 + ,1 + ,3 + ,3 + ,3 + ,5 + ,4 + ,4 + ,4 + ,5 + ,3 + ,3 + ,6 + ,3 + ,1 + ,1 + ,2 + ,4 + ,5 + ,5 + ,5 + ,4 + ,1 + ,1 + ,3 + ,5 + ,4 + ,4 + ,4 + ,5 + ,5 + ,5 + ,6 + ,3 + ,3 + ,3 + ,2 + ,4 + ,4 + ,4 + ,4 + ,6 + ,3 + ,3 + ,5 + ,4 + ,3 + ,3 + ,4 + ,5 + ,6 + ,6 + ,5 + ,4 + ,2 + ,2 + ,2 + ,4 + ,4 + ,4 + ,3 + ,5 + ,5 + ,5 + ,6 + ,4 + ,4 + ,4 + ,3 + ,3 + ,4 + ,4 + ,5 + ,3 + ,3 + ,3 + ,4 + ,6 + ,4 + ,4 + ,5 + ,5 + ,5 + ,5 + ,5 + ,2 + ,3 + ,3 + ,1 + ,4 + ,5 + ,5 + ,4 + ,3 + ,2 + ,2 + ,2 + ,3 + ,6 + ,6 + ,5 + ,5 + ,6 + ,6 + ,7 + ,3 + ,2 + ,2 + ,3 + ,3 + ,3 + ,3 + ,4 + ,4 + ,4 + ,4 + ,4 + ,5 + ,4 + ,4 + ,5 + ,5 + ,4 + ,4 + ,4 + ,5 + ,4 + ,4 + ,1 + ,5 + ,3 + ,3 + ,2 + ,3 + ,7 + ,7 + ,5 + ,4 + ,4 + ,4 + ,2 + ,4 + ,4 + ,4 + ,4 + ,3 + ,4 + ,4 + ,5 + ,4 + ,4 + ,4 + ,6 + ,4 + ,4 + ,4 + ,4 + ,3 + ,3 + ,3 + ,4 + ,2 + ,5 + ,5 + ,6 + ,4 + ,3 + ,3 + ,4 + ,4 + ,5 + ,5 + ,4) + ,dim=c(4 + ,162) + ,dimnames=list(c('Q4' + ,'Q11' + ,'Q18' + ,'Q25') + ,1:162)) > y <- array(NA,dim=c(4,162),dimnames=list(c('Q4','Q11','Q18','Q25'),1:162)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > library(psych) > z <- t(y) > (r <- alpha(z)) In smc, the correlation matrix was not invertible, smc's returned as 1s In smc, the correlation matrix was not invertible, smc's returned as 1s In smc, the correlation matrix was not invertible, smc's returned as 1s In smc, the correlation matrix was not invertible, smc's returned as 1s Reliability analysis Call: alpha(x = z) raw_alpha std.alpha G6(smc) average_r mean sd 0.86 0.86 1 0.6 3.6 1.2 Reliability if an item is dropped: raw_alpha std.alpha G6(smc) average_r Q4 0.89 0.89 1.00 0.73 Q11 0.76 0.76 0.68 0.51 Q18 0.76 0.76 0.68 0.51 Q25 0.85 0.85 1.00 0.65 Item statistics n r r.cor r.drop mean sd Q4 162 0.72 0.72 0.52 3.7 1.3 Q11 162 0.92 0.92 0.85 3.5 1.4 Q18 162 0.92 0.92 0.85 3.5 1.4 Q25 162 0.79 0.79 0.62 3.7 1.4 Non missing response frequency for each item 1 2 3 4 5 6 7 miss Q4 0.05 0.14 0.27 0.27 0.18 0.10 0.00 0 Q11 0.09 0.19 0.18 0.28 0.17 0.09 0.01 0 Q18 0.09 0.19 0.18 0.28 0.17 0.09 0.01 0 Q25 0.08 0.15 0.18 0.29 0.20 0.09 0.01 0 > > #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,'Cronbach Alpha and Related Statistics',5,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Items',1,TRUE) > a<-table.element(a,'Cronbach Alpha',1,TRUE) > a<-table.element(a,'Std. Alpha',1,TRUE) > a<-table.element(a,'G6(smc)',1,TRUE) > a<-table.element(a,'Average R',1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'All itmes',header=TRUE) > a<-table.element(a,round(r$total['raw_alpha'],4)) > a<-table.element(a,round(r$total['std.alpha'],4)) > a<-table.element(a,round(r$total['G6(smc)'],4)) > a<-table.element(a,round(r$total['average_r'],4)) > a<-table.row.end(a) > n = length(z[1,]) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,paste(colnames(z)[i],' excluded',sep=''),header=TRUE) + a<-table.element(a,round(r$alpha.drop[i,'raw_alpha'],4)) + a<-table.element(a,round(r$alpha.drop[i,'std.alpha'],4)) + a<-table.element(a,round(r$alpha.drop[i,'G6(smc)'],4)) + a<-table.element(a,round(r$alpha.drop[i,'average_r'],4)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/1gyta1322663864.tab") > > > > proc.time() user system elapsed 0.378 0.038 0.464