There’s also a multivariate normal approximation: prior_from_idata — pymc_experimental 0.0.15 documentation