library(R2jags)
source("tadpole_BUGS_funs.R")
source("../../TOOLS/misc_funs.R") ## for dropdev
## will this make trouble for sims ???

ReedfrogSizepred <- read.table("../DATA/tadpole.dat",header=TRUE)

set.seed(101)
st1 <- system.time(tfit_jags <- tadpole_BUGS_fit(ReedfrogSizepred))["elapsed"]

tfit_jags_mcmc <- dropdev(as.mcmc(tfit_jags))
results <- list(obj=NULL,
                coef=dropdev(unlist(tfit_jags$BUGSoutput$mean)),
                confint.mcmc=dropdev(HPDinterval(tfit_jags_mcmc)),
                time=c(fit=unname(st1),mcmc=unname(st1)),
                terminfo=NULL,
                convinfo=list(effsize=effectiveSize(tfit_jags_mcmc),
                  geweke=geweke.diag(tfit_jags_mcmc))
)

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