gp function

Set up Gaussian process terms in brms