I have no particular reason, I am using Metropolis in other parts of the code, where it is significantly faster than NUTS. Are you referring to NUTS as the “default sampler”? I have read some problems with the initialization - needed to use ADVI, which is allegedly at its best with large (and multidimensional) data (according to this post ).