Uniform没有logp方法

Question

我一定漏掉了什么。我可以从正态分布计算对数概率，但不能计算均匀分布：

>>> import pymc as pm
>>> with pm.Model():
...     b = pm.Normal('b', 0, 1)
...
>>> b.logp({'b': 2)
array(-2.9189385332046727)

>>> with pm.Model():
...     a = pm.Uniform('a', 0, 1)
...
>>> a.logp({'a': 0.5})
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
AttributeError: 'TransformedRV' object has no attribute 'logp'

Answer 1

PyMC3将有界分布(如均匀分布)变换为无界分布。因此，您可以要求PyMC3不转换变量：

with pm.Model() as model:
    a = pm.Uniform('a', 0, 1, transform=None)
a.logp({'a': 0.5})

或者请求转换后的变量，如下所示：

with pm.Model() as model:
    a = pm.Uniform('a', 0, 1)
model.logp({'a_interval__': 0.5})

您可能需要检查model的unobserved_RVs和free_RVs属性