Skip to contents

Plot metabolism input measurements in separate panels to help identify outliers and other data-quality issues before fitting a model.

Usage

plot_metab_data(
  data,
  cols = c("DO.obs", "DO.sat", "depth", "temp.water", "light")
)

Arguments

data

A data frame or tibble with a solar.time column and the measurement columns named in cols. The measurement columns must be numeric.

cols

Character vector of measurement columns to plot. Any subset of "DO.obs", "DO.sat", "depth", "temp.water", and "light". Defaults to all five. The dissolved oxygen saturation percentage panel is shown when both "DO.obs" and "DO.sat" are selected.

Value

A ggplot object. The dissolved oxygen saturation percentage is calculated as 100 * DO.obs / DO.sat; values where DO.sat is zero are shown as missing.

Examples

hours <- 0:47
data <- tibble::tibble(
  solar.time = as.POSIXct("2024-06-01", tz = "UTC") + hours * 3600,
  DO.obs = 8 + 1.5 * sin((hours - 10) / 24 * 2 * pi),
  DO.sat = 9 - 0.2 * sin((hours - 9) / 24 * 2 * pi),
  depth = 0.4,
  temp.water = 18 + 2 * sin((hours - 9) / 24 * 2 * pi),
  light = pmax(0, sin((hours - 6) / 12 * pi)) * 1500
)
plot_metab_data(data)

plot_metab_data(data, cols = c("DO.obs", "temp.water"))