我想要适合混合效果线性回归。因变量是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,并没有看到任何帖子显示类似错误的情况。
以下是代码和数据集的链接:
答案 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)中的错误:选择了未定义的列