[R] Which.min bug?

Martin Maechler maechler at stat.math.ethz.ch
Sat Feb 24 00:31:59 CET 2001


>>>>> "MM" == Martin Maechler <maechler at stat.math.ethz.ch> writes:

>>>>> "Elliot" == Elliot Williams <ewilliams at ucsd.edu> writes:
    Elliot> Hi, I'm not sure this is a bug, so I thought I'd bounce it off
    Elliot> the help group first.

    Elliot> I had a dataset which I was subsetting, and occasionally I get
    Elliot> an empty subset.  If I don't check for emptiness and go
    Elliot> straight to a which.min call on the subset, the program gets a
    Elliot> big negative number back.

    Elliot> One-line Example:
  
    >>> which.max(NULL)
    Elliot> [1] -2147483647

    MM> (this happens in C code, by "NA_INTEGER + 1" which unfortunately is
    MM> *not* NA ..)

    Elliot> This caused an indexing error which was hard(ish) to track
    Elliot> down.  Should it just return NA or throw an error?  Should I
    Elliot> re-submit to r-bugs?  Or is there some deep reason for this
    Elliot> behavior that I don't understand?

    MM> Yes it is a bug (and I am the culprit), and the same for
    MM> which.max().  As I have been away for almost two weeks, I'm sorry
    MM> that nobody else has confirmed this.

    MM> I believe, for consistency

    MM>   which.min(v) should return the same as which(min(v) == v)

    MM> and hence should return

    MM>   integer(0)

which I have committed for "R-release", to be version 1.2.2 in a few days.

    MM> here.  Returning NA instead is another intuitive possibility.  What
    MM> do others think?

Martin Maechler <maechler at stat.math.ethz.ch>	http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum  LEO D10	Leonhardstr. 27
ETH (Federal Inst. Technology)	8092 Zurich	SWITZERLAND
phone: x-41-1-632-3408		fax: ...-1228			<><
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list