[Experimental]

Calculates the yield of a species for a range of fishing mortalities for that species while the fishing mortalities for the other species are held fixed.

plotYieldVsF(
  params,
  species,
  no_steps = 10,
  F_max = 3,
  F_range = seq(0, F_max, length.out = no_steps)
)

Arguments

params

An object of class MizerParams.

species

Name of the target species

no_steps

The number of steps to use. Only used if F_range is missing.

F_max

The maximum fishing mortality. Used only if F_range is missing.

F_range

A sequence of fishing mortalities at which to evaluate the yield. If missing, it is set to seq(from = 0, to = F_max, length.out = no_steps).

Value

A ggplot object

See also

getYieldVsF

Other plotting functions: animateSpectra(), displayFrames()

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.
plotYieldVsF(params, "Cod")
#> Convergence was achieved in 46.5 years.
#> Convergence was achieved in 21 years.
#> Convergence was achieved in 12 years.
#> Convergence was achieved in 9 years.
#> Convergence was achieved in 7.5 years.
#> Convergence was achieved in 6 years.
#> Convergence was achieved in 6 years.
#> Convergence was achieved in 6 years.
#> Convergence was achieved in 4.5 years.
#> Convergence was achieved in 4.5 years.
# }