这是我的问题 .
样本数据：
library(lme4)
library(ISLR)
data(Auto)
mdl<-lmer(mpg ~ horsepower + I(horsepower^2) + displacement + I(displacement^2) + (1|name) + (1|year),data=Auto)
我想使用这个模型来预测 horsepower 的范围，同时保持 displacement 的平均值 .
horsepower <- min(Auto$horsepower):max(Auto$horsepower)
displacement <- rep(mean(Auto$displacement),185)
data <- data.frame(horsepower,displacement)
# Use predict
yVals <- predict(mdl, newdata = data)
Error in eval(expr, envir, enclos) : object 'name' not found
我认为这个错误正在发生，因为我没有指定 name 这是我在 predict 函数中的随机效果 . 有谁知道如何解决这个错误 .
另外，我想根据模型生成的系数拟合 mpg ~ horsepower 图中的二次函数 .