Pm.sample_posterior_predictive() from an AR1 process

Sorry, I thought you were doing prior prediction.

How is sampling going and what do your parameter estimates look like? Are you getting divergences? I suspect there may be some sampling issues and that you are either just getting your own observed data back or some other degenerate behavior. But that’s just a guess.