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+hatch. > 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-mixed-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] ______________________________________________ [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.
