By the way I refrained from using the new rstanarm function sbc() , as it does not provide the automatic thinning (while guaranteeing a fixed number of sample outputs). On the other hand the rstanarm implementation is parallelised, which I think is a huge advantage given the embarrassingly parallel nature of SBC. I guess sooner or later we combine this…