library(stockassessment)
load("run/model.RData")
Retmat <- cbind(2022-0:4, 3000, 2023-0:4, 2022-0:4)
RETRO<-retro(fit, year=Retmat)
save(RETRO, file="run/retro.RData")