Estou tentando ajustar um modelo de efeitos mistos e depois usá-lo para gerar estimativas em um novo conjunto de dados que pode ter níveis diferentes. Eu esperava que as estimativas em um novo conjunto de dados usassem o valor médio dos parâmetros estimados, mas esse não parece ser o caso. Aqui está um exemplo de trabalho mínimo:

library(lme4)
d = data.frame(x = rep(1:10, times = 3),
               y = NA,
               grp = rep(1:3, each = 10))
d$y[d$grp == 1] = 1:10 + rnorm(10)
d$y[d$grp == 2] = 1:10 * 1.5 + rnorm(10)
d$y[d$grp == 3] = 1:10 * 0.5 + rnorm(10)
fit = lmer(y ~ (1+x)|grp, data = d)
newdata = data.frame(x = 1:10, grp = 4)
predict(fit, newdata = newdata, allow.new.levels = TRUE)

Neste exemplo, estou essencialmente definindo três grupos com diferentes equações de regressão (inclinações de 1, 1,5 e 0,5). No entanto, quando tento prever um novo conjunto de dados com um nível invisível, recebo uma estimativa constante. Eu esperaria que o valor esperado da inclinação e interceptação fossem usados para gerar previsões para esses novos dados. Estou esperando a coisa errada? Ou, o que estou fazendo de errado com meu código?