[R-es] factorización QR
Carlos J. Gil Bellosta
cgb en datanalytics.com
Lun Oct 6 14:54:54 CEST 2014
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
Más información sobre la lista de distribución R-help-es