dat <- t(harborSealWA); dat <- dat[2:4, ]
modlist = list(U="zero", A="zero", x0="zero", Q="unconstrained")
MARSS(dat[,14:22], method="TMB", silent=TRUE, model=modlist)$logLik
MARSS(dat[,14:22], silent=TRUE, model=modlist)$logLik
maxit not getting set for nlminb
a = MARSS(dat[,14:22], method="TMB", model=modlist, control=list(maxit=20000))
false convergence for nlminb but seems ok
modlist = list(U="zero", A="zero", Q="unconstrained", tinitx=1)
a = MARSS(dat[,14:22], method="TMB", model=modlist)
a = MARSS(dat[,14:22], method="BFGS_TMB", model=modlist)
a = MARSS(dat[,14:22], model=modlist)
false convergence for nlminb and no convergence at 5000+ iter for optim,
modlist = list(U="zero", A="zero", Q="unconstrained", x0="zero", tinitx=1)
a = MARSS(dat[,14:22], method="TMB", model=modlist) # reports false convergence but better LL than EM
a = MARSS(dat[,14:22], method="BFGS_TMB", model=modlist) # not remotely converged
a = MARSS(dat[,14:22], method="BFGS", model=modlist) # bad
a = MARSS(dat[,14:22], model=modlist)
a = MARSS(dat[,14:22], method="BFGS", model=modlist, inits=a) # good same as TMB
false/no convergence but LL looks ok
modlist = list(U="zero", A="zero", Q="unconstrained")
a = MARSS(dat[,14:22], method="TMB", model=modlist)
a = MARSS(dat[,14:22], method="BFGS_TMB", model=modlist)
a = MARSS(dat[,14:22], model=modlist)
This is lower with BFGS. Even if. Perhaps has something to do with the NAs in the data. Seems to go away wo the NAs.
maxit not getting set for nlminb
false convergence for nlminb but seems ok
false convergence for nlminb and no convergence at 5000+ iter for optim,
false/no convergence but LL looks ok
Not handling convergence messages from nlminb