Displaying 3 results from an estimated 3 matches for "devlaplace".
2007 Jan 26
0
R crash with modified lmer code
...;- .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("lmer", mer, frame = if (model) mf else
data.frame(),
terms = mt, call = match.call()), weights = weights,
family = family))
}
fixInd <- seq(ncol(X))
cons...
2007 May 14
1
parsing an lmer error with interaction term
...l
STP560JG118S20 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_...
2007 Jan 23
0
New lmer: How to recode random effect ?
...vance.
PS: 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...