ヒストグラムの階級を細かくして正規分布の密度関数を近似的に描写する


#標準正規分布から乱数を取得
dat <- rnorm(1000000)
#レンジを調べる
range(dat)
png("120812_hist.png")
par(mfrow=c(2,2))
hist(dat, breaks=seq(-10,10,10))
hist(dat, breaks=seq(-10,10,1))
hist(dat, breaks=seq(-10,10,0.1))
hist(dat, breaks=seq(-10,10,0.01))
dev.off()

dat <- rnorm(1000000)
png("120812_hist_2.png")
par(mfrow=c(1,2))
hist(dat,breaks=seq(-6,6,0.1))
curve(length(dat)*dnorm(x))
dev.off()