定义法线下阴影的多边形间隔

时间:2018-06-27 20:26:44

标签: r plot

我正在使用一些非常简单的代码来生成以下图像:

mean=1000
sd=200
Score <- round ((seq(-3,3,length=120)*sd + mean), -1)
y <- dnorm(Score,mean,sd)
data <- data.frame
plot(Score, y, type="l")
polygon(c( Score[Score>=1250], 1250 ),  c(y[Score>=1250],0 ), col="skyblue")
percentile <- round(pnorm((Score-mean)/sd)*100,2)

enter image description here

我的问题是多边形偏斜。我知道我的错误在于如何定义多边形,但是我不知道如何将其区域定义为介于1250-1600之间的间隔,而不仅仅是大于或等于1250(即使1600是代码)。

此外,我试图将这段代码转换为ggplot2,但不幸失败了。任何帮助,也将不胜感激。

2 个答案:

答案 0 :(得分:2)

我将0替换为min(y[Score>=1250])。现在看起来不错。

mean=1000
sd=200
Score <- round ((seq(-3,3,length=120)*sd + mean), -1)
y <- dnorm(Score,mean,sd)
plot(Score, y, type="l")
polygon(c( Score[Score>=1250], 1250 ),  c(y[Score>=1250], min(y[Score>=1250])), col="skyblue")
percentile <- round(pnorm((Score-mean)/sd)*100,2)

enter image description here

答案 1 :(得分:2)

我希望您插入一个附加点以显示零以上的正确高度。

certificates

Area under curve