Rで解析：要約値箱ひげ図のプロット「lvplot」パッケージ

公開日 2022-02-05

Rの解析に役に立つ記事

Descriptionでletter value boxplotと記述があり、要約値箱ひげ図という表現が良いのか悩みます。データ分布を表現する一つの手段で便利な表現方法です。

基本プロットだけでなく「ggplot2」パッケージのコマンドを利用することができます。データ分布の確認に便利です。

パッケージバージョンは0.2.0.9000。実行コマンドはwindows 11のR version 4.1.2で確認しています。

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

Rによるベイジアン動的線形モデル (統計ライブラリー) | G.ペトリス, S.ペトローネ, P.カンパニョーリ, 和合肇, 和合肇, 萩原淳一郎

AmazonでG.ペトリス, S.ペトローネ, P.カンパニョーリ, 和合肇, 和合肇, 萩原淳一郎のRによるベイジアン動的線形モデル (統計ライブラリー)。アマゾンならポイント還元本が多数。

www.amazon.co.jp

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

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

#パッケージのインストール
install.packages("devtools")
devtools::install_github("hadley/lvplot")

実行コマンド

詳細はコマンド内を確認ください。

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

#要約値箱ひげ図のプロット:LVboxplotコマンド
par(mfrow = c(4,2), mar = c(3,3,3,3))
for (i in 1:4) {
  x <- rexp(10 ^ (i + 1))
  boxplot(x, col = "lightblue", horizontal = TRUE)
  title(paste("rexp, n = ", length(x)))
  LVboxplot(x, col = "lightblue", xlab = "")
}

#ggplot2でプロット:geom_lvplotコマンド
#tidyverseパッケージがなければインストール
if(!require("tidyverse", quietly = TRUE)){
  install.packages("tidyverse");require("tidyverse")
}
#カラーパレット作製
#scalesパッケージがなければインストール
if(!require("scales", quietly = TRUE)){
  install.packages("scales");require("scales")
}

###プロットデータの作製#####
PlotData <- NULL
for (i in 1:4) {

  EXPData <- cbind(i, rexp(10 ^ (i + 1)))
  PlotData <- rbind(PlotData, EXPData)

}
PlotData <- as.data.frame(PlotData)
#######
ggplot(PlotData, aes(PlotData[, 1], PlotData[, 2])) +
  geom_lv(aes(fill=..LV..)) +
  scale_fill_manual(values = seq_gradient_pal(c("#e1e6ea", "#505457", "#4b61ba", "#a87963",
                                                "#d9bb9c", "#756c6d", "#807765", "#ad8a80"))(seq(0, 1, length = 14))) +
  facet_wrap(~i)