# NOT RUN {
data(dataspm)
plotfishM(dataspm,spsname="Pink Ling")
pars <- log(c(0.264,4740,3064,0.05))
ans <- plotspmmod(pars,dataspm,schaefer=FALSE,plotout=FALSE)
bestSP <- fitSPM(pars,dataspm,schaefer=FALSE,maxiter=1000)
outfit(bestSP)
ans <- plotspmmod(bestSP$estimate,dataspm,schaefer=FALSE,plotout=FALSE)
summspm(ans)
# }
Run the code above in your browser using DataLab