• QR Code (size=2)
• QR Code (size=3)
• QR Code (size=4)
• QR Code (size=5)
• QR Code (size=6)

14	41
18	39
11	30
12	31
16	34
18	35
14	39
14	34
15	36
15	37
17	38
19	36
10	38
16	39
18	33
14	32
14	36
17	38
14	39
16	32
18	32
11	31
14	39
12	37
17	39
9	41
16	36
14	33
15	33
11	34
16	31
13	27
17	37
15	34
14	34
16	32
9	29
15	36
17	29
13	35
15	37
16	34
16	38
12	35
12	38
11	37
15	38
15	33
17	36
13	38
16	32
14	32
11	32
12	34
12	32
15	37
16	39
15	29
12	37
12	35
8	30
13	38
11	34
14	31
15	34
10	35
11	36
12	30
15	39
15	35
14	38
16	31
15	34
15	38
13	34
12	39
17	37
13	34
15	28
13	37
15	33
16	37
15	35
16	37
15	32
14	33
15	38
14	33
13	29
7	33
17	31
13	36
15	35
14	32
13	29
16	39
12	37
14	35
17	37
15	32
17	38
12	37
16	36
11	32
15	33
9	40
16	38
15	41
10	36
10	43
15	30
11	31
13	32
14	32
18	37
16	37
14	33
14	34
14	33
14	38
12	33
14	31
15	38
15	37
15	33
13	31
17	39
17	44
19	33
15	35
13	32
9	28
15	40
15	27
15	37
16	32
11	28
14	34
11	30
15	35
13	31
15	32
16	30
14	30
15	31
16	40
16	32
11	36
12	32
9	35
16	38
13	42
16	34
12	35
9	35
13	33
13	36
14	32
19	33
13	34
12	32
13	34

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')