library(stockassessment)
load("run/model.RData")
#fit2<-fit; fit2$data$logobs[is.na(fit2$data$logobs)]<-fit2$pl$missing
#fit2$conf$initState<-1
#fit2<-runwithout(fit2)
RES<-residuals(fit)
RESP<-procres(fit)
save(RES, RESP, file="run/residuals.RData")