saveGIF({
  for(nr in 2:10){
    print(sonad %>% filter(sonapikkus==nr) %>% 
            group_by(sona, taishaalikuid, sulghaalikuid) %>%
            summarise(kogus=n()) %>% 
            ggplot(aes(taishaalikuid, jitter(sulghaalikuid, factor = 2.5), 
                       label=sona, size=kogus))+  geom_text()+
            xlim(0, 5) + ylim(0, 5)+ggtitle(paste("Sõnas tähti", nr)))
  }
}, movie.name="d:/animatsioon1.gif", interval=4, ani.width=800, ani.height=600 )
#kasutage vaid Kungla rahva sõnu

saveGIF({
  for(nr in 2:8){
    print(sonad %>% filter(sonapikkus==nr, lugu=="kungla") %>% 
            group_by(sona, taishaalikuid, sulghaalikuid) %>%
            summarise(kogus=n()) %>% 
            ggplot(aes(taishaalikuid, jitter(sulghaalikuid, factor = 2.5), 
                       label=sona, size=kogus))+  geom_text()+
            xlim(0, 5) + ylim(0, 5)+ggtitle(paste("Sõnas tähti", nr)))
  }
}, movie.name="d:/animatsioon1.gif", interval=4, ani.width=800, ani.height=600 )
#kasutage vaid Kungla rahva sõnu

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


noodid %>% select(P1:P16) %>% unite("viimased", P13:P16)

#Kuvage kahanevas järjekorras, millist viimaste kombinatsiooni mitu korda esineb

noodid %>% select(P1:P16) %>% na.omit() %>% unite("viimased", P13:P16) %>% 
  group_by(viimased) %>% summarise(kogus=n()) %>% arrange(-kogus)

noodid %>% select(P1:P16) %>% filter(P1=="2d") %>% na.omit() %>% unite("viimased", P13:P16) %>% 
  group_by(viimased) %>% summarise(kogus=n()) %>% arrange(-kogus)


noodid

noodid %>% select(IDviis, P1:P16) %>% unite("algus", P1:P3) %>% unite("viimased", P13:P16)

noodid %>% unite("algus", P1:P3) %>%  filter(IDviis==68) %>% .$algus

koht="P4"
noodid %>% unite("algus", P1 : get("P4") ) %>%  filter(IDviis==68) %>% .$algus

noodid %>% select(P4)


ggplot()+xlim(0, 10)+ylim(0, 10)+geom_point()+annotate("text", 5, 5, label="tere")

#Pange videole algusesse selgitav tekst
saveGIF({
  for(nr in 1:3){
    print(ggplot()+xlim(0, 10)+ylim(0, 10)+geom_point()+
            annotate("text", 5, 5, label="tere", size=7))
  }
  for(nr in 2:8){
    print(sonad %>% filter(sonapikkus==nr, lugu=="kungla") %>% 
            group_by(sona, taishaalikuid, sulghaalikuid) %>%
            summarise(kogus=n()) %>% 
            ggplot(aes(taishaalikuid, jitter(sulghaalikuid, factor = 2.5), 
                       label=sona, size=kogus))+  geom_text()+
            xlim(0, 5) + ylim(0, 5)+ggtitle(paste("Sõnas tähti", nr)))
  }
  #Pange ka lõppu selgitav tekst
}, movie.name="d:/animatsioon1.gif", interval=4, ani.width=800, ani.height=600 )
#kasutage vaid Kungla rahva sõnu

Valimis oli sada sõna, millest 15 osutusid Kungla rahva omaks
prop.test(15, 100)

95 percent confidence interval:
  0.0891491 0.2385308

95% tõenäosusega saame mõõdetud andmete põhjal väita, et 
Kungla rahva sõnade osakaal kogu sõnaloendis on 8,9 kuni 23,9 protsenti

prop.test(1, 10)
prop.test(6, 10)
prop.test(15, 30)
prop.test(15, 30)$conf.int[1]
library(tidyverse)
tabel=tibble(x1=c(0.27, 0.33), x2=c(0.86, 0.67))
tabel %>% ggplot(aes(x1, 0, xend=x2, yend=0))+geom_curve()+
   xlim(0, 1)+ylim(-0.05, 0.05)+
       annotate("text", 0.35, -0.025, label="rida") +
       annotate("text", 0.35, 0.025, label="salm") +
       annotate("segment", x=0.35, y=-0.02, xend=0.27, yend=0, 
                 arrow=arrow( angle=10))
#Lisage salmile ka nool