Заштрихованная область под двумя кривыми с использованием R

Я написал следующий код в R

x=seq(-7,10,length=200)
y1=dnorm(x,mean=0,sd=1)
plot(x,y1,type="l",lwd=2,col="red")
y2=dnorm(x,mean=3,sd=2)
lines(x,y2,type="l",lwd=2,col="blue")

Как я могу затенять область под обоими кривыми (известную как перекрытие между двумя кривыми).

Я высоко ценю любые предложения.

Ответ 1

О, ну, @SachaEpskamp избил меня, но вот мое гораздо менее элегантное решение.

shade_under_curve <- function(fun, xmin, xmax, length=100){
  xvals <- seq(xmin, xmax, length=length)
  dvals <- match.fun(fun)(xvals)
  polygon(c(xvals,rev(xvals)),c(rep(0,length),rev(dvals)),col="gray")
}


y1 <- function(x)sapply(x, function(xt)dnorm(xt,mean=0,sd=1))
y2 <- function(x)sapply(x, function(xt)dnorm(xt,mean=3,sd=2))

my.fun <- function(x){sapply(x, function(xt)min(y1(xt), y2(xt)))}

Изменить, чтобы включить начальный график:

plot(y1, -10, 10, col="red")
curve(y2, add=TRUE, col="blue")
shade_under_curve(my.fun, -10, 10, length=1000)

Ответ 2

Добавьте следующую строку:

polygon(x,pmin(y1,y2),col="gray")

Это в основном работает точно так же, как инструмент пера в Photoshop, где первый вектор x - это координаты x, а второй вектор, pmin(y1,y2) - координаты y. pmin дает вектор с минимальными значениями двух векторов по-разному, что соответствует координатам y вершины перекрытия.

EDIT:

Я предпочитаю использовать curve() (как предложил Андри), который можно использовать для построения функции. Вы можете сохранить его координаты при построении графика, и нас это точно так же:

fun1 <- curve(dnorm(x,mean=0,sd=1),type="l",lwd=2,col="red")
fun2 <- curve(dnorm(x,mean=3,sd=2),type="l",lwd=2,col="blue",add=TRUE)
polygon(fun1$x,pmin(fun1$y,fun2$y),col="gray")

Ответ 3

Вот решение, использующее ggplot2

library(ggplot2)
x  = seq(-7, 10, length = 200)
y1 = dnorm(x, mean = 0,sd = 1)
y2 = dnorm(x, mean = 3,sd = 2)

mydf = data.frame(x, y1, y2)

p0 = ggplot(mydf, aes(x = x)) +                         
     geom_line(aes(y = y1), colour = 'blue') +
     geom_line(aes(y = y2), colour = 'red') +
     geom_area(aes(y = pmin(y1, y2)), fill = 'gray60')