R 的 DPpackage 允许一个层次结构，只要你在函数 DPdensity 中的比例矩阵上建议。您可以在他们的手册或相关插图中查看他们所做的工作，以获得一些想法。让$\Sigma$是协方差矩阵。它设置$\Sigma \sim IW(\nu_1, \Psi_1)$和$\Psi_1 \sim IW(\nu_2, \Psi_2)$在哪里$IW(\nu, \Psi)$是具有自由度的逆Wishart$\nu$和意思$\frac{\Psi^{-1}}{\nu - p - 1}$在哪里$p$是数据的维度。起初这在我看来有点倒退，但如果你使用密度，你会发现它是共轭的。Wishart 密度看起来不太适合进行任何分析$\nu$. 你总是可以穿上任何东西$\nu$并使用 Metropolis-Hastings 步骤。

编辑：我刚刚注意到你正在使用锯齿。很有可能我认为如果您尝试先放任何东西，它会呕吐$\Psi_1$，即使逆Wishart 是共轭的。BUGS 实现对于它们允许的多变量分布可能变化无常，因此它可能不知道如何进行共轭更新。不过我不确定。