5

我有数字向量，例如c(1, 2, 3, 3, 2, 1, 3)or c(1, 4, 1, 4, 4, 1)，我想保留单个元素的位置，但交换/反转值，以便我们c(3, 2, 1, 1, 2, 3, 1)分别得到c(4, 1, 4, 1, 1, 4)。

为了实现这一点，我在下面提出了一个相当粗糙和丑陋的代码，其中包含大量的调试和修补......

blah <- c(1, 4, 1, 4, 4, 1, 3)
blah.uniq <- sort(unique(blah))
blah.uniq.len <- length(blah.uniq)
j <- 1
end <- ceiling(blah.uniq.len / 2)
if(end == 1) {end <- 2} # special case like c(1,4,1), should get c(4,1,4) 
for(i in blah.uniq.len:end) {
  x <- blah == blah.uniq[i]
  y <- blah == blah.uniq[j]
  blah[x] <- blah.uniq[j]
  blah[y] <- blah.uniq[i]
  j = j + 1
}
blah

有没有更简单的方法来做到这一点？

4

3 回答 3

12

我认为您正在尝试反转分数。算法是 (1 + max(x_i)) - x_i

所以...

x <- c(1, 2, 3, 3, 2, 1, 3)
y <- c(1, 4, 1, 4, 4, 1)

(max(x, na.rm=T) + 1) - x
(max(y, na.rm=T) + 1) - y

产生：

> (max(x, na.rm=T) + 1) - x
[1] 3 2 1 1 2 3 1
> (max(y, na.rm=T) + 1) - y
[1] 4 1 4 1 1 4

根据 OP 的评论：

rev.score <- function(x) {
    h <- unique(x)
    a <- seq(min(h, na.rm=T), max(h, na.rm=T))
    b <- rev(a)
    dat <- data.frame(a, b)
    dat[match(x, dat[, 'a']), 2]
}

x <- c(1, 2, 3, 3, 2, 1, 3)
rev.score(x)
y <- c(1, 4, 1, 4, 4, 1)
rev.score(y)
z <- c(1, 5, 10, -3, -5, 2)
rev.score(z)
于 2012-11-14T04:17:04.017 回答
6

恭喜！你可能终于找到了factors 的用途，我还在寻找一个:-)

x <- c(1, 2, 3, 3, 2, 1, 3)
# [1] 1 2 3 3 2 1 3
y <- factor(x)
# [1] 1 2 3 3 2 1 3
# Levels: 1 2 3
levels(y) <- rev(levels(y))
# [1] 3 2 1 1 2 3 1
# Levels: 3 2 1

基于这个想法，这里有一个函数，它返回一个与输入具有相同类的对象：

swap <- function(x) {
    f <- factor(x)
    y <- rev(levels(f))[f]
    class(y) <- class(x)
    return(y)
}
swap(c(1, 2, 3, 3, 2, 1, 3))
# [1] 3 2 1 1 2 3 1
swap(c(1, 4, 1, 4, 4, 1))
# [1] 4 1 4 1 1 4
于 2012-11-14T12:29:31.120 回答
5

一个可能的泛化函数。

revscore <- function(x) {
  rx <- min(x):max(x)
  rev(rx)[sapply(1:length(x), function(y) match(x[y],rx))]
}

x1 <- c(-3,-1,0,-2,3,2,1)
x2 <- c(-1,0,1,2)
x3 <- 1:7

一些测试：

> x1
[1] -3 -1  0 -2  3  2  1
> revscore(x1)
[1]  3  1  0  2 -3 -2 -1

> x2
[1] -1  0  1  2
> revscore(x2)
[1]  2  1  0 -1

> x3
[1] 1 2 3 4 5 6 7
> revscore(x3)
[1] 7 6 5 4 3 2 1
于 2012-11-14T06:08:40.733 回答