6

次のコードを試して、predictが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

x = 7の値を見つけようとすると、次のようになります。

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

私は何が間違っているのですか？

4

1 に答える 1

16

のヘルプを読むとpredict.lm、次のような多くの引数が必要であることがわかります。newdata

newdata-予測する変数を探すためのオプションのデータフレーム。省略した場合は、近似値が使用されます。

predict(mypol, newdata = data.frame(x=7))
于 2012-09-26T00:36:19.290 に答える