[R-es] factorización QR

[Carlos J. Gil Bellosta]

Hola, ¿qué tal?

Mira:

> A<-matrix(c(1,2,3,2,1,3,2,1,1), 3)
> kk <- qr(A)
> qr.Q(kk)
           [,1]       [,2]       [,3]
[1,] -0.2672612  0.7715167 -0.5773503
[2,] -0.5345225 -0.6172134 -0.5773503
[3,] -0.8017837  0.1543033  0.5773503
> qr.Q(kk) %*% t(qr.Q(kk))
              [,1]         [,2]          [,3]
[1,]  1.000000e+00 5.551115e-17 -5.551115e-17
[2,]  5.551115e-17 1.000000e+00  0.000000e+00
[3,] -5.551115e-17 0.000000e+00  1.000000e+00
> qr.R(kk)
          [,1]      [,2]      [,3]
[1,] -3.741657 -3.474396 -1.870829
[2,]  0.000000  1.388730  1.080123
[3,]  0.000000  0.000000 -1.154701

Todo está explicado en ?qr

Un saludo,

Carlos J. Gil Bellosta
http://www.datanalytics.com

El día 6 de octubre de 2014, 14:48, José Miguel Contreras García
<jmcontreras en ugr.es> escribió:
> 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
>
> --
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es