e2ep_run_ycurve {StrathE2EPolar} | R Documentation

## Run a set of models to generate fishery yield curve data for either planktivorous or demersal fish.

### Description

Perform a set of StrathE2E model runs along a sequence of values of either planktivorous or demersal fish harvest ratio multiplier, saving the annual average whole-domain biomass, annual landings and annual discards of all exploitable food web guilds from each run plus the annual average biomasses of all the other living components of the food web.

### Usage

    e2ep_run_ycurve(
      model,
      selection,
      nyears = 50,
      HRvector = c(0, 0.5, 1, 1.5, 2, 2.5, 3),
      HRfixed = 1,
      csv.output = FALSE
    )

### Arguments

`model` | R-list object defining a baseline model and configuration compiled by the e2ep_read() function.
`selection` | Text string from a list identifying the fish guild for which a yield curve is to be generated. Select from: "PLANKTIV", "DEMERSAL", Remember to include the phrase within "" quotes.
`nyears` | Number of years to run the StrathE2E model at each level of harvest ratio. Default = 50.
`HRvector` | A vector of ascending order unique multiplier values to be applied to the baseline fish harvest ratio of the guild to be analysed. The values do not have to be evenly spaced (default = c(0,0.5,1.0,1.5,2.0,2.5,3.0)).
`HRfixed` | Single value of the multiplier to be applied as the alternative (planktivorous or demersal) fish harvest ratio for all runs (default=1.0).
`csv.output` | Logical. If TRUE then enable writing of csv output files (default=FALSE).

### Details

The baseline for the sequence of runs (harvest ratio multiplier = 1.0) is a model name and variant as loaded by the e2ep_read() function.

The planktivorous or demersal fish yield curve can be generated for a given fixed setting of the other (demersal or planktivorous) fish harvest ratio multiplier (default = 1.0). All other conditions are held constant as in the baseline model configuration.

The yield curve represents the catch that would be generated from the stationary state of the model attained with long-term repeating annual cycles of all driving data. Hence it is important that each simulation is run for long enough that the model attains its stationary state, which may be some distance from the baseline model initial conditions. It is recommended that each run is at least 50 years.

The data on annual average biomass and annual integrated catches stored in the returned data object (and optionally a csv file) can subsequently be plotted using the function e2ep_plot_ycurve(). Users can easily plot any of the other saved data using their own plotting code.

If csv output is selected then the resulting files in the current user results folder have names Yield_curve_data_PFHRmult-*.csv or Yield_curve_data_DFHRmult-*.csv, depending on the 'selection' argument, and * represents the model.ident text set in the prior e2ep_read() function call.

### Value

Dataframe of annual average biomss, annual landings and annual discards of all exploitable guilds, plus average biomasses of other food web components, for each level of harvest ratio multiplier.

### See Also

`e2ep_read`, `e2ep_run`, `e2ep_plot_ycurve`

### Examples

    # Load the 2011-2019 version of the Barents Sea model supplied with the package :
        model <- e2ep_read("Barents_Sea", "2011-2019", model.ident="11-19base")
    # In this example csv output is directed to a temporary folder since results.path is not set.

    # In this illustrative example the model is run for only 3 years to enable quick
    # return of results. In a real simulation nyear would be at least 50.
    # This example illustrates that the vector of planktivorous fish harvest ratio multiplers
    # does not have to be evenly spaced.
        hr=c(0,0.5,0.75,1.0,1.25,2.0,3.0)
        pf_yield_data <- e2ep_run_ycurve(model,selection="PLANKTIV", nyears=3, HRvector=hr,
                                        HRfixed=1,csv.output=FALSE)

    # View the column names of the results dataframe:
        names(pf_yield_data)

    # Plotting the results...
    # The planktivorous fish yield curve can be plotted using the function:
        e2ep_plot_ycurve(model, selection="PLANKTIV", results=pf_yield_data,
                        title="Planktivorous yield with baseline demersal fishing")

[Package StrathE2EPolar version 1.0.0 Index]