trendMSY: trendMSY calculates the mean MSY per K class

In haddonm/datalowSA: A Package to Faciliate Application of Data poor Stock Assessments

trendMSY calculates the mean MSY per K class

Description

trendMSY subdivides the range of the successful K values into a set of classes of width 'inc' and for each class calculates the mean (geometric mean) of the MSY for the subset of r and K values. This enables the central value of MSY to be plotted on the scatter of successful points.

Usage

trendMSY(inr, inK, inc = 100)

Arguments

inr	the set of r values to be tested, derives from summarycMSY
inK	the set of K values to be tested, derives from summarycMSY
inc	the class width of the K classes

Value

a matrix of r, K and MSY values

Examples

## 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) 
 summcMSY <- summarycMSY(answer,fish,final=TRUE) 
 trendMSY(summcMSY$r,summcMSY$K,inc=200) 

## End(Not run)

haddonm/datalowSA documentation built on Nov. 5, 2023, 6:40 p.m.