ggplotかlatticeで９つのグラフを一度に表示するコードが書けません。

Question

前提・実現したいこと

Rで、複数のヒストグラムのグラフを一度に表示させたいと思っているのですが、ggplot、latticeでうまくできません。

該当のソースコード

# ９列のデータフレーム　DFを使っている。
par(mfrow=c(3,3)
for (i in 1:9){
    ggplot(DF, aes(DF[,i]))+
    geom_histogram()
  }

par(mfrow=c(3,3)
for (i in 1:9){
    histogram(DF[,i])
  }

試したこと

下記のコードでは成功

par(mfrow=c(3,3)
for (i in 1:9) {
hist(DF[,i])
}

補足情報（FW/ツールのバージョンなど）

上の２つだと表示ができません。
列の指定方法が悪いのか。

ggplotやlatticeなら、forを回さなくてもfacetをうまく使えばできそうな気もするのですが、同じデータフレーム（または行列）の中で列単位でグラフを作成するコードが上手く書けません。

お知恵をください。

Answer 1

 ggplot(DF, aes(DF[,i]))

ここが問題ですね。aesではDFの列名を指定する必要があります。DF[,i]では列名になりません。
fiwaさんの例では、x=valueと決め打ちにしています。しかしこれでは汎用性にかけるしforループで処理するのも困難です。列名は変数で指定したいところです。それにはaes()の代わりにaes_()を使い、項目名はas.name()を介して渡してやる必要があります。

以下に、サンプルデータirisを対象として3つのグラフを同時表示する例を書いてみます。

library(ggplot2)
library(gridExtra)
data("iris")
g<-list()
for(i in 1:ncol(iris)){
  g[[i]] <- ggplot(iris, aes_(x=as.name(colnames(iris)[1])))+geom_histogram()
}
grid.arrange(g[[1]], g[[2]], g[[3]], ncol=ncol(iris), nrow=1)

Comments

[valis]

ありがとうございました。
df[,i]では列名の指定にはならないのですね。
for構文を使わないコードを覚えたかったので、fiwa様の方法を参考に解決しました。
ありがとうございました。

Answer 2

複数の図を並べる方法は使用するライブラリによってお作法が違ってきますので、それに沿って記述するのがよろしいかと。latticeの場合はlayout = c(3, 3)のように、ggplot2の場合はfacet_wrap()やfacet_grid()でncol=3のように指定すると描けます。

グラフィックス参考実例集：ラティスグラフィックス - RjpWiki
ggplotのfacet_grid(), facet_wrap()の使い方 | Memo on the Web

# データフレームの準備
mydata <- as.data.frame(matrix(round(rnorm(100*9,50, 10)), 100))
names(mydata) <- sprintf("Sample%01d", 1:9)

# lattice の場合
library(lattice)
lattice.options(default.args = list(as.table = TRUE))
histogram( ~ values | ind, data = stack(mydata) , layout = c(3, 3))

# ggplot2 の場合
library(ggplot2)
ggplot(stack(mydata), aes(x=values)) + geom_histogram() + facet_wrap(~ind, ncol=3)

###lattice

###ggplot2

Answer 3

#lattice

histogram(~ values | ind, data = stack(df), type='count',beraks=seq(0.5,10.5,1),nint=10,
          panel = function(x,...){
            panel.histogram(x,...)
            panel.abline(v=5.5,col='red')
          },
          layout=c(3,3))

#ggplot

ggplot(data=stack(df), aes(x=values)) +
  geom_histogram(breaks = seq(0.5,10.5,1),fill='grey',col='black') +
  geom_vline(xintercept=5.5,col='red') +
  facet_wrap(~ind,ncol=3)