ggplot2 110129
- 1. ggplot2 パッケージ作成者： Wickam さんの話を聞いて
- 2. 目次 自己紹介
- 3. 統数研でのセミナー
- 4. 今日の話の構成
- 6. インストール方法
- 7. 基本的な使い方
- 8. 具体的にやってみよう
- 9. まとめ
- 10. 参考資料
- 11. 自己紹介 名前 ：簑田 高志 （みのだ たかし）
- 14. 出身地 ：熊本県
- 15. 出身学部 ：法学部
- 16. 仕事 ：一般ユーザー向けの Web サービスの 企 画・運営、アナリスト
- 19. 統数研でのセミナー
- 21. 今回は2人の方が参加されました。 ggplot2のパッケージ開発者 Hadley Wickhamさん
- 22. 黄色い本を書かれた Uwe Liggesさん 概要は以下のURLに公開されています。 http://jasp.ism.ac.jp/~nakanoj/workshop10/2010Rmeeting.htm
- 25. 今日の話の構成
- 27. インストール〜使い方まで。
- 31. きれいなグラフを作りたい。
- 33. ggplot2とはどんなパッケージ？ ggplot2とは？ ggplot2 is a plotting system for R, based on the grammar of graphics, which tries to take the good parts of base and lattice graphics and none of the bad parts. It takes care of many of the fiddly details that make plotting a hassle (like drawing legends) as well as providing a powerful model of graphics that makes it easy to produce complex multi-layered graphics
- 34. 超意訳： グラフ文法をベースにした R の描画機能。簡単にきれいでパワフルなグラフがかけまっせ！ 出典： http://had.co.nz/ggplot2/
- 36. インストール方法
- 38. library(ggplot2) 初めてインストールする場合は、CRANのロケーションを選んでください。 Japan(Hyogo) Japan(Tsukuba) を選択してください。 （画面はMacOS版）
- 39. 基本的な使い方
- 40. 基本的な使い方 ggplot2には2つの構文があります。 ggplot([data] , aes(x = X 軸 , y = Y 軸 )) : データ定義 + geom_◯◯◯◯() : 出力するグラフ
- 41. + geom_◯◯◯◯() : 出力するグラフ 2
- 42. + xlim( 最小値 , 最大値 ) + ylim( 最小値 , 最大値 )
- 43. + xlab(X 軸ラベル )+ ylab(Y 軸ラベル ) などなど… 例） ggplot(movies, aes(x=mpaa, y=rating)) +geom_jitter(aes(colour=rating))+xlab(“対象年令”) + ylab(“レーティング”)
- 44. 基本的な使い方 ggplot2には2つの構文があります。 qplot(x , y , data = [data] , geom = グラフ ) 例） qplot(mpaa, rating, data=movies,geom=c("boxplot","jitter")) ggplot~,qplotも構文が違うだけで同じグラフを作ってくれます。
- 45. 具体的にやってみよう
- 49. carat cut color clarity depth table price x y z
- 50. 1 0.23 Ideal E SI2 61.5 55 326 3.95 3.98 2.43
- 51. 2 0.21 Premium E SI1 59.8 61 326 3.89 3.84 2.31
- 52. 3 0.23 Good E VS1 56.9 65 327 4.05 4.07 2.31
- 53. 4 0.29 Premium I VS2 62.4 58 334 4.20 4.23 2.63
- 54. 5 0.31 Good J SI2 63.3 58 335 4.34 4.35 2.75
- 55. 6 0.24 Very Good J VVS2 62.8 57 336 3.94 3.96 2.48
- 56. 具体的にやってみよう（棒グラフ） ■変数 cut , carat の 棒グラフ と、X軸の間隔を調整してみよう。 ･ qplot(cut, data = diamonds) ･ ggplot(diamonds , aes(cut)) + geom_bar() ･ qplot(carat, data = diamonds) ･ ggplot(diamonds , aes(carat)) + geom_bar() ･ qplot(carat, data = diamonds, binwidth = 1) ･ ggplot(diamonds , aes(x = carat)) + geom_bar(binwidth = 1) ･ qplot(carat, data = diamonds, binwidth = 0.1) ･ ggplot(diamonds , aes(x = carat)) + geom_bar(binwidth = 0.1) ･ qplot(carat, data = diamonds, binwidth = 0.01) ･ ggplot(diamonds , aes(x = carat)) + geom_bar(binwidth = 0.1)
- 57. 具体的にやってみよう（棒グラフ・六角形） ■変数 depth を Cutごとでみてみよう ･ qplot(depth, data = diamonds, binwidth = 0.2, fill = cut) + xlim(55, 70) ･ ggplot(diamonds , aes(depth , fill = cut)) + geom_bar(binwidth = 0.2) + xlim (55 , 70) ■変数 depth を Cutごとでみてみよう ･ qplot(depth, data = diamonds, binwidth = 0.2) + xlim(55, 70) + facet_wrap(~ cut) ･ ggplot(diamonds , aes(depth)) + geom_bar(binwidth = 0.2) + xlim(55 , 70) + facet_wrap(~cut) ■もうちょっとかっこ良く(hexbinパッケージが必要かも) 六角形 ･ qplot(log10(carat), log10(price), data = diamonds,geom = "hex", bins = 10) + facet_wrap(~ cut)
- 58. 具体的にやってみよう（散布図・平均直線・回帰直線） ■変数 price と carat の関係をみてみたい。（ 散布図 ） ･ qplot(price, carat, data = diamonds, colour = I("blue")) ･ ggplot(diamonds , aes(x = price , y = carat)) + geom_jitter(colour = I("blue") ) ■ラインを引いてみよう - 条件つき平均 ･ ggplot(diamonds , aes(price , carat)) + geom_jitter(colour = I("blue") ) + geom_smooth(colour = "red") - 回帰直線 ･ ggplot(diamonds , aes(price , carat)) + geom_jitter() + geom_smooth(method = "lm") このラインの詳細はこちら http://had.co.nz/ggplot2/geom_smooth.html
- 59. 具体的にやってみよう（透過グラフ） ■変数 price と carat の関係の強い所だけ見てみたい。 （グラフを透過的にしてみる） ･ qplot(price, carat, data = diamonds, alpha = I(1/10)) ･ qplot(price, carat, data = diamonds, alpha = I(1/50)) ･ qplot(price, carat, data = diamonds, alpha = I(1/100)) ･ qplot(price, carat, data = diamonds, alpha = I(1/250))
- 60. 具体的にやってみよう（箱ヒゲ図） ■最後にBoxplotとの重ね書き（箱ヒゲ図） ※ データを切り替えます：movie ･ ggplot(movies , aes(x = mpaa , y = rating)) + geom_jitter() ･ ggplot(movies , aes(x = mpaa , y = rating)) + geom_jitter(aes(colour = rating)) ･ ggplot(movies , aes(x = mpaa , y = rating)) + geom_jitter(aes(colour = rating)) + geom_boxplot()
- 61. まとめ
- 62. まとめ ggplot の構文 ggplot(data , aes(x軸, y軸)) + geom_(グラフ種類)~~
- 63. qplot(x軸 , y軸 , data , geom = (グラフ種類))~~ 様々なオプション Xlab , ylab , xlim , ylim, logなど… グラフの重ね書き ggplotの場合：~+ geom_jitter() + geom_boxplot()のように連結していけばよい。
- 64. qplot の場合： geom = c(“jitter” , “boxplot”) のようにc()で追加
- 65. 参考文献
- 66. 今日の話の構成 ggplot2 マニュアル http://had.co.nz/ggplot2/ 利用する時にとても参考になります。 統数研セミナー概要 http://jasp.ism.ac.jp/~nakanoj/workshop10/2010Rmeeting.htm セミナー資料 http://had.co.nz/courses/10-tokyo/ SlideShare(Hadly Wickham) http://www.slideshare.net/hadley/presentations