[R-es] factorización QR

[José Miguel Contreras García]

Hola a todos

Una duda con la función qr
La factorización QR de una matriz A dará dos matrices, una Q ortogonal y 
un R triangular superior. Pero todas las funciones y paquetes que 
encuentro dan una matriz qr que no es ninguna de ellas. ¿En qué me 
equivoco? ¿Cómo consigo las dos matrices?

Por ejemplo:
 > A<-matrix(c(1,2,3,2,1,3,2,1,1), 3)
 > qr(A)$qr->kk


 > t(kk)%*%kk
          [,1]      [,2]     [,3]
[1,] 14.92857 14.009965 6.651530
[2,] 14.00997 14.111441 7.614529
[3,]  6.65153  7.614529 6.000000
 > kk%*%t(kk)
           [,1]       [,2]       [,3]
[1,] 29.571429 -8.8457246 -1.9996028
[2,] -8.845725  3.3809524 -0.3550511
[3,] -1.999603 -0.3550511  2.0876314

No es ortogonal ni triangular



Saludos

--