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

sessionInfo()

setup_admb()
ReedfrogSizepred <- read.table("../DATA/tadpole.dat",header=TRUE)
compile_admb("tadpole")
rdat <- c(list(nobs=nrow(ReedfrogSizepred),
               nexposed=rep(10,nrow(ReedfrogSizepred))),ReedfrogSizepred)

write_dat("tadpole",rdat)
write_pin("tadpole",list(c=0.45,d=13,g=1))

st1 <- system.time(tfit_admb <- tadpole_ADMB_fit())["elapsed"]
st1 <- unname(st1)

## ugh: re-inverting vcov to get hessian
ev <- try(eigen(solve(vcov(tfit_admb)))$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")