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

## 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.

# }