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