Rで解析：視覚でデータマイニング!密度もわかるbeanplotを作成「beanplotパッケージ」

公開日 2014-09-21

密度がわかるbeanplotを描写できるbeanplotパッケージを紹介します。データの密度を確認せずboxplotで済ませてします方も多いかと思います。一度、beanplotパッケージでデータをプロットしてみてください。きっと新たな発見があるはずです。

＜おすすめのRに関する書籍です＞

基礎から学ぶ統計学実践編〜PCを使った計算とカテゴリカルデータ解析

PCで実際に計算しながら進める実践編。クロス集計などカテゴリカルデータの解析を扱う。

www.amazon.co.jp

「beanplot」パッケージの導入

下記コードを実行することで導入することができます。
参考URL: http://cran.r-project.org/web/packages/beanplot/vignettes/beanplot.pdf

install.packages("beanplot")

beanplotの描写

正規分布と一様分布のデータをbeanplotとboxplotで描写するコードです。意外と図に日本語を使う方が多いので文字化け防止にちょっとした工夫をしています。ご参考までに。

注意：beanplot初期設定では群毎の平均を実践バーで描写します。boxplotと合わせるように下記コードではbeanlinesをmedianに変更しています。

#パッケージの読み込み
library("beanplot")

#データの準備
a <- rnorm(100) #正規分布
b <- runif(100) #一様分布

#図の描写

#日本語を使い方向け
#文字化け防止,MACの方はこのまま。Windowsの方はfamilyをJapan1GothicBBBに変更
#文字の種類は完全に好みです。?parで使用できる文字を確認してください。
par(mfrow = c(1,2), family = "HiraMinProN-W3") 

#注意！beanplot初期設定では群毎の平均を実践バーで描写します。
beanplot(data.frame("正規分布" = a, "一様分布" = b), beanlines = "median",
         col = "red", ylim = c(-3.5, 3.5), main = "BeanPlot")
boxplot(data.frame("正規分布" = a, "一様分布" = b),
        col = "blue", ylim = c(-3.5, 3.5), main = "BoxPlot")