epsilon: Epsilon

epsilonR Documentation

Epsilon

Description

epsilon returns a vector of error terms for recruitment, based on equation 4 in Babcock & MacCall (2011).

Usage

epsilon(A = 5, TimeT = 70, CR = 6, FDR, NuR, Rho_R = 0)

Arguments

A

numeric value, the number of total areas in the model. Default value is 5.

TimeT

numeric value, the number of years to run the model total. Default value is 70.

CR

numeric value, the number of control rules to be compared. Default value is 6.

FDR

numeric value, the total number of final density ratios.

NuR

numeric vector, the recruitment random normal variable, pulled from a normal distribution of mean 0 and standard deviation equal to Sigma_R.

Rho_R

numeric value, the recruitment autocorrelation on the interval (-1, 1). Default value is 0.

Value

a numeric vector of recruitment error terms, of dimensions A \* timeT \* CR.

Examples

A = 5; TimeT = 70; CR = 6; FDR = 4
NuR <- array(stats::rnorm(A*TimeT*CR*FDR, 0, 0.5), c(A, TimeT, CR, FDR))
epsilon(A, TimeT, CR, FDR, NuR, Rho_R = 0)

vquennessen/densityratio documentation built on Aug. 28, 2022, 5:36 p.m.