使用mixed()时选择lmerTest和未定义列的计算错误

时间:2017-10-01 00:37:27

标签: r mixed-models

我想要适合混合效果线性回归。因变量是4点评定量表的可接受性判断(完全不可接受,完全不可接受)。这些判断已被赋值为数值(1,2,3,4),并且该向量居中并缩放。

我使用以下代码调用模型:

ln1 = lmer(RatingNorm ~ Group + ProfScore + RegularityInflectedForm + RegRhyme*SimilarityReal + Tense + VerbClass + (1|SubjectID) + (1|Infinitive), data=AJT1) 

然后尝试p值:

混合(ln1,AJT1)

安装模型后没有出现任何错误消息。使用afex包中的mixed()来获取p值会产生一个奇怪的错误消息。

Fitting one lmer() model. [DONE]
Calculating p-values. 
anova from lme4 is returned
some computational error has occurred in lmerTest
Error in `[.data.frame`(anova_table, , c("NumDF", "DenDF", "F.value",  : 
  undefined columns selected

当我使用lmerTest包调用相同的模型时,这已经重复了。我还尝试过只有一种固定效应的简单模型(只有Group或Just Tense,它们是分类的,只有ProfScore,它是连续的),并且只使用两种随机效应中的一种。同样的错误总是重复。但是,我能够使用anova(模型)来查看p值。我想知道为什么在这种情况下我不能成功使用mixed()。我也安装了最新版本的R,并没有看到任何帖子显示类似错误的情况。

以下是代码和数据集的链接:

1 个答案:

答案 0 :(得分:0)

我在尝试在R中使用piecewiseSEM时遇到了类似的问题。如果我删除了DVpopStand模型的混合模型部分(并使用lm),它运行正常。

  

dumpster = list(

     

lmer(DVpopStand~Destage + Site +(1 | fLine / fBowl)),

     

lmer(centroidStand~Dreatment + Site +(1 | fLine / fBowl)),

     

lmer(粉碎〜治疗+网站+质心标准+ DVpopStand +(1 | fLine / fBowl))   )

     

Model.result = sem.fit(dumpster,data = F1culled)

并收到错误消息:

  

返回lme4的摘要

     

lmerTest中发生了一些计算错误

     

[.data.frame(ret,3:4)中的错误:选择了未定义的列