Displaying 1 result from an estimated 1 matches for "liqr".
Did you mean:
libr
2012 Dec 03
1
qr.qy and qr.qty give an error message when y is integer and LAPACK=TRUE
...1:nrow(A)
A.laqr <- qr(A, LAPACK=TRUE)
both qr.qy(A.laqr,y) and qr.qty(A.laqr,y) give the respective error messages
Error in qr.qy(A.laqr, y) : 'b' must be a numeric matrix
Error in qr.qty(A.laqr, y) : 'b' must be a numeric matrix
However when Lapack is not used as in
A.liqr <- qr(A, LAPACK=FALSE)
qr.qy(A.liqr,y) and qr.qty(A.liqr,y) don't issue error messages.
Looking at the source of qr.qy and qr.qty in https://svn.r-project.org/R/trunk/src/library/base/R/qr.R
I see that in the case of Lapack the storage.mode of y is not set to "double" (in contras...