Add natage table to forecast outputs

[HSCole]

Hi,

It would be useful to be able to see how stock numbers at age are changing through each year of the forecast.

Can this be added as an output to the forecast function? Perhaps similar to the catch numbers at age which are added to the output as an attribute [attr(simlist, "caytable") <- caytable)]

Thanks

[fishfollower]

Thanks,

I will consider adding it to the output, but for now, try:

fc <- forecast(fit, fscale=c(1,1,1,1), savesim=TRUE)
N <- ntable(fit)
Nfc <- do.call(rbind, lapply(fc, function(x)exp(colMeans(x$sim[,1:ncol(N)]))))
rownames(Nfc) <- lapply(fc, function(x)x$year)
colnames(Nfc) <- colnames(N)
Nfc

Best wishes,

Anders.