data-raw/larkin.R

# Simulate
larkin <- larkin::sim(
  alpha = 2,
  beta = c(8, 6, 4, 2),
  init = rep(0.1, 8),
  p_bar = c(0.003, 0.917, 0.08),
  sigma = 0.1,
  burn = 100,
  span = 100,
  harvest = 0.2,
  extirp = 1e-06,
  seed = 123
)

# Plot spawner abundance
plot(larkin$spawners, type = "l")

# Write to data/
usethis::use_data(larkin, overwrite = TRUE)
luke-a-rogers/pbsedm documentation built on June 3, 2024, 5:20 a.m.