r - 嵌套狄利克雷过程的演示代码

Question

我的问题是关于如何使用 R 代码实现嵌套的狄利克雷过程（NDP）。

NDP 适用于对分布进行聚类并同时在分布内进行聚类。罗德里格斯等人。(2008) 提供了一个模拟示例来证明 NDP 区分不同分布的能力。我试图通过重现此示例的结果来学习这种方法。但未能这样做，因为我无法很好地理解基本分布与混合物成分的关系。

模拟示例使用正态反伽马分布 NIG(0,0.01,3,1) 作为基本分布。但是四种不同的分布是：

第 4 节 (Rodriguez et al.,2008, p.1135) 中提供的算法用于进行模拟。我很难理解和执行这个算法，尤其是第 5 步：

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

Answer 1

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

在此处输入链接描述

我喜欢这篇论文的好处是它写得很好并且有源代码（R）让我更好地理解这种方法。