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(0.072 + ,0.0213 + ,0.073 + ,0.0218 + ,0.073 + ,0.0290 + ,0.073 + ,0.0263 + ,0.074 + ,0.0267 + ,0.073 + ,0.0181 + ,0.074 + ,0.0133 + ,0.074 + ,0.0088 + ,0.076 + ,0.0128 + ,0.076 + ,0.0126 + ,0.077 + ,0.0126 + ,0.077 + ,0.0129 + ,0.078 + ,0.0110 + ,0.078 + ,0.0137 + ,0.080 + ,0.0121 + ,0.081 + ,0.0174 + ,0.081 + ,0.0176 + ,0.082 + ,0.0148 + ,0.081 + ,0.0104 + ,0.081 + ,0.0162 + ,0.081 + ,0.0149 + ,0.080 + ,0.0179 + ,0.082 + ,0.0180 + ,0.084 + ,0.0158 + ,0.084 + ,0.0186 + ,0.085 + ,0.0174 + ,0.086 + ,0.0159 + ,0.085 + ,0.0126 + ,0.083 + ,0.0113 + ,0.078 + ,0.0192 + ,0.078 + ,0.0261 + ,0.080 + ,0.0226 + ,0.086 + ,0.0241 + ,0.089 + ,0.0226 + ,0.089 + ,0.0203 + ,0.086 + ,0.0286 + ,0.083 + ,0.0255 + ,0.083 + ,0.0227 + ,0.083 + ,0.0226 + ,0.084 + ,0.0257 + ,0.085 + ,0.0307 + ,0.084 + ,0.0276 + ,0.086 + ,0.0251 + ,0.085 + ,0.0287 + ,0.085 + ,0.0314 + ,0.085 + ,0.0311 + ,0.085 + ,0.0316 + ,0.085 + ,0.0247 + ,0.085 + ,0.0257 + ,0.085 + ,0.0289 + ,0.085 + ,0.0263 + ,0.086 + ,0.0238 + ,0.086 + ,0.0169 + ,0.086 + ,0.0196 + ,0.086 + ,0.0219 + ,0.084 + ,0.0187 + ,0.080 + ,0.0160 + ,0.079 + ,0.0163 + ,0.080 + ,0.0122 + ,0.080 + ,0.0121 + ,0.080 + ,0.0149 + ,0.080 + ,0.0164 + ,0.079 + ,0.0166 + ,0.079 + ,0.0177 + ,0.079 + ,0.0182 + ,0.080 + ,0.0178 + ,0.079 + ,0.0128 + ,0.075 + ,0.0129 + ,0.072 + ,0.0137 + ,0.070 + ,0.0112 + ,0.069 + ,0.0151 + ,0.071 + ,0.0224 + ,0.071 + ,0.0294 + ,0.072 + ,0.0309 + ,0.071 + ,0.0346 + ,0.069 + ,0.0364 + ,0.068 + ,0.0439 + ,0.067 + ,0.0415 + ,0.067 + ,0.0521 + ,0.069 + ,0.0580 + ,0.073 + ,0.0591 + ,0.074 + ,0.0539 + ,0.073 + ,0.0546 + ,0.071 + ,0.0472 + ,0.070 + ,0.0314 + ,0.071 + ,0.0263 + ,0.075 + ,0.0232 + ,0.077 + ,0.0193 + ,0.078 + ,0.0062 + ,0.077 + ,0.0060 + ,0.077 + ,-0.0037 + ,0.078 + ,-0.0110 + ,0.080 + ,-0.0168 + ,0.081 + ,-0.0078 + ,0.081 + ,-0.0119 + ,0.080 + ,-0.0097 + ,0.081 + ,-0.0012 + ,0.082 + ,0.0026 + ,0.083 + ,0.0062 + ,0.084 + ,0.0070 + ,0.085 + ,0.0166 + ,0.085 + ,0.0180 + ,0.085 + ,0.0227 + ,0.085 + ,0.0246 + ,0.085 + ,0.0257 + ,0.083 + ,0.0232 + ,0.082 + ,0.0291 + ,0.081 + ,0.0301 + ,0.079 + ,0.0286 + ,0.076 + ,0.0310 + ,0.073 + ,0.0322 + ,0.071 + ,0.0339 + ,0.070 + ,0.0352 + ,0.070 + ,0.0341 + ,0.070 + ,0.0335 + ,0.070 + ,0.0367 + ,0.069 + ,0.0375 + ,0.068 + ,0.0360 + ,0.067 + ,0.0355 + ,0.066 + ,0.0357) + ,dim=c(2 + ,120) + ,dimnames=list(c('Werkloosheidsgraad' + ,'Inlfatie') + ,1:120)) > y <- array(NA,dim=c(2,120),dimnames=list(c('Werkloosheidsgraad','Inlfatie'),1:120)) > 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.79 -1.2 -0.38 -0.38 0.05 0.0062 Reliability if an item is dropped: raw_alpha std.alpha G6(smc) average_r Werkloosheidsgraad -0.38 -0.38 0.15 -0.38 Inlfatie -0.38 -0.38 0.15 -0.38 Item statistics n r r.cor r.drop mean sd Werkloosheidsgraad 120 0.56 NaN -0.38 0.078 0.0059 Inlfatie 120 0.56 NaN -0.38 0.022 0.0133 Warning message: In sqrt(Vtc) : NaNs produced > > #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/10p4g1324416826.tab") > > > > proc.time() user system elapsed 0.301 0.050 0.342