我正在尝试在R中生成正交多项式,但我不断收到错误我不明白
> poly(1:1000, 50)
Error in poly(1:1000, 50) :
'degree' must be less than number of unique points
当然,独特点的数量是1000?这是什么意思?这是一个错误吗?如果有的话,有谁知道我可以解决这个问题吗?
编辑: 这似乎是为了获得程度> 27对于任意数量的点 - 这是一个无证的限制吗?
答案 0 :(得分:12)
数字溢出。如果您查看poly
的代码,您会看到它生成单个多项式项作为中间步骤:
X <- outer(x, seq_len(n) - 1, "^")
当n
(你想要的多项式的次数)为50时,得到的项最多为1e132。