Goodness of fit example on the PainGLM data