[R] Lattice plot with fitted curves
David Scott
d.scott at auckland.ac.nz
Sat Sep 2 00:38:49 CEST 2006
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: d.scott at auckland.ac.nz
More information about the R-help
mailing list