我有一个 控制前后的 设计(见下图)。控制组( 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