conc_corr <- read_csv("data/lang_pairwise_tile_correlations.csv",
                      col_names = c("tile1", "tile2", "corr", "lang1", "lang2")) 
  
conc_corr_ms <- conc_corr %>%
  group_by(tile1, tile2) %>%
  multi_boot_standard(col = "corr") %>%
  ungroup() %>%
  mutate(tile1 = as.factor(tile1))


ggplot(conc_corr_ms, aes(x = tile2, y = mean, color = tile1, 
                         group = tile1)) +
  scale_colour_discrete( name="Decile 2") +
  geom_line() +
  ylab("Cross-linguistic word-pairwise correlation\n (Pearson's r)") +
  geom_pointrange(aes(ymin = ci_lower, max = ci_upper), size = .2)  +
  scale_x_continuous(breaks=1:10) +
  ggtitle("Mean correlation by decile") +
  xlab("Decile 1") +
  theme_classic()