FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_Bayesian_Two_Sample_Test.wasp
Title produced by softwareBayesian Two Sample Test
Date of computationTue, 20 Dec 2016 14:27:11 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2016/Dec/20/t1482240487uezsbuhx4q4axnt.htm/, Retrieved Tue, 06 May 2025 07:20:53 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=301655, Retrieved Tue, 06 May 2025 07:20:53 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact138

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Bayesian Two Sample Test] [blauw en geel] [2016-12-20 13:27:11] [8263efc94e08b372ab727a2b95bd56b1] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
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

Tables (Output of Computation)





Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time46 seconds
R ServerBig Analytics Cloud Computing Center

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time46 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=301655&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]46 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=301655&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=301655&T=0

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time46 seconds
R ServerBig Analytics Cloud Computing Center






Bayesian Two Sample Test
> summary(r, credMass = par4, compValeff = par3)
             mean  median     mode HDI%  HDIlo  HDIup compVal %>compVal
mu1        4.1690  4.1719  4.17368   95  2.796  5.576                  
mu2        4.6094  4.6124  4.67881   95  3.180  6.022                  
muDiff    -0.4404 -0.4403 -0.44331   95 -2.432  1.553       0      32.7
sigma1     2.5758  2.4889  2.30626   95  1.613  3.750                  
sigma2     2.6194  2.5273  2.39963   95  1.649  3.812                  
sigmaDiff -0.0436 -0.0376 -0.00611   95 -1.673  1.610       0      47.9
nu        38.6559 30.0926 13.77224   95  2.165 99.910                  
log10nu    1.4576  1.4785  1.56464   95  0.760  2.099                  
effSz     -0.1725 -0.1714 -0.14369   95 -0.927  0.573       0      32.7


\begin{tabular}{lllllllll}
\hline
Bayesian Two Sample Test \tabularnewline

> summary(r, credMass = par4, compValeff = par3)
             mean  median     mode HDI%  HDIlo  HDIup compVal %>compVal
mu1        4.1690  4.1719  4.17368   95  2.796  5.576                  
mu2        4.6094  4.6124  4.67881   95  3.180  6.022                  
muDiff    -0.4404 -0.4403 -0.44331   95 -2.432  1.553       0      32.7
sigma1     2.5758  2.4889  2.30626   95  1.613  3.750                  
sigma2     2.6194  2.5273  2.39963   95  1.649  3.812                  
sigmaDiff -0.0436 -0.0376 -0.00611   95 -1.673  1.610       0      47.9
nu        38.6559 30.0926 13.77224   95  2.165 99.910                  
log10nu    1.4576  1.4785  1.56464   95  0.760  2.099                  
effSz     -0.1725 -0.1714 -0.14369   95 -0.927  0.573       0      32.7

 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=301655&T=1

[TABLE]
[ROW][C]Bayesian Two Sample Test[/C][/ROW]
[ROW][C]
> summary(r, credMass = par4, compValeff = par3)
             mean  median     mode HDI%  HDIlo  HDIup compVal %>compVal
mu1        4.1690  4.1719  4.17368   95  2.796  5.576                  
mu2        4.6094  4.6124  4.67881   95  3.180  6.022                  
muDiff    -0.4404 -0.4403 -0.44331   95 -2.432  1.553       0      32.7
sigma1     2.5758  2.4889  2.30626   95  1.613  3.750                  
sigma2     2.6194  2.5273  2.39963   95  1.649  3.812                  
sigmaDiff -0.0436 -0.0376 -0.00611   95 -1.673  1.610       0      47.9
nu        38.6559 30.0926 13.77224   95  2.165 99.910                  
log10nu    1.4576  1.4785  1.56464   95  0.760  2.099                  
effSz     -0.1725 -0.1714 -0.14369   95 -0.927  0.573       0      32.7

[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=301655&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=301655&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Bayesian Two Sample Test
> summary(r, credMass = par4, compValeff = par3)
             mean  median     mode HDI%  HDIlo  HDIup compVal %>compVal
mu1        4.1690  4.1719  4.17368   95  2.796  5.576                  
mu2        4.6094  4.6124  4.67881   95  3.180  6.022                  
muDiff    -0.4404 -0.4403 -0.44331   95 -2.432  1.553       0      32.7
sigma1     2.5758  2.4889  2.30626   95  1.613  3.750                  
sigma2     2.6194  2.5273  2.39963   95  1.649  3.812                  
sigmaDiff -0.0436 -0.0376 -0.00611   95 -1.673  1.610       0      47.9
nu        38.6559 30.0926 13.77224   95  2.165 99.910                  
log10nu    1.4576  1.4785  1.56464   95  0.760  2.099                  
effSz     -0.1725 -0.1714 -0.14369   95 -0.927  0.573       0      32.7


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
Parameters (R input):
par1 = 2 ; par2 = 1 ; par3 = 0 ; par4 = 0.95 ; par5 = unpaired ; par6 = no ; par7 = ; par8 = ; par9 = ; par10 = ; par11 = ; par12 = ; par13 = ; par14 = ; par15 = ; par16 = ;
R code (references can be found in the software module):
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('
',RC.texteval('summary(r, credMass=par4, compValeff=par3)'),'
',sep=''))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')