A function to plot results model

This function allows you to plot an overview of the model results in function of the total length of fish

plot_cnp(mod, y, x = "tl", probs = c(0.8, 0.95))

Arguments

mod

Model output from cnp_model_mcmc()

y

Output variable(s) to be plotted. Can be a character or a character vector.

x

Variable to be put on x-axis, "biomass" or "tl"

probs

Width of the confidence

Value

a ggplot object

Examples

library(fishflux)
mod <- cnp_model_mcmc(TL = 5:15, param = list(
          Qc_m = 40, Qn_m = 10, Qp_m = 4, Dn_sd = 0.05))
plot_cnp(mod = mod, y = c("Fp", "Gp", "Wp", "Ip"),
         x = "tl", probs = c(0.5, 0.8))

plot_cnp(mod = mod, y = "Fp", x = "tl",
         probs = c(0.5, 0.8, 0.95))