R help - Apr 2005 - How to get predictions, plots, etc. from lmer{lme4}

Kindly send a cc to me when replying to the list.

I'm having trouble using lmer beyond a first step.

My data:
 > some(exp1B)
     sub   ba amplitude   a   b  c  d
2     1 1.00       1.5  65  63  4  8
41    4 1.15       0.0  92  41  3  4
43    4 1.15       3.0  88  48  2  2
63    6 1.00       3.0  50  72  9  9
77    8 1.15       0.0 112  25  2  1
89   10 1.15       0.0  37  33 36 34
126  13 1.15       1.5  80  50  6  4
140  14 1.15       4.5  12 115  6  7
145  20 1.00       0.0  73  65  0  2
147  20 1.00       3.0  63  72  2  3
etc.

The output:
 > summary(exp1B.both.cont.lmer)
Linear mixed-effects model fit by maximum likelihood
Formula: cbind(b, a) ~ ba + amplitude + (1 | sub)
    Data: exp1B
       AIC      BIC    logLik MLdeviance REMLdeviance
  768.0086 783.2579 -379.0043   758.0086     766.3104
Random effects:
  Groups   Name        Variance Std.Dev.
  sub      (Intercept) 0.18787  0.43344
  Residual             6.36355  2.52261
# of obs: 156, groups: sub, 13

Fixed effects:
               Estimate Std. Error  DF  t value  Pr(>|t|)
(Intercept)   4.118806   0.397780 153  10.3545 < 2.2e-16
ba           -4.205518   0.330010 153 -12.7436 < 2.2e-16
amplitude     0.137958   0.023754 153   5.8076 3.547e-08

This is just what I need. But I also need predicted values, plots, etc, 
and can't figure out how to proceed. Have I overlooked a more extended 
document than the rather terse (for me at least) help page?
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: 	P.O.Box 400400	Charlottesville, VA 22904-4400
Parcels:	Room 102		Gilmer Hall
		McCormick Road	Charlottesville, VA 22903
Office:	B011	+1-434-982-4729
Lab:		B019	+1-434-982-4751
Fax:		+1-434-982-4766
WWW:	http://www.people.virginia.edu/~mk9y/

Michael Kubovy wrote:
> Kindly send a cc to me when replying to the list.
> 
> I'm having trouble using lmer beyond a first step.
> 
> My data:
>  > some(exp1B)
>     sub   ba amplitude   a   b  c  d
> 2     1 1.00       1.5  65  63  4  8
> 41    4 1.15       0.0  92  41  3  4
> 43    4 1.15       3.0  88  48  2  2
> 63    6 1.00       3.0  50  72  9  9
> 77    8 1.15       0.0 112  25  2  1
> 89   10 1.15       0.0  37  33 36 34
> 126  13 1.15       1.5  80  50  6  4
> 140  14 1.15       4.5  12 115  6  7
> 145  20 1.00       0.0  73  65  0  2
> 147  20 1.00       3.0  63  72  2  3
> etc.
> 
> The output:
>  > summary(exp1B.both.cont.lmer)
> Linear mixed-effects model fit by maximum likelihood
> Formula: cbind(b, a) ~ ba + amplitude + (1 | sub)
>    Data: exp1B
>       AIC      BIC    logLik MLdeviance REMLdeviance
>  768.0086 783.2579 -379.0043   758.0086     766.3104
> Random effects:
>  Groups   Name        Variance Std.Dev.
>  sub      (Intercept) 0.18787  0.43344
>  Residual             6.36355  2.52261
> # of obs: 156, groups: sub, 13
> 
> Fixed effects:
>               Estimate Std. Error  DF  t value  Pr(>|t|)
> (Intercept)   4.118806   0.397780 153  10.3545 < 2.2e-16
> ba           -4.205518   0.330010 153 -12.7436 < 2.2e-16
> amplitude     0.137958   0.023754 153   5.8076 3.547e-08
> 
> This is just what I need. But I also need predicted values, plots, etc, 
> and can't figure out how to proceed. Have I overlooked a more extended 
> document than the rather terse (for me at least) help page?
I regret to say no.  Methods for generalized linear mixed models fit by 
lmer are still rather rudimentary.