[R] How to get predictions, plots, etc. from lmer{lme4}

From: Michael Kubovy <kubovy_at_virginia.edu>
Date: Sat 16 Apr 2005 - 12:28:47 EST

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/

______________________________________________
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 Apr 17 08:53:04 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:31:13 EST