R: замена значений NA на среднее значение часа с помощью dplyr

Я изучаю пакет dplyr в R, и мне это очень нравится. Но теперь я имею дело с ценностями NA в моих данных.

Я хотел бы заменить любой NA на среднее значение соответствующего часа, например, на этот очень простой пример:

#create an example
day = c(1, 1, 2, 2, 3, 3)
hour = c(8, 16, 8, 16, 8, 16)
profit = c(100, 200, 50, 60, NA, NA)
shop.data = data.frame(day, hour, profit)

#calculate the average for each hour
library(dplyr)
mean.profit <- shop.data %>%
  group_by(hour) %>%
  summarize(mean=mean(profit, na.rm=TRUE))

> mean.profit
Source: local data frame [2 x 2]

  hour mean
1    8   75
2   16  130

могу ли я использовать команду dplyr transform для замены NA дня 3 в прибыли на 75 (для 8:00) и 130 (для 16:00)?

2 ответов


попробовать

  shop.data %>% 
             group_by(hour) %>% 
             mutate(profit= ifelse(is.na(profit), mean(profit, na.rm=TRUE), profit))

  #   day hour profit
  #1   1    8    100
  #2   1   16    200
  #3   2    8     50
  #4   2   16     60
  #5   3    8     75
  #6   3   16    130

или вы можете использовать replace

  shop.data %>% 
            group_by(hour) %>%
            mutate(profit= replace(profit, is.na(profit), mean(profit, na.rm=TRUE)))

a (менее элегантный) подход с базовыми функциями:

transform(shop.data, 
          profit = ifelse(is.na(profit), 
                          ave(profit, hour, FUN = function(x) mean(x, na.rm = TRUE)), 
                          profit))

#   day hour profit
# 1   1    8    100
# 2   1   16    200
# 3   2    8     50
# 4   2   16     60
# 5   3    8     75
# 6   3   16    130