Your results show that you have a rank 2 matrix. > $rank > [1] 2 Therefore the third column of Q for your x is arbitrary. A better example would be x <- matrix(rnorm(12), 4, 3) qr(x) Rich