Displaying 1 result from an estimated 1 matches for "mcmcout2".
Did you mean:
mcmcout
2017 Jul 06
0
Bayes Factor
...#
tapply(competition$biomass, competition$clipping, mean)
tapply(competition$biomass, competition$clipping, sd)
# Bayesian Procedure for ANOVA
# Calculate Bayes Factors
library(BayesFactor)
competitionBayesOut <- anovaBF(biomass ~ clipping, data = competition)
# Run mcmc iterations
mcmcOut2 <- posterior(competitionBayesOut, iterations = 10000)
# 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)...