Border in R chart

In the example below, how to draw the area between the dashed lines?

Here is the code:

X <- rnorm(100) Y<- X + rnorm(100) reg<-lm(Y~X) CI <- predict(reg,interval="confidence") F <- CI[,1] L <- CI[,2] U <- CI[,3] plot(X,F) abline(reg, col="red") reg2<-lm(L~X) abline(reg2, col="red", lty="dashed") reg3<-lm(U~X) abline(reg3, col="red", lty="dashed") 

displayed here

enter image description here

+4
source share
1 answer

It seems like it should be a duplicate, but I could not find it in the SO search. I know this has appeared many times. But there is a condition for this:

 require(fortunes) fortune("every question") Lastly, by what we could call the 'Liaw-Baron principle', every question that can be asked has in fact already been asked. -- Dirk Eddelbuettel (citing Andy Liaw and Jonathan Baron opinion on unique questions on R-help) R-help (January 2006) polygon(x= c(X[order(X)], rev(X[order(X)]) ), # both x and y go "out" in order; y=c(L[order(X)], rev(U[order(X)])), # and then come "back in reverse order col=rgb(1,0,0,0.3) ) # transparent pink. > plot(X,Y, pch="+") > abline(reg, col="red") > polygon(x= c(X\[order(X)\], rev(X\[order(X)\]) ), y=c(L\[order(X)\],rev(U\[order(X)\])), col=rgb(1,0,0,0.3) ) 

enter image description here

+3
source

All Articles