On 2/10/06, Spencer Graves <[EMAIL PROTECTED]> wrote: > Hi, Doug and Chris: > > I just got the same error message with the "lmList" example in lme4: > > > (fm1 <- lmList(breaks ~ wool | tension, warpbreaks)) > Call: > Error in getResponseFormula(form) : "Form" must be a two sided formula > > sessionInfo() > R version 2.2.1, 2005-12-20, i386-pc-mingw32 > > attached base packages: > [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" > [7] "base" > > other attached packages: > lme4 lattice Matrix > "0.995-2" "0.12-11" "0.995-5"
Hmm. I am unable to reproduce that in either R-devel or R-patched under Linux > library(lme4) Loading required package: Matrix Loading required package: lattice > (fm1 <- lmList(breaks ~ wool | tension, warpbreaks)) Call: lmList(formula = breaks ~ wool | tension, data = warpbreaks) Coefficients: (Intercept) woolB L 44.55556 -16.333333 M 24.00000 4.777778 H 24.55556 -5.777778 Degrees of freedom: 54 total; 48 residual Residual standard error: 10.94028 > sessionInfo() Version 2.3.0 Under development (unstable) (2006-02-11 r37328) i686-pc-linux-gnu attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" [7] "base" other attached packages: lme4 lattice Matrix "0.995-2" "0.13-3" "0.995-5" > library(lme4) Loading required package: Matrix Loading required package: lattice > (fm1 <- lmList(breaks ~ wool | tension, warpbreaks)) Call: lmList(formula = breaks ~ wool | tension, data = warpbreaks) Coefficients: (Intercept) woolB L 44.55556 -16.333333 M 24.00000 4.777778 H 24.55556 -5.777778 Degrees of freedom: 54 total; 48 residual Residual standard error: 10.94028 > sessionInfo() R version 2.2.1, 2005-12-20, i486-pc-linux-gnu attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" [7] "base" other attached packages: lme4 lattice Matrix "0.995-2" "0.12-3" "0.995-5" > > I quit then restarted R and tried the same example in nlme: It > seemed to work fine. > > Chris, have you tried this in nlme? If you absolutely need some > feature of lme4, at least you could do this kind of preliminary work in > nlme, then switch to lme4 (after quitting and restarting R to avoid > potential conflicts between nlme and lme4). > > hope this helps. > spencer graves > > Douglas Bates wrote: > > Please check which packages you have attached when you call lmList. > > That error message looks as if it is coming from the version of lmList > > that is in the nlme package, not the one in lme4. > > > > On 2/6/06, Chris Evans <[EMAIL PROTECTED]> wrote: > > > >>I'm sure I'm being stupid so flame away... > >> > >>R2.2.1 on Windoze (boohoo) latest updates of packages. > >> > >>I'm exploring a dataset (land) with three variables looking at an > >>narrowly unbalanced two group (GROUP) ANCOVA of a randomised > >>controlled trial analysing endpoint score (SFQ.LOCF.ENDPOINT) entering > >>the baseline score (SFQ.BASELINE) as covariate and the following work > >>fine: > >> > >> > >>>res.same <- lm(SFQ.LOCF.ENDPOINT ~ SFQ.BASELINE + GROUP,land) > >>>res.diff <- lm(SFQ.LOCF.ENDPOINT ~ SFQ.BASELINE + GROUP + > >>>SFQ.BASELINE*GROUP,land) > >>>anova(res.same,res.diff) > >> > >>I try: > >> > >> > >>>lmList(SFQ.LOCF.ENDPOINT ~ SFQ.BASELINE | GROUP, land) > >> > >>Call: > >>Error in getResponseFormula(form) : "Form" must be a two sided formula > >> > >>I'm puzzled. That looks like a two sided formula very like the one in > >>the help for lme4 (which had been loaded) and the data look OK: > >> > >> > >>>table(land$SFQ.LOCF.ENDPOINT) > >> > >> 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 23 > >> 1 1 2 4 8 5 16 9 7 14 18 7 16 9 6 8 4 6 2 3 > >> > >>>table(land$SFQ.BASELINE) > >> > >> 3 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 > >> 1 1 3 3 4 11 7 7 10 12 9 16 14 9 8 7 8 6 1 1 > >> > >>>table(land$GROUP) > >> > >> 1 2 > >>87 89 > >> > >>Advice accepted gratefully and flames ruefully! > >> > >>Chris > >> > >>-- > >>Chris Evans <[EMAIL PROTECTED]> > >>Consultant Psychiatrist in Psychotherapy, Rampton Hospital; > >>Research Programmes Director, Nottinghamshire NHS Trust, > >>Hon. Professor of Psychotherapy, Nottingham University, > >>Hon. SL Institute of Psychiatry > >>*** My views are my own and not representative of those institutions *** > >> > >>______________________________________________ > >>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 > >> > > > > > > ______________________________________________ > > 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 > ______________________________________________ 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