qr.coef(qr(X), y)