如何遮蔽曲线下的区域？

Question

我正在尝试创建一个与此类似的图表： http://stats.seandolinar.com/wp-content/uploads/2014/12/Probabiliy-of-Finding-Someone-Taller-than-6-Comparison.png

为此，我正在关注此代码示例：

     curve(dnorm(x,0,1),xlim=c(-3,3),main='Normal Density')
cord.x <- c(-2, seq(-2,-1,.01),-1)
cord.y <- c(0, dnorm(seq(-2,-1,.01)),0)
polygon(cord.x,cord.y,col='skyblue')

然后我尝试修改我自己的数据代码，试图将79和120之间的区域遮蔽：

lbound <- (79-94)/6.8
ubound <- (120-94)/6.8
curve(dnorm(x,94,6.8), xlim=c(70,120), main="normal density")
cord.x <- c(79,seq(79,120,.01),120)
cord.y <- c(0, dnorm(seq(79,120,.01)),0)
polygon(cord.x,cord.y,col='skyblue')

我不确定出了什么问题，我试图用z分数和实际数据替换所有数据点，但没有人正在绘制图表。任何知道如何做的人将不胜感激。

Answer 1

您需要在第二个mean函数中定义sd和dnorm。

curve(dnorm(x, mean=94, sd=6.8), xlim=c(70,120), main="normal density")
cord.x <- c(79, seq(79,120,.01), 120)
cord.y <- c(0, dnorm(seq(79,120,.01), mean=94, sd=6.8), 0)
polygon(cord.x, cord.y, col='skyblue')