Initialzing NUTS with covariance matrix as mass matrix

[aseyboldt]

And by the way, you could also try https://github.com/aseyboldt/covadapt
I’m working on further improving and testing this right now, but even as it is it should work better than initializing with the empirical covariance.