tdata<-read.table("../DATA/weeds.dat",header=TRUE)
t<-tdata$t
y0<-tdata$y
rm(tdata) # cleanup

## from source("lhobbslikrun.R")

require("optimx")
sessionInfo()
source("../R/lhobbs.R")
source("../R/lhobbslik.R")
source("../R/lhobbslikn.R")
cat("Using analytic derivatives\n")
astart<-log(c(2,5,3,1))
tryf<-lhobbs.lik(astart, y=y0)


test<-try(ansR5<-optimx(astart,lhobbs.lik, lhobbs.lg, control=list(all.methods=TRUE), y=y0))



if (class(test) != "try-error") {
    # construct the list 
    results<-with(ansR5, mapply(function(obj, par, time) {list(obj=obj, coef=par, sd=rep(NA, length(par)), 
             time=time)}, 
             fvalues, par, xtimes,SIMPLIFY=FALSE))
    names(results) <- unlist(ansR5$method)
} else {
    # construct a list of NAs
    results <- list(obj=NA, coef=rep(NA,length(astart)),  sd=rep(NA, length(astart)), 
             time=NA)
}

print(results)

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