【发布时间】:2017-10-01 00:37:26
【问题描述】:
我想拟合混合效应线性回归。因变量是 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 或仅 Tense,它们是分类的,以及仅使用 ProfScore,它是连续的),以及仅使用两种随机效果中的一种。同样的错误总是重复。但是,我可以使用 anova(model) 查看 p 值。我想知道为什么在这种情况下我不能成功使用混合()。我还安装了最新版本的 R,但没有看到任何帖子显示这种情况下的类似错误。
这里是代码和数据集的链接:
【问题讨论】:
标签: r mixed-models