It will be faster to use predictive = Predictive(guide, num_samples=800, return_sites=...). The pair model, guide is used when you want to make new predictions for y: under the hood, it will collect samples from guide, then using those samples and model to make predictions for y.