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

39	41	55	56
44	44	54	58
48	44	56	58
45	45	57	58
45	45	59	61
47	45	61	61
47	47	63	63
47	47	62	64
50	47	65	64
50	49	65	64
50	50	64	65
50	50	62	66
50	50	65	66
51	50	65	66
53	50	66	66
51	51	65	67
52	51	66	67
52	51	67	67
53	51	66	68
58	51	68	68
51	52	69	68
52	52	70	68
52	52	71	68
53	52	68	69
53	52	69	69
52	53	66	70
53	53	68	70
53	53	69	70
54	53	69	70
56	53	70	70
53	54	75	70
54	54	69	71
55	54	71	71
56	54	71	71
58	54	74	71
50	55	69	73
53	55	69	73
54	55	74	73
54	55	75	73
55	55	73	74
55	55	70	75
55	55	76	75
55	55	76	75
57	55	76	75
59	55	71	76
59	55	75	76
60	55	76	76
60	55	79	76
54	56	80	76
56	56	76	77
56	56	77	77
56	56	78	77
57	56	80	78
57	56	79	79
166	56	78	80
52	57	80	80
55	57	83	80
56	57	79	81
56	57	81	82
56	57	85	82
57	58	83	83
59	58	85	83
54	59	83	84
57	59	82	85
59	59	84	86
59	59	88	86
59	59	88	86
60	59	89	86
63	59	88	87
63	59	87	89
58	60	84	90
61	60	90	91
61	61	90	91
61	61	88	93
62	61	96	94
62	61	96	95
62	61	97	98
63	61	101	100
62	62	92	101
62	62	103	101
63	62	102	107
64	62	119	124

cat1 <- as.numeric(par1) #
cat2<- as.numeric(par2) #
simulate.p.value=FALSE
if (par3 == 'Exact Pearson Chi Square 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])
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) )
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')