【R】世界の火山図

A weekly social data project in RであるTidy Tuesdayに世界の火山のデータがあるので、世界地図にプロットしてみます。

library(tidyverse)
library(ggplot2)
volcano_raw  <- readr::read_csv("https://raw.githubusercontent.com/rfordatascience/tidytuesday/master/data/2020/2020-05-12/volcano.csv")
volcano_raw %>% count(primary_volcano_type, sort= TRUE)

volcano_df <- volcano_raw %>% mutate_if(is.character, factor)

world <- map_data("world")

ggplot() +
  geom_map(data = world, map = world,
           aes(long, lat, map_id = region),
           color = "white", fill = "gray50", alpha = 0.2) +
  geom_point(data = volcano_df,
             aes(longitude, latitude, color = primary_volcano_type),
             alpba = 0.8)
  

これを見ると、日本は地震が頻発するのに納得します。。。。