Rで解析：なんとも挑戦的な箱ヒゲ図です。「boxplotdbl」パッケージ

公開日 2022-04-19

Rの解析に役に立つ記事

コマンドを実行して出力したプロットを見ると「なんじゃこりゃ」と思わず呟きそうになる2軸のboxplotです。かなり革新的で挑戦的なプロットです。

理解するのに戸惑うプロットですが、立体的に考えると「なるほど」と思う表現です。ただ、本プロットを用いて第３者への説明はかなり注意かつ慎重になる必要があります。プロットをそのまま、平面で考えると戸惑います。

実行コマンドでは各データの要約を掲載しています。是非ともプロットと見比べて謎解きをしていただければと思います。

バージョンは1.3.0。windows 11のR version 4.1.3で確認しています。22.04.19:内容更新

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

R入門毎週ちょっと、統計解析 | 山村吉信 | 数学 | Kindleストア

Amazonで山村吉信のR入門毎週ちょっと、統計解析。アマゾンならポイント還元本が多数。一度購入いただいた電子書籍は、KindleおよびFire端末、スマートフォンやタブレットなど、様々な端末でもお楽しみいただけます。

www.amazon.co.jp

パッケージのインストール

下記コマンドを実行してください。

#パッケージのインストール
install.packages("boxplotdbl")

実行コマンド

詳細はコメント、パッケージヘルプを確認してください。

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

###データ例の作成#####
TestData <- data.frame(Group = c(rep("Group1", 7),
                                 rep("Group2", 3),
                                 rep("Group3", 3)),
                       X_Data = c(100, 1:3, 20, 21,
                                 150, 6, 200, 250,
                                 6, 200, 250),
                       Y_Data = c(1:7, 8:10, 7:9))
########

#データの要約を表示
#X_Data
tapply(TestData[, 2], TestData[, 1], summary)
#$Group1
#Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#1.00    2.50   20.00   42.43   60.50  150.00 
#$Group2
#Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#6     103     200     152     225     250 
#$Group3
#Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#6     103     200     152     225     250 

#Y_Data
tapply(TestData[, 3], TestData[, 1], summary)
#$Group1
#Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#1.0     2.5     4.0     4.0     5.5     7.0 
#$Group2
#Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#8.0     8.5     9.0     9.0     9.5    10.0 
#$Group3
#Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
#7.0     7.5     8.0     8.0     8.5     9.0 

#データのプロット:boxplotdouコマンド
#BoxPlot色を指定:colオプション;c("No1_color", "No2_color")
#BoxPlot枠色を指定:boxcolオプション
#BoxPlotの塗/斜線の指定:blackwhiteオプション:FALSE/TRUE
#斜線の細かさを指定:shadingオプション;c("No1", "No2")
#X,Yの最大最少の枠を描写:boxed.whiskersオプション:TRUE/FALSE
#近いBoxPlotのGroupをまとめてプロット:condenseオプション:TRUE/FALSE
boxplotdou(x = TestData[c(1, 2)], y = TestData[c(1, 3)],
           col = c("red", "blue", "green"), boxcol = "black",
           blackwhite = TRUE, shading = c(10, 20),
           boxed.whiskers = TRUE, condense = TRUE,
           outpch = 17, outcex = 2, outcol = "red", plot = TRUE)