Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
12 38 11 32 14 35 12 33 21 37 12 29 22 31 11 36 10 35 13 38 10 31 8 34 15 35 14 38 10 37 14 33 14 32 11 38 10 38 13 32 7 33 14 31 12 38 14 39 11 32 9 32 11 35 15 37 14 33 13 33 9 28 15 32 10 31 11 37 13 30 8 33 20 31 12 33 10 31 10 33 9 32 14 33 8 32 14 33 11 28 13 35 9 39 11 34 15 38 11 32 10 38 14 30 18 33 14 38 11 32 12 32 13 34 9 34 10 36 15 34 20 28 12 34 12 35 14 35 13 31 11 37 17 35 12 27 13 40 14 37 13 36 15 38 13 39 10 41 11 27 19 30 13 37 17 31 13 31 9 27 11 36 10 38 9 37 12 33 12 34 13 31 13 39 12 34 15 32 22 33 13 36 15 32 13 41 15 28 10 30 11 36 16 35 11 31 11 34 10 36 10 36 16 35 12 37 11 28 16 39 19 32 11 35 16 39 15 35 24 42 14 34 15 33 11 41 15 33 12 34 10 32 14 40 13 40 9 35 15 36 15 37 14 27 11 39 8 38 11 31 11 33 8 32 10 39 11 36 13 33 11 33 20 32 10 37 15 30 12 38 14 29 23 22 14 35 16 35 11 34 12 35 10 34 14 34 12 35 12 23 11 31 12 27 13 36 11 31 19 32 12 39 17 37 9 38 12 39 19 34 18 31 15 32 14 37 11 36 9 32 18 35 16 36
Names of X columns:
depression seperate
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
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')
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