statistical mode function

Warnes, Gregory R <gregory_r_warnes at groton.pfizer.com> writes:

From: ben at zoo.ufl.edu [mailto:ben at zoo.ufl.edu]

[snip]

  How about

mode1 <- function(x) {
  t0 <- table(x)
  as.numeric(names(t0)[t0==max(t0)])
}

e.g.
 mode1(rbinom(100,10,0.4))

[snip]

  If someone else can come up with a one-liner I'd be mildly 
interested.

How about:

as.numeric(names(rev(sort(table(x))))[1])

That's better than:

        as.numeric(names(table(x)[table(x) == max(table(x))])[1])

which I have been using.

Mark

--
Mark Myatt


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._