# NOT RUN {
data(lifestresses)
# Prior specification
# We assign a uniform Dirichlet distribution, that is, we set all
# concentration parameters to 1
a <- rep(1, 18)
x <- lifestresses$stress.freq
factor_levels <- lifestresses$month
# Test the following restricted Hypothesis:
# Hr: month1 > month2 > ... > month18
Hr <- paste0(1:18, collapse=">"); Hr
out <- mult_bf_informed(x=x, Hr=Hr, a=a, factor_levels=factor_levels,
niter=1e3, bf_type = 'BFre', seed = 4)
m1 <- summary(out)
m1
# }
Run the code above in your browser using DataLab