Hi I have no experience with lmer and its plotting method. However If it uses plain (not grid) graphics you maybe shall consult abline and/or segments.
If it uses grid, you shall consult panel.abline from lattice package. BTW. Better to copy your answer every time to the list as somebody from BigBoys can definitely be able to answer questions with more intuition than myself. HTH Petr On 24 Aug 2006 at 16:34, Simon Pickett wrote: Date sent: Thu, 24 Aug 2006 16:34:57 +0100 (BST) Subject: Re: [R] help: trouble using lines() From: "Simon Pickett" <[EMAIL PROTECTED]> To: "Petr Pikal" <[EMAIL PROTECTED]> > Hi, > I have a model with a three way interaction > f2<- lmer(late.growth ~ mtf+year+treat+hatch.day+ hatch.day:year+ > hatch.day:treat+ mtf:treat+ treat:year+ year:treat:mtf+(1|fybrood), > data = A) (one continuous and two factors) and so I wanted to plot the > fitted values of the model for each of the two factor levels. Copying > the technique of the pdf I mentioned "Linear mixed models in R" by > Sřren Feodor Nielsen 20003. > http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-mixed-mod > els.pdf#search=%22Linear%20mixed%20models%20in%20R%22. I plotted the > fitted values against the continuous variable. Now I plan to put lines > on the existing plot for each of the combinations of the factors I > have. where tmp[1] and tmp[2] are coefficients from the model... I was > under the impression that lines(y,exp(tmp[1]+tmp[2])) would give me a > line, but it doesnt work... thanks again > > > > Hi > > > > from lines help page > > > > x, y coordinate vectors of points to join. > > > > and lines or points simply adds lines or points to existing plot. > > What do you want to plot with lines? > > > > HTH > > Petr > > > > > > > > On 24 Aug 2006 at 14:52, Simon Pickett wrote: > > > > Date sent: Thu, 24 Aug 2006 14:52:09 +0100 (BST) > > From: "Simon Pickett" <[EMAIL PROTECTED]> > > To: "Petr Pikal" <[EMAIL PROTECTED]> > > Copies to: [email protected], Simon Pickett > > <[EMAIL PROTECTED]> > > Subject: Re: [R] help: trouble using lines() > > > >> Hi, thanks for replying. > >> No, there arent any NA's in the original data set.... > >> I think I must be mis-interpreting the use of lines()? > >> in the example what exactly is "y"? > >> lines(y,exp(tmp[1]+tmp[2])) > >> In my case tmp[1] and tmp[2] are coeficients from the model so just > >> one number (not a vector) and I havent specified y.... Thanks > >> everyone, Simon > >> > >> > >> > Hi > >> > > >> > not sure but are there some NA values in your data? > >> > > >> > what > >> > > >> > length(mtf) > >> > and > >> > length(fitted(f2)) > >> > > >> > tells you? > >> > > >> > And you need not to use assignment > >> > > >> > graph1 <- plot(....) > >> > > >> > to output a plot on screen. > >> > > >> > > >> > HTH > >> > Petr > >> > > >> > > >> > On 24 Aug 2006 at 13:43, Simon Pickett wrote: > >> > > >> > Date sent: Thu, 24 Aug 2006 13:43:18 +0100 (BST) > >> > From: "Simon Pickett" <[EMAIL PROTECTED]> > >> > To: [email protected] > >> > Subject: [R] help: trouble using lines() > >> > > >> >> Hi R experts, > >> >> I have been using ReML as follows... > >> >> model<-lmer(late.growth~mtf+year+treat+hatch.day+hatch.day:year+ > >> >> hat ch. day:treat+ mtf:treat+ treat:year+ > >> >> year:treat:mtf+(1|fybrood), data > A) then I wanted to plot the > >> >> results of the three way interaction using lines() as follows... > >> >> > >> >> tmp<-as.vector(fixef(model)) > >> >> graph1<-plot(mtf,fitted(f2), xlab=list("Brood Size"), > >> >> ylab=list("Early growth rate"), pch=16, col="darkgrey", > >> >> bg="yellow") lines(y,exp(tmp[1]+tmp[2])) > >> >> > >> >> but no matter what I try I always get the error message > >> >> "Error in xy.coords(x, y) : 'x' and 'y' lengths differ" > >> >> > >> >> Can anyone shed some light please? > >> >> I am basically copying the methods of the pdf entitled "Linear > >> >> mixed models in R" by Sřren Feodor Nielsen 20003. > >> >> http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-mix > >> >> ed- mod els.pdf#search=%22Linear%20mixed%20models%20in%20R%22 > >> >> > >> >> > >> >> > >> >> > >> >> > >> >> Simon Pickett > >> >> PhD student > >> >> Centre For Ecology and Conservation > >> >> Tremough Campus > >> >> University of Exeter in Cornwall > >> >> TR109EZ > >> >> Tel 01326371852 > >> >> > >> >> ______________________________________________ > >> >> [email protected] mailing list > >> >> https://stat.ethz.ch/mailman/listinfo/r-help > >> >> PLEASE do read the posting guide > >> >> http://www.R-project.org/posting-guide.html and provide > >> >> commented, minimal, self-contained, reproducible code. > >> > > >> > Petr Pikal > >> > [EMAIL PROTECTED] > >> > > >> > > >> > >> > >> Simon Pickett > >> PhD student > >> Centre For Ecology and Conservation > >> Tremough Campus > >> University of Exeter in Cornwall > >> TR109EZ > >> Tel 01326371852 > >> > >> ______________________________________________ > >> [email protected] mailing list > >> https://stat.ethz.ch/mailman/listinfo/r-help > >> PLEASE do read the posting guide > >> http://www.R-project.org/posting-guide.html and provide commented, > >> minimal, self-contained, reproducible code. > > > > Petr Pikal > > [EMAIL PROTECTED] > > > > > > > Simon Pickett > PhD student > Centre For Ecology and Conservation > Tremough Campus > University of Exeter in Cornwall > TR109EZ > Tel 01326371852 > Petr Pikal [EMAIL PROTECTED] ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
