Calculates the yield of a species for a range of fishing mortalities for species across all gears at each simulation time step. Note: The function just returns the yield at the last time step of the simulation, which might not have converged, or might oscillate.

getYieldVsFcurve(params, ixSpecies, nSteps = 10, Fmax = 3, bPlot = FALSE)

Arguments

params

An object of class MizerParams.

ixSpecies

The species number to make the calculation over

nSteps

The number of steps in F to calculate

Fmax

The maximum fishing mortality

bPlot

Boolean that indicates whether a plot is to be made

Frange

The range of fishing mortalities to run over

Value

An list with yields and fishing mortalities

Examples

# \dontrun{ params <- newMultispeciesParams(NS_species_params_gears, inter)
#> Note: No h provided for some species, so using f0 and k_vb to calculate it.
#> Note: Because you have n != p, the default value is not very good.
#> Note: No ks column so calculating from critical feeding level.
#> Note: Using z0 = z0pre * w_inf ^ z0exp for missing z0 values.
#> Note: Using f0, h, lambda, kappa and the predation kernel to calculate gamma.
y <- getYieldVsFcurve(params,11,bPlot=TRUE)
# }