[Rd] order(x,y, decreasing = c(FALSE, TRUE)) - how / elegantly?
Martin Maechler
maechler at stat.math.ethz.ch
Thu Aug 21 19:41:51 CEST 2008
I've found the need to compute a version of order(x,y)
where I want the sort order
for *increasing* x and *decresing* y ...
something we could imagine could be provided in the future as
order(x,y, decreasing = c(FALSE, TRUE))
i.e., using a 'vectorized' decreasing argument.
{No, I'm not volunteering right now!}
I've found the following R-level solution and
like to quiz you for more elegant / faster solutions
{but I am not really interested in replacing
order(x) by sort.list(x, method="quick") and similar things;
one thing to consider *is* using an 'na.last = . ' correctly,
and I haven't had the need for that and so not bothered to "do it"}
## Here's a script with my version and a small example
## (if you want speed comparisons, use larger examples) :
orderXuYd <- function(x,y)
{
## Purpose: order(x,y): x up, y down
## ----------------------------------------------------------------------
## Arguments: x,y: vectors of the same length
## ----------------------------------------------------------------------
## Author: Martin Maechler, Date: 21 Aug 2008
ix <- order(x)
xx <- x[ix]
iy <- tapply(y[ix], xx, order, decreasing = TRUE)
## Note: 'SIMPLIFY', 'USE.NAMES', 'use.names' are just for efficiency:
unlist(mapply(`[`, split(ix,xx), iy,
SIMPLIFY = FALSE, USE.NAMES = FALSE),
use.names = FALSE)
}
x <- c(1, 1, 2, 0, 0, 2, 1, 2, 2, 0, 2)
y <- c(27, 21, 45, 11, 13, 58, 35, 74, 95, 16, 122)
ii <- orderXuYd(x,y)
## yes, this is it :
cbind(ii=ii, x=x[ii],y=y[ii])
------------------
Yes, the real reason this goes to R-devel is that it might be
neat to provide this (well, its generalization) via an enhanced
order() function.
Martin Maechler, ETH Zurich
PS: I will be basically offline all day tomorrow, so don't
expect my reactions to your ideas quickly
More information about the R-devel
mailing list