[R] Lattice plot with fitted curves

[David Scott]

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

_________________________________________________________________
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