#Juhusliku eksleja mudel
library(tidyverse)
runif(1, 0, 360)
?runif

startpos= c(0,0)
samm = 100
angle=runif(1, 0, 360)

position = c(startpos[1] + samm*cos(angle), 
             startpos[2] + samm*sin(angle))
position
pikkus = sqrt(
  (position[1] - startpos[1])^2 + 
  (position[2] - startpos[2])^2  
)
pikkus



n=2
position = startpos
for(i in 1:n){
  angle=runif(1, 0, 360)
  position = c(position[1] + samm*cos(angle), 
               position[2] + samm*sin(angle))
  print(position)
}

pikkus = sqrt(
  (position[1] - startpos[1])^2 + 
  (position[2] - startpos[2])^2  
)
pikkus

#1000 korda paar sammu
pikkused = c()
katsed = c()
m=1000
for(j in 1:m){
  n=2
  position = startpos
  for(i in 1:n){
    angle=runif(1, 0, 360)
    position = c(position[1] + samm*cos(angle), 
                 position[2] + samm*sin(angle))
    print(position)
  }
  
  pikkus = sqrt(
    (position[1] - startpos[1])^2 + 
      (position[2] - startpos[2])^2  
  )
  pikkused[j] <- pikkus
  katsed[j] <- j
}
min(pikkused) 
max(pikkused)
mean(pikkused)
median(pikkused)
sd(pikkused)
sterr = sd(pikkused)/sqrt(length(pikkused))
sterr


data = tibble(pikkused, katsed)
data %>% ggplot(aes(katsed, pikkused)) + geom_col() + 
  geom_hline(yintercept = min(pikkused), colour="green") + 
  geom_hline(yintercept=max(pikkused), colour="blue") +
  geom_hline(yintercept=mean(pikkused), colour="yellow") + 
  geom_hline(yintercept=median(pikkused), colour="cyan") +
  geom_hline(yintercept=sd(pikkused), colour="magenta") + 
  geom_hline(yintercept=sterr, colour="red")

################
# 10 katset, kus 100 sammu
#võib ka 10 sammu, 100 katset?

pikkused = c()
katsed = c()
m=10
katseta <- function(m, n){
  for(j in 1:m){
    position = startpos
    for(i in 1:n){
      angle=runif(1, 0, 2*pi)
      position = c(position[1] + samm*cos(angle), 
                   position[2] + samm*sin(angle))
      print(position)
    }
    
    pikkus = sqrt(
      (position[1] - startpos[1])^2 + 
        (position[2] - startpos[2])^2  
    )
    pikkused[j] <- pikkus
    katsed[j] <- j
  }
  
  min(pikkused) 
  max(pikkused)
  mean(pikkused)
  median(pikkused)
  sd(pikkused)
  sterr = sd(pikkused)/sqrt(length(pikkused))
  data = tibble(pikkused, katsed)
  data %>% ggplot(aes(katsed, pikkused)) + geom_col() + 
    geom_hline(yintercept = min(pikkused), colour="green") + 
    geom_hline(yintercept=max(pikkused), colour="blue") +
    geom_hline(yintercept=mean(pikkused), colour="yellow") + 
    geom_hline(yintercept=median(pikkused), colour="cyan") +
    geom_hline(yintercept=mean(pikkused) + sd(pikkused), colour="magenta") + 
    geom_hline(yintercept=mean(pikkused) - sd(pikkused), colour="magenta") +
    geom_hline(yintercept=mean(pikkused) + sterr, colour="red") +
    geom_hline(yintercept=mean(pikkused) - sterr, colour="red")
}

katseta(100,2)
max(pikkused)