为什么`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。