Wessa.net

Send output to:

Data X:
40 1 27 5 26 49 35 45 1 36 4 25 45 34 38 1 25 4 17 54 13 28 1 27 3 37 36 35 2 25 3 35 36 28 38 2 44 3 15 53 32 39 1 50 4 27 46 35 37 1 41 4 36 42 36 30 1 48 5 25 41 27 30 2 43 4 30 45 29 30 2 47 2 27 47 27 26 2 41 3 33 42 28 29 1 44 2 29 45 29 31 2 47 5 30 40 28 27 2 40 3 25 45 30 25 2 46 3 23 40 25 39 1 28 3 26 42 15 35 1 56 3 24 45 33 27 2 49 4 35 47 31 40 2 25 4 39 31 37 34 2 41 4 23 46 37 32 2 26 3 32 34 34 34 1 50 5 29 43 32 38 1 47 4 26 45 21 21 1 52 2 21 42 25 33 2 37 5 35 51 32 27 2 41 3 23 44 28 35 1 45 4 21 47 22 33 2 26 4 28 47 25 36 1 NA 3 30 41 26 1 52 4 21 44 34 37 1 46 2 29 51 34 37 1 58 3 28 46 36 37 1 54 5 19 47 36 32 1 29 3 26 46 26 25 2 50 3 33 38 26 31 1 43 2 34 50 34 33 2 30 3 33 48 33 18 2 47 2 40 36 31 42 1 45 3 24 51 33 26 2 48 1 35 35 22 26 2 48 3 35 49 29 32 2 26 4 32 38 24 31 1 46 5 20 47 37 29 2 NA 3 35 36 32 35 2 50 3 35 47 23 44 1 25 4 21 46 29 35 1 47 2 33 43 35 30 2 47 2 40 53 20 32 1 41 3 22 55 28 24 2 45 2 35 39 26 34 2 41 4 20 55 36 27 2 45 5 28 41 26 31 2 40 3 46 33 33 38 1 29 4 18 52 25 41 2 34 5 22 42 29 40 1 45 5 20 56 32 25 2 52 3 25 46 35 19 2 41 4 31 33 24 33 2 48 3 21 51 31 27 2 45 3 23 46 29 45 1 54 2 26 46 27 27 2 25 3 34 50 29 30 2 26 4 31 46 29 42 1 28 4 23 51 27 21 2 50 4 31 48 34 32 2 48 4 26 44 32 31 2 51 3 36 38 31 36 2 53 3 28 42 31 34 1 37 3 34 39 31 11 1 56 2 25 45 16 35 1 43 3 33 31 25 39 1 34 3 46 29 27 32 1 42 3 24 48 32 28 2 32 3 32 38 28 45 2 31 5 33 55 25 18 1 46 3 42 32 25 35 2 30 5 17 51 36 35 2 47 4 36 53 36 36 2 33 4 40 47 36 34 1 25 4 30 45 27 34 1 25 5 19 33 29 38 2 21 4 33 49 32 28 2 36 5 35 46 29 23 2 50 3 23 42 31 37 2 48 3 15 56 34 29 2 48 2 38 35 27 28 1 25 3 37 40 28 30 1 48 4 23 44 32 24 2 49 5 41 46 33 36 1 27 5 34 46 29 40 1 28 3 38 39 32 37 2 43 2 45 35 35 27 2 48 3 27 48 33 25 2 48 4 46 42 27 22 1 25 1 26 39 16 21 2 49 4 44 39 32 28 1 26 3 36 41 26 34 1 51 3 20 52 32 32 2 25 4 44 45 38 23 1 29 3 27 42 24 29 1 29 4 27 44 26 35 1 43 2 41 33 19 31 2 46 3 30 42 37 36 1 44 3 33 46 25 32 1 25 3 37 45 24 35 1 51 2 30 40 23 45 1 42 5 20 48 28 29 2 53 5 44 32 38 41 1 25 4 20 53 28 36 2 49 2 33 39 28 37 1 51 3 31 45 26 25 2 20 3 23 36 21 36 2 44 3 33 38 35 34 2 38 4 33 49 31 33 1 46 5 32 46 34 32 2 42 4 25 43 30 40 1 29 NA 22 37 30 27 2 46 4 16 48 24 24 2 49 2 36 45 27 26 2 51 3 35 32 26 1 38 3 25 46 30 13 1 41 1 27 20 15 22 2 47 3 32 42 28 29 2 44 3 36 45 34 30 2 47 3 51 29 29 24 2 46 3 30 51 26 1 44 4 20 55 31 26 2 28 3 29 50 28 37 2 47 4 26 44 33 36 2 28 4 20 41 32 38 1 41 5 40 40 33 34 2 45 4 29 47 31 35 2 46 4 32 42 37 32 1 46 4 33 40 27 44 2 22 3 32 51 19 40 2 33 3 34 43 27 24 1 41 4 24 45 31 36 2 47 5 25 41 38 20 1 25 3 41 41 22 28 2 42 3 39 37 35 18 2 47 3 21 46 35 23 2 50 3 38 38 30 28 1 55 5 28 39 41 30 1 21 3 37 45 25 30 1 NA 3 26 46 28 43 1 52 3 30 39 45 20 2 49 4 25 21 21 37 2 46 4 38 31 33 24 1 NA 4 31 35 25 33 2 45 3 31 49 29 43 2 52 3 27 40 31 27 1 NA 3 21 45 29 22 2 40 4 26 46 31 28 2 49 4 37 45 31 18 1 38 5 28 34 25 38 1 32 5 29 41 27 23 2 46 4 33 43 26 38 2 32 3 41 45 26 21 2 41 3 19 48 23 25 2 43 3 37 43 27 1 44 4 36 45 24 30 1 47 5 27 45 35 25 2 28 3 33 34 24 17 1 52 1 29 40 32 26 1 27 2 42 40 24 39 2 45 5 27 55 24 27 1 27 4 47 44 38 33 1 25 4 17 44 36 47 1 28 4 34 48 24 37 1 25 3 32 51 18 34 1 52 4 25 49 34 24 1 44 3 27 33 23 25 2 43 3 37 43 35 20 2 47 4 34 44 22 34 2 52 4 27 44 34 22 2 40 2 37 41 28 39 1 42 3 32 45 34 33 1 45 5 26 44 32 35 1 45 2 29 44 24 26 1 50 5 28 40 34 32 1 49 3 19 48 33 22 1 52 2 46 49 33 39 2 48 3 31 46 29 35 2 51 3 42 49 38 21 2 49 4 33 55 24 27 2 31 4 39 51 25 31 2 43 3 27 46 37 20 2 31 3 35 37 33 28 2 28 4 23 43 30 26 2 43 4 32 41 22 36 2 31 3 22 45 28 16 2 51 3 17 39 24 34 2 58 4 35 38 33 30 2 25 5 34 41 37
Names of X columns:
Intrinsieke_waarden geslacht leeftijd opleiding Neuroticisme Extraversie Openheid
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(40, 1, 27, 5, 26, 49, 35, 45, 1, 36, 4, 25, 45, 34, 38, : argument 29 is empty Calls: array Execution halted