Learn R Programming

DLMtool (version 5.2)

getq3: optimize for catchability (q)

Description

Function optimizes catchability (q, where F=qE) required to get to user-specified stock depletion

Usage

getq3(x, D, SSB0, nareas, maxage, N, pyears, M_ageArray, Mat_age, Asize, Wt_age,
  V, retA, Perr, mov, SRrel, Find, Spat_targ, hs, R0a, SSBpR, aR, bR,
  bounds = c(1e-05, 15), maxF, MPA, useCPP = TRUE)

Arguments

x

Integer, the simulation number

D

A numeric vector nsim long of sampled depletion

SSB0

A numeric vector nsim long of total unfished spawning biomass

nareas

The number of spatial areas

maxage

The maximum age

N

Array of the numbers-at-age in population. Dimensions are nsim, maxage, nyears, nareas. Only values from the first year (i.e N[,,1,]) are used, which is the current N-at-age.

pyears

The number of years to project forward. Equal to 'nyears' for optimizing for q.

M_ageArray

An array (dimensions nsim, maxage, nyears+proyears) with the natural mortality-at-age and year

Mat_age

An array (dimensions nsim, maxage, proyears+nyears) with the proportion mature for each age-class

Asize

A matrix (dimensions nsim, nareas) with size of each area

Wt_age

An array (dimensions nsim, maxage, nyears+proyears) with the weight-at-age and year

V

An array (dimensions nsim, maxage, nyears+proyears) with the vulnerability-at-age and year

retA

An array (dimensions nsim, maxage, nyears+proyears) with the probability retained-at-age and year

Perr

A matrix (dimensions nsim, nyears+proyears) with the recruitment deviations

mov

An array (dimensions nsim, nareas, nareas) with the movement matrix

SRrel

A numeric vector nsim long specifying the recruitment curve to use

Find

A matrix (dimensions nsim, nyears) with the historical fishing effort

Spat_targ

A numeric vector nsim long with the spatial targeting

hs

A numeric vector nsim long with the steepness values for each simulation

R0a

A matrix (dimensions nsim, nareas) with the unfished recruitment by area

SSBpR

A matrix (dimensions nsim, nareas) with the unfished spawning-per-recruit by area

aR

A numeric vector nareas long with the Ricker SRR a values

bR

A numeric vector nareas long with the Ricker SRR b values

bounds

A numeric vector of length 2 with bounds for the optimizer

maxF

A numeric value specifying the maximum fishing mortality for any single age class

MPA

A matrix of spatial closures by year

useCPP

logical - use the CPP code? For testing purposes only