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

Question

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

我通常使用带有lme4功能的lmer包执行混合效果模型。但是，我知道这个函数不允许我在模型中包含负方差组件。

我真的想在R模型中包含负差异。有没有人建议我会使用哪些软件包？或者，新的lme4是否允许它？

Answer 1

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

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

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