【R】United KingdomでUFO

1. はじめに

TidyTuesdayにUFOの目撃データという面白いデータが載っていたので、それを地図上にプロットしてみた。

2. データ

まずはデータの取得。TidyTuesdayからもらってきます。そして、イギリス（gb）のデータだけ抽出します。各町ごとの目撃数を集計しています。後々に地図にプロットするために緯度経度の情報も残しておきます。

library(ggplot2)
library(ggrepel)
library(sf)
library(tidyverse)
library(wesanderson)
library(kableExtra)

ufo_sightings <- readr::read_csv("https://raw.githubusercontent.com/rfordatascience/tidytuesday/master/data/2019/2019-06-25/ufo_sightings.csv")

sightings <- ufo_sightings %>%
  filter(country == "gb") %>%
  group_by(city_area, latitude, longitude) %>%
  summarise(count=n()) %>% 
  ungroup()

次に、地図の準備をします。地図データは、いつものようにDIVA-GISからとってきました。いつもお世話になっております。

map <- read_sf("GBR_adm/GBR_adm1.shp") # GBのシェープファイル

dat <- map %>%
  group_by(NAME_1) %>% 
  mutate(
    centroid = st_centroid(geometry),
    x = st_coordinates(centroid)[, 1],
    y = st_coordinates(centroid)[, 2]
  ) 

最後にプロットしてみます。目撃情報の集計結果を表にするので、そのデータもクリップボードに保持します。

dat %>%
  ggplot() +
  geom_sf(aes(fill = NAME_1)) +  # 人口毎に色分け
  coord_sf(datum = NA) + 
  theme_void()+
  geom_point(data = sightings, aes(x = longitude, y = latitude, size = count), col="yellow") +
  scale_fill_manual(name = "country", values = wes_palette("Moonrise2"))+
  geom_text_repel(aes(x = x, y = y, label = NAME_1), col="red", size = 1.5)


table_df<-sightings %>% 
  data.frame(city_area="city_area", count="count")

library(clipr)
table_df %>% 
  kable(align = "c", row.names=FALSE) %>%
  kable_styling(full_width = F) %>%
  column_spec(1, bold = T) %>%
  collapse_rows(columns = 1, valign = "middle") %>% 
  write_clip

/

4. さいごに

イギリスでもUFOの目撃情報が多いことにビックリ（笑）！