嵌套Dirichlet过程的演示代码

嵌套Dirichlet过程的演示代码,r,cluster-analysis,dirichlet,R,Cluster Analysis,Dirichlet,我的问题是如何用R代码实现嵌套的Dirichlet进程(NDP) NDP适用于分布上的聚类和分布内的同时聚类。Rodriguez等人(2008年)提供了一个模拟示例,以说明NDP区分不同分布的能力。我试图通过重现本例的结果来学习这种方法。但未能做到这一点,因为我无法很好地理解基本分布与混合物成分的关系 模拟示例使用正态逆伽马分布NIG(0,0.01,3,1)作为基础分布。但这四种不同的分布是: 使用第4节(Rodriguez等人,2008年,第1135页)中提供的算法进行模拟。我在理解和执行此

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

NDP适用于分布上的聚类和分布内的同时聚类。Rodriguez等人(2008年)提供了一个模拟示例,以说明NDP区分不同分布的能力。我试图通过重现本例的结果来学习这种方法。但未能做到这一点,因为我无法很好地理解基本分布与混合物成分的关系

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

使用第4节(Rodriguez等人,2008年,第1135页)中提供的算法进行模拟。我在理解和执行此算法时遇到问题,尤其是第5步:


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

我自己还不能进行编码,但我最近发现了一篇论文,它使用精确推理而不是截断近似来进行模拟。我想这可能会帮助像我这样有兴趣的人,所以我在这里发布了这篇文章的链接


我喜欢这篇文章的一个优点是它写得很好,并且有源代码(在R中),以便我更好地理解这种方法

欢迎光临!Stack Overflow是一个问答论坛,不是免费的代码编写服务。如果您想获得有关此问题的帮助,请编辑您的问题,以显示您已经尝试了哪些代码,并在遇到特定步骤或错误消息时询问特定问题。当然可以!还在努力。将很快编辑此问题。您现在有本文的示例代码吗?