I think you want something similar to this discussion? If you want to sample from a fixed distribution I think you have to use a custom step function. You could use pm.Interpolated to convert your posteriors to distributions, then use a custom sampler like this one to draw from them without updating. Is that correct @ricardoV94 ?