Re: [R] Is there a bug in qr(..,LAPACK=T)

2003-07-16 Thread Adelchi Azzalini
On Wednesday 16 July 2003 20:50, Mike Meyer wrote: > Several people have kindly (and gently) pointed out that the ?qr > documentation states that rank detection does not work for the LAPACK case. >  Its my fault for assuming that rank detection did work. --Mike sprictly speaking is your "fault", h

Re: [R] Is there a bug in qr(..,LAPACK=T)

2003-07-16 Thread Mike Meyer
Several people have kindly (and gently) pointed out that the ?qr documentation states that rank detection does not work for the LAPACK case. Its my fault for assuming that rank detection did work. --Mike On Wed, 16 Jul 2003 09:54:39 -0700 Mike Meyer <[EMAIL PROTECTED]> wrote: > The following

Re: [R] Is there a bug in qr(..,LAPACK=T)

2003-07-16 Thread Ravi Varadhan
n the help page for "qr" shows using Hilbert matrix order 9. The rank can change depending on "tolerance" option, which is actually not used if LAPACK = TRUE. Ravi. - Original Message - From: Mike Meyer <[EMAIL PROTECTED]> Date: Wednesday, July 16, 2003 12:54

[R] Is there a bug in qr(..,LAPACK=T)

2003-07-16 Thread Mike Meyer
The following snippet suggests that there is either a bug in qr(,LAPACK=T), or some bug in my understanding. Note that the detected rank is correct (= 2) using the default LINPACK qr, but incorrect (=3) using LAPACK. This is running on Linux Redhat 9.0, using the lapack library that comes wi