parameter_predictions_scatterplots: Hydraulic parameter predictions boxplots (the current version...

In jschap1/WSEw: River Geometry Package for Studying Water Surface Elevation-Width Relationships

Description

Makes box and whisker plots showing predicted hydraulic parameter error

Usage

parameter_predictions_scatterplots(z0.l, z0.sb, z0.nl, z0.nlsb)

Arguments

z	list of arrays containing parameter predictions from each class of model
z.true	true value of the parameter
k	exposure level
lumped	flag, if true, lump the cross sections together and just compare methods
legend	flag for legend
absolute	flag for plotting absolute error
...	other parameters as taken by boxplot()

Details

If lumped = FALSE, plots the results over each cross section separately.

Examples

expo <- seq(0.05, 0.95, length.out = 19)
k <- 10
parameter_predictions_boxplots(z0.l, z0.sb, z0.nl, z0.nlsb, z0.true.ra, k = k, lumped = TRUE, 
                               main = paste("z0 error for each method at", 100*expo[k], "percent channel exposure"), 
                               ylab = "z0 error (m)")
parameter_predictions_boxplots(z0.l, z0.sb, z0.nl, z0.nlsb, z0.true.ra, k = 12, absolute = FALSE, lumped = TRUE)
par(mfrow = c(2,2))
k <- 16
parameter_predictions_boxplots(z0.l, z0.sb, z0.nl, z0.nlsb, z0.true.ra, k = k, lumped = FALSE, absolute = FALSE, 
                               main = paste("z0 error at", 100*expo[k], "percent channel exposure"), 
                               ylab = "z0 error (m)", legend = FALSE)  
Parameter predictions scatterplots (without error)

Plot predicted z0 and A0 from each model for each reach-average cross section
Assumes there is M = 1 replicate. This should be the no error case.
This is an earlier iteration of Figure 5

jschap1/WSEw documentation built on May 12, 2019, 10:50 a.m.