曲线下的阴影区

Question

我正试图在R中弯曲一个曲线下的区域。我不能完全正确，我不知道为什么。该曲线由

定义

# Define the Mean and Stdev
mean=1152
sd=84

# Create x and y to be plotted
# x is a sequence of numbers shifted to the mean with the width of sd.  
# The sequence x includes enough values to show +/-3.5 standard deviations in the data set.
# y is a normal distribution for x
x <- seq(-3.5,3.5,length=100)*sd + mean
y <- dnorm(x,mean,sd)

情节是

# Plot x vs. y as a line graph
plot(x, y, type="l")

我用来尝试在x> = 1250

的曲线下着色的代码

polygon(c( x[x>=1250], max(x) ),  c(y[x==max(x)], y[x>=1250] ), col="red")

但这是我得到的结果 如何正确地为x> = 1250

的曲线下方的部分着色

Answer 1

您需要跟随多边形的曲线的x，y点，然后返回到x = 1250，y = 0的点以完成形状。最终的垂直边缘是自动绘制的，因为多边形通过返回其起点来关闭形状。

polygon(c( x[x>=1250], 1250 ),  c(y[x>=1250],0 ), col="red")