glm scripts based on msteps at q 0.1 and 0.05