Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
7.5 12.9 6 12.2 6.5 12.8 1 7.4 1 6.7 5.5 12.6 8.5 14.8 6.5 13.3 4.5 11.1 2 8.2 5 11.4 0.5 6.4 5 10.6 5 12 2.5 6.3 5 11.3 5.5 11.9 3.5 9.3 3 9.6 4 10 0.5 6.4 6.5 13.8 4.5 10.8 7.5 13.8 5.5 11.7 4 10.9 7.5 16.1 7 13.4 4 9.9 5.5 11.5 2.5 8.3 5.5 11.7 3.5 9 2.5 9.7 4.5 10.8 4.5 10.3 4.5 10.4 6 12.7 2.5 9.3 5 11.8 0 5.9 5 11.4 6.5 13 5 10.8 6 12.3 4.5 11.3 5.5 11.8 1 7.9 7.5 12.7 6 12.3 5 11.6 1 6.7 5 10.9 6.5 12.1 7 13.3 4.5 10.1 0 5.7 8.5 14.3 3.5 8 7.5 13.3 3.5 9.3 6 12.5 1.5 7.6 9 15.9 3.5 9.2 3.5 9.1 4 11.1 6.5 13 7.5 14.5 6 12.2 5 12.3 5.5 11.4 3.5 8.8 7.5 14.6 6.5 12.6 6.5 13 6.5 12.6 7 13.2 3.5 9.9 1.5 7.7 4 10.5 7.5 13.4 4.5 10.9 0 4.3 3.5 10.3 5.5 11.8 5 11.2 4.5 11.4 2.5 8.6 7.5 13.2 7 12.6 0 5.6 4.5 9.9 3 8.8 1.5 7.7 3.5 9 2.5 7.3 5.5 11.4 8 13.6 1 7.9 5 10.7 4.5 10.3 3 8.3 3 9.6 8 14.2 2.5 8.5 7 13.5 0 4.9 1 6.4 3.5 9.6 5.5 11.6 5.5 11.1 0.5 4.35 7.5 12.7 9 18.1 9.5 17.85 8.5 16.6 7 12.6 8 17.1 10 19.1 7 16.1 8.5 13.35 9 18.4 9.5 14.7 4 10.6 6 12.6 8 16.2 5.5 13.6 9.5 18.9 7.5 14.1 7 14.5 7.5 16.15 8 14.75 7 14.8 7 12.45 6 12.65 10 17.35 2.5 8.6 9 18.4 8 16.1 6 11.6 8.5 17.75 6 15.25 9 17.65 8 16.35 9 17.65 5.5 13.6 7 14.35 5.5 14.75 9 18.25 2 9.9 8.5 16 9 18.25 8.5 16.85 9 14.6 7.5 13.85 10 18.95 9 15.6 7.5 14.85 6 11.75 10.5 18.45 8.5 15.9 8 17.1 10 16.1 10.5 19.9 6.5 10.95 9.5 18.45 8.5 15.1 7.5 15 5 11.35 8 15.95 10 18.1 7 14.6 7.5 15.4 7.5 15.4 9.5 17.6 6 13.35 10 19.1 7 15.35 3 7.6 6 13.4 7 13.9 10 19.1 7 15.25 3.5 12.9 8 16.1 10 17.35 5.5 13.15 6 12.15 6.5 12.6 6.5 10.35 8.5 15.4 4 9.6 9.5 18.2 8 13.6 8.5 14.85 5.5 14.75 7 14.1 9 14.9 8 16.25 10 19.25 8 13.6 6 13.6 8 15.65 5 12.75 9 14.6 4.5 9.85 8.5 12.65 9.5 19.2 8.5 16.6 7.5 11.2 7.5 15.25 5 11.9 7 13.2 8 16.35 5.5 12.4 8.5 15.85 9.5 18.15 7 11.15 8 15.65 8.5 17.75 3.5 7.65 6.5 12.35 6.5 15.6 10.5 19.3 8.5 15.2 8 17.1 10 15.6 10 18.4 9.5 19.05 9 18.55 10 19.1 7.5 13.1 4.5 12.85 4.5 9.5 0.5 4.5 6.5 11.85 4.5 13.6 5.5 11.7 5 12.4 6 13.35 4 11.4 8 14.9 10.5 19.9 6.5 11.2 8 14.6 8.5 17.6 5.5 14.05 7 16.1 5 13.35 3.5 11.85 5 11.95 9 14.75 8.5 15.15 5 13.2 9.5 16.85 3 7.85 1.5 7.7 6 12.6 0.5 7.85 6.5 10.95 7.5 12.35 4.5 9.95 8 14.9 9 16.65 7.5 13.4 8.5 13.95 7 15.7 9.5 16.85 6.5 10.95 9.5 15.35 6 12.2 8 15.1 9.5 17.75 8 15.2 8 14.6 9 16.65 5 8.1
Names of X columns:
Ex TOT
Column number of first sample
Column number of second sample
Confidence
Alternative
two.sided
less
greater
Are observations paired?
unpaired
paired
Null Hypothesis
Chart options
Title:
R Code
par1 <- as.numeric(par1) #column number of first sample par2 <- as.numeric(par2) #column number of second sample par3 <- as.numeric(par3) #confidence (= 1 - alpha) if (par5 == 'unpaired') paired <- FALSE else paired <- TRUE par6 <- as.numeric(par6) #H0 z <- t(y) if (par1 == par2) stop('Please, select two different column numbers') if (par1 < 1) stop('Please, select a column number greater than zero for the first sample') if (par2 < 1) stop('Please, select a column number greater than zero for the second sample') if (par1 > length(z[1,])) stop('The column number for the first sample should be smaller') if (par2 > length(z[1,])) stop('The column number for the second sample should be smaller') if (par3 <= 0) stop('The confidence level should be larger than zero') if (par3 >= 1) stop('The confidence level should be smaller than zero') (r.t <- t.test(z[,par1],z[,par2],var.equal=TRUE,alternative=par4,paired=paired,mu=par6,conf.level=par3)) (v.t <- var.test(z[,par1],z[,par2],conf.level=par3)) (r.w <- t.test(z[,par1],z[,par2],var.equal=FALSE,alternative=par4,paired=paired,mu=par6,conf.level=par3)) (w.t <- wilcox.test(z[,par1],z[,par2],alternative=par4,paired=paired,mu=par6,conf.level=par3)) (ks.t <- ks.test(z[,par1],z[,par2],alternative=par4)) m1 <- mean(z[,par1],na.rm=T) m2 <- mean(z[,par2],na.rm=T) mdiff <- m1 - m2 newsam1 <- z[!is.na(z[,par1]),par1] newsam2 <- z[,par2]+mdiff newsam2 <- newsam2[!is.na(newsam2)] (ks1.t <- ks.test(newsam1,newsam2,alternative=par4)) mydf <- data.frame(cbind(z[,par1],z[,par2])) colnames(mydf) <- c('Variable 1','Variable 2') bitmap(file='test1.png') boxplot(mydf, notch=TRUE, ylab='value',main=main) dev.off() bitmap(file='test2.png') qqnorm(z[,par1],main='Normal QQplot - Variable 1') qqline(z[,par1]) dev.off() bitmap(file='test3.png') qqnorm(z[,par2],main='Normal QQplot - Variable 2') qqline(z[,par2]) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,paste('Two Sample t-test (',par5,')',sep=''),2,TRUE) a<-table.row.end(a) if(!paired){ a<-table.row.start(a) a<-table.element(a,'Mean of Sample 1',header=TRUE) a<-table.element(a,r.t$estimate[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Mean of Sample 2',header=TRUE) a<-table.element(a,r.t$estimate[[2]]) a<-table.row.end(a) } else { a<-table.row.start(a) a<-table.element(a,'Difference: Mean1 - Mean2',header=TRUE) a<-table.element(a,r.t$estimate) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a,'t-stat',header=TRUE) a<-table.element(a,r.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'df',header=TRUE) a<-table.element(a,r.t$parameter[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,r.t$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'H0 value',header=TRUE) a<-table.element(a,r.t$null.value[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Alternative',header=TRUE) a<-table.element(a,r.t$alternative) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI Level',header=TRUE) a<-table.element(a,attr(r.t$conf.int,'conf.level')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI',header=TRUE) a<-table.element(a,paste('[',r.t$conf.int[1],',',r.t$conf.int[2],']',sep='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'F-test to compare two variances',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'F-stat',header=TRUE) a<-table.element(a,v.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'df',header=TRUE) a<-table.element(a,v.t$parameter[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,v.t$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'H0 value',header=TRUE) a<-table.element(a,v.t$null.value[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Alternative',header=TRUE) a<-table.element(a,v.t$alternative) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI Level',header=TRUE) a<-table.element(a,attr(v.t$conf.int,'conf.level')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI',header=TRUE) a<-table.element(a,paste('[',v.t$conf.int[1],',',v.t$conf.int[2],']',sep='')) 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,paste('Welch Two Sample t-test (',par5,')',sep=''),2,TRUE) a<-table.row.end(a) if(!paired){ a<-table.row.start(a) a<-table.element(a,'Mean of Sample 1',header=TRUE) a<-table.element(a,r.w$estimate[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Mean of Sample 2',header=TRUE) a<-table.element(a,r.w$estimate[[2]]) a<-table.row.end(a) } else { a<-table.row.start(a) a<-table.element(a,'Difference: Mean1 - Mean2',header=TRUE) a<-table.element(a,r.w$estimate) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a,'t-stat',header=TRUE) a<-table.element(a,r.w$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'df',header=TRUE) a<-table.element(a,r.w$parameter[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,r.w$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'H0 value',header=TRUE) a<-table.element(a,r.w$null.value[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Alternative',header=TRUE) a<-table.element(a,r.w$alternative) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI Level',header=TRUE) a<-table.element(a,attr(r.w$conf.int,'conf.level')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'CI',header=TRUE) a<-table.element(a,paste('[',r.w$conf.int[1],',',r.w$conf.int[2],']',sep='')) 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,paste('Wicoxon rank sum test with continuity correction (',par5,')',sep=''),2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'W',header=TRUE) a<-table.element(a,w.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,w.t$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'H0 value',header=TRUE) a<-table.element(a,w.t$null.value[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Alternative',header=TRUE) a<-table.element(a,w.t$alternative) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Kolmogorov-Smirnov Test to compare <i>Distributions</i> of two Samples',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'KS Statistic',header=TRUE) a<-table.element(a,ks.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,ks.t$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Kolmogorov-Smirnov Test to compare <i>Distributional Shape</i> of two Samples',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'KS Statistic',header=TRUE) a<-table.element(a,ks1.t$statistic[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,ks1.t$p.value) 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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation