stock_sim {TropFishR} | R Documentation |
Stock simulation
Description
This function estimates stock size, biomass and yield of a stock from
fishing mortality per age class or length group. This function is embedded in the
Thompson and Bell model (prediction model: predict_mod
).
Usage
stock_sim(param, age_unit = "year", stock_size_1 = NA, plus_group = NA)
Arguments
param |
a list consisting of following parameters:
|
age_unit |
indicates if the age groups are per month ( |
stock_size_1 |
stock size of smallest age/length group |
plus_group |
indicates age/length group, which should be turned into a plus group (i.e. all groups above are comprised in one group) |
Details
better to treat last group always as a plus group... if stock size 1 not provided assumes 1000 as intital population size make sure that FM is also in same unit as the classes, e.g. when classes in months than also FM has to be provided in 1/months
Value
A list with the input parameters and following list objects:
-
dt: delta t,
-
N: population numbers,
-
dead: number of deaths due to natural mortality,
-
C: catch,
-
Y: yield,
-
B: biomass,
-
V: value,
-
totals: summarised output:
-
totC total catch,
-
totY total yield,
-
totV total value,
-
meanB mean biomass.
,
-
References
Garcia, S. and N.P. van Zalinge, 1982. Shrimp fishing in Kuwait: methodology for a joint analysis of the artisanal and industrial fisheries. pp. 119-142 In: Report on the Workshop on assessment of the shrimp stocks of the west coast of the Gulf between Iran and the Arabian Peninsula. Fisheries development in the Gulf. Rome, FAO, FI:DP/RAB/80/015/1, 163 p.
Millar, R. B., & Holst, R. (1997). Estimation of gillnet and hook selectivity using log-linear models. ICES Journal of Marine Science: Journal du Conseil, 54(3):471-477
Sparre, P., Venema, S.C., 1998. Introduction to tropical fish stock assessment. Part 1. Manual. FAO Fisheries Technical Paper, (306.1, Rev. 2). 407 p.
Examples
# age-based stock simulation
data(shrimps)
# option 1: without plus group
stock_sim(shrimps, age_unit = "month")
# option 2: with plus group
stock_sim(param = shrimps, age_unit = "month", plus_group = 11)
# length-based stock simulation
data(hake)
stock_sim(param = hake, stock_size_1 = 98919.3)