[R] monte carlo simulations/lmer

From: Eduardo Leoni <e.leoni_at_gmail.com>
Date: Sun 14 Aug 2005 - 09:26:05 EST


Hi - I am doing some monte carlo simulations comparing bayesian (using Plummer's jags) and maximum likelihood (using lmer from package lme4 by Bates et al).

I would like to know if there is a way I can flag nonconvergence and exceptions. Currently the simulations just stop and the output reads things like:

Error in optim(.Call("lmer_coef", x, 2, PACKAGE = "Matrix"), fn, gr, method = "L-BFGS-B", :

    L-BFGS-B needs finite values of 'fn' In addition: Warning message:
Leading minor of size 1 of downdated X'X is indefinite

Error in .local(object, ...) : Leading 2 minor of Omega[[1]] not positive definite
In addition: Warning messages:
1: optim or nlminb returned message ERROR: ABNORMAL_TERMINATION_IN_LNSRCH in: "LMEoptimize<-"(`*tmp*`, value = list(maxIter = 200, msMaxIter = 200, 2: optim or nlminb returned message ERROR: ABNORMAL_TERMINATION_IN_LNSRCH in: "LMEoptimize<-"(`*tmp*`, value = list(maxIter = 200, msMaxIter = 200,

thanks for any help.

-eduardo



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sun Aug 14 09:32:06 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:39:48 EST