### Name: QR.Auxiliaries ### Title: Reconstruct the Q, R, or X Matrices from a QR Object ### Aliases: qr.X qr.Q qr.R ### Keywords: algebra array ### ** Examples p <- ncol(x <- LifeCycleSavings[,-1]) # not the 'sr' qrstr <- qr(x) # dim(x) == c(n,p) qrstr $ rank # = 4 = p Q <- qr.Q(qrstr) # dim(Q) == dim(x) R <- qr.R(qrstr) # dim(R) == ncol(x) X <- qr.X(qrstr) # X == x range(X - as.matrix(x))# ~ < 6e-12 ## X == Q %*% R if there has been no pivoting, as here. Q %*% R