R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-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(1,3,2,1,1,1,2,1,1,2,2,1,1,1,1,1,1,2,1,1,1,1,1,1,1,2,1,2,1,1,1,1,1,1,1,2,3,3,3,2,1,3,1,1,1,2,2,3,2,1,2,2,3,2,3,3,3,3,2,2,3,2,3,3,4,4,4,4,4,3,3,3,4,2,1,1,2,3,3,3,2,1,2,3,3,4,1,2,1,1,3,3,3,3,2,1,2,2,2,2,2,3,2,3,2,2,3,2,3,2,2,3,2,2,2,2,2,2,2,2,3,2,2,3,2,2,2,1,1,1,2,2,1,1,1,1,1,1,4,3,2,1,2,1,1,1,3,3,3,3,1,3,2,3,3,3,3,3,3,3,3,3,1,3,2,1,2,1,2,2,2,2,3,2,1,1,2,1,2,1,1,1,2,2,5,3,2,1,3,2,2,2,3,3,3,3,1,2,2,3,3,3,3,3,2,3,3,2,2,2,2,2,2,1,1,1,2,1,1,1,1,1,1,1,2,1,1,1,1,1),dim=c(46,5),dimnames=list(c('Participantes','P1','P2','P3','P4','P5','P6','P7','P8','P9','P10','P11','P12','P13','P14','P15','P16','P17','P18','P19','P20','P21','P22','P23','P24','P25','P26','P27','P28','P29','P30','P31','P32','P33','P34','P35','P36','P37','P38','P39','P40','P41','P42','P43','P44','P45'),1:5)) > y <- array(NA,dim=c(46,5),dimnames=list(c('Participantes','P1','P2','P3','P4','P5','P6','P7','P8','P9','P10','P11','P12','P13','P14','P15','P16','P17','P18','P19','P20','P21','P22','P23','P24','P25','P26','P27','P28','P29','P30','P31','P32','P33','P34','P35','P36','P37','P38','P39','P40','P41','P42','P43','P44','P45'),1:5)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par1 = 'TRUE' > par1 <- 'TRUE' > #'GNU S' R Code compiled by R2WASP v. 1.2.327 () > #Author: root > #To cite this work: Wessa P., 2014, Cronbach alpha (v1.0.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_cronbach.wasp/ > #Source of accompanying publication: > # > library(psych) > par1b <- T > if(par1=='FALSE') par1b <- F > z <- t(y) > (r <- alpha(z,check.keys=par1b)) The determinant of the smoothed correlation was zero. This means the objective function is not defined. Chi square is based upon observed residuals. 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.stats, the correlation matrix is singular, an approximation is used In factor.scores, the correlation matrix is singular, an approximation is used I was unable to calculate the factor score weights, factor loadings used instead 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 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 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 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 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 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 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 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 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 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 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, check.keys = par1b) raw_alpha std.alpha G6(smc) average_r S/N 0.96 0.96 1 27 0.0079 Reliability if an item is dropped: raw_alpha std.alpha G6(smc) average_r S/N Participantes 0.96 0.96 1 0.40 26 P4 0.96 0.96 1 0.39 26 P5 0.96 0.96 1 0.40 26 P6 0.96 0.97 1 0.41 28 P7 0.96 0.96 1 0.40 26 P8 0.96 0.96 1 0.39 25 P9 0.96 0.96 1 0.38 25 P10 0.96 0.96 1 0.39 26 P11 0.96 0.96 1 0.38 25 P12 0.96 0.96 1 0.40 27 P13 0.96 0.96 1 0.39 26 P14 0.96 0.96 1 0.39 25 P15 0.96 0.96 1 0.39 26 P16 0.96 0.96 1 0.38 25 P17 0.96 0.96 1 0.39 26 P18 0.96 0.96 1 0.39 25 P19 0.96 0.96 1 0.38 25 P20 0.96 0.96 1 0.39 25 P21 0.96 0.96 1 0.39 25 P22 0.96 0.96 1 0.39 25 P23 0.96 0.96 1 0.39 25 P24 0.96 0.96 1 0.39 26 P25 0.96 0.96 1 0.40 26 P26 0.96 0.96 1 0.39 25 P27- 0.96 0.97 1 0.41 28 P28 0.96 0.96 1 0.41 27 P29 0.96 0.97 1 0.41 28 P30 0.96 0.96 1 0.39 26 P31 0.96 0.96 1 0.39 26 P32 0.96 0.96 1 0.38 25 P33 0.96 0.96 1 0.39 26 P34 0.96 0.96 1 0.40 26 P35- 0.96 0.96 1 0.40 26 P36- 0.96 0.96 1 0.39 26 P37- 0.96 0.96 1 0.40 27 P38- 0.96 0.96 1 0.41 27 P39 0.96 0.97 1 0.41 28 P40 0.96 0.96 1 0.41 27 P41- 0.96 0.96 1 0.39 26 P44 0.96 0.96 1 0.40 26 P45 0.96 0.97 1 0.41 28 Item statistics r r.cor r.drop Participantes 0.54 0.54 0.482 P4 0.68 0.68 0.670 P5 0.60 0.60 0.573 P6 0.17 0.17 0.145 P7 0.60 0.60 0.573 P8 0.91 0.91 0.907 P9 0.96 0.96 0.953 P10 0.70 0.70 0.701 P11 0.96 0.96 0.951 P12 0.50 0.50 0.463 P13 0.75 0.75 0.736 P14 0.84 0.84 0.814 P15 0.70 0.70 0.685 P16 0.96 0.96 0.951 P17 0.70 0.70 0.701 P18 0.91 0.91 0.909 P19 0.99 0.99 0.987 P20 0.86 0.86 0.852 P21 0.91 0.91 0.909 P22 0.91 0.91 0.909 P23 0.90 0.90 0.900 P24 0.66 0.66 0.637 P25 0.60 0.60 0.590 P26 0.80 0.80 0.788 P27- 0.23 0.23 0.213 P28 0.30 0.30 0.245 P29 0.11 0.11 0.061 P30 0.68 0.68 0.659 P31 0.65 0.65 0.604 P32 0.92 0.92 0.922 P33 0.71 0.71 0.694 P34 0.57 0.57 0.538 P35- 0.60 0.60 0.573 P36- 0.71 0.71 0.678 P37- 0.37 0.37 0.319 P38- 0.31 0.31 0.273 P39 0.24 0.24 0.187 P40 0.29 0.29 0.274 P41- 0.71 0.71 0.678 P44 0.62 0.62 0.601 P45 0.14 0.14 0.122 Non missing response frequency for each item 1 2 3 4 5 miss Participantes 0.2 0.2 0.2 0.2 0.2 0 P4 0.2 0.6 0.2 0.0 0.0 0 P5 0.4 0.6 0.0 0.0 0.0 0 P6 0.2 0.6 0.2 0.0 0.0 0 P7 0.4 0.6 0.0 0.0 0.0 0 P8 0.2 0.2 0.6 0.0 0.0 0 P9 0.0 0.2 0.8 0.0 0.0 0 P10 0.0 0.4 0.6 0.0 0.0 0 P11 0.2 0.0 0.8 0.0 0.0 0 P12 0.6 0.4 0.0 0.0 0.0 0 P13 0.2 0.6 0.2 0.0 0.0 0 P14 0.2 0.4 0.4 0.0 0.0 0 P15 0.2 0.4 0.4 0.0 0.0 0 P16 0.2 0.0 0.8 0.0 0.0 0 P17 0.0 0.4 0.6 0.0 0.0 0 P18 0.2 0.2 0.4 0.2 0.0 0 P19 0.2 0.0 0.6 0.2 0.0 0 P20 0.2 0.4 0.2 0.2 0.0 0 P21 0.2 0.2 0.4 0.2 0.0 0 P22 0.2 0.2 0.4 0.2 0.0 0 P23 0.2 0.4 0.4 0.0 0.0 0 P24 0.4 0.4 0.2 0.0 0.0 0 P25 0.0 0.6 0.4 0.0 0.0 0 P26 0.2 0.6 0.0 0.2 0.0 0 P27 0.2 0.8 0.0 0.0 0.0 0 P28 0.4 0.4 0.2 0.0 0.0 0 P29 0.8 0.2 0.0 0.0 0.0 0 P30 0.4 0.6 0.0 0.0 0.0 0 P31 0.4 0.2 0.4 0.0 0.0 0 P32 0.2 0.6 0.2 0.0 0.0 0 P33 0.4 0.4 0.2 0.0 0.0 0 P34 0.4 0.4 0.2 0.0 0.0 0 P35 0.6 0.4 0.0 0.0 0.0 0 P36 0.6 0.2 0.2 0.0 0.0 0 P37 0.6 0.0 0.4 0.0 0.0 0 P38 0.2 0.4 0.4 0.0 0.0 0 P39 0.4 0.4 0.0 0.2 0.0 0 P40 0.6 0.4 0.0 0.0 0.0 0 P41 0.6 0.2 0.2 0.0 0.0 0 P44 0.6 0.2 0.2 0.0 0.0 0 P45 0.4 0.4 0.2 0.0 0.0 0 There were 50 or more warnings (use warnings() to see the first 50) > > #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/1pew41440406038.tab") > > > > proc.time() user system elapsed 0.403 0.072 0.454