Я хочу соответствовать смешанной модели, но с дисперсионной структурой.
Я создаю свой фрейм данных
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()
Я не знаю, что мне нужно изменить.Спасибо за вашу помощь!