R Tips：ヒストグラム

Rでヒストグラムの使い方は、色々なサイトに載っていますが、最近知った面白い使い方を紹介。

# サンプルデータを作成
# 標準正規分布に従うデータを1000レコード作成
dat <- rnorm(1000, mean=0, sd=1.0)
　
# ヒストグラムを作成
hist(dat)

ここで、hist(dat)を別の変数 "x" に渡すことが出来るようです。
これを使えば、色々な情報をコマンドで得られますね。

# 変数 "x" にヒストグラムを格納
x <- hist(dat)
　
# ヒストグラムの情報を確認
x
　
# 例えば、ヒストグラムの密度を取り出す場合
x$density
　
# ヒストグラムの密度の最大値の位置
which.max(x$density)
　
# ヒストグラムの密度の最大値の値
x$density[which.max(x$density)]