python - PyMC 中的狄利克雷分布

Question

有人可以解释一下在 PyMC 中使用 Dirichlet 发行版（带有工作示例）吗？

我意识到这是微不足道的，但我无法找到所有组件的踪迹。有什么出路吗？

我已经在这里发布了我的代码的相关部分

在 PyMC 中打印跟踪时出现 KeyError

Answer 1

也许困扰你的是，当你定义一个 k 分量 Dirichlet 分布时，pymc 只给出 k-1 个分量。剩余分量假定为 1 减去所有分量的总和。见这里。

顺便说一句，这种行为在 pymc3 中发生了变化，它返回了 k 个分量（显然加到了 1）。