Hi @Martin_Ingram,
I completely agree and would do that if it was computationally doable. The problem is that I am using GPs as priors in a hierarchical model and running it with datasets with 500+ time series with 300+ points each.
Running NUTs with a simple model like this and just one time series takes several hours because of the GP complexity.
Any suggestion on how to improve performance and ensure that the posterior space is free of these pathological shapes?