Pp_check by percentile as test statistic?


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)