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

Question

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