A simpler fix: return rank=NA when Lapack is used.
Keith
> On Jan 23, 2018, at 5:36 AM, Serguei Sokol wrote:
>
> Le 23/01/2018 à 08:47, Martin Maechler a écrit :
>>> Serguei Sokol
>>> on Mon, 22 Jan 2018 17:57:47 +0100 writes:
>> > Le 22/01/2018 à 17:40, Keith O'Hara a écrit :
Le 23/01/2018 à 08:47, Martin Maechler a écrit :
Serguei Sokol
on Mon, 22 Jan 2018 17:57:47 +0100 writes:
> Le 22/01/2018 à 17:40, Keith O'Hara a écrit :
>> This behavior is noted in the qr documentation, no?
>>
>> rank - the rank of x as computed by the decomposition(*
> Serguei Sokol
> on Mon, 22 Jan 2018 17:57:47 +0100 writes:
> Le 22/01/2018 à 17:40, Keith O'Hara a écrit :
>> This behavior is noted in the qr documentation, no?
>>
>> rank - the rank of x as computed by the decomposition(*): always full
rank in the LAPACK case.
I agree the result is a little confusing, but the behavior is in line with the
documentation and so not ‘unexpected’ as such...
I don’t think this is a matter of semantics, more of a ‘return the rank when we
have it for free’ situation—when A is real-valued, qr(A,LAPACK=false) calls a
modified
Le 22/01/2018 à 17:40, Keith O'Hara a écrit :
This behavior is noted in the qr documentation, no?
rank - the rank of x as computed by the decomposition(*): always full rank in
the LAPACK case.
For a me a "full rank matrix" is a matrix the rank of which is indeed
min(nrow(A), ncol(A))
but here
This behavior is noted in the qr documentation, no?
rank - the rank of x as computed by the decomposition(*): always full rank in
the LAPACK case.
> On Jan 22, 2018, at 11:21 AM, Serguei Sokol wrote:
>
> Hi,
>
> I have noticed different rank values calculated by qr() depending on
> LAPACK p
Hi,
I have noticed different rank values calculated by qr() depending on
LAPACK parameter. When it is FALSE (default) a true rank is estimated and
returned.
Unfortunately, when LAPACK is set to TRUE, the min(nrow(A), ncol(A)) is returned
which is only occasionally a true rank.
Would not it be m
