R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. Natural language support but running in an English locale 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(1 + ,1 + ,2 + ,2 + ,1 + ,3 + ,2 + ,1 + ,1 + ,1 + ,3 + ,1 + ,4 + ,2 + ,1 + ,1 + ,2 + ,2 + ,3 + ,3 + ,2 + ,2 + ,1 + ,2 + ,1 + ,2 + ,3 + ,1 + ,2 + ,2 + ,3 + ,3 + ,3 + ,1 + ,3 + ,1 + ,1 + ,2 + ,2 + ,2 + ,4 + ,2 + ,3 + ,4 + ,4 + ,1 + ,3 + ,2 + ,3 + ,1 + ,1 + ,2 + ,1 + ,2 + ,3 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,1 + ,2 + ,1 + ,2 + ,2 + ,3 + ,2 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,3 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,1 + ,2 + ,3 + ,2 + ,3 + ,3 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,2 + ,1 + ,2 + ,1 + ,3 + ,3 + ,2 + ,2 + ,1 + ,2 + ,2 + ,1 + ,3 + ,3 + ,1 + ,1 + ,1 + ,2 + ,1 + ,4 + ,3 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,1 + ,2 + ,1 + ,1 + ,1 + ,3 + ,4 + ,1 + ,1 + ,1 + ,1 + ,0 + ,4 + ,1 + ,1 + ,1 + ,2 + ,2 + ,2 + ,3 + ,1 + ,2 + ,2 + ,2 + ,2 + ,2 + ,4 + ,2 + ,1 + ,1 + ,1 + ,1 + ,1 + ,3 + ,3 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,1 + ,2 + ,1 + ,2 + ,3 + ,3 + ,3 + ,3 + ,1 + ,2 + ,1 + ,1 + ,3 + ,2 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1 + ,1 + ,2 + ,2 + ,3 + ,1 + ,3 + ,3 + ,1 + ,1 + ,2 + ,2 + ,1 + ,4 + ,1 + ,1 + ,1 + ,1 + ,2 + ,1 + ,3 + ,2 + ,1 + ,1 + ,2 + ,1 + ,2 + ,3 + ,3 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,3 + ,2 + ,2 + ,3 + ,2 + ,2 + ,4 + ,2 + ,1 + ,2 + ,1 + ,2 + ,4 + ,2 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,1 + ,1 + ,1 + ,2 + ,1 + ,3 + ,1 + ,2 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,2 + ,1 + ,2 + ,1 + ,3 + ,4 + ,3 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,3 + ,1 + ,2 + ,1 + ,2 + ,2 + ,1 + ,3 + ,2 + ,1 + ,1 + ,1 + ,2 + ,2 + ,3 + ,1 + ,1 + ,1 + ,2 + ,2 + ,1 + ,4 + ,2 + ,1 + ,1 + ,2 + ,1 + ,1 + ,3 + ,1 + ,2 + ,1 + ,2 + ,2 + ,3 + ,3 + ,1 + ,2 + ,2 + ,3 + ,4 + ,3 + ,4 + ,2 + ,1 + ,2 + ,2 + ,2 + ,2 + ,3 + ,2 + ,1 + ,1 + ,1 + ,1 + ,3 + ,4 + ,2 + ,2 + ,1 + ,1 + ,3 + ,2 + ,3 + ,0 + ,1 + ,1 + ,2 + ,3 + ,2 + ,3 + ,1 + ,1 + ,2 + ,1 + ,1 + ,1 + ,4 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,3 + ,1 + ,1 + ,2 + ,3 + ,3 + ,3 + ,2 + ,3 + ,2 + ,3 + ,2 + ,3 + ,2 + ,3 + ,1 + ,1 + ,2 + ,3 + ,2 + ,4 + ,1 + ,2 + ,1 + ,2 + ,3 + ,1 + ,4 + ,1 + ,1 + ,1 + ,2 + ,3 + ,2 + ,2 + ,1 + ,1 + ,1 + ,2 + ,2 + ,2 + ,3 + ,2 + ,1 + ,1 + ,2 + ,1 + ,1 + ,3 + ,2 + ,2 + ,2 + ,2 + ,3 + ,2 + ,2 + ,2 + ,1 + ,1 + ,1 + ,3 + ,1 + ,4 + ,3 + ,1 + ,1 + ,2 + ,3 + ,2 + ,3 + ,1 + ,2 + ,2 + ,2 + ,2 + ,2 + ,3 + ,1 + ,2 + ,1 + ,2 + ,1 + ,2 + ,3 + ,2 + ,1 + ,1 + ,2 + ,2 + ,4 + ,3 + ,2 + ,1 + ,1 + ,1 + ,3 + ,3 + ,3 + ,1 + ,1 + ,1 + ,2 + ,1 + ,1 + ,4 + ,2 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,1 + ,2 + ,2 + ,1 + ,2 + ,3 + ,2 + ,2 + ,2 + ,3 + ,2 + ,2 + ,2 + ,2 + ,2 + ,1 + ,1 + ,2 + ,2 + ,4 + ,3 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1 + ,3 + ,1 + ,2 + ,1 + ,1 + ,4 + ,1 + ,1 + ,1 + ,1 + ,3 + ,1 + ,4 + ,1 + ,1 + ,1 + ,1 + ,2 + ,1 + ,4 + ,1 + ,1 + ,1 + ,2 + ,1 + ,4 + ,4 + ,1 + ,1 + ,1 + ,3 + ,2 + ,2 + ,4 + ,1 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,4 + ,1 + ,1 + ,2 + ,1 + ,3 + ,3 + ,2 + ,2 + ,3 + ,3 + ,3 + ,4 + ,3 + ,4 + ,1 + ,1 + ,2 + ,2 + ,1 + ,4 + ,4 + ,2 + ,1 + ,2 + ,1 + ,3 + ,3 + ,3 + ,2 + ,1 + ,2 + ,1 + ,1 + ,3 + ,3 + ,2 + ,1 + ,2 + ,2 + ,2 + ,2 + ,2 + ,1 + ,1 + ,0 + ,2 + ,1 + ,4 + ,3 + ,1 + ,1 + ,2 + ,1 + ,1 + ,4 + ,3 + ,2 + ,1 + ,2 + ,2 + ,3 + ,3 + ,3 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,1 + ,1 + ,2 + ,2 + ,2 + ,4 + ,1 + ,1 + ,1 + ,2 + ,1 + ,2 + ,4 + ,1 + ,1 + ,1 + ,2 + ,1 + ,2 + ,4 + ,1 + ,2 + ,2 + ,2 + ,3 + ,2 + ,4 + ,3 + ,1 + ,1 + ,1 + ,2 + ,2 + ,3 + ,2 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,4 + ,2 + ,1 + ,2 + ,3 + ,2 + ,4 + ,4 + ,3 + ,2 + ,3 + ,2 + ,3 + ,3 + ,3 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,3 + ,2 + ,2 + ,2 + ,1 + ,2 + ,1 + ,3 + ,3 + ,3 + ,3 + ,4 + ,2 + ,4 + ,2 + ,2 + ,2 + ,2 + ,2 + ,3 + ,1 + ,2 + ,2 + ,2 + ,2 + ,2 + ,3 + ,2 + ,2 + ,1 + ,1 + ,1 + ,1 + ,3 + ,2 + ,3 + ,1 + ,2 + ,1 + ,2 + ,4 + ,4 + ,2 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,1 + ,1 + ,1 + ,1 + ,3 + ,4 + ,1 + ,1 + ,1 + ,2 + ,2 + ,2 + ,2 + ,2 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1 + ,1 + ,1 + ,1 + ,3 + ,4 + ,4 + ,3 + ,2 + ,2 + ,3 + ,3 + ,1 + ,3 + ,1 + ,2 + ,1 + ,1 + ,2 + ,2 + ,3 + ,3 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,3 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,2 + ,1 + ,1 + ,1 + ,1 + ,3 + ,2 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,1 + ,1 + ,1 + ,1 + ,1 + ,1 + ,4 + ,3 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,2 + ,1 + ,1 + ,1 + ,2 + ,3 + ,3 + ,1 + ,2 + ,1 + ,1 + ,1 + ,3 + ,2 + ,3 + ,2 + ,2 + ,3 + ,3 + ,3 + ,4 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,2 + ,1 + ,2 + ,1 + ,3 + ,4 + ,1 + ,1 + ,1 + ,2 + ,2 + ,3 + ,4 + ,3 + ,4 + ,4 + ,4 + ,3 + ,2 + ,2 + ,2 + ,2 + ,1 + ,2 + ,2 + ,2 + ,3 + ,2 + ,2 + ,2 + ,2 + ,2 + ,2 + ,3 + ,3 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,3 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,3 + ,1 + ,1 + ,1 + ,4 + ,1 + ,4 + ,4 + ,1 + ,1 + ,1 + ,1 + ,3 + ,3 + ,2 + ,1 + ,1 + ,2 + ,2 + ,2 + ,3 + ,1 + ,2 + ,1 + ,1 + ,2 + ,2 + ,4 + ,1 + ,1 + ,1 + ,2 + ,1 + ,3 + ,4 + ,2 + ,1 + ,1 + ,1 + ,1 + ,3 + ,3 + ,3 + ,2 + ,2 + ,3 + ,1 + ,1 + ,2 + ,4 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,3 + ,2 + ,3 + ,3 + ,2 + ,2 + ,1 + ,2 + ,2 + ,1 + ,2 + ,1 + ,3 + ,3 + ,3 + ,2 + ,1 + ,1 + ,2 + ,1 + ,3 + ,4 + ,1 + ,1 + ,1 + ,2 + ,2 + ,4 + ,2 + ,1 + ,1 + ,1 + ,1 + ,2 + ,4 + ,1) + ,dim=c(7 + ,142) + ,dimnames=list(c('Q31_1' + ,'Q31_2' + ,'Q31_3' + ,'Q31_4' + ,'Q31_5' + ,'Q31_6' + ,'Q31_7') + ,1:142)) > y <- array(NA,dim=c(7,142),dimnames=list(c('Q31_1','Q31_2','Q31_3','Q31_4','Q31_5','Q31_6','Q31_7'),1:142)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Dr. Ian E. Holliday > #To cite this work: Ian E. Holliday, 2009, YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: > #Technical description: > library(psych) > z <- t(y) > (r <- alpha(z)) Reliability analysis Call: alpha(x = z) raw_alpha std.alpha G6(smc) average_r mean sd 0.38 0.4 0.54 0.088 1.9 0.35 Reliability if an item is dropped: raw_alpha std.alpha G6(smc) average_r Q31_1 0.22 0.22 0.40 0.045 Q31_2 0.22 0.19 0.36 0.037 Q31_3 0.21 0.20 0.35 0.040 Q31_4 0.32 0.34 0.52 0.079 Q31_5 0.29 0.32 0.50 0.072 Q31_6 0.64 0.68 0.68 0.262 Q31_7 0.32 0.35 0.52 0.081 Item statistics n r r.cor r.drop mean sd Q31_1 142 0.66 0.65 0.41 1.5 0.66 Q31_2 142 0.70 0.74 0.47 1.3 0.57 Q31_3 142 0.69 0.74 0.42 1.6 0.71 Q31_4 142 0.51 0.31 0.22 1.8 0.76 Q31_5 142 0.54 0.37 0.26 1.9 0.81 Q31_6 142 -0.33 -0.71 -0.53 3.4 0.70 Q31_7 142 0.50 0.30 0.21 2.0 0.99 Non missing response frequency for each item 0 1 2 3 4 miss Q31_1 0.00 0.63 0.29 0.07 0.01 0 Q31_2 0.00 0.78 0.18 0.02 0.01 0 Q31_3 0.01 0.46 0.44 0.08 0.01 0 Q31_4 0.00 0.39 0.44 0.16 0.01 0 Q31_5 0.01 0.31 0.47 0.18 0.04 0 Q31_6 0.00 0.01 0.08 0.40 0.50 0 Q31_7 0.01 0.35 0.32 0.23 0.09 0 > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/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/www/html/freestat/rcomp/tmp/1ei2k1291457847.tab") > > > > proc.time() user system elapsed 0.731 0.025 0.747