if (FALSE) {
# fit a model with non-uniform priors
fit <- brm(rating ~ treat + period + carry,
data = inhaler, family = sratio(),
prior = set_prior("normal(0, 0.5)"))
summary(fit)
# The following code requires the 'priorsense' package to be installed:
library(priorsense)
# perform power-scaling of the prior
powerscale(fit, alpha = 1.5, component = "prior")
# perform power-scaling sensitivity checks
powerscale_sensitivity(fit)
# create power-scaling sensitivity plots (for one variable)
powerscale_plot_dens(fit, variable = "b_treat")
}
Run the code above in your browser using DataLab