Skip to contents

ggplot2-wrapper to summarise and plot data (of rOstluft::format_rolf() format) for diurnal time course-stats

Source: R/ggperiodic.R
ggdiurnal.Rd

ggplot2-wrapper to summarise and plot data (of rOstluft::format_rolf() format) for diurnal time course-stats

Usage

ggdiurnal(
  data,
  y = "value",
  fun = list(middle = "mean", middle2 = "median", ymin = "percentile05", lower =
    "percentile25", upper = "percentile75", ymax = "percentile95"),
  nmin = 3,
  ribbon_color = "gray60",
  middle_color = "gray30",
  expand_xscale = 0.025
)

Arguments

data

a data.frame or tibble containing the data

y

a character string specifying the taget variable to be summarised, usually "value". outcomes of cut_timeseries_periodic().

fun

function or list of functions for summary, can be named (then the outut stat is named after the function's name); Strings matching the regular expression /^percentile([0-9]){1,2}$/i are converted into the respective function. "percentile95" => function(x, ...) quantile(x, 95 / 100, ...)

nmin

numeric, minimum number of values for fun, if n < nmin: NA is returned

ribbon_color

character string for specifying the ribbon color.

middle_color

character string for specifying the middle's color.

expand_xscale

multiplicative range expansion factor

Value

ggplot2::ggplot() object

Examples

library(ggplot2)
fn <- rOstluft.data::f("Zch_Stampfenbachstrasse_2010-2014.csv")

data_min30 <-
  rOstluft::read_airmo_csv(fn) %>%
  dplyr::filter(parameter %in% c("NOx", "O3", "PM10"))

data_h1 <- rOstluft::resample(data_min30, new_interval = "h1")

ggdiurnal(data_h1)


# change scale_expands: less margin on x, no expand lower y and set to zero
ggdiurnal(data_h1, expand_xscale = 0.01) +
  scale_y_continuous(limits = c(0, NA), expand = expansion(mult = c(0, 0.05))) +
  theme(panel.spacing.y = unit(12, "pt"))


# 30min resolution, add title, x/y label (must set theme elements)
ggdiurnal(data_min30) +
  ggtitle("Diurnal Zuerich Stampfenbachstrasse 2010-2014 ") +
  ylab("parameter") +
  xlab("time") +
  theme(
    axis.title = element_text(),
    axis.title.x = element_text()
  )