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

統計学備忘録 since2016

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

正規分布の重ね描き

text関数でグラフに文字の挿入

curve(dnorm(x, -2, 4), from=-10, to=10, ylim=c(0,0.4),ylab ="")
text(-5, 0.1, "N(-2,4)")
par(new=T)
curve(dnorm(x,3, 1), from=-10, to=10, ylim=c(0,0.4),ylab ="")
text(1.5, 0.3, "N(3,1)")
par(new=T)
curve(dnorm(x, 3, 4), from=-10, to=10, ylim=c(0,0.4),ylab ="")
text(6, 0.1, "N(3,4)")
dev.off()

f:id:yoshida931:20180419170138p:plain:w600

一部分を塗りつぶし

curve(dnorm(x,3, 1), from=-10, to=10, ylim=c(0,0.4),ylab ="")
# X4~10を100個の多角形(台形)に等分割
xvals <- seq(4, 10, length=100)  
# 対応するグラフの高さ
dvals <- dnorm(xvals,3, 1)      
# 塗りつぶす
polygon(c(xvals,rev(xvals)),c(rep(0,100),rev(dvals)),col="2")     
#x軸に4を追加
axis(side=1,at=c(4))
dev.off()

f:id:yoshida931:20180419171823p:plain:w400