QRforX function

QR decomposition of linear model design matrices