search for: randomslop

...eaction ~ Days + (1| Subject), sleepstudy), lmer(Reaction ~ Days + (Days | Subject), sleepstudy)) # # models is an unnamed list, do.call works (although with warning): do.call(anova, models) # # after labeling the models, do.call gives an error: names(models) <- c("randomIC", "randomSlope") do.call(anova, models) # # without do.call, anova works fine: anova(models[[1]], models[[2]]) Is there a possibility to use do.call(anova, models) on named lmer-models? Regards, Denes

Dear lmer-ers, My thanks for all of you who are sharing your trials and tribulations publicly. I was hoping to elicit some feedback on my thoughts on denominator degrees of freedom for F ratios in mixed models. These thoughts and practices result from my reading of previous postings by Doug Bates and others. - I start by assuming that the appropriate denominator degrees lies between n

...ply greatful. My R script: library (nlme) #Datei laden randomInterceptDIQAGQ <- lme(NoteD ~ IQ + AGQ, data = Gind, random = ~1|Klnr, method = "ML", na.action = na.exclude) summary (randomInterceptDIQAGQ) intervals (randomInterceptDIQAGQ) my Output: Final model, : 2 predictors, no RandomSlope > randomInterceptDIQAGQ <- lme(NoteD ~ IQ + AGQ, data = Gind, random = ~1|Klnr, method = "ML", na.action = na.exclude) > summary (randomInterceptDIQAGQ) Linear mixed-effects model fit by maximum likelihood Data: Gind AIC BIC logLik 943.9653 964.7289 -466....