[R] Help with plot(augPred()) plot(comparePred()) in nlme

[Michael Kubovy]

Dear r-helpers,

When I issue the commands
plot(augPred(mcc.lme, primary = ~ mcc$age, length.out = 2, level = c 
(0, 1)))
or
plot(comparePred(mcc1.lis, mcc.lme, primary = ~ mcc$age, length.out =  
2), layout = c(4, 3)))
no lines are drawn.

But
fm1 <- lme(Orthodont)
plot(augPred(fm1, level = 0:1, length.out = 2))
does plot lines.

augPred() and comparePredict() seem to be working correctly, as shown  
below

augPred(mcc.lme, primary = ~ mcc$age, length.out = 2, level = c(0, 1))
   mcc$age .groups score    .type
1      30       1   108 original
2      30       2   103 original
3      30       3    96 original
4      30       4    84 original
5      30       5   118 original
6      30       6   110 original
...
91      30      10   101 predict.subject
92      48      10   101 predict.subject
93      30      11   109 predict.subject
94      48      11   109 predict.subject
95      30      12   122 predict.subject
96      48      12   122 predict.subject

comparePred(mcc1.lis, mcc.lme, primary = ~ mcc$age, length.out = 2))
   mcc$age subject score    .type
1      30       1   108 original
2      30       2    96 original
3      30       3   110 original
4      30       4   122 original
5      30       5   103 original
6      30       6   117 original
...
671      30      10   101 mcc.lme
681      48      10   101 mcc.lme
691      30      11   109 mcc.lme
701      48      11   109 mcc.lme
711      30      12   122 mcc.lme
721      48      12   122 mcc.lme


_____________________________
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/