Displaying 1 result from an estimated 1 matches for "bestmcmc".
2017 Jul 06
0
Bayes Factor
...# boxplot of the posteriors for the groups
boxplot(as.matrix(mcmcOut2[,2:6]))
# Show the HDIs
summary(mcmcOut2)
# Calculate the Bayes Factor
competitionBayesOut
# Pairwise "post hoc" tests
library(rjags)
library(BEST)
# competitionare 'r5' vs. 'control'
plot(BESTmcmc(competition[competition$clpping=="r5",2],
competition[competition$clipping=="control",2]))
#
plot(BESTmcmc(competition[competition$clpping=="r10",2],
competition[competition$clipping=="control",2]))
# competitionare 'n10' vs. 'control'
plo...