R中包lme4中多级模型的规范

时间:2018-03-05 22:03:42

标签: r lme4

我有一个 控制前后的 设计(见下图)。控制组( c )经过预先测试(pre.c)并经过后期测试(pos.c)。类似地,治疗组( t )经过预先测试(pre.t)和后测试(pos.t)。因此,我们在两个时间点(时间因素)测试了两组(组因子)。

此设计的数据显示在下方(R代码中)。

问题:我想知道下面我的R代码是否正确捕获了两个因素(即时间)以及此相互作用设计?

这是我的数据以及我尝试的内容(没有成功):

pre.c = rnorm(100)

pos.c = rnorm(100)

pre.t = rnorm(100)

pos.t = rnorm(100, 5)


data = data.frame(subjects = rep(1:200, each = 2), 
              y = c(pre.c, pos.c, pre.t, pos.t), 
              time = rep(0:1, 200), 
              group = rep(c(0, 1), each = 200))

   library(lme4)

   m1 <- lmer(y ~ time + (time | subjects), data = data)         # Gives Error

   m2 <- lmer(y ~ time * group + (time | subjects), data = data) # Gives Error

enter image description here

0 个答案:

没有答案