library(tidyverse)
noodid=read_csv("http://www.tlu.ee/~jaagup/andmed/muu/muusika/regiviisid.txt")
noodid %>% group_by(P1) %>% summarise(kogus=n()) %>% arrange(-kogus)
noodid %>% filter(tugiheli=="g") %>%  group_by(P16) %>% summarise(kogus=n()) %>% arrange(-kogus)
noodid %>% filter(tugiheli=="g") %>% unite("kolmik", P1, P2, P3) %>%
    group_by(kolmik) %>% summarise(kogus=n()) %>% arrange(-kogus) %>% head(7)


noodid %>% filter(tugiheli=="g") %>% unite("kolmik", P1, P2, P3) %>%
  group_by(kolmik) %>% summarise(kogus=n()) %>% arrange(-kogus) %>% head(7) %>%
  ggplot(aes(kolmik, kogus)) + geom_col()

noodid %>% filter(tugiheli=="g") %>% unite("kolmik", P1, P2, P3) %>%
  group_by(kolmik) %>% summarise(kogus=n()) %>% arrange(-kogus) %>% head(7) %>%
  .$kolmik -> kolmikud

noodid %>% filter(tugiheli=="g") %>% unite("kolmik", P1, P2, P3) %>%
  group_by(kolmik, P16) %>% summarise(kogus=n()) %>% arrange(-kogus)  %>%
  ggplot(aes(kolmik, kogus)) + geom_col() + facet_wrap(~P16)


noodid %>% filter(tugiheli=="g") %>% unite("kolmik", P1, P2, P3) %>%
  filter(kolmik %in% kolmikud, P16 %in% c("g", "h", "d", "a")) %>%
  group_by(kolmik, P16) %>% summarise(kogus=n()) %>% arrange(-kogus)  %>%
  ggplot(aes(kolmik, kogus)) + geom_col() + facet_wrap(~P16)


noodid %>% filter(tugiheli=="g") %>% mutate(kolmik=paste(P1, P2, P3, sep="_")) %>%
  filter(kolmik %in% kolmikud, P16 %in% c("g", "h", "d", "a")) %>%
  group_by(kolmik, P16) %>% summarise(kogus=n()) %>% arrange(-kogus)  %>%
  ggplot(aes(kolmik, kogus)) + geom_col() + facet_wrap(~P16)

noodid=read_csv("http://www.tlu.ee/~jaagup/andmed/muu/muusika/regiviisid.txt")
metaandmed=read_csv("http://www.tlu.ee/~jaagup/andmed/muu/muusika/regiviisidmeta.txt")


koos=metaandmed %>% inner_join(noodid, by=c("ID"="FKey"))

#metaandmed %>% View()

koos %>% group_by(kihelkond) %>% summarise(kogus=n()) %>% arrange(-kogus)

koos %>% mutate(kolmik=paste(P1, P2, P3, sep="_"))%>% group_by(kihelkond, kolmik) %>% summarise(kogus=n()) %>% arrange(-kogus)


koos %>% mutate(kolmik=paste(P1, P2, P3, sep="_"))%>% 
  filter(kolmik %in% kolmikud) %>%
  group_by(kihelkond, kolmik) %>% summarise(kogus=n()) %>% arrange(-kogus)

kolmikukogused=koos %>% mutate(kolmik=paste(P1, P2, P3, sep="_"))%>% 
  filter(kolmik %in% kolmikud) %>%
  group_by(kihelkond, kolmik) %>% summarise(kogus=n()) %>% arrange(-kogus)

kihelkonnakogused=koos %>% group_by(kihelkond) %>% summarise(kihelkonnakogus=n()) %>% arrange(-kogus)

kolmikukogused %>% inner_join(kihelkonnakogused, by=c("kihelkond"="kihelkond"))

kolmikukogused %>% inner_join(
  kihelkonnakogused %>% arrange(-kihelkonnakogus) %>% head(5), by=c("kihelkond"="kihelkond"))


kolmikukogused %>% inner_join(
  kihelkonnakogused %>% arrange(-kihelkonnakogus) %>% head(5), by=c("kihelkond"="kihelkond")) %>%
    mutate(osakaal=kogus/kihelkonnakogus)


kolmikukogused %>% inner_join(
  kihelkonnakogused %>% arrange(-kihelkonnakogus) %>% head(5), by=c("kihelkond"="kihelkond")) %>%
  mutate(osakaal=kogus/kihelkonnakogus) %>%
  ggplot(aes(kolmik, osakaal))+geom_col() + facet_wrap(~kihelkond) + coord_flip()




###########

metaandmed=read_csv("http://www.tlu.ee/~jaagup/andmed/muu/muusika/regiviisidmeta.txt")
koos=metaandmed %>% inner_join(noodid, by=c("ID"="FKey"))

noodid %>% filter(tugiheli=="g") %>% unite("kolmik", P1, P2, P3) %>%
  group_by(kolmik) %>% summarise(kogus=n()) %>% arrange(-kogus) %>% head(7) %>%
  .$kolmik -> kolmikud

kolmikukogused=koos %>% mutate(kolmik=paste(P1, P2, P3, sep="_"))%>% 
  filter(kolmik %in% kolmikud) %>%
  group_by(kihelkond, kolmik) %>% summarise(kogus=n()) %>% arrange(-kogus)

kihelkonnakogused=koos %>% group_by(kihelkond) %>% summarise(kihelkonnakogus=n()) %>%   
  arrange(-kogus)

print(kihelkonnakogused)
print(kolmikukogused)