我对R中predictInterval包中的merTools函数有疑问。这是一个例子
数据
library(lme4)
library(ISLR)
data(Auto)
mdl<-lmer(mpg ~ horsepower + I(horsepower^2) + displacement + I(displacement^2) + (1|name) + (1|year),data=Auto)
使用此模型，我想预测horesepower的值，同时保持displacement的平均值。
我创建了一个名为auto1的新数据框。
auto1<-Auto
auto1$displacement <- mean(auto1$displacement)
PI <- predictInterval(merMod = mdl, newdata = auto1,level = 0.95, n.sims = 1000, stat = "median",include.resid.var = TRUE)
我的问题是：我是否还必须指定displacement二次项的平均值
I(displacement^2)或者R内部使用auto1$displacement来计算预测函数中的I(displacement^2)。
horsepower相同。我是否需要单独指定二次值horsepower
或者R在内部使用horsepower来处理I(horsepower^2)。
由于