Hi, My current models is as follows; with pm.Model() as model: l = pm.Gamma("l", alpha=2, beta=1) offset = pm.Gamma("offset", alpha=2, beta=1) nu = pm.HalfCauchy("nu", beta=1) cov = nu ** 2 * pm.gp.cov.Polynomial(X.shape[1], l, 2, offset) gp = pm.gp.Marginal(cov_func=cov) sigma = p…

So what you’re doing here is drawing full samples from the GP, then calculating the pointwise mean and std: pred_samples = pm.sample_posterior_predictive(map_trace, vars=[f_pred], samples=2000) y_pred_custom, uncer = pred_samples['f_pred'].mean(axis=0), pred_samples['f_pred'].std(axis=0) This is…

Predicting on new data with gp.conditional

Questions

BioGoertz November 10, 2020, 7:57pm 8

That’s good to know! I always forget the power of shared variables…

Topic		Replies	Views
GP Predict Point Questions	9	1033	March 6, 2020
Improving design of conditional / predictive sampling from GPs Development	1	586	May 4, 2018
Problem when using shared variable and linear mean function with GP conditional Questions	10	810	December 16, 2020
Predictions from Latent GP in hierarchical model appear not conditioned on anything Questions	6	832	July 6, 2020
GP conditional prediction when GP was construced inside user defined function v5 modeling	3	201	November 27, 2023

Predicting on new data with gp.conditional

Related topics