x <- array(list('HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'HI' ,'LO' ,'HI' ,'LO' ,'LO' ,'HI' ,'LO' ,'HI' ,'LO' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI' ,'HI' ,'HI' ,'LO' ,'LO' ,'LO' ,'HI') ,dim=c(6 ,162) ,dimnames=list(c('Connected' ,'Separate' ,'Learning' ,'Software' ,'Happiness' ,'Depression') ,1:162)) y <- array(NA,dim=c(6,162),dimnames=list(c('Connected','Separate','Learning','Software','Happiness','Depression'),1:162)) for (i in 1:dim(x)[1]) { for (j in 1:dim(x)[2]) { y[i,j] <- as.numeric(x[i,j]) } } par3 = 'Pearson Chi-Squared' par2 = '6' par1 = '5' main = 'Association Plot' #'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(vcd) cat1 <- as.numeric(par1) # cat2<- as.numeric(par2) # simulate.p.value=FALSE if (par3 == 'Exact Pearson Chi-Squared by Simulation') simulate.p.value=TRUE x <- t(x) (z <- array(unlist(x),dim=c(length(x[,1]),length(x[1,])))) (table1 <- table(z[,cat1],z[,cat2])) (V1<-dimnames(y)[[1]][cat1]) (V2<-dimnames(y)[[1]][cat2]) postscript(file="/var/www/html/freestat/rcomp/tmp/1tb9n1289920633.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) assoc(ftable(z[,cat1],z[,cat2],row.vars=1,dnn=c(V1,V2)),shade=T) dev.off() #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,'Tabulation of Results',ncol(table1)+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste(V1,' x ', V2),ncol(table1)+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, ' ', 1,TRUE) for(nc in 1:ncol(table1)){ a<-table.element(a, colnames(table1)[nc], 1, TRUE) } a<-table.row.end(a) for(nr in 1:nrow(table1) ){ a<-table.element(a, rownames(table1)[nr], 1, TRUE) for(nc in 1:ncol(table1) ){ a<-table.element(a, table1[nr, nc], 1, FALSE) } a<-table.row.end(a) } a<-table.end(a) table.save(a,file="/var/www/html/freestat/rcomp/tmp/2ebqs1289920633.tab") (cst<-chisq.test(table1, simulate.p.value=simulate.p.value) ) if (par3 == 'McNemar Chi-Squared') { (cst <- mcnemar.test(table1)) } if (par3 != 'McNemar Chi-Squared') { a<-table.start() a<-table.row.start(a) a<-table.element(a,'Tabulation of Expected Results',ncol(table1)+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste(V1,' x ', V2),ncol(table1)+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, ' ', 1,TRUE) for(nc in 1:ncol(table1)){ a<-table.element(a, colnames(table1)[nc], 1, TRUE) } a<-table.row.end(a) for(nr in 1:nrow(table1) ){ a<-table.element(a, rownames(table1)[nr], 1, TRUE) for(nc in 1:ncol(table1) ){ a<-table.element(a, round(cst$expected[nr, nc], digits=2), 1, FALSE) } a<-table.row.end(a) } a<-table.end(a) table.save(a,file="/var/www/html/freestat/rcomp/tmp/30cpy1289920633.tab") } a<-table.start() a<-table.row.start(a) a<-table.element(a,'Statistical Results',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, cst$method, 2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Chi Square Statistic', 1, TRUE) a<-table.element(a, round(cst$statistic, digits=2), 1,FALSE) a<-table.row.end(a) if(!simulate.p.value){ a<-table.row.start(a) a<-table.element(a, 'Degrees of Freedom', 1, TRUE) a<-table.element(a, cst$parameter, 1,FALSE) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a, 'P value', 1, TRUE) a<-table.element(a, round(cst$p.value, digits=2), 1,FALSE) a<-table.row.end(a) a<-table.end(a) table.save(a,file="/var/www/html/freestat/rcomp/tmp/4al6j1289920633.tab") try(system("convert tmp/1tb9n1289920633.ps tmp/1tb9n1289920633.png",intern=TRUE))