On Fri, 1 Sep 2006, hadley wickham 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.
I haven't looked at ggplot yet so it might take a while to explore this solution, especially since Deepayan's worked for me. Seems like ggplot is interesting though. Thanks David > > This is very easy to do with ggplot: > > install.packages("ggplot") > library(ggplot) > > p <- ggplot(temp, . ~ LocCode, aes=list(y=Density, x=Year)) > (p <- ggpoint(p)) > ggsmooth(p, method=lm, formula=y~poly(x,3)) > > You can also do it with fitted values in another column. > > temp$fitted <- fitted(paraslm1) > > p <- ggplot(temp, . ~ LocCode, aes=list(y=Density, x=Year)) > (p <- ggpoint(p)) > ggline(p, aes=list(y=fitted)) > > Hadley > _________________________________________________________________ David Scott Visiting (July 06 to January 07) Department of Probability and Statistics The University of Sheffield The Hicks Building Hounsfield Road Sheffield S3 7RH United Kingdom Phone: +44 114 222 3908 Email: [EMAIL PROTECTED] ______________________________________________ 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.