library(stockassessment)

load("run/model.RData")

rY <- t(outer(c(2024,2025,2025,2025,2024), 0:4, "-"))
#       f1   f2   f3   f4   f5 
#[1,] 2024 2025 2025 2025 2024  run 1 (exclude from and including)
#[2,] 2023 2024 2024 2024 2023  run 2
#[3,] 2022 2023 2023 2023 2022  run 3 
#[4,] 2021 2022 2022 2022 2021  run 4
#[5,] 2020 2021 2021 2021 2020  run 5

RETRO<-retro(fit, year=rY)
save(RETRO, file="run/retro.RData")