You could technically use the MAP to draw samples from scipy distributions, but I think what you want is to sample from the posterior predictive distribution. Take a look at this talk:
Talk Abstract
PyMC3 is great for inferring parameter values in a model given some observations, but sometimes we also want to generate random samples from the model as predictions given what we already inferred from the observed data. This kind of sampling is called posterior predictive sampling, and it can be very hard. The typical problems that show up are related to shape mismatches in hierarchical models, latent categorical values that aren’t correctly re-sampled or changing the shape of th…
and/or this notebook: https://docs.pymc.io/pymc-examples/examples/diagnostics_and_criticism/posterior_predictive.html
2 Likes