我在嵌套模型中遇到了很大的问题,我试图适应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
不幸的是,我甚至不知道从哪里开始查看问题所在。 有任何想法吗?请拜托好吗? =) 干杯!
答案 0 :(得分:0)
这可能更适合here,但是:condition
试图告诉您一个或多个随机效果与剩余方差混淆。正如您所描述的数据,我不明白为什么:您应该有2 * 46 * 32 = 2944总观察值,2 * 46 = 92 person
和measure
的组合,以及46 * 32 = person
和lf <- 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
)