search for: pqlpar

...;binomial", "poisson")) mer at devComp[8] <- -1 mer at status["glmm"] <- as.integer(switch(method, PQL = 1, Laplace = 2, AGQ = 3)) GSpt <- .Call(glmer_init, environment(), fltype) if (cv$usePQL) { .Call(glmer_PQL, GSpt) PQLpars <- c(fixef(mer), .Call(mer_coef, mer, 2)) } else { PQLpars <- c(coef(glmFit), .Call(mer_coef, mer, 2)) } if (method == "PQL") { .Call(glmer_devLaplace, PQLpars, GSpt) .Call(glmer_finalize, GSpt) return(new("glmer", new(&quot...

...18S20 no g int And here is the model: lmer_int_ageXdisc <- lmer(Concord ~ Age*Disc + (1|Subject), family="binomial", data=d, method="Laplace") In the actual data set the AgexDisc interaction should have 16 levels. This produces the error: Error in devLaplace(PQLpars) : Leading minor of order 7 in downdated X'X is not positive definite What does this mean? This error is not thrown for other interactions I specify from my data. Here is my session info: R version 2.4.1 (2006-12-18) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=E...

...l; lmer(a~b + (1|c), family=poisson, data=tab), What I want to do is to see if number of recruits (a) is dependent on the brood sex ratio (b) including brood identity (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 nlm...

...ace: > model4b<-lmer(RESPONSE~ D_TO_FORAL + + I((DIST_GREEN-300)*(DIST_GREEN<300))+ + I((DIST_WATER-200)*(DIST_WATER<200)) + + I((DIST_VILL-900)*(DIST_VILL<900)) + + I((DIST_HOUSE-200)*(DIST_HOUSE<200)) + + (1|TIER), family=binomial, method="Laplace") Fehler in optim(PQLpars, obj, method = "L-BFGS-B", lower = ifelse(const, : non-finite finite-difference value [7] [[alternative HTML version deleted]]

...I tested fm2<-lmer(data=NGud,family=poisson, seed~hab*seedtray +(1|site)+(1|site/hab)) which also made R crashing. And, fm2<-lmer(data=NGud,family=poisson, seed~hab*seedtray +(1|site/hab)) and for this, I obtained : CHOLMOD warning: matrix not positive definite Error in devLaplace(PQLpars) : Cholmod error `matrix not positive definite' at file:../Supernodal/t_cholmod_super_numeric.c, line 614 J?r?me Lema?tre ?tudiant au doctorat Chaire de recherche industrielle CRSNG - Universit? Laval en sylviculture et faune & D?partement de biologie, Facult? des sciences et de g?ni...