Wessa.net

Send output to:

Data X:
68.95 70.18 70.84 69.94 70.26 69.61 70.93 71.39 70.71 72.54 70.73 70.48 70.20 70.48 69.76 71.68 71.46 70.59 70.31 71.61 69.24 70.56 70.94 72.74 68.98 72.63 71.16 71.58 69.96 70.14 69.66 71.22 70.20 70.72 69.99 71.37 71.97 69.73 68.29 71.29 67.67 73.19 70.07 70.64 70.89 69.30 69.66 71.76 71.12 70.28 70.60 69.77 70.09 70.37 70.55 71.70 71.63 70.85 69.30 71.07
Names of X columns:
Drug Placebo
Alternative


Column number for the first variable


Column number for the second variable


Test to Perform


Are observations paired?


Value of Difference of Means to Test


Confidence Level to test (?)


Display a Boxplot


Chart options


Title:
Label y-axis:
Label x-axis:

R Code
par2 <- as.numeric(par2) par3 <- as.numeric(par3) par4 <- as.character(par4) par5 <- as.character(par5) par6 <- as.numeric(par6) par7 <- as.numeric(par7) par8 <- as.logical(par8) if ( par5 == 'unpaired') paired <- FALSE else paired <- TRUE x <- t(y) if(par8){ bitmap(file='test1.png') (r<-boxplot(x ,xlab=xlab,ylab=ylab,main=main,notch=FALSE,col=2)) dev.off() } load(file='createtable') if( par4 == 'Wilcoxon-Mann_Whitney'){ a<-table.start() a <- table.row.start(a) a <- table.element(a,'Wilcoxon Test',3,TRUE) a <- table.row.end(a) a <- table.row.start(a) a <- table.element(a,'',1,TRUE) a <- table.element(a,'Statistic',1,TRUE) a <- table.element(a,'P-value',1,TRUE) a <- table.row.end(a) W <- wilcox.test(x[,par2],x[,par3],alternative=par1, paired = paired) a<-table.row.start(a) a<-table.element(a,'Wilcoxon Test',1,TRUE) a<-table.element(a,W$statistic[[1]]) a<-table.element(a,round(W$p.value, digits=5) ) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab') } if( par4 == 'T-Test') { T <- t.test(x[,par2],x[,par3],alternative=par1, paired=paired, mu=par6, conf.level=par7) a<-table.start() a <- table.row.start(a) a <- table.element(a,'T-Test',3,TRUE) a <- table.row.end(a) if(paired){ a <- table.row.start(a) a <- table.element(a,'Difference: Mean1 - Mean2',1,TRUE) a<-table.element(a,round(T$estimate, digits=5) ) a <- table.row.end(a) } if(!paired){ a <- table.row.start(a) a <- table.element(a,'Mean1',1,TRUE) a<-table.element(a,round(T$estimate[1], digits=5) ) a <- table.row.end(a) a <- table.row.start(a) a <- table.element(a,'Mean2',1,TRUE) a<-table.element(a,round(T$estimate[2], digits=5) ) a <- table.row.end(a) } a <- table.row.start(a) a <- table.element(a,'T Statistic',1,TRUE) a<-table.element(a,round(T$statistic, digits=5) ) a <- table.row.end(a) a <- table.row.start(a) a <- table.element(a,'P-value',1,TRUE) a<-table.element(a,round(T$p.value, digits=5) ) a <- table.row.end(a) a <- table.row.start(a) a <- table.element(a,'Lower Confidence Limit',1,TRUE) a<-table.element(a,round(T$conf.int[1], digits=5) ) a <- table.row.end(a) a<-table.row.start(a) a <- table.element(a,'Upper Confidence Limit',1,TRUE) a<-table.element(a,round(T$conf.int[2], digits=5) ) a <- table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab') } a<-table.start() a <- table.row.start(a) a <- table.element(a,'Standard Deviations',3,TRUE) a <- table.row.end(a) a <- table.row.start(a) a <- table.element(a,'Variable 1',1,TRUE) a<-table.element(a,round(sd(x[,par2]), digits=5) ) a <- table.row.end(a) a <- table.row.start(a) a <- table.element(a,'Variable 2',1,TRUE) a<-table.element(a,round(sd(x[,par3]), digits=5) ) a <- table.row.end(a) a<-table.end(a) table.save(a,file='mytable1.tab')

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	Big Analytics Cloud Computing Center

T-Test
Mean1	70.204
Mean2	70.98033
T Statistic	-3.08443
P-value	0.00312
Lower Confidence Limit	-1.28016
Upper Confidence Limit	-0.27251

Standard Deviations
Variable 1	0.96628
Variable 2	0.98326