probability - 在 Stan 或 PyMC3 中是否可以组合 2 个多项分布

Question

我有 2 个不同的多项式分布，我想将它们相乘以得到一个矩阵。该矩阵与我观察到的数据相匹配。有没有办法使用 PyMC3 或 Stan 分解这些数据？有什么好的例子吗？这似乎类似于非负矩阵分解的贝叶斯版本。