[R] Shading in prediction intervals

[David Winsemius]

On Sep 12, 2012, at 10:38 PM, Jonathan Zhang wrote:

> 
> I have the following code for the minimum and maximum of my prediction interval
> 
>> y.down=lines(x[x.order], set1.pred[,2][x.order], col=109)
>> y.up=lines(x[x.order], set1.pred[,3][x.order], col=109)

You do not understand what the lines function is doing:

> y <- lines(1:3, 1:3)
> y
NULL


> 
> domain=min(x):max(x)
> 
> polygon(c(domain,rev(domain)),c(y.up,rev(y.down)),col=109)
> 
> It doesnt seem to shade the right region, it gives me a trapezoid.
> 
> Any help? Thanks!
> 		 	   		  
> 	[[alternative HTML version deleted]]

David Winsemius, MD
Alameda, CA, USA