Working R code with naive poisson, linear, and quadratic