library(stockassessment)
load("run/model.RData")

dat<-fit$data


datm1<-reduce(dat,fleet=2, conf=fit$conf)
confm1<-attr(datm1,"conf")
parm1<-defpar(datm1,confm1)
parm1$logSdLogN=c(0,log(0.001))
fitm1<-sam.fit(datm1,confm1,parm1, map=list(logSdLogN=factor(c(1,NA))))

datm2<-reduce(dat,fleet=3, conf=fit$conf)
confm2<-attr(datm2,"conf")
parm2<-defpar(datm2,confm2)
parm2$logSdLogN=c(0,log(0.001))
fitm2<-sam.fit(datm2,confm2,parm2, map=list(logSdLogN=factor(c(1,NA))))

runs<-c(fitm1,fitm2)

names(runs) <- paste0("w.o. ", lapply(2:3, function(x) paste(attr(fit$data, 
        "fleetNames")[x], collapse = " and ")))
    attr(runs, "fit") <- fit
    class(runs) <- "samset"
LO<-runs


save(LO, file="run/leaveout.RData")