R help - Oct 2006 - Adding locfit confidence intervals in trelis xyplot

Dear all,
I am trying to include confidence intervals in a xyplot.

This is what I am doing: 
xyplot(x ~ y|z, alpha = 1,band = "global",panel = panel.locfit)
(more specifically, in my case x is a binary response from a logistic
regression model)

The output plot was fine but it did not include the confidence intervals
Anyone knows how to do it? (xYplot did not work either) 

many thanks
- juan

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

On 10/18/06, juan f poyatos <jpoyatos at cnio.es>
wrote:
> Dear all,
> I am trying to include confidence intervals in a xyplot.
>
> This is what I am doing:
> xyplot(x ~ y|z, alpha = 1,band = "global",panel = panel.locfit)
> (more specifically, in my case x is a binary response from a logistic
> regression model)
>
> The output plot was fine but it did not include the confidence intervals
> Anyone knows how to do it? (xYplot did not work either)
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.

-Deepayan

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