Re: [R] Adding locfit confidence intervals in trelis xyplot

From: juan f poyatos <jpoyatos_at_cnio.es>
Date: Thu 19 Oct 2006 - 08:24:29 GMT

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!!

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@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}}



______________________________________________ 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.

Received on Thu Oct 19 18:35:11 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 19 Oct 2006 - 18:30:10 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.