[R-sig-ME] lmer nonconvergent: care to run and explain?

Paul Johnson pauljohn32 at gmail.com
Wed Oct 21 17:53:27 CEST 2015


Thanks to everybody for looking at the example. If that code can be
re-used in any way that helps lme4 development, I give permission to
re-use or edit it and put it to use. I'm happy to let everybody who
actually understands this debate it.  I don't (yet).

I need to explain to users why we have these warnings with lmer but
not SAS or Stata.  In the output I pasted in to the original email, it
reports convergence in a few steps of EM.  But lmer is going for a lot
more iterations.  How to explain that difference to students?

I'm reading through the papers that Doug has written in the last 10
years or so explaining the estimation process in PLS.   Bates and
Debroy makes this clear for LMM.  In comparison, the mainstream HLM
folks treated MLM a a GLS problem. Raudenbush & Bryk, for example, or
Snidjers & Bosker, describe calculation of predictions for the b's as
a posterior calculation, rather than an element of the optimization.

It appears to me Stata is written that GLS way.  Stata has a parameter
vector with fixed effects and variances of random effects (Beta,
Sigma).   In contrast, lmer i optimising over (Beta, Sigma, b).

Am I just making up a story here?

pj

On Fri, Oct 16, 2015 at 1:35 PM, Douglas Bates <bates at stat.wisc.edu> wrote:
> For those who may be interested, these are the results of timing the fits of
> two models on these simulated data.  For consistency within the timings I
> have renamed the grouping factor Mind to G and named the three continuous
> covariates as S, T and U.  The optimizers whose names start with LN_ are
> timings from the Julia MixedModels package using the NLopt package for
> optimization.  Those whose names start with NLOPT_LN_ are the same optimizer
> code accessed through the nloptr package for R.  The others are from the
> optimx package, bobyqa from the minqa package (the default for lmer) and the
> build-in Nelder_Mead optimizer, which is generally pretty bad and I can say
> that because I wrote it.
>
> dsname = "paulsim"
> form = Formula: Y ~ 1 + S + T + U + (1 | G) + ((0 + S) | G)
> -2log(likelihood) time(s) feval geval optimizer
>    143232.6341     1.5120   606     0 bobyqa
>    143564.1597     0.2770    70     0 Nelder_Mead
>    143232.9465     0.2680    66     0 NLOPT_LN_BOBYQA
>    143272.7444     0.2430    53     0 NLOPT_LN_COBYLA
>    143803.9823     0.3420    40     0 NLOPT_LN_NELDERMEAD
>    143232.6341     0.4570   147     0 NLOPT_LN_SBPLX
>    143232.6582     0.6320    58     0 optimx:L-BFGS-B
>    143232.6341     0.5480   104     0 optimx:nlminb
>    143232.6341     6.7930    NA     0 optimx:spg
>    143232.6341     1.6930    NA     0 optimx:bobyqa
>    143232.6341     0.0489   107     0 LN_BOBYQA
>    143232.6382     1.9885 69711     0 LN_COBYLA
>    143803.9823     0.0474    56     0 LN_NELDERMEAD
>    143232.6341     0.0527   147     0 LN_SBPLX
> form = Formula: Y ~ 1 + S + T + U + ((0 + S) | G)
> -2log(likelihood) time(s) feval geval optimizer
>    143232.6341     0.1400    41     0 bobyqa
>    143232.6341     0.1510    49     0 Nelder_Mead
>    143232.6343     0.1360    36     0 NLOPT_LN_BOBYQA
>    143232.6503     0.1170    24     0 NLOPT_LN_COBYLA
>    143232.6341     0.1540    48     0 NLOPT_LN_NELDERMEAD
>    143232.6341     0.1900    74     0 NLOPT_LN_SBPLX
>    143232.6368     0.3560    70     0 optimx:L-BFGS-B
>    143232.6341     0.2470    29     0 optimx:nlminb
>    143232.6341     0.3660    NA     0 optimx:spg
>    143232.6341     0.2650    NA     0 optimx:bobyqa
>    143232.6341     0.0240    43     0 LN_BOBYQA
>    143232.6341     0.0240    34     0 LN_COBYLA
>    143232.6341     0.0242    52     0 LN_NELDERMEAD
>    143232.6341     0.0246    81     0 LN_SBPLX
>



-- 
Paul E. Johnson
Professor, Political Science        Director
1541 Lilac Lane, Room 504      Center for Research Methods
University of Kansas                 University of Kansas
http://pj.freefaculty.org              http://crmda.ku.edu



More information about the R-sig-mixed-models mailing list