search for: lmeopt

...09e-01 4.4632e-01 637 -1.8083 0.07103 . --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 There were 11 warnings (use warnings() to see them) > warnings() Warning messages: 1: optim or nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv) 2: optim or nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv) 3: optim or nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv) 4: optim or nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv)...

...tity (c) as random factor. However, I get the following error message; lmer(a~b + (1|c), family=poisson, data=tab) Error in devAGQ(PQLpars, 1) : Unable to invert singular factor of downdated X'X In addition: Warning messages: 1: optim or nlminb returned message singular convergence (7) in: LMEopt(x = mer, value = cv) 2: optim or nlminb returned message singular convergence (7) in: LMEopt(x = mer, value = cv) 3: optim or nlminb returned message singular convergence (7) in: LMEopt(x = mer, value = cv) 4: optim or nlminb returned message singular convergence (7) in: LMEopt(x = mer, va...

...imate Std. Error t value (Intercept) -0.2431 0.6619 -0.3672 In.Out 1.6004 0.5645 2.8349 Correlation of Fixed Effects: (Intr) In.Out -0.975 There were 30 warnings (use warnings() to see them) > warnings() Warning messages: 1: Estimated variance-covariance for factor ?Replicate? is singular in: LMEopt(x = mer, value = cv) 2: nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv) I have made sure I have the latest lme4 version, any help would be much appreciated. Cheers, Brendan Connors Behavioural Ecology Research Group

...imate Std. Error t value (Intercept) -0.2431 0.6619 -0.3672 In.Out 1.6004 0.5645 2.8349 Correlation of Fixed Effects: (Intr) In.Out -0.975 There were 30 warnings (use warnings() to see them) > warnings() Warning messages: 1: Estimated variance-covariance for factor ?Replicate? is singular in: LMEopt(x = mer, value = cv) 2: nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv) but as Mr. Bates pointed out, this is inappropriate b/c I am trying to use 3 distinct replicates to estimate 3 variance-covariance parameters. "It won't work. Notice that the estimated c...

...s? Here is the function and the warning message: > model.growth.mcas5 <- lmer(response ~ monthElapsed + (monthElapsed|studentID), + data= mcas5, family=binomial(link="logit"), method='ML') Warning messages: 1: nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv) 2: nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv) 3: nlminb returned message false convergence (8) in: LMEopt(x = mer, value = cv) I tried to increase maximum iteration limit by adding control=list(msMaxIter=600, maxIter=300, PQLmaxIt=1...

Hi all! My question is why, and what I can do about that I sometimes, but not always, get warning-messages like nlminb returned message singular convergence (7) in: LMEopt(x = mer, value = cv) or IRLS iterations for PQL did not converge when trying to fit a model looking like this: lmer<-(cbind(Diseased,Healthy)~Fungus+(1|Family)+(1|Fungus:Family), family="binomial") to four similar data sets? All four data sets consists of four columns; Fungus (1...

I am getting the following error message using the lmer function for mixed models with method="Laplace": "nlminb returned message false convergence (8) in: LMEopt(x=mer,value=cv)" Could anyone explain what this means, and how I might overcome (or track down) the problem? Bill Shipley [[alternative HTML version deleted]]

...lace") Both model 1 and 2 run fine. For model 3, I got error message: ---------------------------------- Error in fn(par, ...) : Unable to invert singular factor of downdated X'X In addition: Warning messages: 1: optim or nlminb returned message ERROR: ABNORMAL_TERMINATION_IN_LNSRCH in: LMEopt(x = mer, value = cv) 2: Leading minor of size 1 of downdated X'X is indefinite ---------------------------------- What is going on here? Any workarounds? Thanks! Best, Shige

...1 Item (Intercept) 0.509 0.713 number of obs: 322, groups: Subject, 46; Item, 7 B Subject (Intercept) 3.317 1.821 Item (Intercept) 0.725 0.852 number of obs: 322, groups: Subject, 46; Item, 7 "ORIGINAL" (PQL) A 1: Estimated variance for factor Item is effectively zero in: LMEopt(x = mer, value = cv) 2: Estimated variance for factors Subject, Item is effectively zero in: LMEopt(x = mer, value = cv) B Estimated variance for factors Subject, Item is effectively zero in: LMEopt(x = mer, value = cv) I understand that the between-Item variance is low, and probably it is no...

...)1 -0.334 fctr(gcmp)j -0.417 0.066 factor(sx)m -0.505 -0.002 -0.173 ResCondCorp -0.309 -0.010 0.302 -0.032 ### Here is the error message: ### Warning message: Estimated variance for factor 'factor(cdrgsaou2$ids)' is effectively zero in: LMEopt(x = mer, value = cv) ### Thank you very much by advance for any help. J?r?me Lema?tre Ph.D. student Silviculture-wildlife research chair in irregular boreal forests & D?partment of biology, Faculty of Sciences and Engineering Alexandre-Vachon building University Laval Quebec, QC G1K 7P4...

...fl <- FL$fl if (fltype < 0) { mer <- .Call(mer_create, fl, Zt, X, as.double(Y), method == "REML", nc, cnames) if (!is.null(start)) mer <- setOmega(mer, start) .Call(mer_ECMEsteps, mer, cv$niterEM, cv$EMverbose) LMEoptimize(mer) <- cv return(new("lmer", mer, frame = if (model) fr$mf else data.frame(), terms = mt, call = mc)) } if (method %in% c("ML", "REML")) method <- "Laplace" if (method == "AGQ") stop(&qu...

I am fitting models to the responses to a questionnaire that has seven yes/no questions (Item). For each combination of Subject and Item, the variable Response is coded as 0 or 1. I want to include random effects for both Subject and Item. While I understand that the datasets are fairly small, and there are a lot of invariant subjects, I do not understand something that is happening

Dear R-users, Half a year ago we put out the R package "glmmADMB" for fitting overdispersed count data. http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html Several people who used this package have requested additional features. We now have a new version ready. The major new feature is that glmmADMB allows Bernoulli responses with logistic and probit links. In addition there

...-0.346605 0.026666 -12.9979 < 2e-16 *** sigma 2.608 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Warning message: optim or nlminb returned message ERROR: ABNORMAL_TERMINATION_IN_LNSRCH in: LMEopt(x = mer, value = cv) The R routine correctly identifies the treatment effect as not significant. However the parameter estimates are poor. Likelihood Ratio Testing Accurate calculation of the log-likelihood value is desirable so that hypothesis testing can be carried out using l...