使用lmerTest无法获得lmer的p值

时间:2014-07-15 07:15:23

标签: p-value lmer

我使用lmerTest并使用lme4

运行了以下模型
model2 = lmer(log(RT)~Group*A*B*C+(1|item)+(1+A+B+C|subject),data=dt)

使用lmerTest 键入summary()命令时出现以下错误:

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

我发现这对其他用户来说已经是一个问题,并且一个用户能够绕过运行lsmeans()的问题。 当我尝试lsmeans时,我收到了错误:

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

在研究协方差矩阵时,我没有看到任何NA。 请注意,如果我只是反转组因子中的对比,我就可以运行此模型。 我很难理解为什么会这样。

当我使用lme4而不是lmerTest运行相同的模型时,我能够获得summary()的所有输出但没有p值(如预期的那样)。 pvals.fnc在lme4中停止,我还没有找到替代方案。另外,对于model2和我成功使用lmerTest的其他模型,以相同的方式估算p值会很好。

有人知道我此刻应该做些什么吗?任何帮助将不胜感激!

1 个答案:

答案 0 :(得分:0)

如果A或B或C是因素,那么您可能会收到错误 - lmerTest包尚不支持此类模型(我们会在帮助页面中将警告消息与此类模型的限制一起提供)