podaci = dataN$x
n = length(podaci)
# uzoracka medijana (ocena)
Tn = median(podaci)

B = 1000
T_B = rep(0, B)
for (i in 1:B) {
  # reuzorkujemo i racunamo vrednost ocene
  x_B = sample(podaci, n, replace = T)
  T_B[i] = median(x_B)
}
(se = sqrt(sum((T_B - mean(T_B))^2) / B))
# tj.
(se = sqrt( (B - 1) / B * var(T_B) ))