[R] R equivalent to `estimate' in SAS proc mixed

Dieter Menne dieter.menne at menne-biomed.de
Fri Aug 19 08:54:40 CEST 2005


Randy Johnson <rjohnson <at> ncifcrf.gov> writes:

>     > model <- lmer(response ~ time * trt * bio + (time|id), data = dat)
..
> 
> and I would like to obtain an estimate (with standard error) of the change
> in response over time for individuals in the treatment group with the
> biological factor.

Greg Warnes' gmodels package has "estimable", which works for lme (but probably
not for lmer). You model should work the same with lme with slight syntax
changes, or you could ask Greg to make estimable lmer-aware.

Dieter




More information about the R-help mailing list