e2ep_calculate_hrscale {StrathE2EPolar} | R Documentation

## Calculate initial values of the scaling parameters for the fishing fleet model which link effort to harvest ratios.

### Description

The scaling parameters convert the effort applied by each gear to each living resource guild in the model, into a harvest ratio (proportion of mass captured per day). Effort is the product of activity and the relative power of each gear with respect to each guild.

### Usage

    e2ep_calculate_hrscale(model)

### Arguments

`model` | R-list object generated by the e2ep_read() function which defined the model configuration

### Details

In order to estimate these scaling parameters, data are needed for a 'calibration' time period when activity, catch and harvest ratio are all known.

The function assumes that the relevant model configuration has already been loaded using the e2ep_read() function. It is expected that this has loaded the array of power parameters (catch per unit activity by gear and guild during the calibration period) given in the file /Param/fishing_power_*.csv, the corresponding activity in the file /Param/fishing_activity_*.csv, and the known harvest ratios of each guild in the file /Target/region_harvest_r_*.csv

The function returns rough estimates of the scaling parameters. These are only rough because they assume that the effort and resource biomass are distributed uniformly over the model domain. Hence they should be used as initial estimates to be refined by optimization.

The scaling parameters are dispayed on the screen and saved in a dataframe. The values then need to be manually edited into rows 12-21 (excluding the header row) of the file /Param/fishing_fleet_*.csv

### Value

Dataframe of the effort-harvest ratio scaling parameter values for each guild

### See Also

`e2ep_ls` , `e2ep_read`

### Examples

    # Load the 2011-2019 version of the Barents Sea model supplied with the package and calculate
    # scaling parameter values:
        model <- e2ep_read("Barents_Sea", "2011-2019", quiet=FALSE)
        scale_values <- e2ep_calculate_hrscale(model)
        scale_values

[Package StrathE2EPolar version 2.1.0 Index]