データを概観するためには,記述統計量を求めるとともに,グラフを作成して視覚的にデータの特性を捉えることも重要である。一変数の分布を概観するにはヒストグラムや箱ひげ図,変数間の関係を概観するには散布図を描いてみるのが有効。
R本体でもさまざまなグラフを作成することができるが,ggplot2というパッケージを用いればより簡単にモダンなグラフを作成することができる。ggplot2は,tidyverseに含まれているので,tidyverseを読み込んでおく。
── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
✔ dplyr     1.1.4     ✔ readr     2.1.5
✔ forcats   1.0.0     ✔ stringr   1.5.1
✔ ggplot2   3.5.2     ✔ tibble    3.2.1
✔ lubridate 1.9.4     ✔ tidyr     1.3.1
✔ purrr     1.0.4     
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
 
 
また,練習用データとしてirisを用いる。
ヒストグラム
まず,150のアヤメの萼片の長さについて,ヒストグラムを作成してみよう。
ggplot2でグラフを作成するには,まずグラフ・オブジェクトを作成する。
g <- iris |>
  ggplot(aes(x = Sepal.Length))
 
ここでは,gというオブジェクトにグラフの素地を作成している。ggplot()関数の振舞いを理解するのは簡単ではないが,まずは決まったやり方を覚えよう。ggplot()関数を使うときには,まずデータフレームと,x軸,y軸にとる変数などをaes()で指定しグラフの素地を作成する。ここではggplot()関数にパイプでデータフレームirisを渡し,aes()でグラフのx軸に用いる変数Sepal.Lengthを指定している。ヒストグラムの場合に最低限必要なのは,x軸の変数(ヒストグラムを作成する変数)だけである。
次にヒストグラムを追加する。
g <- g + geom_histogram()
 
このように,ggplot2では作成した素地に,グラフやさまざまなパーツを足していくことで,好みのグラフを作成する。
グラフを表示させるにはplot()関数を用いる。
`stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
 
 
とくに何も指定しなくても,Rが適当に区間を分割してヒストグラムを描いてくれる。
このままだと見にくいので,ヒストグラムのビンの幅を変更し,枠線や塗り色を指定してみよう。ヒストグラムの細かい設定は,geom_histogram()の引数に指定すれば良い。
g <- iris |>
  ggplot(aes(x = Sepal.Length))+
  geom_histogram(
    binwidth = 0.25,
    color = "black",
    fill = "blue"
  )
plot(g)
 
次に,グラフの軸のラベルを変更して,グラフ・タイトルをつけてみよう。下のように,グラフ・オブジェクトにlabs()で軸ラベルやタイトルの設定を足す。
g <- iris |>
  ggplot(aes(x = Sepal.Length)) +
  geom_histogram(
    binwidth = 0.25,
    color = "black",
    fill = "blue"
  ) +
  labs(
    x = "萼片長",
    y = "頻度",
    title = "アヤメの萼片長の分布"
  )
plot(g)
 
箱ひげ図
一変数の分布だけに関心があるときにはヒストグラムが有用だが,複数の変数の分布を比較したいときには箱ひげ図が便利である。
まず,アヤメの萼片長の箱ひげ図を描いてみよう。
bp <- iris |>
  ggplot(aes(y = Sepal.Length)) +
  geom_boxplot()
plot(bp)
 
箱ひげ図の場合には,aes()でy軸に変数を指定してやるだけで良い。ここで,yではなくxを指定すれば,横向きの箱ひげ図を描くことができる。
次に,アヤメの種類ごとに箱ひげ図を描いて,萼片長の分布を比較しよう。
bp <- iris |>
  ggplot(
    aes(
      y = Sepal.Length, 
      fill = Species
    )
  ) +
  geom_boxplot()
plot(bp)
 
ggplot関数では,aes()にfillでファクター型変数を指定すれば,指定された変数のレベルごとに塗り色を分けてグラフを描くことができる。ここでは,Speciesを指定して,アヤメの種類ごとに塗り分けてグラフを描いている。
散布図
次に,変数の間の関係を分析するために,散布図を作成してみよう。ここでは,アヤメの萼片の長さを横軸に,萼片の幅を縦軸にとって,この二変数の関係を表す散布図を作成する。
scatter <- iris |>
  ggplot(
    aes(
      x = Sepal.Length,
      y = Sepal.Width
    )
  ) +
  geom_point()
plot(scatter)
 
散布図の場合は,aes()でx軸の変数とy軸の変数を指定する。あとは,geom_point()を足してやれば,単純な散布図を描くことができる。
散布図も,アヤメの種類ごとに色分けして描くことができる。色分けを行うには,aes()のcolorでファクター型変数を指定すれば良い。この場合は,Speciesによって点の色を変えるので,以下のように書く。
scatter <- iris |>
  ggplot(
    aes(
      x = Sepal.Length,
      y = Sepal.Width,
      color = Species
    )
  ) +
  geom_point() +
  labs(
    x = "萼片長",
    y = "萼片幅",
    title = "萼片の長さと幅の関係"
  )
plot(scatter)
 
ここでは,あわせてlabs()を足してやることにより,グラフタイトルや軸のラベルも指定している。
ggplot2は,使い方を覚えるというよりは,グラフを描く必要があるときに,どうすれば目的のグラフを描くことができるかをググりながら使えば良い。インターネット上にいくらでも丁寧な解説がある。