On Fri, 1 Sep 2006, Deepayan Sarkar wrote:
> 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
>
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: [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.
