没有用于优化的掉头采样器

Question

我希望解决以下性质的优化问题：

                     theta' ~ argmax {p(theta|data)}

这可以通过马尔可夫链蒙特卡罗采样和模拟退火来完成。

是否可以通过HMC变体NUTS方法和模拟退火使用PyMC3库优化参数集theta？

到目前为止，在PyMC3库的教程中，我只看到了使用NUTS进行采样的示例。这对于计算难以处理的集成很有用，但我不知道如何用它来解决优化问题。