robustASPM | R Documentation |
robustASPM conducts a robustness test on the quality of fit of an ASPM. This is done by using the original optimal model parameters or the original guessed parameter values, add random variation to each of them, and re-fit the model. This process needs to be repeated multiple times. This should enable an analysis of the stability of the modelling outcomes. If the optimum parameters are used then add more variation, if initial guesses are used you may need to select different starting points so that the random variation covers the parameter space reasonably well.
robustASPM(
inpar,
fish,
glb,
props,
N = 10,
scaler = 15,
Hrange = c(0.01, 0.45, 0.005),
numyrs = 50,
console = TRUE
)
inpar |
the parameter set to begin the trials with |
fish |
the fisheries data: at least year, catch, and cpue |
glb |
the global variables containing the biological information |
props |
the properties calculated from the globals |
N |
the number of random trials to run; defaults to 10, which is too few |
scaler |
the divisor that sets the degree of normal random variation to add to the parameter values; default = 15 the smaller the value the more variable the outcome |
Hrange |
a vector of three numbers denoting the range of harvest rates to use when characterizing the productivity implied by each fitted parameter set. defaults to c(0.01, 0.45, 0.005) |
numyrs |
the number of years used to drive each harvest rate to an equilibrium population structure |
console |
print summary statistics to the screen? default = TRUE |
a list of results from each run, the range of values across runs, and the median values.
## Not run:
data(dataspm)
fish <- dataspm$fish
glb <- dataspm$glb
props <- dataspm$props
pars <- c(14,0.19,0.6)
out <- robustASPM(pars,fish,glb,props)
str(out)
print(out$results)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.