x <- array(list(4
,7
,5
,7
,5
,5
,5
,5
,4
,6
,5
,5
,3
,6
,6
,6
,6
,6
,6
,6
,5
,5
,5
,7
,5
,7
,5
,6
,1
,6
,6
,7
,4
,7
,7
,7
,5
,6
,6
,6
,6
,7
,7
,7
,7
,7
,6
,7
,5
,6
,6
,5
,6
,6
,6
,6
,5
,7
,6
,7
,4
,6
,6
,6
,7
,7
,7
,7
,7
,7
,7
,7
,6
,7
,7
,7
,6
,7
,6
,6
,2
,5
,5
,4
,7
,7
,7
,7
,5
,5
,5
,6
,4
,6
,7
,7
,7
,7
,6
,7
,1
,6
,2
,5
,1
,6
,2
,5
,7
,6
,7
,7
,4
,6
,4
,6
,5
,6
,4
,6
,5
,6
,5
,6
,5
,5
,5
,6
,4
,3
,4
,1
,5
,3
,7
,5
,5
,5
,7
,7
,4
,6
,4
,5
,7
,7
,7
,7
,7
,7
,4
,5
,7
,7
,6
,7
,4
,6
,6
,6
,7
,7
,6
,7
,7
,7
,6
,6
,4
,6
,6
,6
,3
,6
,7
,5
,2
,5
,6
,6
,6
,6
,6
,6
,4
,6
,4
,5
,5
,7
,7
,7
,4
,6
,6
,6
,7
,7
,7
,7
,4
,6
,6
,5
,6
,5
,7
,5
,7
,7
,7
,7
,1
,4
,2
,3
,5
,7
,6
,6
,4
,5
,6
,6
,5
,6
,5
,5
,5
,7
,6
,6
,5
,7
,6
,6
,5
,7
,6
,7
,5
,7
,6
,7
,5
,6
,6
,6
,6
,6
,4
,5
,3
,6
,3
,4
,4
,4
,4
,5
,6
,6
,7
,6
,6
,6
,6
,6
,3
,7
,6
,7
,5
,7
,7
,6
,2
,7
,7
,7
,7
,7
,7
,7
,7
,7
,5
,7
,4
,6
,6
,6
,6
,5
,3
,5
,6
,7
,6
,7
,3
,5
,5
,5
,3
,5
,5
,5
,6
,7
,6
,6
,7
,7
,6
,7
,3
,2
,4
,3
,1
,6
,4
,4
,5
,6
,5
,5
,5
,6
,6
,6
,5
,6
,5
,6
,5
,5
,5
,6
,6
,6
,6
,6
,6
,7
,7
,7
,5
,5
,5
,7
,7
,7
,7
,7
,6
,6
,7
,7
,1
,5
,2
,2
,3
,6
,5
,5
,5
,7
,4
,6
,1
,7
,7
,6
,6
,6
,6
,6
,4
,7
,7
,7
,5
,6
,5
,5
,5
,5
,5
,5
,6
,5
,4
,6
,5
,6
,7
,6
,5
,6
,5
,7
,4
,6
,4
,4
,6
,6
,6
,6
,6
,6
,6
,6
,4
,5
,5
,5
,5
,5
,5
,5
,5
,5
,3
,5
,2
,7
,5
,6
,7
,7
,6
,7
,5
,6
,6
,6
,5
,7
,7
,7
,2
,7
,7
,7
,3
,4
,5
,5
,5
,5
,5
,6
,5
,7
,6
,7
,5
,7
,4
,7
,6
,7
,7
,7
,6
,6
,5
,7
,4
,5
,6
,5
,6
,7
,7
,7
,3
,4
,7
,7
,6
,6
,6
,6
,4
,6
,6
,5
,3
,3
,4
,5
,4
,7
,7
,7
,6
,6
,6
,7
,4
,7
,5
,6
,6
,6
,6
,6
,5
,5
,5
,6
,5
,6
,5
,6
,5
,6
,6
,6
,3
,7
,6
,6
,5
,6
,4
,5
,4
,7
,6
,6
,5
,6
,7
,7
,5
,6
,6
,6
,1
,7
,7
,7
,4
,7
,7
,7
,7
,7
,6
,6
,4
,7
,5
,6
,6
,6
,5
,7
,7
,7
,5
,5
,6
,7
,7
,5
,5
,6
,3
,6
,6
,7
,6
,7
,5
,6
,5
,5
,5
,7
,7
,7
,6
,6
,4
,5
,5
,5
,6
,6
,3
,4
,4
,5
,6
,7
,6
,7
,7
,7
,7
,7
,4
,5
,4
,5
,5
,7
,6
,7
,4
,6
,7
,6
,5
,5
,5
,6
,2
,6
,6
,6
,7
,7
,4
,7
,5
,7
,4
,6
,4
,6
,6
,5
,2
,5
,4
,5
,4
,6
,3
,5)
,dim=c(4
,162)
,dimnames=list(c('Q1'
,'Q8'
,'Q15'
,'Q22')
,1:162))
 y <- array(NA,dim=c(4,162),dimnames=list(c('Q1','Q8','Q15','Q22'),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))

#Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/www/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/rcomp/tmp/1pq5n1318954792.tab")