R version 3.1.0 (2014-04-10)
lmer package version 1.1-6
lmerTest package version 2.0-6
我目前正在使用lmer和lmerTest进行分析 . 每次我向随机结构添加效果时,运行summary()时会出现以下错误:
#Fitting a mixed model:
TRT5ToVerb.lmer3 = lmer(TRT5ToVerb ~ Group + Condition + (1+Condition|Participant) + (1|Trial), data=AllData, REML=FALSE, na.action=na.omit)
summary(TRT5ToVerb.lmer3)
Error in `colnames<-`(`*tmp*`, value = c("Estimate", "Std. Error", "df", : length of 'dimnames' [2] not equal to array extent
如果我离开这样的结构:
TRT5ToVerb.lmer2 = lmer(TRT5ToVerb ~ Group + Condition + (1|Participant) + (1|Trial), data=AllData, REML=FALSE, na.action=na.omit)
没有错误运行摘要(TRT5ToVerb.lmer2),返回AIC,BIC,logLik偏差,随机效应的估计,固定效应的估计及其相应的p值等等 .
所以,当我运行lmerTest时,显然会发生一些事情,尽管TRT5ToVerb.lmer3对象就在那里 . 两者之间的唯一区别是随机结构:(1条件|参与者)与(1 |参与者)
我的数据的一些特征:
-
条件和组都是分类变量:条件包含3个级别,第2组
-
因变量(TRT5ToVerb)是连续的:它对应于以ms为单位的读取时间
-
这是一项重复测量实验,每位参与者有48次观察(参与者= 28次)
我读了this threat,但我看不出明确的解决方案 . 我必须将我的数据帧转换为长格式吗?如果是这样,那么我如何在lmer中使用它呢?我希望不是那样的 .
谢谢!
免责声明:我既不是R的专家,也不是统计学的专家,所以请耐心等待 .
1 回答
(应该是评论,但是太长/代码格式化等)
这个假的例子似乎与lmerTest 2.0-6和lme4的开发版本(1.1-8;但我不希望在这个例子中与1.1-6有任何相关的差异......)工作正常
生产环境 :