control = list(maxit = 1000, trace = 0)
minqa_control <- list(iprint = control$trace, maxfun = 20*control$maxit)
mod <- uobyqa(param, nll, control=minqa_control)