R - 将 SD 区域作为虚线添加到正态分布

Question

我有以下 R 脚本，它生成以下图像。我想将 SD 区域添加到这条正常曲线(正负 3 SD)。换句话说，两条线表示 68% 的区域，依此类推。您可以看到我正在谈论的行的示例 here .

g = c(1,2,3,4,5,6,7,8,9,2,3,3,3,4,4,4,5,5,6,6,6,6,6,7,8)
h<-hist(g, breaks=10, density=10, col="lightgray", xlab="Accuracy", main="Overall") 
xfit<-seq(min(g),max(g),length=40) 
yfit<-dnorm(xfit,mean=mean(g),sd=sd(g)) 
yfit <- yfit*diff(h$mids[1:2])*length(g) 
lines(xfit, yfit, col="black", lwd=2)

Answer 1

您需要使用 ?segments .试试这个:

x.vals <- c(mean(g)-(1:3*sd(g)), mean(g), mean(g)+(1:3*sd(g)))
y.vals <- dnorm(x.vals, mean=mean(g),sd=sd(g))*diff(h$mids[1:2])*length(g)
segments(x0=x.vals, y0=0, x1=x.vals, y1=y.vals)

请注意，并非所有 SD 都显示在图中。要使原始图更宽，您可以使用 xlim=c(-3,13)在调用 hist() (并且您还需要相应地更改 xfit 的范围)。