Michael Kubovy
2006-Nov-09 14:57 UTC
[R] Help with plot(augPred()) plot(comparePred()) in nlme
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/