具有负方差的混合效应模型

Question

我知道这是一个老问题，但我只是想知道现在是否有解决方案。

我通常使用lme4 包和lmer 函数来执行混合效果模型。但是，我知道这个函数不允许我在模型中包含负方差分量。

我真的很想在模型中包含 R 的负方差。有人对我将使用哪些包有任何建议吗？或者，新的lme4 允许吗？

Answer 1

新的lme4 不允许，nlme 也不允许。看起来ASReml 可能会这样做，如果您按照here 的描述设置 IGU 参数——但 ASReml 是商业的，所以您必须购买许可证。

处理这种情况的统计上更合理的方法通常是拟合复合对称方差结构，这将允许组内的负相关。您可以通过nlme 或在几个更具实验性的框架（例如lme4 的“flexLambda”开发分支）中完成此操作。

更多讨论来自邮件列表here。