http://rosettacode.org/wiki/QR_decomposition ---------------------- R <http://rosettacode.org/wiki/Category:R>
# R has QR decomposition built-in (using LAPACK or LINPACK) a <- matrix(c(12, -51, 4, 6, 167, -68, -4, 24, -41), nrow=3, ncol=3, byrow=T) d <- qr(a) qr.Q(d) qr.R(d) # now fitting a polynomial x <- 0:10 y <- 3*x^2 + 2*x + 1 # using QR decomposition directly a <- cbind(1, x, x^2) qr.coef(qr(a), y) # using least squares a <- cbind(x, x^2) lsfit(a, y)$coefficients # using a linear model xx <- x*x m <- lm(y ~ x + xx) coef(m) ---------------------- El 6 de octubre de 2014, 14:48, José Miguel Contreras García < [email protected]> 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 > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-help-es > -- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
