最近３０日間の地震分布

データの配布を行なっている http://earthquake.usgs.gov/earthquakes/feed/ で、"Past 30 days" のデータが提供されるようになった。
さらに、更新は１５分ごとと、かなりキメ細かい対応がなされている。
今回は、その「最近３０日間」に起こった地震のマグニチュードと震源を Google Map にプロットした。

library(ggmap)

## Loading required package: ggplot2

## Find out what's changed in ggplot2 with
## news(Version == "0.9.1", package = "ggplot2")

(map = ggmap(get_googlemap(center = c(137.5, 35.5), zoom = 5)) + 
    geom_point(data = read.csv("http://earthquake.usgs.gov/earthquakes/feed/csv/2.5/month"), 
        aes(Lon, Lat, size = Magnitude, colour = Depth)))

## Warning message: Removed 1010 rows containing missing values (geom_point).