理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

正規分布の色塗り

plot(dnorm, -4, 4, xaxt="n")
xvals <- seq(-4, -1.96, length=10)        # -4以上-1.33以下 領域をx軸方向に10個の多角形(台形)に等分割
dvals <- dnorm(xvals)    # 対応するグラフの高さ
polygon(c(xvals,rev(xvals)),
        c(rep(0,10),rev(dvals)),col=5)       # 塗りつぶす

xvals2 <- seq(1.96, 4, length=10)        
dvals2 <- dnorm(xvals2)   
polygon(c(xvals2,rev(xvals2)),
        c(rep(0,10),rev(dvals2)),col=5) 

name<-c("-1.96","0","1.96")
axis(side=1,at=c(-1.96,0,1.96),labels=name)  

[f:id:yoshida931:20180719112607p:plain:500]