library(ggplot2)
library(tidyverse)
tibble(xmin = rep(1, 100), xmax = rep(10, 100), ymin = 1:100, ymax = ymin + 1, id = factor(c(rep(1, 10), rep(2, 30), rep(3, 60)))) %>% 
    ggplot(aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax, fill = id)) +
    geom_rect() -> plt
print(plt)

library(tidyverse)
expand.grid(1:100, 1:100) %>% 
  as_tibble() %>% 
  rowwise() %>% 
  mutate(stat_diff = chisq.test(c(Var1, Var2))$p.value < 0.05,
         value = abs(Var1-Var2)) %>% 
  ggplot(aes(Var1, Var2, fill = stat_diff, label = round(value, 2))) +
  geom_tile()+
  ggtitle("χ2 < 0.05")