R에서 신뢰구간을 정규분포에 색칠하기

아래는 코드입니다. 95% 신뢰구간입니다. 

#표준정규분포 그래프
x=seq(-4,4,by=0.01)
y=dnorm(x)
plot(x,y,type="l",ann=FALSE,axes=FALSE)

#임계값 설정(기각역)
#99%로 할 경우 0.995 입력
margin=qnorm(0.975)

#축 및 값 추가, pos=c(0,0) 으로 설정해야 축과 그래프 사이 간격 없어짐
#mgp 는 tick 과 tick label 사이의 간격임 mgp=c(3, 2, 0)
#padj=1 이 나은듯. 윗선이 맞춰짐
#99%로 할 경우 2.58로 수정
axis(1,at=c(-4,-margin,0,margin,4),
     pos=c(0,0),
     label=c("",expression(-1.96* frac(sigma, sqrt(n))),expression(mu),expression(1.96* frac(sigma, sqrt(n))),""),
     padj=1)

#색칠을 위한 값 설정
x_margin=seq(-margin,margin,by=0.01)
y_margin=dnorm(x_margin)

#색칠
polygon(c(x_margin,rev(x_margin)),c(y_margin,rep(0,length(y_margin))),col="red")

#테두리
box("outer",col="gray")

 

아래는 결과입니다.