Re: [R] lme4: Error in getResponseFormula(form) : "Form" must be a two sided formula

From: Douglas Bates <bates_at_stat.wisc.edu>
Date: Sun 12 Feb 2006 - 03:37:19 EST

On 2/10/06, Spencer Graves <spencer.graves@pdf.com> 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 <stats@psyctc.org> 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 <chris@psyctc.org>
> >>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 Received on Sun Feb 12 03:41:24 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:27 EST