R не распознает квадратичную c функцию в формуле модели?

Question

Действительно простой вопрос, я работаю через Пинейро и Бейтса (nlme) на стр. 42, используя набор данных Pixel.

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 ...

Кто-нибудь еще с такой же проблемой?

Я что-то печатал неправильно?

Благодарен за вашу помощь,

Эри c

Answer 1

Вам нужно поместить ваши преобразования в элемент I(), например:

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