Hatena::ブログ(Diary)

盆栽日記

2016-02-12

カーネル密度推定とヒートマップ

サーモグラフィみたいなヒートマップをggplot2で描きたい。

たとえば以下のirisデータのプロット

ggplot(iris, aes(x = Sepal.Length, y = Sepal.Width))+ 
  geom_point(color = "black")

f:id:dichika:20160212175625p:image:w360

点が密集しているところを赤く、薄いところを緑で塗りたい。

こういう時はstat_density2dを使う。

ggplot(iris, aes(x = Sepal.Length, y = Sepal.Width))+ 
  stat_density2d(aes(fill=..level..), geom="polygon") +
  scale_fill_continuous(low = "green", high = "red") +
  geom_point(color = "black")

f:id:dichika:20160212175622p:image:w360

カーネル密度推定した結果を色分けしてくれる。

はてなユーザーのみコメントできます。はてなへログインもしくは新規登録をおこなってください。

トラックバック - http://d.hatena.ne.jp/dichika/20160212/p1
Connection: close