我希望适合混合模型但具有方差结构。
我创建了数据框
DTF <- data.frame(Y = rnorm(7*4, 0, 1),
TIME = rep(1:7, 4),
GROUP = c(rep(1,7), rep(2,7), rep(3,7), rep(4,7)),
X = rep(rnorm(7, 20, 5), 4))
而不是我运行此命令
model <- lmer(Y ~ TIME:GROUP + X + (1 | GROUP), DTF, REML = FALSE,
varWeights(varFixed(~ TIME)))
我有
Error in lmer(...): 'control' is not a list; use lmerControl()
我不知道我要改变什么。 谢谢你的帮助!