Как преобразовать ежедневный набор данных в месячный набор данных с 95% квантильными значениями в ежемесячном формате?

Question

У меня есть набор дневных данных за 10 лет. Я хочу преобразовать его в месячный набор данных за 10 лет, однако месячные значения будут в квантилях, а не в среднем за месяц. Как это сделать с помощью R?

Answer 1

Вы можете выполнить эту задачу с пакетом tsibble .Я использую данные о погоде из пакета nycflights13 R, то есть почасовой набор данных.

library(nycflights13)
library(tsibble)

nycflights13::weather %>% 
  mutate(time_month = yearmonth(time_hour)) %>% 
  group_by(origin, time_month) %>% 
  summarise(temp_q25 = quantile(temp, probs = .25, na.rm=T))

   origin time_month temp_q25
   <chr>       <mth>    <dbl>
 1 EWR      2013 Jan     28.9
 2 EWR      2013 Feb     28.9
 3 EWR      2013 Mär     35.1
 4 EWR      2013 Apr     46.9
 5 EWR      2013 Mai     55.9
 6 EWR      2013 Jun     66.9
 7 EWR      2013 Jul     75.9
 8 EWR      2013 Aug     71.1
 9 EWR      2013 Sep     60.1
10 EWR      2013 Okt     53.1