example.R

library(devtools)
install_github("wstolte/DelwaqR")

library(DelwaqR)

arr <- his2arr(filename = "extdata/NZBLOOM.HIS", timestamp = F, begintime = "2003-01-01 00:00:00")
dimnames(arr)
submod <- c("Chlfa", "OXY")
locmod <- c("NZR6NW020", "NZR9TS010")
df <- arr2df(arr, locmod=locmod, submod=submod)
df$value[df$variable == "fResptot"] <- -df$value[df$variable == "fResptot"]
library(ggplot2)
plot <- ggplot(df, aes(time, value))
plot +
  geom_line(aes(color = variable), size = 1) +
  geom_point(aes(color = variable), fill = "white",  shape = 21, size = 4) +
  facet_grid((. ~ location))

plot +
  geom_line(aes(), size = 1) +
  geom_point(aes(), fill = "white",  shape = 21, size = 4) +
  facet_grid((variable ~ location)) +
  theme_bw() +
  theme(text = element_text(size = 14),
        axis.text.x = element_text(angle=90, vjust=0.5))

ggsave("example.png", width = 10, height = 7)

limmod = c("Limit e", "Limit nit", "Limit pho", "Limit sil")
DelwaqEcoplot(arr = arr, locmod = locmod, submod = submod, limmod = limmod, plottype = 1)
wstolte/DelwaqR documentation built on June 20, 2021, 12:03 p.m.