Calculates the AR1 autocorrelation coefficient (phi) for recruitment deviations.
Usage
get_rho(first_yr = 1931, last_yr = 2022, rdev)
Arguments
- first_yr
First model year.
- last_yr
Last model year.
- rdev
Vector of recruitment deviations.
Value
Estimated autocorrelation.
Examples
first_yr <- 1931
last_yr <- 2022
N <- length(first_yr:last_yr)
rdev <- arima.sim(list(order = c(1, 0, 0), ar = 0.5), n = N)
get_rho(first_yr, last_yr, rdev)
#> [1] 0.3793482
