?gls in package nlme. It's like lme but with no random effects. But you can still model the variance-covariance properties of the data.
Simon. On Tue, 2007-06-05 at 19:11 -0700, [EMAIL PROTECTED] wrote: > Can lme or lmer fit a plain regular fixed effects anova? Ie a model without a > random effect, or have there be at least one random effect in order for these > functions to work? > > Trying to run such, (1) without specifying a random effect produces an error, > (2) specifying that there is no random effect does not produce the same > output > as an anova run in lm(); (2b) specifying that there is no random effect in > lmer > crashed R (division by zero, I think). > > Just trying to see the connection of fixed and random effects anova in R. > STATA > gives me same results for both models up to the point where they differ. > > Best Toby > > > > > > dt1 = > as.data.frame(cbind(c(28,35,27,21,21,36,25,18,26,38,27,17,16,25,22,18),c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4))) > > summary(a1 <- lm(V1~factor(V2)-1, dt1)) > anova(a1) > > summary(a1 <- lm(V1~factor(V2), dt1)) > anova(a1) > > dt1$f = factor(dt1$V2) > > summary(a2 <- lme(V1~f, dt1)) #1a > > summary(a2 <- lme(V1~f, dt1, ~-1|f)) #2a > anova(a2) > > lmer(V1~f, dt1) #1b > > lmer(V1~f+(-1|f), dt1) #2b > > ______________________________________________ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- Simon Blomberg, BSc (Hons), PhD, MAppStat. Lecturer and Consultant Statistician Faculty of Biological and Chemical Sciences The University of Queensland St. Lucia Queensland 4072 Australia Room 320, Goddard Building (8) T: +61 7 3365 2506 email: S.Blomberg1_at_uq.edu.au The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. - John Tukey. ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.