R

时间:2015-05-25 19:30:48

标签: r statistics nested mixed lmer

我在嵌套模型中遇到了很大的问题,我试图适应R。

我有2个条件的响应时间实验,每个条件有46个人,每个条件有32个。我希望措施能够嵌套在人和人群中嵌套的条件下,但我无法让它发挥作用。

我认为应该有意义的代码是:

nestedmodel <- lmer(responsetime ~ 1 + condition + 
 (1|condition:person) + (1|person:measure), data=dat) 

但是,我得到的只是一个错误:

Error in checkNlevels(reTrms$flist, n = n, control) : 
  number of levels of each grouping factor must be < number of observations

不幸的是,我甚至不知道从哪里开始查看问题所在。 有任何想法吗?请拜托好吗? =) 干杯!

1 个答案:

答案 0 :(得分:0)

这可能更适合here,但是:condition试图告诉您一个或多个随机效果与剩余方差混淆。正如您所描述的数据,我不明白为什么:您应该有2 * 46 * 32 = 2944总观察值,2 * 46 = 92 personmeasure的组合,以及46 * 32 = personlf <- lFormula(responsetime ~ 1 + condition + (1|condition:person) + (1|person:measure), data=dat) 的1472种组合。

如果你这样做

lapply(lf$reTrms$Ztlist,dim)

然后

lmerControl(check.nobs.vs.nRE="ignore")

查看每个术语的转置随机效应设计矩阵,你会得到什么?你应该(根据你对数据的描述)看到这些矩阵分别是1472乘2944和92乘2944。

正如@MrFlick所说,一个可重复的例子会很好。您可以向我们展示的其他事项包括:

  • 无论如何都适合模型,使用with(dat,table(table(interaction(condition,person)))忽略测试,并向我们显示结果(尤其是随机效应差异和组数量的声明)
  • 向我们展示measure的结果,以提供有关每个组合的重复次数的信息(同样适用于pkg-config is required for building PyAV; aborting! Complete output from command python setup.py egg_info: pkg-config is required for building PyAV; aborting! ---------------------------------------- Command "python setup.py egg_info" failed with error code 1 in /private/var/folders/lk/fv_597y15lz7wqz2c4n3hb200000gn/T/pip-build-z1a_1u_k/av