def METHOD_NAME(data: Tuple[AnyNDArray, AnyNDArray]) -> GPR:
    kernel = gpflow.kernels.Matern52(lengthscales=0.3)
    meanf = gpflow.mean_functions.Linear(1.0, 0.0)
    model = GPR(data, kernel, meanf, noise_variance=0.01)
    for p in model.parameters:
        p.prior = Gamma(to_default_float(1.0), to_default_float(1.0))
    return model