Problem with Gaussian processes in v5

It seems like you’re running into API differences between v3 and v5. The main one is that sampling functions return arviz InferenceData objects by default now. Here’s a thread discussing some of these changes that might be helpful.