On 9/1/06, David Scott <[EMAIL PROTECTED]> wrote: > > I have some data which consists of time series for a number of sites. It > appears that there is not much autocorrelation in the data and I have > fitted a cubic for each site using lm. I would like to obtain a lattice > plot with one panel for each site and showing the original data, and the > fitted cubic. > > The closest I have got to doing what I want is: > > fit <- fitted(paraslm1) > temp <- cbind(paras, fit) > xyplot(Density ~ Year | LocCode, data = temp, > panel = function(x, y, ...){ > panel.xyplot(x,y) > panel.xyplot(x[LocCode],fit[LocCode],type="l") > }) > > This doesn't give an error (most of my other attempts did), and draws the > panels correctly with the original data, but doesn't draw the fitted > lines. > > paraslm1 is my fitted linear model with cubics for each location.
I'm guessing since I don't have actual data to work with, but the following is probably what you are trying to do: xyplot(Density + fit ~ Year | LocCode, data = temp, panel = panel.superpose.2, type = c('p', 'l')) -Deepayan ______________________________________________ R-help@stat.math.ethz.ch 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.