Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
9 23 11 16 12 33 12 32 7 37 12 14 12 52 12 75 10 72 15 15 10 29 15 13 10 40 15 19 9 24 15 121 12 93 13 36 12 23 12 85 8 41 9 46 15 18 12 35 12 17 15 4 11 28 12 44 6 10 14 38 12 57 12 23 12 36 11 22 12 40 12 31 12 11 12 38 8 24 8 37 12 37 12 22 11 15 10 2 11 43 12 31 13 29 12 45 12 25 10 4 10 31 11 -4 8 66 12 61 9 32 12 31 9 39 11 19 15 31 8 36 8 42 11 21 11 21 11 25 13 32 7 26 12 28 8 32 8 41 4 29 11 33 10 17 7 13 12 32 11 30 9 34 10 59 8 13 8 23 11 10 12 5 10 31 10 19 12 32 8 30 11 25 8 48 10 35 14 67 9 15 9 22 10 18 13 33 12 46 13 24 8 14 3 12 8 38 12 12 11 28 9 41 12 12 12 31 12 33 10 34 13 21 9 20 12 44 11 52 14 7 11 29 9 11 12 26 8 24 15 7 12 60 14 13 12 20 9 52 9 28 13 25 13 39 15 9 11 19 7 13 10 60 11 19 14 34 14 14 13 17 12 45 8 66 13 48 9 29 12 -2 13 51 11 2 11 24 13 40 12 20 12 19 10 16 9 20 10 40 13 27 13 25 9 49 11 39 12 61 8 19 12 67 12 45 12 30 9 8 12 19 12 52 11 22 12 17 6 33 7 34 10 22 12 30 10 25 12 38 9 26
Names of X columns:
CONFSOFTTOT PRH
Factor 1
Factor 2
Type of test to use
(?)
Pearson Chi-Squared
Pearson Chi-Squared
Exact Pearson Chi-Squared by Simulation
McNemar Chi-Squared
Fisher Exact Test
Chart options
Title:
R Code
par3 <- 'Pearson Chi-Squared' par2 <- '2' par1 <- '1' 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=='Fisher Exact Test') { (cst <- fisher.test(table1)) } if ((par3 != 'McNemar Chi-Squared') & (par3 != 'Fisher Exact Test')) { 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) if (par3=='Pearson Chi-Squared') a<-table.element(a, 'Pearson Chi Square Statistic', 1, TRUE) if (par3=='Exact Pearson Chi-Squared by Simulation') a<-table.element(a, 'Exact Pearson Chi Square Statistic', 1, TRUE) if (par3=='McNemar Chi-Squared') a<-table.element(a, 'McNemar Chi Square Statistic', 1, TRUE) if (par3=='Fisher Exact Test') a<-table.element(a, 'Odds Ratio', 1, TRUE) if (par3=='Fisher Exact Test') { if ((ncol(table1) == 2) & (nrow(table1) == 2)) { a<-table.element(a, round(cst$estimate, digits=2), 1,FALSE) } else { a<-table.element(a, '--', 1,FALSE) } } else { a<-table.element(a, round(cst$statistic, digits=2), 1,FALSE) } a<-table.row.end(a) if(!simulate.p.value){ if(par3!='Fisher Exact Test') { 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')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation