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

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