library(tidyverse)
sonad=read_csv("http://www.tlu.ee/~jaagup/andmed/keel/kunglarahvas_lambipirn_pikkused_haalikud.txt")

sapply(1:10, function(x){
  sonad %>% filter(lugu=="kungla") %>% sample_n(10) %>% summarise(k=mean(sonapikkus)) %>% .$k
}) %>% hist()

estonia=read_csv("http://www.tlu.ee/~jaagup/andmed/muu/estonia_reisijad.txt")

vanused=estonia %>% sample_n(100) %>% .$Age
shapiro.test(vanused)
hist(vanused)


vanused=estonia %>% sample_n(100) %>% .$Age
ks.test(vanused, "gamma", 3, 2)
hist(vanused)

kogus=100
ooteajad=runif(kogus, 0, max=10)
soiduajad=rnorm(kogus, 8, 2)
koguajad=ooteajad+soiduajad
print(summary(koguajad))
hist(koguajad)
print(sum(koguajad>18))
t.test(koguajad)


jalgsiajad=rnorm(kogus, 14, 2)
sum(jalgsiajad>18)
t.test(jalgsiajad)


t1=tibble(reanr=1:length(koguajad), ajad=sort(koguajad, decreasing=TRUE))
t1 %>% ggplot(aes(reanr, ajad))+geom_col()
t1 %>% head(10) %>% ggplot(aes(reanr, ajad))+geom_col()

t2=tibble(koguaeg=koguajad, ooteaeg=ooteajad, soiduaeg=soiduajad)
t2=t2 %>% arrange(-koguaeg) %>% head(10) %>% mutate(reanr=row_number())

t2 %>% select(reanr, ooteaeg, soiduaeg) %>% gather(tunnus, vaartus, -reanr) %>%
   ggplot(aes(reanr, vaartus, fill=tunnus)) + geom_col()