if (FALSE) {
library(brms)
if ("log_prior_draws.brmsfit" %in% methods(log_prior_draws) &&
("log_lik_draws.brmsfit" %in% methods(log_lik_draws))) {
fit <- brm(
yield ~ N * P * K,
data = npk,
prior = prior(normal(0, 1), class = "b"),
refresh = 0
)
powerscale_sensitivity(
fit,
variable = "_pred",
prediction = function(x) predictions_as_draws(
x, brms::posterior_epred
)
)
}
}
Run the code above in your browser using DataLab