两条曲线之间阴影填充并指定颜色

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我有两条曲线x1，x2，想在两条曲线之间填充阴影，并且x1大于x2时阴影为蓝色，x1小于于x2时阴影为绿色。
有没有大神知道怎么用ggplot做出来呢？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：怎么用 有没有

GAILY2014 发表于 2019-2-26 18:39
我有两条曲线x1，x2，想在两条曲线之间填充阴影，并且x1大于x2时阴影为蓝色，x1小于于x2时阴影为绿色。
有 ...

如果时间精力有限，摸索代码耗时过长，以下“偏方”，仅供参考:<br>
导出图像(.eps)格式，然后在AI界面内，右键解除合并，再进一步针对目标区域进行编辑。

使用geom_ribbon函数

给你个例子参考下吧

1. x=seq(1.1,10,len=100)
   
2. y1=log(x)
   
3. y2=sin(x)+1.4
   
4. 
   
5. plot(x,y1,type = 'l',lwd=4,ylim = c(0,3.3),ylab = 'y')
   
6. lines(x,y2,lwd=4,col=2)
   
7. #find change point of abs(y1/y2)
   
8. end_id = function(id1){
   
9.   if(max(id1)!=length(x)){
   
10.     return(c(max(id1),length(x)))
    
11.   }else{
    
12.     return(length(x))
    
13.   }
    
14. }
    
15. id1=which(abs(y1/y2)<=1)
    
16. id2=which(diff(id1)>1)
    
17. id = c(1,id1[id2:(id2+1)],end_id(id1))
    
18. for (i in 1:(length(id)-1)) {
    
19.   y1x = y1[id[i]:id[i+1]]
    
20.   y2x = y2[id[i]:id[i+1]]
    
21.   xx = x[id[i]:id[i+1]]
    
22.   #判断每个区域中间的点是否满足条件
    
23.   if(all(y1x[2:4]>y2x[2:4])){
    
24.     polygon(c(xx,rev(xx)),c(y2x,rev(y1x)),col = 3)
    
25.   }else{
    
26.     polygon(c(xx,rev(xx)),c(y1x,rev(y2x)),col = 4)
    
27.   }
    
28. }
    
29. legend('topright',c('y1=log(x)','y2=sin(x)+1.4','y1/y2>1','y1/y2<1'),
    
30.        col = 1:4, lty = c(1,1,NA,NA),pch = c(NA,NA,15,15))

复制代码

ggplot版本

1. x=seq(1.1,10,len=100)
   
2. y1=function(x) log(x)
   
3. y2=function(x) sin(x)+1.4
   
4. #找y1,y2的焦点
   
5. 
   
6. q=rootSolve::uniroot.all(function(x) y1(x)-y2(x),range(x))
   
7. p <- y1(q)
   
8. 
   
9. p1=ggplot() +
   
10.   stat_function(aes(x, color = "y1"), fun = y1, size = 1) +
    
11.   stat_function(aes(x, color = "y2"), fun = y2, size = 1)
    
12. 
    
13. z1 <- seq(1.1, q[1], 0.01)
    
14. z2 <- seq(q[1], q[2], 0.01)
    
15. z3 <- seq(q[2], q[3], 0.01)
    
16. z4 <- seq(q[3], 10, 0.01)
    
17. p1 +
    
18.   geom_ribbon(aes(x = z1, ymin = y1(z1), ymax = y2(z1), fill = "y1/y2<1"),
    
19.               alpha = 0.5)+
    
20.   geom_ribbon(aes(x = z2, ymin = y2(z2), ymax = y1(z2), fill = "y1/y2>1"),
    
21.               alpha = 0.5)+
    
22.   geom_ribbon(aes(x = z3, ymin = y1(z3), ymax = y2(z3), fill = "y1/y2<1"),
    
23.               alpha = 0.5)+
    
24.   geom_ribbon(aes(x = z4, ymin = y2(z4), ymax = y1(z4), fill = "y1/y2>1"),
    
25.               alpha = 0.5)+
    
26.   theme_bw() +
    
27.   scale_y_continuous(limits = c(0,3.3))+
    
28.   theme(panel.grid = element_blank(),
    
29.         legend.position = c(0.95, 0.99),
    
30.         legend.justification = c(1, 1),
    
31.         legend.spacing = unit(0, "cm"),
    
32.         legend.margin = margin(0, 0, 0, 0, "cm")) +
    
33.   scale_fill_manual(values = c("blue", "orange"), name = "") +
    
34.   scale_color_manual(values = c("black", "red"), name = "")
    
35.   

复制代码

zhou1_20 发表于 2019-7-14 19:11
ggplot版本

感谢