Многомерный HPD Интервал в R

Coda и другие пакеты R позволяют рассчитывать интервалы HPD для нескольких задних цепей. Однако интервалы HPD всегда рассчитываются незначительно для каждой цепи.

Есть ли способ в R найти совместный интервал HPD с учетом корреляционной структуры цепей?