We use hdi function in our analysis. The present helper function is needed for the Bayesian parameter estimation robustness test.

mode_HDI(scale, density, crit_width = 0.95, n_samples = 1e+05)

Arguments

scale

numeric, scale of the sample

density

numeric, density

crit_width

numeric, the critical width

n_samples

numeric, number of iteration