I am definitely thinking of a hierarchical structure. Would this be achieved using another GP to model the mean of the one below?
That’s right. And it’s true that there isn’t a ton online. It’s still in progress, but example 3 in this PR has some working code you could try: Add HSGP reference and example nb by bwengals · Pull Request #647 · pymc-devs/pymc-examples · GitHub