Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
12.9 50 12.2 68 12.8 62 7.4 54 6.7 71 12.6 54 14.8 65 13.3 73 11.1 52 8.2 84 11.4 42 6.4 66 10.6 65 12 78 6.3 73 11.3 75 11.9 72 9.3 66 9.6 70 10 61 6.4 81 13.8 71 10.8 69 13.8 71 11.7 72 10.9 68 16.1 70 13.4 68 9.9 61 11.5 67 8.3 76 11.7 70 9 60 9.7 77 10.8 72 10.3 69 10.4 71 12.7 62 9.3 70 11.8 64 5.9 58 11.4 76 13 52 10.8 59 12.3 68 11.3 76 11.8 65 7.9 67 12.7 59 12.3 69 11.6 76 6.7 63 10.9 75 12.1 63 13.3 60 10.1 73 5.7 63 14.3 70 8 75 13.3 66 9.3 63 12.5 63 7.6 64 15.9 70 9.2 75 9.1 61 11.1 60 13 62 14.5 73 12.2 61 12.3 66 11.4 64 8.8 59 14.6 64 12.6 60 13 66 12.6 78 13.2 53 9.9 67 7.7 59 10.5 66 13.4 68 10.9 71 4.3 66 10.3 73 11.8 72 11.2 71 11.4 59 8.6 64 13.2 66 12.6 78 5.6 68 9.9 73 8.8 62 7.7 65 9 68 7.3 65 11.4 60 13.6 71 7.9 65 10.7 68 10.3 64 8.3 74 9.6 69 14.2 76 8.5 68 13.5 72 4.9 67 6.4 63 9.6 59 11.6 73 11.1 66 4.35 62 12.7 69 18.1 66 17.85 51 16.6 56 12.6 67 17.1 69 19.1 57 16.1 56 13.35 55 18.4 63 14.7 67 10.6 65 12.6 47 16.2 76 13.6 64 18.9 68 14.1 64 14.5 65 16.15 71 14.75 63 14.8 60 12.45 68 12.65 72 17.35 70 8.6 61 18.4 61 16.1 62 11.6 71 17.75 71 15.25 51 17.65 56 16.35 70 17.65 73 13.6 76 14.35 59 14.75 68 18.25 48 9.9 52 16 59 18.25 60 16.85 59 14.6 57 13.85 79 18.95 60 15.6 60 14.85 59 11.75 62 18.45 59 15.9 61 17.1 71 16.1 57 19.9 66 10.95 63 18.45 69 15.1 58 15 59 11.35 48 15.95 66 18.1 73 14.6 67 15.4 61 15.4 68 17.6 75 13.35 62 19.1 69 15.35 58 7.6 60 13.4 74 13.9 55 19.1 62 15.25 63 12.9 69 16.1 58 17.35 58 13.15 68 12.15 72 12.6 62 10.35 62 15.4 65 9.6 69 18.2 66 13.6 72 14.85 62 14.75 75 14.1 58 14.9 66 16.25 55 19.25 47 13.6 72 13.6 62 15.65 64 12.75 64 14.6 19 9.85 50 12.65 68 19.2 70 16.6 79 11.2 69 15.25 71 11.9 48 13.2 66 16.35 73 12.4 74 15.85 66 18.15 71 11.15 74 15.65 78 17.75 75 7.65 53 12.35 60 15.6 50 19.3 70 15.2 69 17.1 65 15.6 78 18.4 78 19.05 59 18.55 72 19.1 70 13.1 63 12.85 63 9.5 71 4.5 74 11.85 67 13.6 66 11.7 62 12.4 80 13.35 73 11.4 67 14.9 61 19.9 73 11.2 74 14.6 32 17.6 69 14.05 69 16.1 84 13.35 64 11.85 58 11.95 60 14.75 59 15.15 78 13.2 57 16.85 60 7.85 68 7.7 68 12.6 73 7.85 69 10.95 67 12.35 60 9.95 65 14.9 66 16.65 74 13.4 81 13.95 72 15.7 55 16.85 49 10.95 74 15.35 53 12.2 64 15.1 65 17.75 57 15.2 51 14.6 80 16.65 67 8.1 70
Names of X columns:
TOT AMS.E
Column number of first sample
Column number of second sample
Confidence
Alternative
two.sided
two.sided
less
greater
Are observations paired?
unpaired
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation