Some useful functions about bayes object parsing

stan_tidy_tools.R

get_bimodality = function(fit) {
fit %>%
    tidybayes::gather_draws(exposure_rate[S]) %>% summarise(
        bimodal = diptest::dip.test(`.value`) %$% `p.value`
       ) 
}