Embora sua dúvida seja muito ampla, o pacote "languageR" faz intervalos HPD. Sobre o uso da técnica, procure um estatístico local que entenda de inferência bayesiana.
(S,f,P)
Allaman
\begin{signature}
<<>>=
Prof. Dr. Ivan Bezerra Allaman
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e Tecnológicas
Ilhéus/BA - Brasil
Fone: +55 73 3680-5596
E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com
@
\end{signature}