Rewrote linreg.[ch], regression.q and glm.q to use the new