f_mortality: Fishing mortality

In vquennessen/densityratio: Using Density Ratios to Help Manage Fisheries In And Around Marine Reserves

Fishing mortality

Description

f_mortality returns a numeric vector that represents the fishing mortality for all ages, on the interval (0, 1) in a specific area, at a specific time step, under a certain control rule, with a certain estimate of natural mortality

Usage

f_mortality(t, cr, fdr, FM, A, Fb, E, S)

Arguments

t	temporary numeric value, the current time step.
cr	temporary numeric value, the current control rule.
fdr	temporary numeric value, the current final target density ratio.
FM	numeric array, the values of fishing mortality for all ages, areas, timesteps, control rules, and natural mortality estimates.
A	numeric value, the number of total areas in the model. Default value is 5.
Fb	numeric value, the historical fishing effort for the fished species.
E	numeric matrix, the relative fishing effort displayed in each area, at each time step, under each control rule, and for each natural mortality estimate.
S	numeric vector, the selectivities at age from age at recruitment to maximum age, on the interval (0, 1).

Value

a numeric vector that corresponds to the fishing mortality at each age in a certain area, at a certain timestep, under a certain control rule, with a certain estimate of natural mortality.

Examples

n = 34; A = 5; TimeT = 70; CR = 6; FDR = 4
FM <- array(rep(0, n*A*TimeT*CR*FDR), c(n, A, TimeT, CR, FDR))
E <- array(rep(0.2, A*TimeT*CR*FDR), c(A, TimeT, CR, FDR))
L <- length_age(Rec_age = 2, Max_age = 35, A1 = 5, L1 = 32.21, A2 = 15,
   L2 = 47.95, K = 0.2022, All_ages = FALSE)
S <- selectivity(Rec_age = 2, Max_age = 35, A1 = 5, L1 = 32.21, A2 = 15,
   L2 = 47.95, K = 0.2022, Fleets = c('sport', 'hook', 'trawl'),
   A50_up = c(2, 5, 10), A50_down = c(6, 16, 35), Alpha = c(0.33, 0.6, 0.64),
   F_fin = c(0.25, 0.06, 1), Beta = c(1.2, 0.6, 0), Cf = c(0.71, 0.28, 0.01))
f_mortality(t = 1, cr = 1, fdr = 1, FM, A = 5, Fb = 0.2, E, S)

vquennessen/densityratio documentation built on Aug. 28, 2022, 5:36 p.m.