[R] Adding locfit confidence intervals in trelis xyplot

[juan f poyatos]

But it is *within* a trellis plot when I cannot plot the bands!

Building up on Dieter's example:

x <- rnorm(100)
y <- dnorm(x) + rnorm(100) / 5
plot(locfit(y~x), band="global")

Plots two nice confidence bands. 

But try this
z <- dnorm(y) + rnorm(100) / 3
Z <- equal.count(z, number = 4, overlap = .1)
xyplot(z ~ y|Z, alpha = 1,band = "global",panel = panel.locfit)

and you will not see the bands,

indeed band is completely ignored as
xyplot(z ~ y|Z, alpha = 1,band = "rubbish",panel = panel.locfit)

works fine but without the bands!!

- juan


On Thu, 2006-10-19 at 07:49 +0000, Dieter Menne wrote:
> Deepayan Sarkar <deepayan.sarkar <at> gmail.com> writes:
> 
> > 
> > On 10/18/06, juan f poyatos <jpoyatos <at> cnio.es> wrote:
> > > Dear all,
> > > I am trying to include confidence intervals in a xyplot.
> > >
> ...
>  
> > Well, panel.locfit doesn't have any options to draw confidence bands,
> > so you'll have to write a panel function that does. Shouldn't be hard
> > to extend panel.locfit if you know how to extract that information
> > from a locfit object.
> 
> plot.locfit has bands. I think Juan mixed trellis and standard plot 
> documentation.
> 
> 
> x <- rnorm(100)
> y <- dnorm(x) + rnorm(100) / 5
> plot(locfit(y~x), band="global")
> 
> ______________________________________________
> R-help at 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.
-- 
Juan F. Poyatos
Structural and Computational Biology Programme
Spanish National Cancer Centre (CNIO)
Melchor Fernandez Almagro, 3/E-28029 Madrid SPAIN
Phone:+34 912 246 900/Fax: +34 912 246 980

http://bioinfo.cnio.es/~jpoyatos/




**NOTA DE CONFIDENCIALIDAD** Este correo electrónico, y en s...{{dropped}}