Rで解析：タイルグリッドでアメリカ合衆国のヒートマップ「minimap」パッケージ

公開日 2023-01-20

Rの解析に役に立つ記事

タイルグリッドで州を示し、アメリカ合衆国のヒートマップを作成するパッケージです。

実行コマンドでは、指定した色をプロットする方法と、データの値を判断して付与した色をプロットする方法を紹介します。

パッケージバージョンは0.1.0。実行コマンドはR version 4.2.2で確認しています。

＜おすすめのRに関する書籍です＞

RStudioではじめるRプログラミング入門 | ギャレット・グロ-ルマンド, 大橋真也, 長尾高弘

Amazonでギャレット・グロ-ルマンド, 大橋真也, 長尾高弘のRStudioではじめるRプログラミング入門。アマゾンならポイント還元本が多数。

www.amazon.co.jp

パッケージのインストール

下記、コマンドを実行してください。

#パッケージのインストール
install.packages("devtools")
devtools::install_github("seankross/minimap")

実行コマンド

詳細はコメント、パッケージのヘルプを確認してください。

#パッケージの読み込み
library("minimap")

###データ例の作成#####
#パッケージ付属のssmデータから州名を抽出
data("ssm")
StateNames <- unique(ssm$State)
#作成
TestData <- data.frame(StateNames = StateNames,
                       Data = sample(1:6, 51, replace = TRUE),
                       StateColors = I(sample(c("#4b61ba", "#a87963", "#deb7a0"),
                                              51, replace = TRUE)))
########

##指定した色でプロットする方法
#プロット:miniusaコマンド
#各州名データ:statesオプション
#塗り色:state_colorsオプション
#州名の表示設定:state_namesオプション
#州名の色:state_name_colorsオプション
#州名のサイズ:state_name_cexオプション
miniusa(states = TestData[, 1], state_colors = TestData[, 3],
        state_names = TRUE, state_name_colors = rep("white", 51),
        state_name_cex = 1)

##値で色を判断して塗り分け
#判断のfunctionを作成
ColorJudge <- function(x){
  if(x <= 2){
    "#90736a"

    } else if(x >= 2 & x <= 4) {

      "#28231e"

  } else {

    "#4b61ba"
    }
}

#処理
DataColor <- sapply(TestData[, 2], ColorJudge)
#プロット
miniusa(states = TestData[, 1], state_colors = DataColor,
        state_names = TRUE, state_name_colors = rep("white", 51),
        state_name_cex = 1)

#カラーバーを追加
#fieldsパッケージがなければインストール
if(!require("fields", quietly = TRUE)){
  install.packages("fields");require("fields")
}
image.plot(z = 1:5, nlevel = 6, col = c("#90736a", "#28231e", "#4b61ba"),
           legend.only = TRUE)