Wessa.net

Send output to:

Data X:
26 NA 51 NA 57 NA 37 NA 67 NA 43 NA 52 NA 52 NA 43 NA 84 NA 67 NA 49 NA 70 NA 52 NA 58 NA 68 NA NA 62 43 NA 56 NA NA 56 74 NA 65 NA 63 NA 58 NA 57 NA 63 NA 53 NA NA 57 NA 51 64 NA 53 NA 29 NA 54 NA 51 NA 58 NA 43 NA 51 NA 53 NA 54 NA NA 56 61 NA 47 NA 39 NA 48 NA 50 NA 35 NA NA 30 68 NA 49 NA NA 61 67 NA NA 47 NA 56 NA 50 43 NA NA 67 62 NA 57 NA NA 41 54 NA NA 45 NA 48 61 NA 56 NA 41 NA 43 NA 53 NA NA 44 66 NA 58 NA 46 NA NA 37 51 NA 51 NA NA 56 NA 66 45 NA 37 NA 59 NA 42 NA NA 38 66 NA NA 34 53 NA NA 49 NA 55 NA 49 NA 59 NA 40 NA 58 NA 60 NA 63 NA 56 NA 54 NA 52 NA 34 NA 69 NA 32 NA 48 NA 67 NA 58 NA 57 NA 42 NA 64 NA 58 NA 66 NA 26 NA 61 NA 52 NA 51 NA 55 NA 50 NA 60 NA 56 NA 63 NA 61 52 NA 16 NA 46 NA 56 NA NA 52 NA 55 50 NA 59 NA 60 NA 52 NA 44 NA 67 NA 52 NA 55 NA 37 NA 54 NA NA 72 51 NA 48 NA 60 NA 50 NA 63 NA 33 NA 67 NA 46 NA 54 NA 59 NA 61 NA NA 33 47 NA 69 NA 52 NA 55 NA 55 NA 41 NA 73 NA 51 NA 52 NA 50 NA 51 NA 60 NA 56 NA 56 NA 29 NA NA 66 NA 66 73 NA 55 NA NA 64 NA 40 NA 46 NA 58 43 NA 61 NA NA 51 NA 50 NA 52 NA 54 NA 66 NA 61 NA 80 NA 51 NA 56 56 NA 56 NA NA 53 47 NA 25 NA NA 47 46 NA NA 50 NA 39 51 NA NA 58 NA 35 NA 58 NA 60 NA 62 NA 63 NA 53 NA 46 NA 67 NA 59 NA 64 NA 38 NA 50 48 NA NA 48 NA 47 NA 66 47 NA NA 63 58 NA NA 44 51 NA NA 43 55 NA NA 38 NA 56 NA 45 NA 50 NA 54 57 NA 60 NA NA 55 56 NA 49 NA NA 37 43 NA 59 NA NA 46 NA 51 58 NA NA 64 53 NA 48 NA 51 NA NA 47 59 NA NA 62 62 NA 51 NA 64 NA 52 NA NA 67 50 NA 54 NA 58 NA NA 56 63 NA 31 NA NA 65 71 NA NA 50 NA 57 NA 47 54 NA NA 47 NA 57 43 NA 41 NA 63 NA 63 NA 56 NA 51 NA NA 50 NA 22 41 NA NA 59 NA 56 66 NA NA 53 NA 42 NA 52 NA 54 NA 44 NA 62 NA 53 NA 50 NA 36 NA 76 NA 66 NA 62 NA 59 NA 47 NA 55 NA 58 NA 60 44 NA NA 57 NA 45 58 NA 51 NA 57 NA 30 NA 46 NA 51 NA 56 NA 58 NA 44 NA 14 NA 53 NA 42 NA NA 49 44 NA NA 62 30 NA 46 NA NA 56 50 NA 54 NA 48 NA 55 NA 35 NA 55 NA 41 NA 59 NA 54 NA 66 NA 55 NA 45 NA 51 NA 47 NA 42 NA 53 NA 53 NA 41 NA 55 NA 55 NA 46 NA 63 NA 43 NA 65 NA 59 NA 39 NA 44 NA NA 60 57 NA NA 67 NA 52 NA 52 69 NA 46 NA 46 NA NA 53 40 NA 70 NA 54 NA 77 NA NA 45 60 NA NA 47 50 NA 66 NA 60 NA NA 41 NA 53 NA 34 51 NA 69 NA 60 NA NA 45 58 NA 39 NA 51 NA 52 NA 49 NA 63 NA NA 44 51 NA 52 NA NA 60 NA 53 NA 53 52 NA 31 NA NA 51 NA 65 NA 51 NA 49 61 NA NA 58 NA 62 54 NA NA 52 72 NA NA 50 65 NA NA 53 56 NA 63 NA NA 62 NA 66 NA 50 45 NA NA 58 52 NA NA 53 68 NA NA 59 NA 58 NA 52 45 NA NA 58 70 NA 69 NA NA 71 46 NA NA 58 39 NA NA 46 NA 64 NA 67 NA 44 54 NA 41 NA 68 NA 63 NA 57 NA 61 NA 39 NA NA 69 NA 64 NA 38 59 NA 51 NA NA 59 51 NA 65 NA NA 47 50 NA NA 57 21 NA 47 NA NA 51 37 NA NA 67 NA 43 58 NA 51 NA 40 NA NA 41 NA 58 NA 64 64 NA 58 NA NA 50 NA 59 NA 55 NA 59 NA 58 NA 41 56 NA 63 NA NA 77 60 NA NA 58 64 NA 47 NA 46 NA NA 62 NA 60 50 NA 46 NA 44 NA 58 NA NA 56 NA 43 NA 54 NA 54 NA 56 NA 65 NA 66 NA 62 58 NA NA 67 25 NA 56 NA NA 53 56 NA 59 NA 46 NA NA 49 NA 56 NA 76 NA 33 49 NA 53 NA 58 NA NA 72 NA 51 NA 42 NA 69 NA 51 54 NA 52 NA 59 NA NA 51 NA 67 NA 64 NA 58 NA NA NA NA NA NA NA NA 53 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
Names of X columns:
AMS.I_(schakel) AMS.I_(bachelor)
Column number of first sample


Column number of second sample


Confidence


Alternative


Are observations paired?


Null Hypothesis


Chart options


Title:

R Code
par6 <- '0.0' par5 <- 'unpaired' par4 <- 'two.sided' par3 <- '0.95' par2 <- '2' par1 <- '1' 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