Displaying 1 result from an estimated 1 matches for "spoptim".
2010 Jun 26
4
optim() not finding optimal values
...10.5,99.1,84.5,95.7,74.1,70.2,63.1,66.4,60.5,89.9,117.0,93.0,116.6,90.0,105.1))
pars <- c(800000,1000000,0.0001,0.17) # put all parameters into one vector
names(pars) <- c("B0","K","q","r") # name the parameters
( SPoptim <- optim(pars,SPsse,B=d$catch,CPE=d$cpe) ) # run optim()
This produces parameter estimates, however, that are not at the minimum value of the SPsse function. For example, these parameter estimates produce a smaller SPsse,
parsbox <- c(732506,1160771,0.0001484,0.4049)
names(parsbox) <...