Can you provide some examples of multivariate gaussian process models in pymc3

When the input variable is more than two, the function y_ = gp.marginal_likelihood(“y”, X=t[:,None], y=y, noise=σ)often reports error, I want to find an example for reference. Thanks