pymc - 从 PyMC3 中用户提供的目标密度采样

Question

是否可以以简单的方式从 PyMC3 中用户提供的目标度量中进行采样？即我希望能够提供黑盒功能logposterior(theta)，grad_logposterior(theta)并从中采样，而不是在 PyMC3s 建模语言中指定模型。

Answer 1

这有点笨拙。您需要创建一个新的 Theano Op。这里有几个例子：https ://github.com/Theano/Theano/blob/master/theano/tensor/slinalg.py#L32

然后，您需要创建一个通过 new 评估 logp 的分发类，Op例如：https ://github.com/pymc-devs/pymc3/blob/master/pymc3/distributions/continuous.py#L70