我最近开始在贝叶斯框架中拟合回归混合模型，使用 MCMC 算法（实际上是 R 中的函数 MCMCglmm）。

我相信我已经了解如何诊断估计过程的收敛性（迹线、geweke 图、自相关、后验分布......）。

在贝叶斯框架中让我印象深刻的一件事是，似乎投入了大量精力来进行这些诊断，而在检查拟合模型的残差方面似乎做得很少。例如，在 MCMCglmm 中，residual.mcmc() 函数确实存在，但实际上尚未实现（即返回：“MCMCglmm 对象尚未实现残差”；predict.mcmc() 的情况相同）。其他软件包似乎也缺少它，而且在我发现的文献中更普遍地很少讨论（除了 DIC，它也被大量讨论）。

谁能指出一些有用的参考资料，以及我可以使用或修改的理想 R 代码？

非常感谢。