search for: jpvt

...utine dqeqp3 Execution halted As an R novice, I am puzzled that the call to dgeqp3 has more arguments than dgeqp3 itself. The call has an extra argument "rwork". But this could explain the strange error code! The call: F77_CALL(dgeqp3)(&m, &n, REAL(A), &m, INTEGER(jpvt), REAL(tau), &tmp, &lwork, rwork, &info); The subroutine: SUBROUTINE DGEQP3( M, N, A, LDA, JPVT, TAU, WORK, LWORK, INFO ) Thanks for R! Regards Richard E. Roger Dr. R. E. Roger NSW Agriculture Spatial Information Officer Systems...

Hi, I think there is a problem with the LAPACK version of qr() in version 1.7.0. (version below). 1. The documentation states that LAPACK=TRUE is the default, but the code has LAPACK=FALSE. 2. With LAPACK=TRUE qr() is never pivoting, even in cases where it very clearly should be. e.g. set.seed(0) X<-matrix(rnorm(40),10,4);X[,1]<-X[,2] qrx<-qr(X,LAPACK=TRUE) qrx$pivot # note, no

Dear all, I’m not sure if it is O.K. to ask this question here. But where can I find the code for the function ‘C_Cdqrls’ which is called by the R function ‘lsfit‘. Thank you all. Sorry for being naïve if so. -------------------- Ziqiang Zhao 2013-04-23 [[alternative HTML version deleted]]

...A. This seems somewhat inconvenient. So I'd really like R to return this: > x y zero 0 one 7 More on the scrambling; I see in dqrsl.f that: c if pivoting was requested in dqrdc, the j-th c component of b will be associated with column jpvt(j) c of the original matrix x that was input into dqrdc.) which I think is referring to the helpful un-permuting, but I think qr.coef() in R needs to correspondingly un-permute the dimnames as well? Code from qr.coef: if (!is.null(nam <- colnames(qr$qr))) rownames...

...lly > like > R to return this: > > > x > y > zero 0 > one 7 > > More on the scrambling; I see in dqrsl.f that: > > c if pivoting was requested in dqrdc, the j-th > c component of b will be associated with column jpvt(j) > c of the original matrix x that was input into dqrdc.) > > which I think is referring to the helpful un-permuting, but I think > qr.coef() in R needs to correspondingly un-permute the dimnames as > well? > > Code from qr.coef: > if (!is.null(nam <-...