Yes, I meant to write data minus shift. But this gives the can’t resolve inputs error. It works for everything else though (i.e. I can fit the model and get meaningful results) — just sampling from the PP gives the error. Not sure if that’s a bug or expected behavior.
That’s for your further suggestions! I’ll try that.