require("optimx")
source("../R/lhobbs.R")
source("../R/lhobbslik.R")
source("../R/lhobbslikn.R")
cat("Using analytic derivatives\n")
test<-try(ansR5<-optimx(log(c(2,5,3,1)),lhobbs.lik, lhobbs.lg, control=list(all.methods=TRUE)))
if (class(test) != "try-error") {
    print(ansR5)
} else {
    cat("ML attempt with optimx on scaled likelihood and analytic gradients failed for some method\n")
}

cat("Using numerical derivative approximations\n")
test<-try(ansR5n<-optimx(log(c(2,5,3,1)),lhobbs.lik, control=list(all.methods=TRUE)))
if (class(test) != "try-error") {
    print(ansR5n)
} else {
    cat("ML attempt with optimx on scaled likelihood and numerical gradients failed for some method\n")
}