Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
NA 21 NA 26 22 NA NA 22 18 NA 23 NA 12 NA NA 20 22 NA 21 NA 19 NA 22 NA 15 NA 20 NA NA 19 NA 18 NA 15 20 NA NA 21 21 NA NA 15 16 NA 23 NA NA 21 18 NA 25 NA 9 NA 30 NA NA 20 23 NA NA 16 NA 16 NA 19 25 NA 25 NA 18 NA 23 NA 21 NA NA 10 14 NA 22 NA NA 26 23 NA 23 NA 24 NA 24 NA 18 NA NA 23 15 NA 19 NA NA 16 25 NA 23 NA 17 NA 19 NA 21 NA 18 NA 27 NA NA 21 13 NA NA 8 29 NA 28 NA NA 23 NA 21 19 NA NA 19 20 NA NA 18 19 NA 17 NA NA 19 NA 25 NA 19 NA 22 23 NA 26 NA NA 14 28 NA NA 16 24 NA NA 20 NA 12 24 NA NA 22 NA 12 NA 22 20 NA NA 10 23 NA 17 NA NA 22 NA 24 NA 18 21 NA 20 NA 20 NA NA 22 19 NA NA 20 26 NA 23 NA 24 NA 21 NA 21 NA NA 19 8 NA 17 NA 20 NA NA 11 NA 8 NA 15 NA 18 NA 18 NA 19 19 NA 23 NA 22 NA 21 NA 25 NA NA 30 17 NA 27 NA NA 23 23 NA NA 18 NA 18 23 NA 19 NA 15 NA 20 NA 16 NA 24 NA 25 NA 25 NA NA 19 19 NA 16 NA 19 NA 19 NA 23 NA 21 NA NA 22 19 NA 20 NA 20 NA 3 NA 23 NA NA 14 NA 23 NA 20 15 NA NA 13 NA 16 NA 7 24 NA NA 17 24 NA 24 NA NA 19 25 NA 20 NA 28 NA NA 23 NA 27 NA 18 NA 28 21 NA NA 19 23 NA NA 27 22 NA NA 28 25 NA NA 21 NA 22 28 NA NA 20 29 NA 25 NA 25 NA 20 NA 20 NA NA 16 20 NA NA 20 NA 23 NA 18 25 NA NA 18 19 NA NA 25 NA 25 NA 25 NA 24 19 NA 26 NA 10 NA 17 NA NA 13 NA 17 30 NA NA 25 NA 4 NA 16 NA 21 23 NA 22 NA NA 17 NA 20 20 NA NA 22 16 NA 23 NA 16 NA NA 0 18 NA 25 NA 23 NA NA 12 NA 18 NA 24 11 NA 18 NA NA 14 23 NA 24 NA NA 29 NA 18 NA 15 29 NA 16 NA NA 19 NA 22 NA 16 23 NA 23 NA NA 19 NA 4 NA 20 24 NA 20 NA 4 NA 24 NA NA 22 16 NA 3 NA 15 NA NA 24 NA 17 20 NA NA 27 23 NA 26 NA 23 NA NA 17 20 NA NA 22 19 NA 24 NA NA 19 23 NA NA 15 27 NA NA 26 22 NA NA 22 NA 18 15 NA 22 NA NA 27 10 NA 20 NA NA 17 23 NA NA 19 NA 13 27 NA 23 NA NA 16 25 NA NA 2 NA 26 20 NA NA 23 NA 22 24 NA 22 NA 17 NA NA 23 23 NA 28 NA 29 NA 21 NA NA 24 20 NA NA 7 NA 19 28 NA NA 18 26 NA NA 21 NA 19 20 NA 23 NA 24 NA NA 16 NA 19 NA 24 21 NA NA 16 16 NA 21 NA 28 NA NA 16 23 NA NA 26 29 NA NA 18 NA 19 NA 19 NA 16 NA 16 NA 16 18 NA 22 NA 14 NA NA 20 NA 15 NA 22 NA 24 NA 16 19 NA 24 NA 19 NA 15 NA NA 11 15 NA NA 17 20 NA 21 NA NA 16 17 NA NA 20 NA 15 NA 21 NA 16 NA 18 NA 25 21 NA NA 21 NA 16 20 NA 24 NA 28 NA 27 NA NA 22 20 NA 27 NA NA 17 NA 22 NA 23 NA 15 22 NA NA 13 NA 21 NA 18 NA 22 NA 19 NA 15 20 NA 17 NA 21 NA NA 23 NA 20 18 NA NA 22 24 NA 24 NA 18 NA 27 NA 19 NA NA 20 NA 15 NA 20 NA 27 NA 20 20 NA NA 13 NA 21 23 NA NA 26 NA 24 25 NA NA 18 21 NA 23 NA NA 16 19 NA NA 20 25 NA NA 22 20 NA 25 NA 27 NA NA 20 18 NA 26 NA NA 26 24 NA 27 NA 16 NA 15 NA NA 25 27 NA NA 18 NA 16 18 NA NA 23 21 NA 21 NA NA 14 NA 24 18 NA 16 NA 25 NA 22 NA NA 13 20 NA 17 NA 23 NA 22 NA NA 23 NA 22 23 NA NA 10 18 NA 25 NA NA 26 14 NA NA 23 22 NA NA 23 NA 19 14 NA 26 NA 24 NA 21 NA NA 17 NA 16 15 NA NA 11 NA 19 21 NA 20 NA 16 NA 19 NA 16 NA 11 NA 22 NA 20 NA NA 26 26 NA NA 20 NA 24 20 NA 15 NA 23 NA 25 NA 27 NA 23 NA 20 NA NA 25 24 NA 22 NA 27 NA NA 20 17 NA 22 NA 26 NA 19 NA NA 19 24 NA 22 NA NA 16 NA 22 23 NA 19 NA 20 NA 16 NA NA 19 20 NA NA 15 22 NA 22 NA NA 12 NA 15 21 NA 26 NA 27 NA 23 NA 21 NA NA 22 26 NA 24 NA 27 NA NA 18 NA 18 25 NA NA 12 19 NA NA 24 NA 17 NA 22 NA 15 NA 20 24 NA 17 NA NA 24 NA 15 NA 20 NA 20 NA 17 NA 11 NA 21 NA 28 NA 14 NA 13 12 NA NA 21 NA 13 NA 19 23 NA NA 27 NA 25 22 NA NA 27 NA 16 NA 20 NA 18 NA 19 NA 17 NA 10 NA 11 16 NA NA 13 NA 14 NA 12 NA 15 19 NA NA 15 NA 14 NA 14 NA 10 NA 13 21 NA NA 11 NA 14 NA 20 NA 7 22 NA NA 24 NA 16 NA 22 NA 25 NA 5 NA 19 23 NA 13 NA NA 10 NA 12 21 NA NA 22 NA 20 NA 17 20 NA NA 13 NA 9 22 NA 15 NA NA 12 25 NA NA 14 NA 14 NA 17 NA 9 10 NA 15 NA 15 NA 15 NA NA 14 NA 21 NA 13 18 NA NA 20 NA 16 NA 28 NA 12 NA 20 26 NA NA 18 NA 21 NA 23 NA 13 NA 22 NA 14 NA 23 NA 16 NA 14 NA 22 NA 19 NA 23 NA 16 NA 20 NA 8 16 NA NA 11 NA 16 NA 10 NA 17 NA 16 17 NA NA 10 NA 15 NA 13 NA 19 NA 14 NA 18 NA 25 NA 10 NA 22 NA 15 NA 18 NA 22 18 NA NA 15 NA 20 NA 18 6 NA NA 17 NA 12 NA 12 19 NA NA 23 NA 26 NA 28 NA 19 NA 16 NA 3 NA 11 NA 15 NA 22 NA 12 NA 21 NA 25 12 NA NA 14 NA 24 NA 12 NA 13 NA 15 NA 17 NA 12 28 NA 25 NA NA 14 21 NA 18 NA NA 23 NA 16 NA 15 NA 5 NA 19 NA 22 NA 19 NA 12 22 NA NA 18 NA 24 19 NA NA 4 NA 20 NA 24 NA 26 NA 22 19 NA NA 9 NA 22 NA 18 NA 16 NA 19 NA 20 NA 21 NA 17 NA 9 26 NA NA 28 NA 13 NA 16 NA 22 18 NA 21 NA NA 10 15 NA 15 NA NA 13 NA 10 NA 23 NA 21 NA 14 NA 17 NA 15 NA 15 NA 17 NA 26 NA 12 NA 14 26 NA NA 18 NA 17 NA 20 NA 16 NA 19 NA 12 NA 20 NA 19 NA 25 NA 19 NA 15
Names of X columns:
numeracy_mannen numeracy_vrouwen
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