R 从lme（nlme包）计算HPD间隔

我想知道是否有人知道如何创建mcmc采样并计算线性混合模型的HPD间隔，该模型由

lme

从

nlme

软件包制作而成。我熟悉如何使用

lme4

包来实现这一点，但我需要

nlme

中的时间自相关属性

这并不容易：我可能会尝试通过破解

lme

并使用

mcmcmpack

中的

MCMCmetrop1R

函数来展示如何做到这一点，但除此之外，我想不出一个简单的方法。如果这里没有答案，您也可以尝试r-sig-mixed-models@r-project.org.Hi，很抱歉回复太晚，感谢您的快速回复。但是，在当时，在我的例子中，自相关性似乎不足以让人担心。但出于好奇（可能是未来的情况），也可以为nlme计算HPD间隔。