# NOT RUN {
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))
# }
Run the code above in your browser using DataLab