この記事を参考に、ごく基本的な地図作成のスクリプト。 データソースはUSGSの地震データで、過去１時間から１日、７日、３０日、また、マグニチュードの下限についても、いくつか選択できる。 ここでは、Past 7 Days > All Earthquakes を用いる。 ダウンロ…

政府統計 e-Stat では、統計データがエクセル形式などのファイルのほかにデータベースとして提供されているものもある。後者のデータにAPI機能でアクセスするときに、estatapi パッケージを使うとRで比較的簡単に必要とするデータを抽出・取得できる。 e-Sta…

sf パッケージを用いたコロプレス・マップのプロット。 今回は千葉県市区町村の高齢化率を例にとる。 まずは、千葉県行政区分のシェイプファイルが必要だが、これは ESRI ジャパンが無償配布している「全国市区町村界データ」から、千葉県の部分を取り出して…

これの第２章に、アカウミガメの移動軌跡の例が載っている。 発信機をつけて人工衛星で追跡したら、メキシコから日本まで、太平洋横断の旅が観察された。その経過をプロットするというもの。 データファイルとRスクリプトは上記サイトからダウンロードできる…

以前の記事「関東甲信・梅雨入り日の推移」の続き。 気象庁のウェブページにある表： これをスクレイピングして、各年の梅雨の期間（日数）を求め、これと降水量の平年比との関係を調べる。 library(rvest)library(dplyr)library(stringr)library(ggplot2) d…