#interval poevrenja za m, kad je sigma=1 poznato animacija

sigma =1
m=0
n=100
N_CI= 50
beta=0.95
brojac=0

interval_poverenja_m <- function(...) {
xsr=mean(rnorm(n,m, sigma))
L <- xsr - qnorm((1+beta)/2)/ sqrt(n)
U <- xsr + qnorm((1+beta)/2)/ sqrt(n)
return(c(L, U))
}

CIS <- lapply(1:N_CI,interval_poverenja_m)



brojac=0

library(animation)
ani.options(interval = .25)
saveGIF({
  for (i in 1:N_CI) {
    plot(
      0,
      xlim = c(1, N_CI),
      ylim = c(-0.6, 0.6),
      xlab = "Uzorak",
      ylab = "Interval poverenja za m",
      main = "n"
    )
    abline(h = m, col = "gray", lwd = 2)
    x = 1:i
    y1 = unlist(lapply(CIS[1:i], function(ci) ci[[1]]))
    y2 = unlist(lapply(CIS[1:i], function(ci) ci[[2]]))
    color = unlist(lapply(CIS[1:i], function(ci) {
      if (all(ci > m) || all(ci < m)) {brojac=brojac+1;
        return("red") 
      }
      return("gray")
    }))
    segments(x, y1, x, y2, color, lwd = 3)
  }
})