Various R functions
Plots uncertainty of a spict-estimated production curve by sampling from the joint posterior distribution of the parameters. Written by Cóilín Minto and Paul Bouch.
Plotting of empirical surplus-production inspired by Hilborn (2001).
library(mvtnorm)
library(spict)
source("sp_uncertainty.R")
data(pol)
rep <- fit.spict(pol$hake)
sp_uncertainty(rep)