sandbox/examples/arrests-lra-polygon.r

In corybrunson/ordr: A Tidyverse Extension for Ordinations and Biplots

# state abbreviations
state <- data.frame(
  .name = state.name,
  .abb = state.abb
)
# Log-ratio analysis of 1973 violent crime arrests in the United States
USArrests %>%
  subset(select = -UrbanPop) %>%
  lra() %>%
  print() -> arrests_lra
arrests_lra %>%
  as_tbl_ord() %>%
  augment() %>%
  left_join_rows(state, by = ".name") %>%
  print() -> arrests_lra
# Adapt Exhibit 7.1 in Greenacre (2010)
arrests_lra %>%
  confer_inertia(1) %>%
  ggbiplot(sec.axes = "cols", scale.factor = 1/20) +
  theme_bw() +
  geom_rows_text(
    aes(label = .abb),
    size = 3, color = "darkgreen", alpha = .5
  ) +
  geom_cols_polygon(fill = NA, linetype = "dashed", color = "brown4") +
  geom_cols_text(
    aes(label = .name),
    color = "brown4", fontface = "bold"
  ) +
  ggtitle(
    "Log-ratio analysis of violent crime arrest rates",
    "United States, 1973"
  ) +
  guides(color = "none", size = "none")