Usando prever para encontrar valores do modelo não linear
Eu estou tentando o próximo código para tentar ver se prever pode me ajudar a encontrar os valores da variável dependente para um polinômio de ordem 2, neste caso, é óbvio 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
Se eu tentar encontrar o valor de x = 7, recebo isto:
> predict(mypol, 7)
Error in eval(predvars, data, env) : not that many frames on the stack
O que estou fazendo de errado?