Wessa.net

Send output to:

Data X:
'B' 'HI' 'A' 'LO' 'D' 'HI' 'D' 'LO' 'A' 'HI' 'A' 'LO' 'B' 'LO' 'B' 'HI' 'B' 'HI' 'B' 'HI' 'A' 'LO' 'A' 'HI' 'D' 'HI' 'A' 'HI' 'A' 'HI' 'B' 'LO' 'B' 'LO' 'A' 'HI' 'B' 'HI' 'A' 'LO' 'A' 'LO' 'D' 'LO' 'B' 'HI' 'D' 'HI' 'A' 'LO' 'D' 'LO' 'A' 'HI' 'B' 'HI' 'B' 'LO' 'D' 'LO' 'A' 'LO' 'C' 'LO' 'A' 'LO' 'B' 'HI' 'B' 'LO' 'A' 'LO' 'D' 'LO' 'B' 'LO' 'A' 'LO' 'C' 'LO' 'B' 'LO' 'A' 'LO' 'A' 'LO' 'D' 'LO' 'D' 'LO' 'D' 'HI' 'B' 'HI' 'B' 'HI' 'A' 'HI' 'C' 'LO' 'A' 'HI' 'B' 'LO' 'D' 'LO' 'D' 'HI' 'D' 'LO' 'B' 'LO' 'A' 'HI' 'B' 'HI' 'D' 'HI' 'D' 'HI' 'D' 'LO' 'C' 'HI' 'D' 'HI' 'B' 'HI' 'B' 'LO' 'D' 'HI' 'D' 'HI' 'D' 'LO' 'B' 'HI' 'B' 'HI' 'B' 'HI' 'A' 'HI' 'B' 'HI' 'B' 'HI' 'C' 'LO' 'D' 'LO' 'A' 'HI' 'C' 'LO' 'B' 'LO' 'C' 'LO' 'B' 'HI' 'A' 'HI' 'B' 'HI' 'A' 'LO' 'B' 'HI' 'B' 'LO' 'B' 'HI' 'B' 'HI' 'C' 'LO' 'D' 'LO' 'A' 'HI' 'C' 'LO' 'B' 'HI' 'B' 'LO' 'C' 'LO' 'A' 'HI' 'D' 'HI' 'B' 'LO' 'A' 'HI' 'B' 'HI' 'A' 'HI' 'D' 'HI' 'A' 'HI' 'D' 'LO' 'B' 'HI' 'D' 'LO' 'A' 'HI' 'B' 'HI' 'D' 'HI' 'D' 'HI' 'B' 'HI' 'D' 'LO' 'C' 'HI' 'B' 'LO' 'A' 'HI' 'A' 'LO' 'B' 'HI' 'B' 'HI' 'B' 'HI' 'B' 'HI' 'D' 'HI' 'B' 'LO' 'B' 'HI' 'B' 'HI' 'B' 'LO' 'C' 'LO' 'A' 'LO' 'A' 'HI' 'A' 'HI' 'B' 'LO' 'C' 'LO' 'D' 'LO' 'B' 'HI' 'B' 'LO' 'B' 'HI' 'A' 'LO' 'D' 'LO' 'B' 'HI' 'D' 'HI' 'B' 'HI' 'C' 'HI' 'B' 'HI' 'A' 'HI' 'B' 'HI' 'B' 'LO' 'A' 'LO' 'A' 'LO' 'D' 'HI' 'D' 'LO' 'D' 'LO' 'A' 'HI' 'C' 'HI' 'A' 'HI' 'D' 'HI' 'D' 'HI' 'C' 'LO' 'C' 'LO' 'B' 'HI' 'A' 'HI' 'C' 'LO' 'D' 'HI' 'C' 'HI'
Names of X columns:
Happiness Separate
Factor 1


Factor 2


Type of test to use (?)


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

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

Tabulation of Results
Happiness x Separate
	HI	LO
A	25	18
B	40	20
C	5	14
D	20	20

Tabulation of Expected Results
Happiness x Separate
	HI	LO
A	23.89	19.11
B	33.33	26.67
C	10.56	8.44
D	22.22	17.78

Statistical Results
Pearson's Chi-squared test
Chi Square Statistic	10.2
Degrees of Freedom	3
P value	0.02