PyMC3：如何從轉換后的隨機變量中采樣？

Question

我有一個非常基本的PYMC3問題：如何從轉換后的RV采樣？ 我主要是出於調試目的。

例如：

import pymc3 as pm

with pm.Model():
   A = pm.Normal('A')

   B = pm.Deterministic('B', A + 1)
   # or
   B = A + 1

我可以使用A.random()從A生成樣本，但未為B定義該函數。 + 1是微不足道的變換，但如果涉及更多，我想從變形的A ，即B抽取一些樣本，以確保一切正常。

Answer 1

就是這樣PyMC3要求您首先構建一個已經完成的模型，然后使用pm.sample從后驗中采樣（通常在有數據的情況下！）。

您要查找的語法是：

with pm.Model():
    A = pm.Normal('A')
    B = pm.Deterministic('B', A + 1)
    trace = pm.sample(1000)

(trace['B'] == trace['A'] + 1).all()  # True