library(R2admb)  ## requires version >= 0.7.2
library(coda)
source("skate_ADMB_funs.R")

sessionInfo()

setup_admb()
compile_admb("skate", re=TRUE)
#compile_admb("skate", re=FALSE) # I think steve commented out the above so it didn't recompile every time.
st1 <- system.time(tfit_admb_mcmc<- skate_ADMB_mcmc())["elapsed"]
st1 <- unname(st1)


ev <- try(eigen(solve(vcov(tfit_admb_mcmc)))$value)
eratio <- if (inherits(ev,"try-error")) NA else min(ev)/max(ev)


  
#results <- list(obj=-logLik(tfit_admb),  ## objective function/Neg log likelihood (NULL for MCMC)
#                coef=coef(tfit_admb),  ## coefficients
#                sd=stdEr(tfit_admb),
#                confint.quad=confint(tfit_admb,method="quad"),
#                time=c(fit=st1),
#                terminfo=c(maxgrad=tfit_admb$maxgrad,eratio=eratio) ## MLE termination info
 #               )
#save("results",file="fit.RData")

save.image()
save("tfit_admb_mcmc",file="fit2_mcmc.RData")