Gram-Schmidt orthogonalization
Either classical or modified algorithms. The modified algorithm is the more accurate.
gs_orth_modified(X) gs_orth(X)
X |
a numerical matrix with ncol(X) <= nrow(X) |
A list with two components, Q, R, as usually defined.
set.seed(1234) X <- matrix(rnorm(10*7), 10) gs_orth_modified(X) all.equal(gs_orth(X), gs_orth_modified(X)) all.equal(gs_orth_modified(X), givens_orth(X))
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.