Delete comment from: One R Tip A Day
Paolo said...
This should fit your need:
p = seq(-8.5,4.9,.001)
x1 <- dnorm(p,-1.5,1.5)
x2 <- dnorm(p,0,1.5)
plot(range(p),range(x1,x2),type="n")
lines(p,x1,col="red",lwd=4)
lines(p,x2,col="blue",lwd=4,lty=2)
line <-- 1.13
abline(v=line)
xt1 <- c(p[1], p[p < line], line)
yt1 <- c(x1[1], x1[p < line],0)
polygon(xt1, yt1, col="grey")
Sep 11, 2008, 3:50:00 PM