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

Question

我试图在R中生成正交多项式,但我一直得到一个我不明白的错误

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

当然,独特点的数量是1000？这是什么意思？这是一个错误吗？如果有的话,有谁知道我的工作？

编辑:对于任意数量的点数,这似乎是> 27的程度 - 这是一个无证件限制吗？

Answer 1

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

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

当n(你想要的多项式的次数)为50时,结果项达到1e132.