PacifichakeMSE/Simple MSE/plotUncertainty.R

plotUncertainty <- function(est, data){

yl = ylimits(est$name, data)
year <- 1:length(est$name)

plot(year,est$name, ylim = yl, type = 'l', col = 'red', lwd = 2,
     ylab = '', xlab = '')
points(year,data, lwd =2)
polygon(c(year, rev(year)), c(est$max, rev(est$min)), 
        border = NA, col = "#FF000050")


}
nissandjac/commspectrum documentation built on June 1, 2019, 5:02 a.m.