Estoy buscando un buen tutorial sobre la agrupación de datos al R
usar el proceso de dirichlet jerárquico (HDP) (uno de los métodos bayesianos no paramétricos recientes y populares).
Hay DPpackage
(en mi humilde opinión, el más completo de todos los disponibles) en R
el análisis bayesiano no paramétrico. Pero no puedo entender los ejemplos proporcionados en R News
o en el manual de referencia del paquete lo suficientemente bien como para codificar HDP.
Cualquier ayuda o puntero es apreciada.
Una implementación de C ++ de HDP para modelado de temas está disponible aquí (mire la parte inferior para ver el código de C ++)
fuente
Estos dos enlaces proporcionan algunos ejemplos de código R (y C) para implementar una mezcla normal DP:
Encontró otra referencia. El Capítulo 15 tiene el código DP winbugs:
fuente