Как добавить среднее в интервалах на участке

Question

У меня есть набор данных, состоящий из дат, начиная с 2018-04-22T11:48:53 и заканчивая 2018-04-22T12:03:24.Каждая отметка времени имеет свою ценность.Мне нужно создать функцию, которая принимает data.frame и на основе этого рисует график с реальными значениями и средними значениями для каждого интервала.

Answer 1

ОП запросил создать функцию, которая принимает data.frame и на основе этого рисует график с реальными значениями и средними значениями для каждого интервала.

Это то, что I будет соответствовать этому требованию, используя data.table для агрегирования, scales для создания "хороших" интервалов, ggplot2 для построения графиков и checkmate для проверки входных параметров:

plot_avg <- function(DF, x_arg, y_arg, interval_width, show_avg = NULL) {
  # check arguments
  checkmate::assert_data_frame(DF, min.rows = 1L, min.cols = 2L, 
                               col.names = "strict")
  checkmate::assert_string(x_arg, min.chars = 1L)
  checkmate::assert_string(y_arg, min.chars = 1L)
  checkmate::assert_subset(c(x_arg, y_arg), names(DF))
  checkmate::assert_number(interval_width, lower = .Machine$double.xmin)
  checkmate::assert_character(show_avg, null.ok = TRUE)
  checkmate::assert_subset(show_avg, c("segm", "step", ""))
  # load required packages
  library(data.table)
  library(ggplot2)
  # compute averages
  breaks <- scales::fullseq(range(DF[[x_arg]]), interval_width)
  aggDT <- as.data.table(DF)[
    , .(avg = mean(get(y_arg))), 
    by = .(start  = breaks[cut(get(x_arg), breaks, right = FALSE, labels = FALSE)])]
  # start plotting
  g <- ggplot(DF) + aes_string(x_arg, y_arg) +
    geom_point(color = "blue") + 
    {if (length(breaks) < 30) scale_x_continuous(breaks = breaks)} +
    theme_bw()
  if ("segm" %in% show_avg) 
    g <- g + geom_segment(aes(x = start, xend = start + interval_width, 
                              y = avg, yend = avg), aggDT)
  if ("step" %in% show_avg) 
    g <- g + geom_step(aes(start, avg), aggDT, linetype = "dashed")
  # return plot object
  return(g)
}

С образцом данных, созданным, как описано ниже, мы можем создавать различные графики.

Только точки данных:

plot_avg(DT, "sec", "value", 60L)

Со средними в виде горизонтальных отрезков:

plot_avg(DT, "sec", "value", 60L, "segm")

Со средними в виде горизонтальных отрезков плюс ступени

plot_avg(DT, "sec", "value", 60L, c("segm", "step"))

Ширина интервала может варьироваться:

plot_avg(DT, "sec", "value", 20L, c("segm", "step"))

Данные

ОП не предоставил общедоступных образцов данных.Итак, я должен был создать свой собственный набор данных:

library(data.table)
secs <- seq(lubridate::ymd_hms("2018-04-22T11:48:53"), 
            lubridate::ymd_hms("2018-04-22T12:03:24"), 
            by = "sec")
n_secs <- length(secs)
n_row <- as.integer(n_secs / 10)
set.seed(0)
DT <- data.table(times = sort(sample(secs, n_row)))
DT[, sec := as.integer(times - min(times))]
f <- 2*pi/n_secs
DT[, value := cos(f*sec) + sin(2*f*sec) + 0.1 * rnorm(.N)]