Hi, first a happy new year ;) I would like to plot the regression lines of a model like lme <- lme(conc ~ name/time - 1, random=conc~time|nr,method="ML",data=measurements) Ideally a seperate plot for each regression line (respectively name) and seperate colors for the measurements (nr). How would you do that? thx Christof
Christof Klu? <ckluss <at> email.uni-kiel.de> writes:> I would like to plot the regression lines of a model like > > lme <- lme(conc ~ name/time - 1, > random=conc~time|nr,method="ML",data=measurements) > > Ideally a seperate plot for each regression line (respectively name) and > seperate colors for the measurements (nr). > > How would you do that? > > thx > Christof >see plot.augPred in the nlme package -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neurosciences 18 avenue du Doyen L?pine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http://www.sbri.fr/members/kenneth-knoblauch.html
Am 02-01-2012 10:54, schrieb ken knoblauch:> Christof Klu?<ckluss<at> email.uni-kiel.de> writes: >> >> lme<- lme(conc ~ name/time - 1, >> random=conc~time|nr,method="ML",data=measurements)> see plot.augPred in the nlme packagethx, but how to set "primary"? I always get the error plot(augPred(lme)) augPred.lme without "primary" can only be used with fits of groupedData objects greetings Christof
Christof Klu? <ckluss <at> email.uni-kiel.de> writes:> Am 02-01-2012 10:54, schrieb ken knoblauch: > > Christof Klu?<ckluss<at> email.uni-kiel.de> writes: > >> lme<- lme(conc ~ name/time - 1, > >> random=conc~time|nr,method="ML",data=measurements) > > see plot.augPred in the nlme package > thx, but how to set "primary"? I always get the error > plot(augPred(lme)) > augPred.lme without "primary" can only be used with fits of > groupedData objects > > greetings > Christofsee ?augPred for explanation of primary -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neurosciences 18 avenue du Doyen L?pine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http://www.sbri.fr/members/kenneth-knoblauch.html