Подтвердить что ты не робот

Извлечь максимальное значение в каждой группе в кадре данных

У меня есть кадр данных с переменной группировки ( "Gene" ) и переменной значения ( "Значение" ):

Gene   Value
A      12
A      10
B      3
B      5
B      6
C      1
D      3
D      4

Для каждого уровня моей переменной группировки я хочу извлечь максимальное значение. Таким образом, результатом должен быть кадр данных с одной строкой на уровне переменной группировки:

Gene   Value
A      12
B      6
C      1
D      4

Может ли aggregate сделать трюк?

4b9b3361

Ответ 1

Есть много возможностей сделать это в R. Вот некоторые из них:

df <- read.table(header = TRUE, text = 'Gene   Value
A      12
A      10
B      3
B      5
B      6
C      1
D      3
D      4')

# aggregate
aggregate(df$Value, by = list(df$Gene), max)
aggregate(Value ~ Gene, data = df, max)

# tapply
tapply(df$Value, df$Gene, max)

# split + lapply
lapply(split(df, df$Gene), function(y) max(y$Value))

# plyr
require(plyr)
ddply(df, .(Gene), summarise, Value = max(Value))

# dplyr
require(dplyr)
df %>% group_by(Gene) %>% summarise(Value = max(Value))

# data.table
require(data.table)
dt <- data.table(df)
dt[ , max(Value), by = Gene]

# doBy
require(doBy)
summaryBy(Value~Gene, data = df, FUN = max)

# sqldf
require(sqldf)
sqldf("select Gene, max(Value) as Value from df group by Gene", drv = 'SQLite')

# ave
df[as.logical(ave(df$Value, df$Gene, FUN = function(x) x == max(x))),]

Ответ 2

Использование sqldf и стандартного sql для получения максимальных значений, сгруппированных по другой переменной

https://cran.r-project.org/web/packages/sqldf/sqldf.pdf

library(sqldf)
sqldf("select max(Value),Gene from df1 group by Gene")

или

Использование превосходного пакета Hmisc для группы с применением функции (max) https://www.rdocumentation.org/packages/Hmisc/versions/4.0-3/topics/summarize

library(Hmisc)
summarize(df1$Value,df1$Gene,max)

Ответ 3

df$Gene <- as.factor(df$Gene)
do.call(rbind, lapply(split(df,df$Gene), function(x) {return(x[which.max(x$Value),])}))

Просто используя базу R