compareToGLM function

Compare renewals fit to glm models fit