You can use nutpie for posterior sampling, and then use the trace with the standard PyMC model for sampe_posterior_predictive, but perhaps nutpie is not returning all the info that sample_posterior_predictive needs (constant/observed data?)
Or you’re not updating the data on the regular pymc model?