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

sessionInfo()

setup_admb()
compile_admb("weed")

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

results <- list(obj=-logLik(tfit_admb),  ## objective function/Neg log likelihood (NULL for MCMC)
                coef=coef(tfit_admb),  ## coefficients
                sd=rep(NA,length(coef)),
                time=c(fit=st1) )

save("results",file="fit.RData")