Oi Flávia, Extendi uma função que eu tinha para calcular intervalos HPD de funções com expressão analítica para aplicar a densidade empíricas neste link: http://leg.ufpr.br/~elias/diversos/hpd-to-density-function.R Att. Elias.