Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
12.2 12.9 7.4 12.8 13.3 6.7 6.3 12.6 11.3 14.8 11.9 11.1 9.6 8.2 6.4 11.4 13.8 6.4 9.9 10.6 8.3 12.0 11.7 9.3 12.7 10.0 5.9 13.8 11.8 10.8 12.3 11.7 14.3 10.9 13.3 16.1 7.6 13.4 15.9 11.5 9.1 9.0 13.0 9.7 12.3 10.8 11.4 10.3 8.8 10.4 14.6 9.3 12.6 11.8 9.9 11.4 7.7 13.0 13.4 10.8 10.9 12.3 4.3 11.3 11.8 7.9 8.6 12.7 13.2 11.6 12.6 6.7 7.7 10.9 7.3 12.1 8.3 13.3 13.5 10.1 4.9 5.7 6.4 8.0 9.6 9.3 11.6 12.5 11.1 9.2 12.6 11.1 16.1 14.5 18.4 12.2 14.7 12.6 14.75 13.0 16.1 13.2 16.35 10.5 17.65 10.3 13.6 11.2 14.35 11.4 14.75 5.6 18.25 9.9 16 8.8 14.6 9.0 14.85 11.4 11.75 13.6 18.45 7.9 15.9 10.7 16.1 10.3 10.95 9.6 15.1 14.2 11.35 8.5 15.95 4.35 14.6 12.7 15.35 18.1 13.4 17.85 13.9 16.6 19.1 17.1 12.9 19.1 17.35 13.35 13.15 10.6 12.15 12.6 12.6 16.2 13.6 13.6 14.85 18.9 14.1 14.1 14.9 14.5 16.25 16.15 19.25 14.8 15.65 12.45 12.75 12.65 9.85 17.35 16.6 8.6 13.2 18.4 16.35 11.6 12.4 17.75 18.15 15.25 17.75 17.65 7.65 9.9 12.35 18.25 15.2 16.85 17.1 13.85 15.6 18.95 18.55 15.6 19.1 17.1 13.1 19.9 13.6 18.45 11.4 15 14.9 18.1 14.05 15.4 13.35 15.4 14.75 17.6 13.2 13.35 7.85 19.1 12.6 7.6 7.85 15.25 9.95 16.1 13.4 10.35 15.7 15.4 16.85 9.6 12.2 18.2 17.75 14.75 15.2 13.6 16.65 13.6 8.1 14.6 12.9 12.65 12.8 19.2 6.7 11.2 12.6 15.25 14.8 11.9 11.1 15.85 8.2 11.15 11.4 15.65 6.4 15.6 10.6 19.3 12.0 18.4 9.3 19.05 10.0 12.85 13.8 9.5 10.8 4.5 11.7 11.85 10.9 11.7 16.1 12.4 13.4 13.35 11.5 19.9 9.0 11.2 9.7 14.6 10.8 17.6 10.3 16.1 10.4 11.85 9.3 11.95 11.8 15.15 11.4 16.85 13.0 7.7 10.8 10.95 12.3 12.35 11.3 14.9 7.9 16.65 12.7 13.95 11.6 10.95 6.7 15.35 10.9 15.1 12.1 14.6 13.3 NA 10.1 NA 5.7 NA 8.0 NA 9.3 NA 12.5 NA 9.2 NA 11.1 NA 14.5 NA 12.2 NA 12.6 NA 13.0 NA 13.2 NA 10.5 NA 10.3 NA 11.2 NA 11.4 NA 5.6 NA 9.9 NA 8.8 NA 9.0 NA 11.4 NA 13.6 NA 7.9 NA 10.7 NA 10.3 NA 9.6 NA 14.2 NA 8.5 NA 4.35 NA 12.7 NA 18.1 NA 17.85 NA 16.6 NA 17.1 NA 19.1 NA 13.35 NA 10.6 NA 12.6 NA 16.2 NA 13.6 NA 18.9 NA 14.1 NA 14.5 NA 16.15 NA 14.8 NA 12.45 NA 12.65 NA 17.35 NA 8.6 NA 18.4 NA 11.6 NA 17.75 NA 15.25 NA 17.65 NA 9.9 NA 18.25 NA 16.85 NA 13.85 NA 18.95 NA 15.6 NA 17.1 NA 19.9 NA 18.45 NA 15 NA 18.1 NA 15.4 NA 15.4 NA 17.6 NA 13.35 NA 19.1 NA 7.6 NA 15.25 NA 16.1 NA 10.35 NA 15.4 NA 9.6 NA 18.2 NA 14.75 NA 13.6 NA 13.6 NA 14.6 NA 12.65 NA 19.2 NA 11.2 NA 15.25 NA 11.9 NA 15.85 NA 11.15 NA 15.65 NA 15.6 NA 19.3 NA 18.4 NA 19.05 NA 12.85 NA 9.5 NA 4.5 NA 11.85 NA 11.7 NA 12.4 NA 13.35 NA 19.9 NA 11.2 NA 14.6 NA 17.6 NA 16.1 NA 11.85 NA 11.95 NA 15.15 NA 16.85 NA 7.7 NA 10.95 NA 12.35 NA 14.9 NA 16.65 NA 13.95 NA 10.95 NA 15.35 NA 15.1 NA 14.6 NA
Names of X columns:
Vrouwen Mannen
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation