A de Finetti diagram showing Hardy-Weinberg Equilibrium

HWE_deFinetti.R

# de Finetti diagram
library(dplyr)
library(ggplot2)
library(ggtern)
library(magrittr)

HWEfrequencies <- data.frame(
  p = seq(0, 1, 0.01)
) %>%
  mutate(q = 1 - p) %>%
  mutate(AA = p^2, Aa = 2 * p * q, aa = q^2)

ggtern(
  data = HWEfrequencies,
  aes(x = AA, y = Aa, z = aa)
) + geom_line(size = 1) +
  labs(
    x = "aa (%)", y = "Aa (%)", z = "AA (%)",
    title = "Hardy-Weinberg Equilibrium",
    subtitle = "de Finetti diagram"
  ) +
  theme_bvbw()