R不能识别模型公式中的二次函数

Question

非常简单的问题，我正在使用Pixel数据集研究第42页的Pinheiro和Bates(nlme(。

library (nlme)
fm1Pixel <- lme(pixel ~ day + day^2, data=Pixel, random = list(Dog = ~ day, Side = ~ 1))

看起来很简单，但输出根本无法识别day^2:的存在

intervals (fm1Pixel)
Approximate 95% confidence intervals
Fixed effects
lower         est.        upper
(Intercept) 1071.415261 1093.2153217 1115.0153825
day           -1.126045   -0.1486644    0.8287158
attr(,"label")
[1] "Fixed effects:"...

如果我手动创建一个^2变量，它工作得很好，但模型本身(也有简单的lm()(无法识别^2。。。

还有其他人有同样的问题吗？

我打错什么了吗？

感谢你的帮助，

Eric

Answer 1

您需要将转换放入I()元素中，如下所示：

lme (pixel ~ day + I(day^2), data=Pixel, random = list(Dog = ~ day, Side = ~ 1))