source('tools.R')
sessionInfo()
file.copy('../DATA/min.dat', './min.dat')
system('admb min')
time<-unname(system.time(system('./min -nox'))["elapsed"])

fit<-read.fit('min')

idx<-1:4
coef=fit$est[idx]
sd=fit$std[idx]
ci.q=coef+sd%o%c(-2,2)
obj=fit$nlogl
maxgrad=fit$maxgrad
ev<-eigen(solve(fit$cov[idx,idx]))$value
eratio <- min(ev)/max(ev)

results <- list(obj=obj,                     ## negative log L
                coef=coef,                   ## coefficients (MLE/mode)
                coef.median=NULL,
                sd=sd,                       ## standard deviations 
                confint.quad=ci.q,           ## Nx2 - matrix 95% ci 
                confint.profile=NULL,        
                confint.mcmc=NULL,           ## 0.025 0.975 quantiles
                time=c(fit=time,profile=NULL),  
                terminfo=c(maxgrad=maxgrad,eratio=eratio), ## MLE termination info
                convinfo=NULL,
		fit=fit
                )

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

rmlist<-c("admodel.cov","admodel.hes","admodel.dep","fmin.log",
          "eigv.rpt","variance",paste('min',c(".bar",".cpp",
          ".eva",".log", ".htp",".o","",".cor",".par",".std",
          ".dat"),sep=""))
file.remove(rmlist)