嵌套Dirichlet过程的演示代码

Question

我的问题是如何用R代码实现嵌套Dirichlet进程(NDP)。

NDP适用于对分布进行聚类，同时在分布中进行聚类。Rodriguez等人(2008)提供了一个模拟例子，证明NDP区分不同分布的能力。我试图通过复制这个例子的结果来学习这个方法。但没有做到这一点，因为我不能很好地理解基本分布是如何与混合物的成分相关的。

模拟算例采用正态逆伽马分布NIG(0,0.01,3,1)作为基分布.但这四种不同的发行版本是：

​

​

第4节提供的算法(Rodriguez等人，2008年，第1135页)被用来进行模拟。我有问题要理解和执行这个算法，特别是步骤5：

​

​

您能提供一个示例代码来演示这个算法吗？非常感谢您的帮助！

Answer 1

我无法自己进行编码，但我发现最近的一篇论文用精确推理而不是截断近似来进行模拟。我认为这可能会帮助其他人谁有兴趣，就像我，所以我张贴链接到该论文在这里。

在这里输入链接描述

我喜欢这篇文章的好处是它写得很好，并且有源代码(在R中)，让我更好地理解这个方法。