Я пробую следующий код, чтобы попытаться выяснить, может ли предикат помочь мне найти значения зависимой переменной для полинома порядка 2, в этом случае очевидно, что y = x ^ 2:

x <- c(1, 2, 3, 4, 5 , 6)
y <- c(1, 4, 9, 16, 25, 36)
mypol <- lm(y ~ poly(x, 2, raw=TRUE))

> mypol

Call:
lm(formula = y ~ poly(x, 2, raw = TRUE))

Coefficients:
            (Intercept)  poly(x, 2, raw = TRUE)1  poly(x, 2, raw = TRUE)2  
                      0                        0                        1  

Если я пытаюсь найти значение х = 7, я получаю это:

> predict(mypol, 7)
Error in eval(predvars, data, env) : not that many frames on the stack

Что я делаю неправильно?