randJack: randJack

In laurieKell/biodyn: Biomass Dynamic Model

Description

Simulates a biodyn object for a catch series, given the parameter estimates in the param slot and variance covariance matrix http://young.physics.ucsc.edu/jackboot.pdf

Usage

randJack(n, object, sim, ...)

Arguments

n	numeric with number of simulations
object	biodyn
...	other arguments

Value

biodyn with estimates of stock based on catch time series

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))


bdRnd  =randJack(100,bd,bdJK)

plot(randJack(100,stock(bd)[,40:45],stock(bdJK)[,40:45]))

## End(Not run)

laurieKell/biodyn documentation built on May 20, 2019, 7:58 p.m.