Hi,
How would we do pp_check by percentiles for brms() ? i.e. for specified test statistics other than min, max, mean, median?
Thank you!
Hi,
How would we do pp_check by percentiles for brms() ? i.e. for specified test statistics other than min, max, mean, median?
Thank you!
I think you should be able to pass the stat
argument via ...
to bayesplot::ppc_stat()
, which is the function called behind the scenes for the plots of statistics. stat
can be any function that takes a vector input and returns a scalar (statistic). For example, stat = function(x) quantile(x, probs=0.9)