Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
26 12.9 51 12.2 57 12.8 37 7.4 67 6.7 43 12.6 52 14.8 52 13.3 43 11.1 84 8.2 67 11.4 49 6.4 70 10.6 52 12 58 6.3 68 11.3 62 11.9 43 9.3 56 9.6 56 10 74 6.4 65 13.8 63 10.8 58 13.8 57 11.7 63 10.9 53 16.1 57 13.4 51 9.9 64 11.5 53 8.3 29 11.7 54 9 51 9.7 58 10.8 43 10.3 51 10.4 53 12.7 54 9.3 56 11.8 61 5.9 47 11.4 39 13 48 10.8 50 12.3 35 11.3 30 11.8 68 7.9 49 12.7 61 12.3 67 11.6 47 6.7 56 10.9 50 12.1 43 13.3 67 10.1 62 5.7 57 14.3 41 8 54 13.3 45 9.3 48 12.5 61 7.6 56 15.9 41 9.2 43 9.1 53 11.1 44 13 66 14.5 58 12.2 46 12.3 37 11.4 51 8.8 51 14.6 56 12.6 45 13 37 12.6 59 13.2 42 9.9 38 7.7 66 10.5 34 13.4 53 10.9 49 4.3 55 10.3 49 11.8 59 11.2 40 11.4 58 8.6 60 13.2 63 12.6 56 5.6 54 9.9 52 8.8 34 7.7 69 9 32 7.3 48 11.4 67 13.6 58 7.9 57 10.7 42 10.3 64 8.3 58 9.6 66 14.2 26 8.5 61 13.5 52 4.9 51 6.4 55 9.6 50 11.6 60 11.1 56 4.35 63 12.7 61 18.1 52 17.85 16 16.6 46 12.6 56 17.1 52 19.1 55 16.1 50 13.35 59 18.4 60 14.7 52 10.6 44 12.6 67 16.2 52 13.6 55 18.9 37 14.1 54 14.5 72 16.15 51 14.75 48 14.8 60 12.45 50 12.65 63 17.35 33 8.6 67 18.4 46 16.1 54 11.6 59 17.75 61 15.25 33 17.65 47 16.35 69 17.65 52 13.6 55 14.35 55 14.75 41 18.25 73 9.9 51 16 52 18.25 50 16.85 51 14.6 60 13.85 56 18.95 56 15.6 29 14.85 66 11.75 66 18.45 73 15.9 55 17.1 64 16.1 40 19.9 46 10.95 58 18.45 43 15.1 61 15 51 11.35 50 15.95 52 18.1 54 14.6 66 15.4 61 15.4 80 17.6 51 13.35 56 19.1 56 15.35 56 7.6 53 13.4 47 13.9 25 19.1 47 15.25 46 12.9 50 16.1 39 17.35 51 13.15 58 12.15 35 12.6 58 10.35 60 15.4 62 9.6 63 18.2 53 13.6 46 14.85 67 14.75 59 14.1 64 14.9 38 16.25 50 19.25 48 13.6 48 13.6 47 15.65 66 12.75 47 14.6 63 9.85 58 12.65 44 19.2 51 16.6 43 11.2 55 15.25 38 11.9 56 13.2 45 16.35 50 12.4 54 15.85 57 18.15 60 11.15 55 15.65 56 17.75 49 7.65 37 12.35 43 15.6 59 19.3 46 15.2 51 17.1 58 15.6 64 18.4 53 19.05 48 18.55 51 19.1 47 13.1 59 12.85 62 9.5 62 4.5 51 11.85 64 13.6 52 11.7 67 12.4 50 13.35 54 11.4 58 14.9 56 19.9 63 11.2 31 14.6 65 17.6 71 14.05 50 16.1 57 13.35 47 11.85 54 11.95 47 14.75 57 15.15 43 13.2 41 16.85 63 7.85 63 7.7 56 12.6 51 7.85 50 10.95 22 12.35 41 9.95 59 14.9 56 16.65 66 13.4 53 13.95 42 15.7 52 16.85 54 10.95 44 15.35 62 12.2 53 15.1 50 17.75 36 15.2 76 14.6 66 16.65 62 8.1
Names of X columns:
AMS.I TOT
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