Как получить звания без пробелов, когда есть связи между ценностями?

Когда в исходных данных есть связи, существует ли способ создать рейтинг без пробелов в рядах (последовательные, целочисленные значения рангов)? Пусть:

x <-  c(10, 10, 10, 5, 5, 20, 20)
rank(x)
# [1] 4.0 4.0 4.0 1.5 1.5 6.5 6.5

В этом случае желаемым результатом будет:

my_rank(x)
[1] 2 2 2 1 1 3 3

Я играл со всеми опциями для опции ties.method (average, max, min, random), ни одна из которых не предназначена для обеспечения желаемого результата.

Возможно ли это выполнить с помощью функции rank()?

Ответ 1

Я могу придумать быструю функцию для этого. Это не оптимально с циклом for, но он работает:)

x=c(1,1,2,3,4,5,8,8)

foo <- function(x){
    su=sort(unique(x))
    for (i in 1:length(su)) x[x==su[i]] = i
    return(x)
}

foo(x)

[1] 1 1 2 3 4 5 6 6

Ответ 2

Изменено решение crayola, но вместо merge используется match:

x_unique <- unique(x)
x_ranks <- rank(x_unique)
x_ranks[match(x,x_unique)]

редактировать

или в одну строку, как в комментарии @hadley:

match(x, sort(unique(x)))

Ответ 3

"Без петлевого" способа сделать это - просто обработать вектор как упорядоченный множитель, а затем преобразовать его в числовой:

> as.numeric( ordered( c( 10,10,10,10, 5,5,5, 10, 10 ) ) )
[1] 2 2 2 2 1 1 1 2 2
> as.numeric( ordered( c(0.5,0.56,0.76,0.23,0.33,0.4) ))
[1] 4 5 6 1 2 3
> as.numeric( ordered( c(1,1,2,3,4,5,8,8) ))
[1] 1 1 2 3 4 5 6 6

Обновление. Другим способом, который кажется быстрее, является использование findInterval и sort(unique()):

> x <- c( 10, 10, 10, 10, 5,5,5, 10, 10)
> findInterval( x, sort(unique(x)))
[1] 2 2 2 2 1 1 1 2 2

> x <- round( abs( rnorm(1000000)*10))
> system.time( z <- as.numeric( ordered( x )))
   user  system elapsed 
  0.996   0.025   1.021 
> system.time( z <- findInterval( x, sort(unique(x))))
   user  system elapsed 
  0.077   0.003   0.080

Ответ 4

Другая функция, которая делает это, но кажется неэффективной. Нет цикла for, но я сомневаюсь, что он более эффективен, чем предложение Sacha!

x=c(1,1,2,3,4,5,8,8)
fancy.rank <- function(x) {
    x.unique <- unique(x)
    d1 <- data.frame(x=x)
    d2 <- data.frame(x=x.unique, rank(x.unique))
    merge(d1, d2, by="x")[,2]
}

fancy.rank(x)

[1] 1 1 2 3 4 5 6 6

Ответ 5

попытайтесь подумать о другом пути

x <-  c(10,10,10,5,5,20,20)
as.numeric(as.factor(x))
[1] 2 2 2 1 1 3 3

Ответ 6

Если вы не против оставить Base-R:

library(data.table)
frank(x, ties.method = "dense")
[1] 2 2 2 1 1 3 3

данные:

x <-  c(10, 10, 10, 5, 5, 20, 20)

Ответ 7

Для тех, кто любит использовать dplyr:

dense_rank(x)

[1] 2 2 2 1 1 3 3

Ответ 8

Как насчет sort()?

x <- c(1,1,2,3,4,5)
sort(x)

> sort(x) 
[1] 1 1 2 3 4 5