Hi, I am performing deviance minimization on a biological model of reproductive incompatibility due to symbiotic bacteria, I have up to 15 parameter estimates, but even with 8 parameters. Using simulated data, I observe my nlm procedure do not converge to minimum value, how can I improve and make sure it converges? thanks, St??phane -- St??phane Dupas IRD c/o CNRS Laboratoire Populations G??n??tique et Evolution 1 av de la Terrasse 91198 Gif sur Yvette + 33 1 69 82 37 04 http://www.cnrs-gif.fr/pge/index.html http://www.cnrs-gif.fr/pge/index.php?lang=en
when I tried the example in glmmPQL I got an error
library(nlme)
summary(glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID,
family = binomial, data = bacteria))
iteration 1
iteration 2
iteration 3
iteration 4
iteration 5
iteration 6
Error: No slot of name "reStruct" for this object of class
"lme"
Error in logLik(object at reStruct) : Unable to find the argument
"object" in
selecting a method for function "logLik">
I thinks I understand the porblem, you can not use glmmPQL if you have open lme4 On Wed, 28 Apr 2004, Liliana Forzani wrote:> when I tried the example in glmmPQL I got an error > > library(nlme) > summary(glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID, > family = binomial, data = bacteria)) > > iteration 1 > iteration 2 > iteration 3 > iteration 4 > iteration 5 > iteration 6 > Error: No slot of name "reStruct" for this object of class "lme" > Error in logLik(object at reStruct) : Unable to find the argument "object" in > selecting a method for function "logLik" > > > > ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://www.stat.math.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html >