未能使用lmerTest获取lmer的p值

Question

我使用lmerTest和lme4运行了以下模型

model2 = lmer(log(RT)~Group*A*B*C+(1|item)+(1+A+B+C|subject),data=dt)

使用lmerTest 输入命令时会出现以下错误：

> summary(model1)
Error in `colnames<-`(`*tmp*`, value = c("Estimate", "Std. Error", "df",  : 
  length of 'dimnames' [2] not equal to array extent

我看到，对于其他用户来说，这已经是一个问题了，而且有一个用户能够绕过运行lsmeans()的问题。当我尝试的时候，我发现了一个错误：

Error in asMethod(object) : not a positive definite matrix. 

在观察协方差矩阵时，我没有看到任何NAs。请注意，如果我只是逆组因子中的对比，我就能够运行这个模型。我很难理解为什么会这样。

当我使用lme4而不是lmerTest运行同一个模型时，我能够获得汇总()的所有输出，但没有p值(如预期的那样)。pvals.fnc在lme4中已经停用了，我还没有找到替代方案。另外，model2的p值估计方法与我成功地使用lmerTest的其他模型相同，这将是很好的。

现在有人知道我该怎么做吗？任何帮助都将不胜感激！

Answer 1

如果A、B或C是因素，那么您可能会出错--这些模型还没有得到lmerTest包的支持(我们将在帮助页面中将警告消息连同对此类模型的限制放在一起)