[R] Lattice plot with fitted curves

David Scott d.scott at auckland.ac.nz
Sat Sep 2 00:34:03 CEST 2006

On Fri, 1 Sep 2006, Deepayan Sarkar wrote:

> On 9/1/06, David Scott <d.scott at auckland.ac.nz> 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

Your guess was spot on Deepayan. That gave exactly what I wanted.
I was not even close to this solution.

Many thanks


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