title: High (or very high) order polynomial regression in R (or alternatives?)
body: <p>I would like to fit a (very) high order regression to a set of data in R, however the <code>poly()</code> function has a limit of order 25.</p>
  
  <p>For this application I need an order on the range of 100 to 120.</p>
  
  <pre><code>model &lt;- lm(noisy.y ~ poly(q,50))
  # Error in poly(q, 50) : 'degree' must be less than number of unique points
  model &lt;- lm(noisy.y ~ poly(q,30))
  # Error in poly(q, 30) : 'degree' must be less than number of unique points
  model &lt;- lm(noisy.y ~ poly(q,25))
  # OK
  </code></pre>