output <-glm(formula = breaks ~ wool+tension, data = warpbreaks,

   family = poisson)

print(summary(output))

​