4

以下の例では、破線の間の領域をどのようにペイントできますか？

コードは次のとおりです。

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")

これが出力です

ここに画像の説明を入力してください

4

1 に答える 1

3

重複しているようですが、SO検索では見つかりませんでした。rhelpに何度も登場していることは知っています。しかし、それには幸運があります。

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's and Jonathan Baron's 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) )

ここに画像の説明を入力

于 2013-01-14T09:32:28.080 に答える