我正在使用一些非常简单的代码来生成以下图像:
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)
我的问题是多边形偏斜。我知道我的错误在于如何定义多边形,但是我不知道如何将其区域定义为介于1250-1600之间的间隔,而不仅仅是大于或等于1250(即使1600是代码)。
此外,我试图将这段代码转换为ggplot2,但不幸失败了。任何帮助,也将不胜感激。
答案 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)
答案 1 :(得分:2)