# NOT RUN {
data(journals)
# Prior specification
# We assign a uniform Beta distribution on each binomial probability
a <- rep(1, 8)
b <- rep(1, 8)
x <- journals$errors
n <- journals$nr_NHST
factor_levels <- levels(journals$journal)
# restricted hypothesis
Hr1 <- c('JAP , PS , JCCP , PLOS , DP , FP , JEPG < JPSP')
out <- binom_bf_informed(x=x, n=n, Hr=Hr1, a=a, b=b,
factor_levels=factor_levels, niter = 2e3)
summary(out)
# }
Run the code above in your browser using DataLab