Демонстрационный код для вложенного процесса Дирихле

Question

У меня вопрос о том, как реализовать вложенный процесс Дирихле (NDP) с кодом R.

NDP подходит для кластеризации по распределениям и одновременной кластеризации внутри распределения.Родригес и соавт.(2008) предоставили пример моделирования, чтобы продемонстрировать способность NDP различать различные распределения.Я пытаюсь изучить этот подход, воспроизводя результаты для этого примера.Но не смог этого сделать, потому что я не могу хорошо понять, как базовое распределение связано с компонентами смеси.

В примере моделирования использовалось нормальное распределение обратной гаммы, NIG (0,0.01,3,1), в качествебазовое распределение.Но четыре различных распределения:

Алгоритм, представленный в Разделе 4 (Rodriguez et al., 2008, p.1135), был использован для моделирования.У меня есть проблема, чтобы понять и выполнить этот алгоритм, особенно шаг 5:

Не могли бы вы предоставить пример кода для демонстрации этого алгоритма?Ваша помощь высоко ценится!

Answer 1

Я не смог самостоятельно выполнить кодирование, но нашел недавнюю работу, в которой моделирование выполняется с использованием точного вывода вместо усечения.Я думаю, что это может помочь кому-то, у кого такой же интерес, как и у меня, поэтому я публикую ссылку на эту статью здесь.

введите описание ссылки здесь

В этой статье мне нравится то, что она хорошо написана и имеет исходный код (на R), чтобы я мог понять эту методологиюлучше.