Skip to content
Prev 874 / 7420 Next

Plot a polynomial regression

Hi,
Doing that with xyplot is a bit of a PITA, but
try:
fits<-fitted(model)
xyplot((GPP)~(Discharge), groups=Site, data=nighttimeds1.dat, fits=fits,
       panel=function(x,y,fits,...){
             panel.xyplot(x,y,...)
             panel.lines(x,fits)
        },col=c("red","blue","green","yellow"), type="p")

by the way, polynomials may have some properties that maybe undesirable
in your 
model. Try 
library(splines)
model<-lm(GPP ~ poly(Discharge, 2), data=nighttimeds1.dat)
model2<-lm(GPP ~ ns(Discharge, k=3), data=nighttimeds1.dat)
anova(model, model2) ## likelihood ratio test

Nicholas