Wessa.net

Send output to:

Data X:
1 5 4 5 6 1 1 7 8 2 3 4 8 1 7 2 7 2 4 5 6 3 2 7 5 6 2 9 5 4
Names of X columns:
blauw geel
Column number of first sample


Column number of second sample


Null Hypothesis H0


Confidence


Are observations paired?


Use informative priors?


prior value of mu1 (?)


prior value of mu2 (?)


prior value of sigma 1 of the population mean (?)


prior value of sigma 2 of the population mean (?)


prior value of the mode of sigma 1 (?)


prior value of the mode of sigma 2 (?)


prior value of SD of sigma 1 (?)


prior value of SD of sigma 2 (?)


prior value of normality parameter (?)


prior value of SD of normality parameter (?)


Chart options



R Code
par16 <- '' par15 <- '' par14 <- '' par13 <- '' par12 <- '' par11 <- '' par10 <- '' par9 <- '' par8 <- '' par7 <- '' par6 <- 'no' par5 <- 'unpaired' par4 <- '' par3 <- '0' par2 <- '1' par1 <- '2' library(BEST) par1 <- as.numeric(par1) #column number of first sample par2 <- as.numeric(par2) #column number of second sample par3 <- as.numeric(par3) par4 <- as.numeric(par4) if(par6=='yes') { par7 <- as.numeric(par7) par8 <- as.numeric(par8) par9 <- as.numeric(par9) par10 <- as.numeric(par10) par11 <- as.numeric(par11) par12 <- as.numeric(par12) par13 <- as.numeric(par13) par14 <- as.numeric(par14) par15 <- as.numeric(par15) par16 <- as.numeric(par16) } 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(par6=='no') { if(par5=='unpaired') { (r <- BESTmcmc(z[,par1],z[,par2], parallel=F)) } if(par5=='paired') { (r <- BESTmcmc(z[,par1]-z[,par2], parallel=F)) } } else { yy <- cbind(z[1,],z[2,]) if(par5=='unpaired') { (r <- BESTmcmc(z[,par1],z[,par2], priors=list(muM = c(par7,par8), muSD = c(par9,par10), sigmaMode = c(par11,par12), sigmaSD = c(par13,par14), nuMean = par15, nuSD = par16), parallel=F)) } if(par5=='paired') { (r <- BESTmcmc(z[,par1]-z[,par2], priors=list(muM = c(par7,par8), muSD = c(par9,par10), sigmaMode = c(par11,par12), sigmaSD = c(par13,par14), nuMean = par15, nuSD = par16), parallel=F)) } } bitmap(file='test2.png') plot(r, credMass=par4, compVal=par3) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Bayesian Two Sample Test',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('<pre>',RC.texteval('summary(r, credMass=par4, compValeff=par3)'),'</pre>',sep='')) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.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