r - lme（nlme-package）からHPD間隔を計算する

Question

誰かがmcmcサンプリングを作成しlme、パッケージから作成された線形混合モデルのHPD間隔を計算する方法を知っているかどうか疑問に思いましたnlme。私はlme4パッケージでこれを行うことに精通していますが、からの一時的な自己相関プロパティが必要nlmeです。