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