search

n-kall / priorsense

priorsense: an R package for prior diagnostics and sensitivity
https://n-kall.github.io/priorsense/
GNU General Public License v3.0
53 stars 5 forks source link

powerscale_sensitivity fails, if the draws object has only one variable #46

Closed avehtari closed 4 months ago

avehtari commented 4 months ago 
data_bern <- data.frame(y = c(1, 1, 1, 0, 1, 1, 1, 0, 1, 0))
fit_bern <- brm(y ~ 1, family = bernoulli(), data = data_bern,
                prior = prior(student_t(7, 0, 1.5), class='Intercept'),
                seed = SEED, refresh = 0)
draws <- as_draws_df(fit_bern)

This works:

powerscale_sensitivity(draws, fit=fit_bern, variable="theta")

This gives an error Error in apply(has_prefix, 1, any) : dim(X) must have a positive length

powerscale_sensitivity(subset_draws(draws, variable="theta"), fit=fit_bern)