Skip to content

Instantly share code, notes, and snippets.

## Replication of https://archive.org/details/ost-geography-completegeograph00sauluoft/page/n574/mode/1up
# Inspired by @PaoloCrosetto's question and @geokaramanis's implementation
# https://twitter.com/PaoloCrosetto/status/1448594997251657731
# https://twitter.com/geokaramanis/status/1447913599712825349
library(tidyverse)
countries <- c("Belgium", "England", "Japan", "Italy", "China", "Germany", "France", "India", "Spain", "Phil. Is'ds", "Russia", "Cuba", "U. States", "Mexico", "Hawaii", "C.Colony", "Brazil", "Argentina", "Canada", "Australia")
densities <- c(563,500,284,280,270,250,186,184,89,72,54,36,20,16,16,5,4.5,3,2,1.23)
df <- tibble(country=factor(countries, levels=countries), # Trick to keep the order