[R] Anomaly in MIN (MAX) ?

Martin Maechler maechler at stat.math.ethz.ch
Mon Jul 5 15:20:54 CEST 1999


    Athula> I have observed the following differences in the use of min/max
    Athula> functions between R and Splus (4.5)

    Athula> R (0.64.1):

    Athula> min(c(NA,NA,NA),na.rm=T)

    Athula> returns the following

    Athula> Warning in min(..., na.rm = na.rm) : no finite arguments to
    Athula> min/max; returning extreme.  
    Athula> [1] 2147483647


    Athula> S-Plus 4.5( Win32 edition)

    Athula> min(c(NA,NA,NA),na.rm=T)

    Athula> returns
    Athula> [1] NA


    Athula> The S-PLUS behaviour seems to be the right one

Some of us disagree.
One property of min()/max() should be the following:

x , y : any (vector like) objects :

  min(x1,y,.., x2,...) == min(x1, min(y),.. min(x2),..)

This property is fulfilled by R, but not by S-plus.
The idea is that the logically entails

    min({empty}) = +Inf
    max({empty}) = -Inf

and then  min({empty},x) == min(x) = min(x, min({empty})
etc.

Only problem being that "+/-Inf" for integers have to be coded as
"+/- Max.Integer".
---------

The only thing I could think of changing
  --  with some extra effort, unfortunately -- 
is the default type from "integer" to "numeric" ("double").

Reason why this seems to  need quite an extra effort:

Think of  
      min(i1,i2,i3,i4,  d1,...)
where i* are integer (vectors) and d1 is a double one.

The current implementation starts with integer type (and hence +Inf = +MaxInt)
and only switches to double when 'd1' is encountered.

    Athula> Thanks,
    Athula> Athula.


    Athula> -- Athula Herath, PhD Oxford GlycoSciences,10, The Quadrant
    Athula> Abingdon Science Park, Abingdon, OXON OX14 3YS
    Athula> Email:Athula.Herath at ogs.co.uk Phone:+44 [0]1235 544202
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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