I want to plot ,export, and save the result of predictive of my observations，I know how to do this in v3.I use with model_reg: # pred_samples = pm.sample_posterior_predictive(f_pred) pred_samples = pm.sample_posterior_predictive(trace_reg, var_names=['f_pred']) to get the prediction, then I…

You have only partial code, and the plotting code you posted will not make the plot you posted (there is nothing to get the HDI, draw sample trajectories, etc). Here is a simple model and a plot of a conditional distribution. Perhaps it will be helpful? # Generate fake data rng = np.random.default…

What were you doing in v3 to make plots? The arviz API hasn’t changed much (at all?), so my guess is you just need some pointers related to working with xarrays. It would be good to see your starting point.

I want to plot the predicted samples f_pred at X_new points and error bar, in v3,I used X_new = np.linspace(np.floor(x.min()), np.ceil(x.max()), 100)[:,None] with model_reg: f_pred = gp.conditional('f_pred', X_new with model_reg: pred_samples = pm.sample_posterior_predictive(f_pred) _, ax =…

Thanks very much, it’s helpful. Besides, when I use gaussian process, it takes a very long time to finish the process, if there any method can make it faster?

How to plot the predict results in v5

Questions v5

jessegrabowski December 27, 2023, 1:53pm 6

You can try an HSGP approximation, see this thread and the linked PyMCon videos for all the details.

Topic		Replies	Views
Problem with Gaussian processes in v5 v5 gaussian_process , modeling	11	847	January 2, 2024
PyMC plot pm.gp.util.plot_gp_dist version agnostic gaussian_process	1	909	October 11, 2022
Plotting Sample_Posterior_Predictive output with Arviz v5 arviz	6	2019	January 26, 2023
GP Predict Point Questions	9	1099	March 6, 2020
Plot posterior predictive results scatterplot in specific format using arviz version agnostic arviz	3	608	July 11, 2023

How to plot the predict results in v5

Related topics