RICselect: RICselect

In PAskey/SPDT: Stocking Product Development Team analyses (SPDT)

RICselect

Description

A function to correct for gillnet selectivity

A function to convert sby_code (brood year) to age (Jan 1st birthday convention)

Usage

RICselect(FLengths_mm, Millar_model = FALSE, meshSizes_in = NULL)

sby2age(Spp_code = NULL, sby_code = NULL, Year = NULL)

Arguments

Millar_model	a TRUE/FALSE to indicate whether Millar model should be used for predictions
meshSizes_in	a vecotr of mesh sizes in inches. Only applicable to Millar model and defaults to full RIC net
Spp_code	an character string or vector describing a standard species code.
sby_code	an integer value or vecotr descibing the brood year
Year	an integer value or vecotr descibing the year of sample (to be aged)
FLengths	an integer or vector of fork lengths in mm for which to calculate relative probability of capture by RIC gillnet.

Details

This function estimates the relative size-dependent vulnerability of trout to the standard RIC gillnets. The function is used to multiply observed catches by a multiplier proportional to the fish fork length. In the future the selectivity function can be made more robust to account for non-standard combinations of gillnet panels. Using a simple multiplier to correct for selectivity only works with good sample sizes, because 0 * anything = 0.

This function calculates the age of fish based on the species-specific spawn period and the brood year.

Examples

Must be connected to VPN if working remotely


#Create a vector of fish lengths from 50 to 650mm
Fish_lengths = c(50:650)

#Estimate the relative probability of capture for each of the fish lengths
pvals = RICselect(Fish_lengths)

#Plot the selectivity function
plot(pvals~Fish_lengths)

#However, you can access this selectivity function directly by simply typing
select_lookup

Must be connected to VPN if working remotely

sby_code = 2005
Year = 2010
sby2age("RB", sby_code, Year)
sby2age("KO", sby_code, Year)

PAskey/SPDT documentation built on May 31, 2024, 12:21 a.m.