| 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])
 bitmap(file='pic1.png')
 assoc(ftable(z[,cat1],z[,cat2],row.vars=1,dnn=c(V1,V2)),shade=T)
 dev.off()
 load(file='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='mytable.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='mytable1.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='mytable2.tab')
 
 |