plotJack: plotJack

In laurieKell/mpb: Management Procedures

plotJack

Description

Create a ggplot2 plot based on a jack knifed biodyn and plots time series of biomass and harvest rate. The basic object can then be modified by adding ggpot2 layers.

Usage

plotJack(x, ncol = 1)

Arguments

x	an object of class biodyn that has been jack knifed, i.e. by providing a jack knifed CPUE series to fit
ncol	number of colums in plot panel

Details

plotJack

Value

an ggplot2 object

Examples

## Not run: 
#simulate an object with known properties
bd=sim()
bd=window(bd,end=49)

#simulate a proxy for stock abundance
cpue=(stock(bd)[,-dims(bd)$year]+stock(bd)[,-1])/2
cpue=rlnorm(1,log(cpue),.2)

#set parameters
setParams(bd) =cpue
setControl(bd)=params(bd)
control(bd)[3:4,"phase"]=-1

#fit
bd=fit(bd,cpue)

bdJK=fit(bd,jackknife(cpue))
plotJack(bdJK)
bd  =randJack(100,bd)

## End(Not run)

laurieKell/mpb documentation built on Sept. 9, 2023, 9:47 p.m.