为什么`poly`抱怨程度小于独特点的数量?

时间:2013-08-01 15:56:48

标签: r polynomial-math

我正在尝试在R中生成正交多项式,但我不断收到错误我不明白

> poly(1:1000, 50)
Error in poly(1:1000, 50) : 
  'degree' must be less than number of unique points

当然,独特点的数量是1000?这是什么意思?这是一个错误吗?如果有的话,有谁知道我可以解决这个问题吗?

编辑: 这似乎是为了获得程度> 27对于任意数量的点 - 这是一个无证的限制吗?

1 个答案:

答案 0 :(得分:12)

数字溢出。如果您查看poly的代码,您会看到它生成单个多项式项作为中间步骤:

X <- outer(x, seq_len(n) - 1, "^")

n(你想要的多项式的次数)为50时,得到的项最多为1e132。