[R-sig-ME] Purpose of .nloptr object?
Ben Bolker
bbo|ker @end|ng |rom gm@||@com
Wed Jun 9 03:58:19 CEST 2021
Yes, that's probably it.
I should look into updating that vignette anyway (it would be good to
run the machinery and cache some timings to be included in the text).
cheers
Ben
On 6/8/21 7:47 PM, Phillip Alday wrote:
> This seems to be the docs getting slightly out of date in some sections
> (as happens with any big piece of software). The actual nloptwrap
> function in lme4 (now the default optimizer, IIRC) doesn't include that:
>
> https://github.com/lme4/lme4/blob/c3d0c643d444785fad5e554380c028706f14e274/R/utilities.R#L1151-L1169
>
> A wild guess on my part (until Ben Bolker has a few): that code was
> there to debug writing a wrapper / show how you can further examine the
> return value from nloptr to get more information and got left in, even
> once the further debugging/examination was removed from the example.
>
>
> Phillip
>
> On 08/06/2021 14:19, Michael Chirico via R-sig-mixed-models wrote:
>> I see this code snippet suggested in the performance vignette:
>>
>> https://cran.r-project.org/web/packages/lme4/vignettes/lmerperf.html
>>
>> nlopt <- function(par, fn, lower, upper, control) {
>> *.nloptr <<-* res <- nloptr(par, fn, lb = lower, ub = upper,
>> opts = list(algorithm = "NLOPT_LN_BOBYQA", print_level = 1,
>> maxeval = 1000, xtol_abs = 1e-6, ftol_abs = 1e-6))
>> list(par = res$solution,
>> fval = res$objective,
>> conv = if (res$status > 0) 0 else res$status,
>> message = res$message
>> )
>> }
>>
>> That's writing an .nloptr object to .GlobalEnv (usually), but there's no
>> mention of why and I don't see any usage of this object anywhere else in
>> the code base:
>>
>> grep -Fr ".nloptr" lme4/R
>>
>> Am I missing something?
>> Mike C
>>
>> [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> R-sig-mixed-models using r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
More information about the R-sig-mixed-models
mailing list