Regression using RBF kernel

Thanks for the suggestion. I implemented that, but now I don’t understand how to make the predictions because the number of predictions are always equal to the number of samples used for training.

I used the same rbf function and this is the PyMC3 implementation.

with pm.Model() as model:
    cov = tf.cholesky(rbf(X.get_value().astype(np.float32)))
    sess = tf.Session()
    with sess.as_default():
        cov = cov.eval()
    z = pm.MvNormal("z", mu=0, cov=cov.astype(np.float64), shape=100)
    sigma = pm.HalfCauchy("sigma", beta=10)
 
    likelihood = pm.Normal("y", mu=z, sd=sigma, observed=y)

Can you please help me with figuring out how to make the predictions?