Нумерація рядків у групах у кадрі даних


163

Робота з фреймом даних, подібним до цього:

set.seed(100)  
df <- data.frame(cat = c(rep("aaa", 5), rep("bbb", 5), rep("ccc", 5)), val = runif(15))             
df <- df[order(df$cat, df$val), ]  
df  

   cat        val  
1  aaa 0.05638315  
2  aaa 0.25767250  
3  aaa 0.30776611  
4  aaa 0.46854928  
5  aaa 0.55232243  
6  bbb 0.17026205  
7  bbb 0.37032054  
8  bbb 0.48377074  
9  bbb 0.54655860  
10 bbb 0.81240262  
11 ccc 0.28035384  
12 ccc 0.39848790  
13 ccc 0.62499648  
14 ccc 0.76255108  
15 ccc 0.88216552 

Я намагаюся додати стовпчик з нумерацією всередині кожної групи. Зробити це таким чином, очевидно, не використовується повноваження R:

 df$num <- 1  
 for (i in 2:(length(df[,1]))) {  
   if (df[i,"cat"]==df[(i-1),"cat"]) {  
     df[i,"num"]<-df[i-1,"num"]+1  
     }  
 }  
 df  

   cat        val num  
1  aaa 0.05638315   1  
2  aaa 0.25767250   2  
3  aaa 0.30776611   3  
4  aaa 0.46854928   4  
5  aaa 0.55232243   5  
6  bbb 0.17026205   1  
7  bbb 0.37032054   2  
8  bbb 0.48377074   3  
9  bbb 0.54655860   4  
10 bbb 0.81240262   5  
11 ccc 0.28035384   1  
12 ccc 0.39848790   2  
13 ccc 0.62499648   3  
14 ccc 0.76255108   4  
15 ccc 0.88216552   5  

Який був би гарний спосіб зробити це?


1
Я б запропонував додати в назву питання щось на кшталт "seq по рівнях" або "підрахунок реплік", оскільки саме так я знайшов це питання, і саме це я шукав
crazysantaclaus

2
@crazysantaclaus Якби це заголовок, я б не знайшов того, що шукав :-( Я буквально шукав "як
нумерувати

Відповіді:


280

Використовуйте ave, ddply, dplyrабо data.table:

df$num <- ave(df$val, df$cat, FUN = seq_along)

або:

library(plyr)
ddply(df, .(cat), mutate, id = seq_along(val))

або:

library(dplyr)
df %>% group_by(cat) %>% mutate(id = row_number())

або (найбільш ефективна пам'ять, як вона призначається за посиланням всередині DT):

library(data.table)
DT <- data.table(df)

DT[, id := seq_len(.N), by = cat]
DT[, id := rowid(cat)]

2
Можливо, варто згадати, що aveдає float замість int тут. Також можна змінити df$valна seq_len(nrow(df)). Я щойно натрапив на це тут: stackoverflow.com/questions/42796857/…
Френк

1
Цікаво, що це data.tableрішення здається швидшим, ніж використання frank: library(microbenchmark); microbenchmark(a = DT[, .(val ,num = frank(val)), by = list(cat)] ,b =DT[, .(val , id = seq_len(.N)), by = list(cat)] , times = 1000L)
hannes101

4
Дякую! dplyrРішення добре. Але якщо, як і я, ви намагалися отримувати дивні помилки при спробі такого підходу, переконайтеся, що у вас не виникають конфлікти між ними, plyrі dplyrяк це пояснено в цьому дописі, цього можна уникнути, прямо зателефонувавшиdplyr::mutate(...)
EcologyTom

2
інший data.tableметодsetDT(df)[, id:=rleid(val), by=.(cat)]
chinsoon12

Як модифікувати library(plyr)та library(dplyr)відповідати, щоб скласти рейтинг val-колонки у порядку зменшення?
Перемислав Ремін

26

Для виготовлення цього питання більш повне, базова R альтернатива з sequenceі rle:

df$num <- sequence(rle(df$cat)$lengths)

що дає намічений результат:

> df
   cat        val num
4  aaa 0.05638315   1
2  aaa 0.25767250   2
1  aaa 0.30776611   3
5  aaa 0.46854928   4
3  aaa 0.55232243   5
10 bbb 0.17026205   1
8  bbb 0.37032054   2
6  bbb 0.48377074   3
9  bbb 0.54655860   4
7  bbb 0.81240262   5
13 ccc 0.28035384   1
14 ccc 0.39848790   2
11 ccc 0.62499648   3
15 ccc 0.76255108   4
12 ccc 0.88216552   5

Якщо df$catє коефіцієнтом змінної, його потрібно спочатку обговорити as.character:

df$num <- sequence(rle(as.character(df$cat))$lengths)

Щойно помічено, що для цього рішення потрібно catсортувати стовпчик?
zx8754

@ zx8754 так, якщо ви не хочете нумерувати послідовні випадкиcat
Jaap

9

Ось варіант, що використовує forцикл за групами, а не за рядками (як це робив OP)

for (i in unique(df$cat)) df$num[df$cat == i] <- seq_len(sum(df$cat == i))

9

Ось невеликий трюк вдосконалення, який дозволяє сортувати 'val' всередині груп:

# 1. Data set
set.seed(100)
df <- data.frame(
  cat = c(rep("aaa", 5), rep("ccc", 5), rep("bbb", 5)), 
  val = runif(15))             

# 2. 'dplyr' approach
df %>% 
  arrange(cat, val) %>% 
  group_by(cat) %>% 
  mutate(id = row_number())

Ви не можете сортувати за group_by?
zcoleman

6

Я хотів би додати data.tableваріант, використовуючи rank()функцію, яка надає додаткову можливість змінити замовлення і, таким чином, робить його трохи більш гнучким, ніж seq_len()рішення, і досить схожий на функції row_number в RDBMS.

# Variant with ascending ordering
library(data.table)
dt <- data.table(df)
dt[, .( val
   , num = rank(val))
    , by = list(cat)][order(cat, num),]

    cat        val num
 1: aaa 0.05638315   1
 2: aaa 0.25767250   2
 3: aaa 0.30776611   3
 4: aaa 0.46854928   4
 5: aaa 0.55232243   5
 6: bbb 0.17026205   1
 7: bbb 0.37032054   2
 8: bbb 0.48377074   3
 9: bbb 0.54655860   4
10: bbb 0.81240262   5
11: ccc 0.28035384   1
12: ccc 0.39848790   2
13: ccc 0.62499648   3
14: ccc 0.76255108   4

# Variant with descending ordering
dt[, .( val
   , num = rank(-val))
    , by = list(cat)][order(cat, num),]

5

Іншою dplyrможливістю може бути:

df %>%
 group_by(cat) %>%
 mutate(num = 1:n())

   cat      val   num
   <fct>  <dbl> <int>
 1 aaa   0.0564     1
 2 aaa   0.258      2
 3 aaa   0.308      3
 4 aaa   0.469      4
 5 aaa   0.552      5
 6 bbb   0.170      1
 7 bbb   0.370      2
 8 bbb   0.484      3
 9 bbb   0.547      4
10 bbb   0.812      5
11 ccc   0.280      1
12 ccc   0.398      2
13 ccc   0.625      3
14 ccc   0.763      4
15 ccc   0.882      5

3
У деяких випадках замість 1:n()використання seq_len(n())безпечніше, якщо у вашій послідовності операцій виникає ситуація, коли n()можна повернутись 0, тому що 1:0дає вам вектор довжини два, в той час як seq_len(0)дає нульовий вектор довжини, тим самим уникаючи помилки невідповідності довжини mutate().
Брайан Стампер

0

Використання rowid()функції в data.table:

> set.seed(100)  
> df <- data.frame(cat = c(rep("aaa", 5), rep("bbb", 5), rep("ccc", 5)), val = runif(15))
> df <- df[order(df$cat, df$val), ]  
> df$num <- data.table::rowid(df$cat)
> df
   cat        val num
4  aaa 0.05638315   1
2  aaa 0.25767250   2
1  aaa 0.30776611   3
5  aaa 0.46854928   4
3  aaa 0.55232243   5
10 bbb 0.17026205   1
8  bbb 0.37032054   2
6  bbb 0.48377074   3
9  bbb 0.54655860   4
7  bbb 0.81240262   5
13 ccc 0.28035384   1
14 ccc 0.39848790   2
11 ccc 0.62499648   3
15 ccc 0.76255108   4
12 ccc 0.88216552   5

1
Дякую за вашу відповідь, але, здається, вже висвітлено останню пропозицію у відповіді @ mnel
eli-k
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.