如何填充R中某些特定区域的颜色？

Question

这里的问题是：

x<-seq(0,10,length.out = 1000) 
y1<-dnorm(x,mean = 2,sd=1) 
y2<-dnorm(x,mean = 6,sd=1) 
plot(x,y1,type="l") 
lines(x,y2) 
abline(v=x[380]) 

图表如下所示。如何在垂直线的两侧填充2种不同的颜色，例如红色和蓝色，但仍低于两个正常密度函数。我以为我可以使用多边形，但失败了。

这是无填充颜色的图形：

Answer 1

这里有一种方法：

首先，我们会得到你的密度的并行最低 - 这是对一个向量顶部y我们的多边形坐标。

y = pmin(y1, y2) 

# set up your plot as in the question  
plot(x, y1, type="l") 
lines(x, y2) 

# define a re-usable variable for the vertical line placement 
x_vert = 380 
abline(v = x[x_vert]) 

# Now we'll draw 2 polygons, one for the left side, one for the right. 
# The first (x,y) pairs of the polygon are just the (x,y) coords of the 
# density we're filling to, until the vertical line 
# Then we need to connect the "bottom" points, which have coordinates 
# (x[x_vert], 0) and (x[1], 0)  
polygon(x = c(x[1:x_vert], x[x_vert], x[1]), 
     y = c(y[1:x_vert], 0, 0), 
     col = "blue") 
# similar for the right hand polygon, but now going from x_vert to length(x) 
polygon(x = c(x[x_vert:length(x)], x[length(x)], x[x_vert]), 
     y = c(y[x_vert:length(x)], 0, 0), 
     col = "red") 

瞧！