В R: как агрегировать целые строки на основе максимального значения в определенном столбце за день?

ниже приведен пример моего набора данных:

date                O3  NOX SO2
01/01/1994 06:00    7   14  29
01/01/1994 07:00    5   18  30
01/01/1994 08:00    3   18  29
01/01/1994 09:00    6   24  35
01/01/1994 10:00    20  42  52
01/02/1994 06:00    19  7   11
01/02/1994 07:00    19  6   16
01/02/1994 08:00    42  10  4
01/02/1994 09:00    25  17  26
01/02/1994 10:00    36  29  13
01/03/1994 06:00    10  1   6
01/03/1994 07:00    14  2   5
01/03/1994 08:00    24  4   4
01/03/1994 09:00    57  11  6
01/03/1994 10:00    39  34  11

Используя R, я хочу агрегировать или подмножать целые строки на основе максимального значения в данном столбце по дням, в данном случае выбирая всю строку, содержащую максимальное значение O3 за каждый день.

Мой код объединяет только столбец O3 и дату, но удаляет час, в котором он был записан.

с совокупностью:

agg.df <- aggregate(df["O3"]~., format(df["date"],"%Y-%j"),max, na.rm = TRUE)

или с пакетом dplyr:

agg.df<-df %>% collapse_by("daily") %>% group_by(date) %>% summarise(O3 = max(var))

Я хочу получить что-то вроде этого:

date             O3 NOX SO2
01/01/1994 10:00 20 42  52
02/01/1994 08:00 42 10  4 
03/01/1994 09:00 57 11  6
...

Где окончательное подмножество включает целые строки с максимальным значением в O3 в день.

Большое спасибо.

10.04.2019 00:04

Стоит ли изучать PHP в 2023-2024 годах?

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

Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией

В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.

Приемы CSS-макетирования - floats и Flexbox

Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...

Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest

В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...

Концепция локализации и ее применение в приложениях React ⚡️

Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...

Пользовательский скаляр GraphQL

Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...

Перейти к ответу Данный вопрос помечен как решенный

Ответы 2

Ответ принят как подходящий

Примечание. При этом будут сохранены все строки, в которых O3 равно максимальному значению этого дня. Если вы хотите только 1 в строке, вы можете использовать slice или top_n для достижения этой цели.

library(tidyverse)

dat %>%
  group_by(day = lubridate::date(date)) %>%
  filter(O3 == max(O3)) %>% 
  ungroup %>%
  select(date, O3, NOX, SO2, -day)

# A tibble: 3 x 4
  date                   O3   NOX   SO2
  <dttm>              <int> <int> <int>
1 1994-01-01 10:00:00    20    42    52
2 1994-02-01 08:00:00    42    10     4
3 1994-03-01 09:00:00    57    11     6

данные

dat <- structure(list(O3 = c(7L, 5L, 3L, 6L, 20L, 19L, 19L, 42L, 25L, 
                             36L, 10L, 14L, 24L, 57L, 39L), NOX = c(14L, 18L, 18L, 24L, 42L, 
                                                                    7L, 6L, 10L, 17L, 29L, 1L, 2L, 4L, 11L, 34L), SO2 = c(29L, 30L, 
                                                                                                                          29L, 35L, 52L, 11L, 16L, 4L, 26L, 13L, 6L, 5L, 4L, 6L, 11L), 
                      date = structure(c(757404000, 757407600, 757411200, 757414800, 
                                         757418400, 760082400, 760086000, 760089600, 760093200, 760096800, 
                                         762501600, 762505200, 762508800, 762512400, 762516000), tzone = "UTC", class = c("POSIXct", 
                                                                                                                          "POSIXt"))), class = "data.frame", row.names = c(NA, -15L
                                                                                                                          ), .Names = c("O3", "NOX", "SO2", "date"))

10.04.2019 01:21

Подумайте о том, чтобы рассчитать макс. O3 с помощью ave, а затем соответственно разделить на подмножества:

dat$Max_O3 <- with(dat, ave(O3, as.Date(date), FUN=max))

subset(dat, O3 == Max_O3)
#    O3 NOX SO2                date Max_O3
# 5  20  42  52 1994-01-01 10:00:00     20
# 8  42  10   4 1994-02-01 08:00:00     42
# 14 57  11   6 1994-03-01 09:00:00     57

Рад слышать! Не забудьте способ StackOverflow сказать Спасибо!

— 14.04.2019 18:27

10.04.2019 02:20