Hello I am using GuassianProcessess.jl, I have been looking at their tutorial (Simple GP Regression · GaussianProcesses.jl) and implement it. Likewise, I want to implement a kernel with summing its value. Now I would like run a chain to optimize the hyper parameters. But I am running into this error, can anybody help me?
using Random Random.seed!(20140430) # Training data n=10; #number of training points x = 2π * rand(n); #predictors y = sin.(x) + 0.05*randn(n); #Select mean and covariance function #here values are log scaled. mZero = MeanZero() #Zero mean function kern1 = Periodic(log(1.0), log(1.0), log(1.5)) #Sqaured exponential kernel (note that hyperparameters are on the log scale) kern2 = Noise(log(0.01)) kern = kern1+kern2 # log standard deviation of observation noise (this is optional) gp = GP(x,y,mZero,kern, 0) μ, σ² = predict_y(gp,range(0,stop=2π,length=100)); mcmc(gp)