polygon如何在两个函数曲线间填充 R

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

目标：已有

1. x = runif(n, min=-2, max=2)
   
2. y = x^3 + rnorm(n)

复制代码

1. q<-qnorm(0.9)
   
2. yy1<-x^3+q
   
3. yy2<-x^3-q

复制代码

要在yy1,yy2两个函数之间的区域填上灰色（Confidence band ）

请问如何做？

ps：我的问题代码
plot(x,yy1,col="blue",type = "p")  #,xlim = c(0,2)
points(x,yy2)
#polygon(c(x,rev(x)),c(x^3+q,rev(x^3-q)),col = "gray")
#polygon(c(x),c(x^3+q,rev(x)^3-q),col = "gray")
polygon(c(xx),c(yy1,yy2),col = "gray")

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

以及

1. xx<-rep(rev(x),each=2)

复制代码

你的x需要排一下序，这样应该就可以了。

1. n <- 100
   
2. x = runif(n, min=-2, max=2)
   
3. x <- sort(x)
   
4. y = x^3 + rnorm(n)
   
5. 
   
6. q <- qnorm(0.9)
   
7. yy1 <- x^3+q
   
8. yy2 <- x^3-q
   
9. 
   
10. plot(x,yy1,col="blue",type = "p")  #,xlim = c(0,2)
    
11. points(x,yy2)
    
12. polygon(c(x,rev(x)),c(x^3+q,rev(x^3-q)),col = "gray")

复制代码

你的x向量需要排一下序，这样应该就可以了。

1. n <- 100
   
2. x <- runif(n, min=-2, max=2)
   
3. x <- sort(x)
   
4. y <- x^3 + rnorm(n)
   
5. 
   
6. q <- qnorm(0.9)
   
7. yy1 <- x^3+q
   
8. yy2 <- x^3-q
   
9. 
   
10. plot(x,yy1,col="blue",type = "p")  #,xlim = c(0,2)
    
11. points(x,yy2)
    
12. polygon(c(x,rev(x)),c(x^3+q,rev(x^3-q)),col = "gray")

复制代码

奇怪了，发了好几次都没发出来。你的向量x需要排一下序，增加一行代码就可以了。

1. x <- runif(n, min=-2, max=2)
   
2. x <- sort(x)
   
3. y <- x^3 + rnorm(n)

复制代码