use lm.fit rather than .Call(stats:::C_Cdqrls)