## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>")

## -----------------------------------------------------------------------------
library(islandcodes)

is_sids(c("Aruba", "Curacao", "Bonaire", "Brazil"))
is_snij(c("Aruba", "Curacao", "Bonaire", "Brazil"))

## -----------------------------------------------------------------------------
is_sids(c("AW", "CW", "BQ-BO", "AX", "BR"))

## -----------------------------------------------------------------------------
df <- data.frame(
  country  = c("Aruba", "Curacao", "Bonaire", "Sint Maarten", "Brazil"),
  variable = c(3.5, 3.1, 0.5, 1.2, 1900)
)

add_island_cols(df, "country")

## -----------------------------------------------------------------------------
library(countrycode)

df$iso2 <- countrycode(df$country, "country.name", "iso2c")
df$iso2  # note Bonaire collapses to NA in countrycode

# islandcodes recovers the BES cases via direct hyphenated lookup
add_island_cols(df, "country",
                cols = c("iso_code", "is_sids", "is_snij"))

## -----------------------------------------------------------------------------
nrow(small_islands(sids_only = TRUE))
nrow(small_islands(snij_only = TRUE))

head(small_islands(criteria = c(small = TRUE, island = TRUE, sovereign = TRUE)),
     8)

## -----------------------------------------------------------------------------
island_coords(c("Aruba", "Curacao", "Bonaire", "Sint Maarten"))

## -----------------------------------------------------------------------------
island_coords(c("PF", "US"))
island_coords(c("PF", "US"), which = "capital")

## -----------------------------------------------------------------------------
round(island_distance(c("AW", "CW", "BQ-BO", "SX", "BQ-SA", "BQ-SE")))

## -----------------------------------------------------------------------------
island_distance("Aruba", c("Curacao", "Bonaire", "Sint Maarten"))

## -----------------------------------------------------------------------------
island_distance("AW", "CW", which = "capital", unit = "nmi")

## -----------------------------------------------------------------------------
abc <- island_distance(c("AW", "CW", "BQ-BO"))
diag(abc) <- NA
apply(abc, 1, min, na.rm = TRUE)  # km to nearest of the other two