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

Chris Evans stats at psyctc.org
Mon Feb 6 23:29:59 CET 2006


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 at 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 ***




More information about the R-help mailing list