Wessa.net

Send output to:

Data X:
15 21 21 26 30 22 20 22 14 18 18 23 19 12 25 20 23 22 17 21 21 19 21 22 8 15 29 20 20 19 19 18 22 20 23 21 24 15 12 16 22 23 12 21 22 18 20 25 10 9 23 23 17 16 22 16 24 19 18 25 21 25 20 18 20 23 22 21 19 10 20 22 26 26 23 23 24 23 21 24 21 24 19 23 8 15 17 16 20 19 11 18 8 27 15 13 18 28 18 23 19 21 19 19 30 19 17 18 24 19 20 17 25 25 20 19 27 26 18 14 28 28 21 16 27 20 22 24 28 23 25 22 21 21 22 25 28 27 20 23 29 23 20 18 20 18 23 23 18 19 18 15 19 20 25 16 25 25 25 25 24 19 19 19 26 16 10 19 17 19 13 23 17 21 30 22 4 19 16 20 21 3 22 23 20 14 22 23 23 20 16 15 0 13 18 16 25 7 18 24 18 17 24 24 29 24 15 19 22 28 23 23 24 19 22 23 15 25 17 25 20 20 27 16 26 20 23 25 23 25 15 23 26 17 22 20 18 16 15 23 22 12 27 24 10 11 20 14 17 23 23 18 19 29 13 16 27 19 23 16 16 23 25 19 2 4 26 20 20 20 22 4 24 24 18 16 21 3 24 24 19 23 24 17 19 20 17 22 20 19 21 24 21 19 21 27 16 22 27 23 15 22 21 17 18 23 22 23 20 28 17 29 21 21 23 24 18 20 22 7 24 19 27 28 20 26 27 19 20 20 20 23 21 24 26 16 25 19 18 24 21 21 16 16 25 16 20 21 27 28 20 16 18 23 26 26 18 29 16 18 18 19 21 19 18 16 25 16 20 16 23 18 22 22 10 14 18 20 25 15 23 22 22 16 23 15 19 11 14 15 26 20 17 21 15 16 21 17 20 15 22 16 20 18 26 25 26 20 20 24 24 28 20 22 15 20 25 27 20 17 27 22 20 23 17 22 22 13 24 19 22 15 16 20 22 24 23 18 19 19 20 15 15 20 22 13 24 23 22 24 18 23 28 19 20 22 25 26 24 27 16 15 25 27 23 21 14 24 16 22 13 17 23 22 23 26 14 24 21 16 11 19 16 19 16 11 23 27 23 25 24 22 26 19 19 19 20 16 26 24 17 22 15 20 17 24 15 20 20 17 11 21 28 14 13 12 21 13 19 23 27 25 22 27 16 20 18 19 17 10 11 16 13 14 12 15 19 15 14 14 10 13 21 11 14 20 7 24 16 22 25 5 19 23 13 10 12 21 22 20 17 20 13 9 22 15 12 25 14 14 17 9 10 15 15 15 14 21 13 20 16 12 20 26 18 21 23 13 22 14 23 16 14 22 19 23 16 20 8 16 11 16 10 17 16 17 10 15 13 19 14 18 25 10 22 15 18 22 18 15 20 18 6 17 12 12 19 23 26 28 19 16 3 11 15 22 12 21 25 12 14 24 12 13 15 17 12 28 25 14 21 18 23 16 15 5 19 22 19 12 22 18 24 19 4 20 24 26 22 19 9 22 18 16 19 20 21 17 9 26 28 13 16 22 18 21 10 15 15 13 10 23 21 14 17 15 15 17 26 12 14 26 18 17 20 16 19 12 20 19 25 19 15 12 32
Names of X columns:
Bachelor Schakel
Column number of first sample


Column number of second sample


Confidence


Alternative


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

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
R Engine error message	Error in list(15, 21, 21, 26, 30, 22, 20, 22, 14, 18, 18, 23, 19, 12, : argument 455 is empty Calls: array Execution halted