Nothing
## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.align = 'center'
)
old <- options(useFancyQuotes = FALSE)
library("hydroroute")
## -----------------------------------------------------------------------------
Q_path <- system.file("testdata", "Q.csv", package = "hydroroute")
Q <- read.csv(Q_path)
head(Q)
## -----------------------------------------------------------------------------
relation_path <- system.file("testdata", "relation.csv", package = "hydroroute")
relation <- read.csv(relation_path)
relation
## -----------------------------------------------------------------------------
Sx <- system.file("testdata", "Events", "100000_2_2014-01-01_2014-02-28.csv",
package = "hydroroute")
Sx <- read.csv(Sx)
head(Sx)
## -----------------------------------------------------------------------------
Q_file <- system.file("testdata", "Q.csv", package = "hydroroute")
relation <- system.file("testdata", "relation.csv", package = "hydroroute")
(get_lag_file(Q_file, relation, inputsep = ",", format = "%Y-%m-%d %H:%M",
save = FALSE, overwrite = TRUE))
## -----------------------------------------------------------------------------
Q_file <- system.file("testdata", "Q.csv", package = "hydroroute")
relations_path <- file.path(tempdir(), "testdata")
dir.create(relations_path)
file.copy(list.files(system.file("testdata", package = "hydroroute"),
full.name = TRUE),
relations_path, recursive = FALSE, copy.mode = TRUE)
(get_lag_dir(Q_file, relations_path, inputsep = ",",
inputdec = ".", format = "%Y-%m-%d %H:%M", overwrite = TRUE))
## -----------------------------------------------------------------------------
# file paths
Sx <- system.file("testdata", "Events", "100000_2_2014-01-01_2014-02-28.csv",
package = "hydroroute")
Sy <- system.file("testdata", "Events", "200000_2_2014-01-01_2014-02-28.csv",
package = "hydroroute")
relation <- system.file("testdata", "relation.csv", package = "hydroroute")
# read data
Sx <- utils::read.csv(Sx)
Sy <- utils::read.csv(Sy)
relation <- utils::read.csv(relation)
relation <- relation[1:2, ]
# estimate AE, exact time matches
results <- estimate_AE(Sx, Sy, relation, timeLag = c(0, 1, 0))
## -----------------------------------------------------------------------------
results$settings
## ---- fig.width = 6, fig.height = 3-------------------------------------------
results$plot_threshold
## -----------------------------------------------------------------------------
head(results$real_AE)
## -----------------------------------------------------------------------------
initial_values_path <- system.file("testdata", "initial_value_routing.csv",
package = "hydroroute")
initials <- read.csv(initial_values_path)
initials
## -----------------------------------------------------------------------------
relation_path <- system.file("testdata", "relation.csv", package = "hydroroute")
events_path <- system.file("testdata", "Events", package = "hydroroute")
initial_values_path <- system.file("testdata", "initial_value_routing.csv",
package = "hydroroute")
res <- peaktrace(relation_path, events_path, initial_values_path)
## -----------------------------------------------------------------------------
res$`2`$settings
## ---- fig.width = 8, fig.height = 3-------------------------------------------
grid::grid.draw(res$`2`$plot_threshold)
## -----------------------------------------------------------------------------
head(res$`2`$real_AE)
## ---- fig.width = 8, fig.height = 15------------------------------------------
grid::grid.draw(res$`2`$plot_scatter)
## -----------------------------------------------------------------------------
res$`2`$models
## ---- fig.width = 8, fig.height = 10------------------------------------------
gridExtra::grid.arrange(grobs = res$`2`$plot_predict$grobs, nrow = 3, ncol = 2)
## -----------------------------------------------------------------------------
initials[1, ]
res$`2`$models[1, ]
## -----------------------------------------------------------------------------
res$`2`$models[6, ]
## -----------------------------------------------------------------------------
res$`2`$models[11, ]
## -----------------------------------------------------------------------------
res$`4`$settings
## -----------------------------------------------------------------------------
head(res$`4`$real_AE)
## ---- fig.width = 8, fig.height = 3-------------------------------------------
grid::grid.draw(res$`4`$plot_threshold)
## ---- fig.width = 8, fig.height = 15------------------------------------------
grid::grid.draw(res$`4`$plot_scatter)
## -----------------------------------------------------------------------------
res$`4`$models
## ---- fig.width = 8, fig.height = 10------------------------------------------
gridExtra::grid.arrange(grobs = res$`4`$plot_predict$grobs, nrow = 3, ncol = 2)
## -----------------------------------------------------------------------------
relation_path <- system.file("testdata", "relation.csv", package = "hydroroute")
events_path <- system.file("testdata", "Events", package = "hydroroute")
settings_path <- system.file("testdata", "Q_event_2_AMP-LAG_aut_settings.csv",
package = "hydroroute")
initials_path <- system.file("testdata", "initial_value_routing.csv",
package = "hydroroute")
real_AE <- extract_AE(relation_path, events_path, settings_path)
head(real_AE)
## -----------------------------------------------------------------------------
relation <- utils::read.csv(relation_path)
initials <- utils::read.csv(initials_path)
res <- routing(real_AE, initials, relation)
## ---- fig.width = 8, fig.height = 15------------------------------------------
grid::grid.draw(res$plot_scatter)
## -----------------------------------------------------------------------------
res$models
## ---- fig.width = 8, fig.height = 10------------------------------------------
gridExtra::grid.arrange(grobs = res$plot_predict$grobs, nrow = 3, ncol = 2)
## -----------------------------------------------------------------------------
settings_path <- system.file("testdata", package = "hydroroute")
res <- peaktrace(relation_path, events_path, initial_values_path,
settings_path)
res$`2`$settings
## ---- include = FALSE---------------------------------------------------------
options(old)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.