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
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
