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)