create a map of linguae francae of Dagestan for Nina Dobrushina

code.R

# 30.05.2022, George Moroz
# GPL-3.0 license

library(tidyverse)
library(lingtypology)

# load data ---------------------------------------------------------------
values <- read_tsv("values.tsv")
individual_points <- read_csv("individual_points.csv")
# from Moroz, Verhees (2020), see https://github.com/sverhees/master_villages
data_for_extrapolation <- read_tsv("https://raw.githubusercontent.com/sverhees/master_villages/master/data/villages.tsv")

# merge and clean for lingtypology ----------------------------------------
data_for_extrapolation %>% 
  filter(!kutans) %>% 
  inner_join(values) %>% 
  mutate(lang = case_when(lang == "Azerbaijani" ~ "Northern Azerbaijani",
                          lang == "Dargwa" ~ "Nuclear North Dargwa",
                          lang == "Khwarshi" ~ "Khwarshi-Inkhoqwari",
                          TRUE ~ lang)) ->
  for_map

# create a map ------------------------------------------------------------
map.feature(languages = for_map$lang,
            latitude = for_map$lat,
            longitude = for_map$lon,
            density.estimation = for_map$value,
            density.estimation.color = c("#1f77b4", "#ff7f0e","#2ca02c"),
            density.method = "kernal density estimation",
            density.points = FALSE, 
            density.legend.position = "bottomright",
            tile="Stamen.TonerLite", 
            minimap = TRUE, 
            minimap.position = "topright") %>%   
  map.feature(languages = individual_points$lang,
              latitude = individual_points$lat,
              longitude = individual_points$lon,
              label = individual_points$lang,
              pipe.data = .,
              tile="Stamen.TonerLite")

# Mehtods are hidden, but retrievable from ?map.feature, look at 
# ?lingtypology::polygon.points_kde()
# ?lingtypology::polygon.points_fd()
# If you have in mind something better, I would love to include it in lingtypology.

# you can see data used in extropolation ----------------------------------
map.feature(languages = for_map$lang,
            features = for_map$value,
            latitude = for_map$lat,
            longitude = for_map$lon,
            density.estimation = for_map$value,
            density.estimation.color = c("#1f77b4", "#ff7f0e","#2ca02c"),
            density.method = "kernal density estimation",
            density.points = TRUE, 
            density.legend.position = "bottomright",
            label = for_map$village,
            tile="Stamen.TonerLite", 
            minimap = TRUE, 
            minimap.position = "topright",
            width = 2)

# Afterwards I used some tricks with repeating some observations in order to 
# get more smooth picture you see in my tweet

individual_points.csv

values.tsv