library(stockassessment)
setwd("run")
load("data.RData")
conf<-loadConf(dat,"../conf/model.cfg", patch=TRUE)
par<-defpar(dat,conf)
par$itrans_rho <- 4
#mymap <- list( itrans_rho = factor(NA) )
##upper = list(itrans_rho=3)
fit<-sam.fit(dat,conf,par) #,map=mymap)
if(fit$opt$convergence!=0) stop("Model did not converge.")
save(fit, file="model.RData")