inst/doc/hake-vignette.R

In smartell/CatchMSY: Catch-MSY method for estimating MSY-based reference points

## ---- echo = FALSE-------------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>")

## ------------------------------------------------------------------------
	library(catchMSY)
	data(NamibianHake)
	names(hake)
	head(hake$data)

## ------------------------------------------------------------------------
	hake <- runModel(hake)
	str(hake)

## ---- echo=FALSE, results='asis'-----------------------------------------
library(catchMSY)
data(NamibianHake)
knitr::kable(head(hake$data, 10))

## ------------------------------------------------------------------------
library(catchMSY)
library(ggplot2)


# load the hake sID object
data(NamibianHake)

# run the age structured model.
hake <- runModel(hake)

# plot the biomass
qplot(hake$year,hake$bt,geom="line")