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