title: Polynomial data and R's glm()
body: <p>How can you get R's glm() to match polynomial data? I've tried several iterations of 'family=AAA(link="BBB")' but I can't seem to get trivial predictions to match.</p>
  
  <p>For example, please help with R's glm to match polynomial data</p>
  
  <pre><code>x=seq(-6,6,2)
  y=x*x
  parabola=data.frame(x,y)
  plot(parabola)
  model=glm(y~x,dat=parabola)
  test=data.frame(x=seq(-5,5,2))
  test$y=predict(model,test)
  plot(test)
  </code></pre>
  
  <p>The plot(parabola) looks as expected, but I can find the incantation of glm() that will make plot(test) look parabolic.</p>