doproject | R Documentation |
after running the run_cMSY analysis this function projects each of the accepted biomass trajectories
doproject(intraj, constC, projn = 5, sigpR = 0.025)
intraj |
the successful biomass trajectories obtained from gettraject |
constC |
the constant catch that is to be applied as a projection |
projn |
the number of years to be projected. |
sigpR |
the process error, the same as used in run_cMSY |
the same biomass trajecotry matrix as input excpet the empty years will have been filled with projections of the surplus production dynamics made under the constant catch level.
## Not run:
data(invert)
fish <- invert$fish
glb <- invert$glb
reps <- 5000 # one would run at least 20000, preferably more
answer <- run_cMSY(fish,glb,n=reps,sigpR=0.04)
traject <- gettraject(answer$R1,projn=5)
newtraj <- doproject(traject,constC=100)
head(newtraj)
trajdepl <- makedeplet(newtraj)
plotconstC(trajdepl,endyear=2017,constC=100,target=0.40)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.