0

我已经使用lmerTest和 using运行了以下模型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 中已停产，我还没有找到替代方案。另外，最好以与我成功使用 lmerTest 的其他模型相同的方式估计模型 2 的 p 值。

有谁知道我此时应该做什么？任何帮助将非常感激！

4

1 回答 1

0

If A or B or C are factors then you might get errors - such models are not yet supported by the lmerTest package (we will put the warning message together with the restrictions for such models in the help page)

于 2014-07-15T16:27:29.967 回答