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

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

「beanplot」パッケージの導入

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

install.packages("beanplot")

beanplotの描写

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

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

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#35501;&#12415;&#36796;&#12415;
library("beanplot")

#&#12487;&#12540;&#12479;&#12398;&#28310;&#20633;
a <- rnorm(100) #&#27491;&#35215;&#20998;&#24067;
b <- runif(100) #&#19968;&#27096;&#20998;&#24067;

#&#22259;&#12398;&#25551;&#20889;

#&#26085;&#26412;&#35486;&#12434;&#20351;&#12356;&#26041;&#21521;&#12369;
#&#25991;&#23383;&#21270;&#12369;&#38450;&#27490;,MAC&#12398;&#26041;&#12399;&#12371;&#12398;&#12414;&#12414;&#12290;Windows&#12398;&#26041;&#12399;family&#12434;Japan1GothicBBB&#12395;&#22793;&#26356;
#&#25991;&#23383;&#12398;&#31278;&#39006;&#12399;&#23436;&#20840;&#12395;&#22909;&#12415;&#12391;&#12377;&#12290;?par&#12391;&#20351;&#29992;&#12391;&#12365;&#12427;&#25991;&#23383;&#12434;&#30906;&#35469;&#12375;&#12390;&#12367;&#12384;&#12373;&#12356;&#12290;
par(mfrow = c(1,2), family = "HiraMinProN-W3") 

#&#27880;&#24847;&#65281;beanplot&#21021;&#26399;&#35373;&#23450;&#12391;&#12399;&#32676;&#27598;&#12398;&#24179;&#22343;&#12434;&#23455;&#36341;&#12496;&#12540;&#12391;&#25551;&#20889;&#12375;&#12414;&#12377;&#12290;
beanplot(data.frame("&#27491;&#35215;&#20998;&#24067;" = a, "&#19968;&#27096;&#20998;&#24067;" = b), beanlines = "median",
         col = "red", ylim = c(-3.5, 3.5), main = "BeanPlot")
boxplot(data.frame("&#27491;&#35215;&#20998;&#24067;" = a, "&#19968;&#27096;&#20998;&#24067;" = b),
        col = "blue", ylim = c(-3.5, 3.5), main = "BoxPlot")

出力例

少しでも、あなたの解析が楽になりますように！！