library(tidyverse)
 noodid=read_csv("http://www.tlu.ee/~jaagup/andmed/muu/muusika/regiviisid.txt")
 meta=read_csv("http://www.tlu.ee/~jaagup/andmed/muu/muusika/regiviisidmeta.txt")
 noodid %>% group_by(P1) %>% summarise(kogus=n()) %>% arrange(desc(kogus))
 noodid %>% filter(tugiheli=="g") %>% group_by(P1, P2, P3) %>%
     summarise(kogus=n()) %>% arrange(desc(kogus)) %>% head(7)
 
 noodid %>% filter(tugiheli=="g") %>% group_by(P1, P2, P3) %>%
   summarise(kogus=n()) %>% arrange(desc(kogus)) %>% head(7) %>%
     unite("kolmik", P1:P3) %>% 
      mutate(kolmik=factor(kolmik, levels=kolmik[order(kogus)])) %>%
      ggplot(aes(kolmik, kogus)) + geom_col() + coord_flip()
 
 noodid %>% filter(tugiheli=="g") %>% group_by(P16) %>% 
   summarise(kogus=n()) %>% arrange(desc(kogus)) %>% head(4) %>%
    .$P16 -> lopud
 
for(lopp in lopud){ 
 joonis=noodid %>% filter(tugiheli=="g", P16==lopp) %>% group_by(P1, P2, P3) %>%
   summarise(kogus=n()) %>% arrange(desc(kogus)) %>% head(7) %>%
   unite("kolmik", P1:P3) %>% 
   mutate(kolmik=factor(kolmik, levels=kolmik[order(kogus)])) %>%
   ggplot(aes(kolmik, kogus)) + geom_col() + ggtitle(lopp)+ coord_flip() 
 print(joonis)
}
 
meta %>% group_by(kihelkond) %>% summarise(kogus=n()) %>%
   arrange(desc(kogus)) %>% head(5) %>% .$kihelkond -> suuredkihelkonnad

noodid %>% filter(tugiheli=="g") %>% group_by(P1, P2, P3) %>%
  summarise(kogus=n()) %>% arrange(desc(kogus)) %>% head(7) %>%
  unite("kolmik", P1:P3) %>% .$kolmik -> suuredkolmikud
 
noodid %>% filter(tugiheli=="g") %>% 
  inner_join(meta, by=c("FKey"="ID")) %>% filter(kihelkond %in% suuredkihelkonnad) %>%
  unite("kolmik", P1:P3) %>% filter(kolmik %in% suuredkolmikud) %>%
  group_by(kolmik, kihelkond) %>% summarise(kogus=n()) %>% arrange(desc(kogus)) -> tabel2


tabel2 %>% mutate(kolmik=factor(kolmik, levels=rev(suuredkolmikud))) %>%
  ggplot(aes(kolmik, kogus)) + geom_col() + coord_flip()+ facet_wrap(~kihelkond)

tabel2 %>% pivot_wider(names_from=kihelkond, values_from=kogus, values_fill=0)