<R> 全ての列のヒストグラムを一括で出力する方法に関して

お世話になっております。R初心者のYukiです。

#全ての列のヒストグラムを一括で出力する方法に関して

お伺いしたいのですが、

df <- data.frame(
  df01 = c(1, 2, 3, 4, 5,1, 2, 3, 4, 5,1, 2, 3, 4, 5),
  df02 = c(200, 203, 193, 193, 187,192, 211, 223, 232, 243, 198, 200, 201, 204, 203),
  df03 = c(4.2, 3.1, 3.2, 3.6, 3.5,4.1, 3.5, 3.6, 4.0, 4.3,4.2, 4.4, 4.3, 4.3, 4.1),
  df04=c(2.3,NA,6.0,5.2,NA,2.8,6.2,4.5,NA,2,3,1,4.2,NA,NA))
df

というデータフレームがあるとして、df01-04すべての列のヒストグラムを”一括で”出したいと考えています。

a<-paste0("df",sprintf("%02d",seq(1,4)))

で列名を取得して、

for (i in a){
    ggplot(df, aes(x = i))+geom_histogram(stat='count',binwidth = 40)
}

のfor ループで出力しようとしたのですがうまくいきません。
https://data-science.gr.jp/implementation/ida_r_histogram.html
などを参考に作成したのですが、進まず困っています。

どなたかご教授いただけないでしょうか。よろしくお願いします。　Yuki

行動規範の内容に同意します

回答1件

ベストアンサー

まず「縦持ち」データを作ります。

Ref:
【tidyr】gather?, spread? もう古い。時代はpivot - Qiita

df <- data.frame(
df01 = c(1, 2, 3, 4, 5,1, 2, 3, 4, 5,1, 2, 3, 4, 5),
df02 = c(200, 203, 193, 193, 187,192, 211, 223, 232, 243, 198, 200, 201, 204, 203),
df03 = c(4.2, 3.1, 3.2, 3.6, 3.5,4.1, 3.5, 3.6, 4.0, 4.3,4.2, 4.4, 4.3, 4.3, 4.1))

library(tidyverse)
df2 <- pivot_longer(df, names_to="Sample", cols=1:ncol(df))

head(df2)
# A tibble: 6 x 2
  Sample value
  <chr>  <dbl>
1 df01     1  
2 df02   200  
3 df03     4.2
4 df01     2  
5 df02   203  
6 df03     3.1

このようにグラフの分割に使う「サンプル名」を独立した列に持たせておき、これをaesに利用します。

ggplot(df2, aes(x=value, fill=Sample)) + geom_histogram(binwidth=40) + scale_fill_manual(values=c("df01"="red", "df02"="orange", "df03"="purple"))  + facet_grid(Sample~.)

facet_grid()でグラフの分割。もちろん一枚のグラフの上に全部をオーバーラップさせることもできます。scale_full_widthは、こういう風に色を変えることも出来ますという例です。
この辺はまあお好みで。

投稿2020/08/25 08:27