R で正規分布の確率密度函数を描いてみよう

できるだけ無駄の少ない PDF 図が欲しいときには、R が良いですね。

R で正規分布の確率密度函数を描いてみましょう。

x <- seq(-4,4, length=100)
plot(dnorm(x)~x, type="l")
abline(v=qnorm(0.025))
abline(v=qnorm(0.975))
abline(h=0)

できあがりは下図です。

qnorm(0.025) で、正規分布時の 97.5 パーセンタイルが求まります。

縦の線2本は、2.5 パーセンタイルと、97.5 パーセンタイルです。

y=0 を追加して描いてあります。