QR分解を行う

> mx1 <- matrix(c(1, 2, 3, 4, 5, 6, 7, 8, 9), 3, 3)
> mx2 <- qr(mx1)
> qr.Q(mx2)
           [,1]       [,2]       [,3]
[1,] -0.2672612  0.8728716  0.4082483
[2,] -0.5345225  0.2182179 -0.8164966
[3,] -0.8017837 -0.4364358  0.4082483
> qr.R(mx2)
          [,1]      [,2]          [,3]
[1,] -3.741657 -8.552360 -1.336306e+01
[2,]  0.000000  1.963961  3.927922e+00
[3,]  0.000000  0.000000  3.443426e-16